Operational calculus

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.

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.

The key element of the operational calculus is to consider differentiation as an operator p = ⁠d/dt⁠ acting on functions. Linear differential equations can then be recast in the form of "functions" F(p) of the operator p acting on the unknown function equaling the known function. Here, F is defining something that takes in an operator p and returns another operator F(p).

Solutions are then obtained by making the inverse operator of F act on the known function. The operational calculus generally is typified by two symbols: the operator p, and the unit function 1. The operator in its use probably is more mathematical than physical, the unit function more physical than mathematical. The operator p in the Heaviside calculus initially is to represent the time differentiator ⁠d/dt⁠. Further, it is desired for this operator to bear the reciprocal relation such that p⁻¹ 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 H(t), such that H(t) = 0 if t < 0, and H(t) = 1 if t > 0.

The simplest example of application of the operational calculus is to solve: p y = H(t), which gives y = p − 1 ⁡ H = ∫ 0 t H ( u ) d u = t H ( t ) . {\displaystyle y=\operatorname {p} ^{-1}H=\int _{0}^{t}H(u)\,\mathrm {d} u=t\,H(t).}

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

Continuing to treat p as if it were a variable, p p − a H ( t ) = 1 1 − a p H ( t ) , {\displaystyle {\frac {\operatorname {p} }{\operatorname {p} -a}}H(t)={\frac {1}{1-{\frac {a}{\operatorname {p} }}}}\,H(t),} which can be rewritten by using a geometric series expansion: 1 1 − a p H ( t ) = ∑ n = 0 ∞ a n p − n ⁡ H ( t ) = ∑ n = 0 ∞ a n t n n ! H ( t ) = e a t H ( t ) . {\displaystyle {\frac {1}{1-{\frac {a}{\operatorname {p} }}}}H(t)=\sum _{n=0}^{\infty }a^{n}\operatorname {p} ^{-n}H(t)=\sum _{n=0}^{\infty }{\frac {a^{n}t^{n}}{n!}}H(t)=e^{at}H(t).}

Using partial fraction decomposition, one can define any fraction in the operator p and compute its action on H(t).

Moreover, if the function 1/F(p) has a series expansion of the form 1 F ( p ) = ∑ n = 0 ∞ a n p − n , {\displaystyle {\frac {1}{F(\operatorname {p} )}}=\sum _{n=0}^{\infty }a_{n}\operatorname {p} ^{-n},} it is straightforward to find 1 F ( p ) H ( t ) = ∑ n = 0 ∞ a n t n n ! H ( t ) . {\displaystyle {\frac {1}{F(\operatorname {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 p, thus establishing a connection between operational calculus and fractional calculus.

Using the Taylor expansion, one can also verify the Lagrange–Boole translation formula, e^ap f(t) = f(t + a), so the operational calculus is also applicable to finite-difference equations and to electrical engineering problems with delayed signals.

• Calculus of finite differences
• Umbral calculus

• IV Lindell HEAVISIDE OPERATIONAL RULES APPLICABLE TO ELECTROMAGNETIC PROBLEMS
• Ron Doerfler Heaviside's Calculus
• Jack Crenshaw essay showing use of operators More On the Rosetta Stone