The Old Laplace Transform 2.4 2001/07/06 2002/07/16 Doug Daniels rainking@alumni.rice.edu Steven Cox cox@rice.edu Doug Daniels rainking@alumni.rice.edu Steven Cox cox@rice.edu characteristic polynomial determinant eigenvalues Laplace transform MATLAB symbolic toolbox This module examines the Laplace Transform, an analytical tool that produces exact solutions for small, closed-form, tractable systems. We use the Laplace transform to move toward a solution for the nerve fiber potentials modeled by the dynamic Strang Quartet in the earlier module of the same name. The Laplace Transform is typically credited with taking dynamical problems into static problems. Recall that the Laplace Transform of the function h is ℒ h s t 0 s t h t . MATLAB is very adept at such things. For example: The Laplace Transform in MATLAB >> syms t >> laplace(exp(t)) ans = 1/(s-1) >> laplace(t*(exp(-t)) ans = 1/(s+1)^2 The Laplace Transform of a matrix of functions is simply the matrix of Laplace transforms of the individual elements. Laplace Transform of a matrix of functions ℒ t t t 1 s 1 1 s 1 2 Now, in preparing to apply the Laplace transform to our equation from the dynamic strang quartet module: x B x g , we write it as ℒ t x ℒ B x g and so must determine how ℒ acts on derivatives and sums. With respect to the latter it follows directly from the definition that ℒ B x g ℒ B x ℒ g B ℒ x ℒ g . Regarding its effect on the derivative we find, on integrating by parts, that ℒ t x t 0 s t t x t 0 x t s t s t 0 s t x t . Supposing that x and s are such that x t s t 0 as t we arrive at ℒ t x s ℒ x x 0 . Now, upon substituting and into we find s ℒ x x 0 B ℒ x ℒ g , which is easily recognized to be a linear system for ℒ x , namely s I B ℒ x ℒ g x 0 . The only thing that distinguishes this system from those encountered since our first brush with these systems is the presence of the complex variable s. This complicates the mechanical steps of Gaussian Elimination or the Gauss-Jordan Method but the methods indeed apply without change. Taking up the latter method, we write ℒ x s I B ℒ g x 0 . The matrix s I B is typically called the transfer function or resolvent, associated with B, at s. We turn to MATLAB for its symbolic calculation. (for more information, see the tutorial on MATLAB's symbolic toolbox). For example, >> B = [2 -1; -1 2] >> R = inv(s*eye(2)-B) R = [ (s-2)/(s*s-4*s+3), -1/(s*s-4*s+3)] [ -1/(s*s-4*s+3), (s-2)/(s*s-4*s+3)] We note that s I B is well defined except at the roots of the quadratic, s 2 4 s 3 . This quadratic is the determinant of s I B and is often referred to as the characteristic polynomial of B. Its roots are called the eigenvalues of B. As a second example let us take the B matrix of the dynamic Strang quartet module with the parameter choices specified in fib3.m, namely B -0.1350.1250 0.5-1.010.5 00.5-0.51 The associated s I B is a bit bulky (please run fib3.m) so we display here only the denominator of each term, i.e., s 3 1.655 s 2 0.4078 s 0.0039 . Assuming a current stimulus of the form i 0 t t 3 t 6 10000 and E m 0 brings ℒ g s 0.191 s 1 6 4 0 0 and so persists in ℒ x s I B ℒ g 0.191 s 1 6 4 s 3 1.655 s 2 0.4078 s 0.0039 s 2 1.5 s 0.27 0.5 s 0.26 0.2497 Now comes the rub. A simple linear solve (or inversion) has left us with the Laplace transform of x. The accursed No Free Lunch Theorem We shall have to do some work in order to recover x from ℒ x . confronts us. We shall face it down in the Inverse Laplace module.