Operational calculus

Template:Short description Operational calculus, also known as operational analysis, is a technique by which problems in analysis, in particular differential equations, are transformed into algebraic problems, usually the problem of solving a polynomial equation.

History

The idea of representing the processes of calculus, differentiation and integration, as operators has a long history that goes back to Gottfried Wilhelm Leibniz. The mathematician Louis François Antoine Arbogast was one of the first to manipulate these symbols independently of the function to which they were applied.^[1]

This approach was further developed by Francois-Joseph Servois who developed convenient notations.^[2] Servois was followed by a school of British and Irish mathematicians including Charles James Hargreave, George Boole, Bownin, Carmichael, Doukin, Graves, Murphy, William Spottiswoode and Sylvester.

Treatises describing the application of operator methods to ordinary and partial differential equations were written by Robert Bell Carmichael in 1855^[3] and by Boole in 1859.^[4]

This technique was fully developed by the physicist Oliver Heaviside in 1893, in connection with his work in telegraphy.

Guided greatly by intuition and his wealth of knowledge on the physics behind his circuit studies, [Heaviside] developed the operational calculus now ascribed to his name.^[5]

At the time, Heaviside's methods were not rigorous, and his work was not further developed by mathematicians. Operational calculus first found applications in electrical engineering problems, for the calculation of transients in linear circuits after 1910, under the impulse of Ernst Julius Berg, John Renshaw Carson and Vannevar Bush.

A rigorous mathematical justification of Heaviside's operational methods came only after the work of Bromwich that related operational calculus with Laplace transformation methods (see the books by Jeffreys, by Carslaw or by MacLachlan for a detailed exposition). Other ways of justifying the operational methods of Heaviside were introduced in the mid-1920s using integral equation techniques (as done by Carson) or Fourier transformation (as done by Norbert Wiener).

A different approach to operational calculus was developed in the 1930s by Polish mathematician Jan Mikusiński, using algebraic reasoning.

Norbert Wiener laid the foundations for operator theory in his review of the existential status of the operational calculus in 1926:^[6]

The brilliant work of Heaviside is purely heuristic, devoid of even the pretense to mathematical rigor. Its operators apply to electric voltages and currents, which may be discontinuous and certainly need not be analytic. For example, the favorite corpus vile on which he tries out his operators is a function which vanishes to the left of the origin and is 1 to the right. This excludes any direct application of the methods of Pincherle…

Although Heaviside’s developments have not been justified by the present state of the purely mathematical theory of operators, there is a great deal of what we may call experimental evidence of their validity, and they are very valuable to the electrical engineers. There are cases, however, where they lead to ambiguous or contradictory results.

Principle

The key element of the operational calculus is to consider differentiation as an operator Template:Math acting on functions. Linear differential equations can then be recast in the form of "functions" Template:Math of the operator Template:Math acting on the unknown function equaling the known function. Here, Template:Math is defining something that takes in an operator Template:Math and returns another operator Template:Math.

Solutions are then obtained by making the inverse operator of Template:Mvar act on the known function. The operational calculus generally is typified by two symbols: the operator Template:Math, and the unit function Template:Math. The operator in its use probably is more mathematical than physical, the unit function more physical than mathematical. The operator Template:Math in the Heaviside calculus initially is to represent the time differentiator Template:Math. Further, it is desired for this operator to bear the reciprocal relation such that Template:Math denotes the operation of integration.^[5]

In electrical circuit theory, one is trying to determine the response of an electrical circuit to an impulse. Due to linearity, it is enough to consider a unit step function Template:Math, such that Template:Math if Template:Math, and Template:Math if Template:Math.

The simplest example of application of the operational calculus is to solve: Template:Math, which gives $y = p^{- 1} H = \int_{0}^{t} H (u) d u = t H (t) .$

From this example, one sees that $p^{- 1}$ represents integration. Furthermore Template:Mvar iterated integrations is represented by $p^{- n},$ so that $p^{- n} H (t) = \frac{t^{n}}{n!} H (t) .$

Continuing to treat Template:Math as if it were a variable, $\frac{p}{p - a} H (t) = \frac{1}{1 - \frac{a}{p}} H (t),$ which can be rewritten by using a geometric series expansion: $\frac{1}{1 - \frac{a}{p}} H (t) = \sum_{n = 0}^{\infty} a^{n} p^{- n} H (t) = \sum_{n = 0}^{\infty} \frac{a^{n} t^{n}}{n!} H (t) = e^{a t} H (t) .$

Using partial fraction decomposition, one can define any fraction in the operator Template:Math and compute its action on Template:Math.

Moreover, if the function Template:Math has a series expansion of the form $\frac{1}{F (p)} = \sum_{n = 0}^{\infty} a_{n} p^{- n},$ it is straightforward to find $\frac{1}{F (p)} H (t) = \sum_{n = 0}^{\infty} a_{n} \frac{t^{n}}{n!} H (t) .$

Applying this rule, solving any linear differential equation is reduced to a purely algebraic problem.

Heaviside went further and defined fractional power of Template:Math, thus establishing a connection between operational calculus and fractional calculus.

Using the Taylor expansion, one can also verify the Lagrange–Boole translation formula, Template:Math, so the operational calculus is also applicable to finite-difference equations and to electrical engineering problems with delayed signals.