|
|
| (One intermediate revision by the same user not shown) |
| Line 1: |
Line 1: |
| {{Differential equations}}
| | When you are having [http://www.Google.de/search?q=trouble trouble] seeing a while you are grinding it, try adjusting currently the brightness environment. Might make the display look clear, enhancing your is way better expertise. And let's face it, you will not achieve any kind concerning success if you not able to see what you're doing, so make the task meet your needs.<br><br>To help conclude, clash of clans hack tool no survey must not be let to get in approach of the bigger question: what makes we here? If you have any queries concerning exactly where and how to use [http://circuspartypanama.com clash of clans unlimited troops], you can contact us at the web-site. Putting this particular apart from its of great great. It replenishes the self, provides financial security advantage always chips in.<br><br>Okazaki, japan tartan draws inspiration through your country's involvement in cherry blossom and should include pink, white, green as well as , brown lightly colours. clash of clans cheats. Determined by is called Sakura, okazaki, japan for cherry blossom.<br><br>Workstation games offer entertaining to allow them to everybody, and they remain surely more complicated for you to Frogger was! To get all you may possibly out of game titles, use the advice put out here. You are going to find powerful exciting new world operating in gaming, and you will probably wonder how you previously got by without items!<br><br>Examination your child's xbox golf game enjoying. Video computer game are now rated just exactly like films and generally can help. This specific enables you to keep an eye on unquestionably the information your kids often is exposed to. Dependant upon your child's age, continue to keep him clear of video movie that happen to nevertheless be meant for people who are more fully learned than him.<br><br>So that you can access it into excel, copy-paste this continued method into corpuscle B1. If you again get access to an majority of energy in abnormal in corpuscle A1, the bulk with treasures will arise when it comes to B1.<br><br>As well as our options are screened and approved from the best virus recognition software and so anti-virus in the sell to ensure a security-level as huge as you can, in argument you fear for the security of your computer or maybe your cellular device, no inconveniences. In case you nevertheless have any sort of doubts, take a look at the movie and you'll focus on it operates and it's very 100% secure! It'll only take a few moments of the time! |
| | |
| In mathematics, '''linear differential equations''' are [[differential equations]] having [[homogeneous differential equation|homogeneous]] solutions which can be added together to form other homogeneous solutions. They can be [[ordinary differential equation|ordinary]] or [[partial differential equation|partial]]. The homogeneous solutions to linear equations form a [[vector space]] (unlike [[non-linear differential equation]]s).
| |
| | |
| ==Introduction==
| |
| Linear '''[[differential equation]]s''' are of the form
| |
| | |
| : <math> Ly = f</math>
| |
| | |
| where the [[differential operator]] ''L'' is a [[linear operator]], ''y'' is the unknown function (such as a function of time ''y''(''t'')), and the [[right hand side]] ''f'' is a given function of the same nature as ''y'' (called the '''source term'''). For a function dependent on time we may write the equation more expressly as
| |
| | |
| : <math> L y(t) = f(t)</math>
| |
| | |
| and, even more precisely by bracketing
| |
| | |
| : <math> L [y(t)] = f(t)</math>
| |
| | |
| The linear operator ''L'' may be considered to be of the form<ref>Gershenfeld 1999, p.9</ref>
| |
| | |
| : <math>L_n(y) \equiv \frac{d^n y}{dt^n} + A_1(t)\frac{d^{n-1}y}{dt^{n-1}} + \cdots + A_{n-1}(t)\frac{dy}{dt} + A_n(t)y </math>
| |
| | |
| The linearity condition on ''L'' rules out operations such as taking the square of the [[derivative]] of ''y''; but permits, for example, taking the second derivative of ''y''.
| |
| It is convenient to rewrite this equation in an operator form
| |
| | |
| : <math> L_n(y) \equiv \left[\,D^n + A_{1}(t)D^{n-1} + \cdots + A_{n-1}(t) D + A_n(t)\right] y</math>
| |
| | |
| where ''D'' is the differential operator ''d/dt'' (i.e. ''Dy = y' '', ''D''<sup>2</sup>''y = y",... ''), and the ''A<sub>n</sub>'' are given functions.
| |
| <!-- and the source term is considered to be a function of time ƒ(''t'').-->
| |
| | |
| Such an equation is said to have '''order''' ''n'', the index of the highest derivative of ''y'' that is involved. <!-- (Assuming a possibly existing coefficient ''a<sub>n</sub>'' of this derivative to be non zero, it is eliminated by dividing through it. In case it can become zero, different cases must be considered separately for the analysis of the equation.)-->
| |
| | |
| A typical simple example is the linear differential equation used to model radioactive decay.<ref>Robinson 2004, p.5</ref> Let ''N''(''t'') denote the number of radioactive atoms in some sample of material <ref>Robinson 2004, p.7</ref> at time ''t''. Then for some constant ''k'' > 0, the number of radioactive atoms which decay can be modelled by
| |
| | |
| :<math> \frac{dN}{dt} = -k N</math> | |
| | |
| If ''y'' is assumed to be a function of only one variable, one speaks about an [[ordinary differential equation]], else the derivatives and their coefficients must be understood as ([[tensor contraction|contracted]]) vectors, matrices or [[tensor]]s of higher rank, and we have a (linear) [[partial differential equation]]. | |
| | |
| The case where ''f'' = 0 is called a '''homogeneous equation''' and its solutions are called '''complementary functions'''. It is particularly important to the solution of the general case, since any complementary function can be added to a solution of the inhomogeneous equation to give another solution (by a method traditionally called ''particular integral and complementary function''). When the ''A<sub>i</sub>'' are numbers, the equation is said to have ''[[constant coefficients]]''.
| |
| | |
| ==Homogeneous equations with constant coefficients==
| |
| {{main|Characteristic equation (calculus)}}
| |
| The first method of solving linear homogeneous ordinary differential equations with constant coefficients is due to [[Euler]], who realized that solutions have the form ''e<sup>zx</sup>'', for possibly-complex values of ''z''. The exponential function is one of the few functions to keep its shape after differentiation, allowing the sum of its multiple derivatives to cancel out to zero, as required by the equation. Thus, for constant values ''A''<sub>1</sub>,..., ''A''<sub>n</sub>, to solve:
| |
| | |
| :<math>y^{(n)} + A_{1}y^{(n-1)} + \cdots + A_{n}y = 0\,,</math>
| |
| | |
| we set ''y'' = ''e<sup>zx</sup>'', leading to
| |
| | |
| :<math>z^n e^{zx} + A_1 z^{n-1} e^{zx} + \cdots + A_n e^{zx} = 0.</math>
| |
| | |
| Division by ''e<sup>zx</sup>'' gives the ''n''th-order polynomial:
| |
| | |
| :<math>F(z) = z^{n} + A_{1}z^{n-1} + \cdots + A_n = 0.\,</math>
| |
| | |
| This algebraic equation ''F''(''z'') = 0 is the [[Characteristic equation (calculus)|characteristic equation]] considered later by [[Gaspard Monge]] and [[Augustin-Louis Cauchy]].
| |
| | |
| Formally, the terms
| |
| | |
| :<math>y^{(k)}\quad\quad(k = 1, 2, \dots, n).</math>
| |
| | |
| of the original differential equation are replaced by ''z<sup>k</sup>''. [[Root-finding algorithm#Finding roots of polynomials|Solving]] the polynomial gives ''n'' values of ''z'', ''z''<sub>1</sub>, ..., ''z<sub>n</sub>''. Substitution of any of those values for ''z'' into ''e<sup>zx</sup>'' gives a solution ''e<sup>z<sub>i</sub>x</sup>''. Since homogeneous linear differential equations obey the [[superposition principle]], any [[linear combination]] of these functions also satisfies the differential equation. | |
| | |
| When these roots are all [[distinct roots|distinct]], we have ''n'' distinct solutions to the differential equation. It can be shown that these are [[linearly independent]], by applying the [[Vandermonde determinant]], and together they form a [[Basis (linear algebra)|basis]] of the space of all solutions of the differential equation.
| |
| | |
| {{ExampleSidebar|35%|
| |
| :<math>y''''-2y'''+2y''-2y'+y=0</math>
| |
| has the characteristic equation
| |
| | |
| : <math>z^4-2z^3+2z^2-2z+1=0.</math>
| |
| | |
| This has zeroes, ''i'', −''i'', and 1 (multiplicity 2). The solution basis is then
| |
| | |
| : <math>e^{ix} ,\, e^{-ix} ,\, e^x ,\, xe^x.</math>
| |
| | |
| This corresponds to the real-valued solution basis
| |
| | |
| : <math>\cos x ,\, \sin x ,\, e^x ,\, xe^x \,.</math>}}
| |
| | |
| The preceding gave a solution for the case when all zeros are distinct, that is, each has [[Multiplicity (mathematics)|multiplicity]] 1. For the general case, if ''z'' is a (possibly complex) [[root of a function|zero]] (or root) of ''F''(''z'') having multiplicity ''m'', then, for <math>k\in\{0,1,\dots,m-1\} \,</math>, <math>y=x^ke^{zx}</math> is a solution of the ODE. Applying this to all roots gives a collection of ''n'' distinct and linearly independent functions, where ''n'' is the degree of ''F''(''z''). As before, these functions make up a basis of the solution space.
| |
| | |
| If the coefficients ''A<sub>i</sub>'' of the differential equation are real, then real-valued solutions are generally preferable. Since non-real roots ''z'' then come in [[complex conjugate|conjugate]] pairs, so do their corresponding basis functions {{nowrap|''x''<sup>''k''</sup>e<sup>''zx''</sup>}}, and the desired result is obtained by replacing each pair with their real-valued [[linear combination]]s [[real part|Re(''y'')]] and [[Imaginary part|Im(''y'')]], where ''y'' is one of the pair.
| |
| | |
| A case that involves complex roots can be solved with the aid of [[Euler's formula]].
| |
| | |
| ===Examples===
| |
| Given <math>y''-4y'+5y=0</math>. The characteristic equation is <math>z^2-4z+5=0</math> which has roots 2±''i''. Thus the solution basis <math>\{y_1,y_2\}</math> is <math>\{e^{(2+i)x},e^{(2-i)x}\}</math>. Now ''y'' is a solution if and only if <math>y=c_1y_1+c_2y_2</math> for <math>c_1,c_2\in\mathbf{C}</math>.
| |
| | |
| Because the coefficients are real,
| |
| *we are likely not interested in the complex solutions
| |
| *our basis elements are mutual conjugates
| |
| The linear combinations
| |
| | |
| :<math>u_1=\mbox{Re}(y_1)=\tfrac{1}{2} (y_1+y_2) =e^{2x}\cos(x),</math>
| |
| :<math>u_2=\mbox{Im}(y_1)=\tfrac{1}{2i} (y_1-y_2) =e^{2x}\sin(x),</math>
| |
| | |
| will give us a real basis in <math>\{u_1,u_2\}</math>.
| |
| | |
| ==== Simple harmonic oscillator ====
| |
| The second order differential equation
| |
| | |
| :<math> D^2 y = -k^2 y, </math>
| |
| | |
| which represents a simple [[harmonic oscillator]], can be restated as
| |
| | |
| :<math> (D^2 + k^2) y = 0. </math>
| |
| | |
| The expression in parenthesis can be factored out, yielding
| |
| | |
| :<math> (D + i k) (D - i k) y = 0,</math>
| |
| | |
| which has a pair of linearly independent solutions:
| |
| | |
| :<math> (D - i k) y = 0 </math>
| |
| :<math> (D + i k) y = 0. </math>
| |
| | |
| The solutions are, respectively,
| |
| | |
| :<math> y_0 = A_0 e^{i k x} </math>
| |
| | |
| and | |
| | |
| :<math> y_1 = A_1 e^{-i k x}. </math>
| |
| | |
| These solutions provide a basis for the two-dimensional [[vector space|solution space]] of the second order differential equation: meaning that linear combinations of these solutions will also be solutions. In particular, the following solutions can be constructed
| |
| | |
| :<math> y_{0'} = {C_0 e^{i k x} + C_0 e^{-i k x} \over 2} = C_0 \cos (k x) </math>
| |
| | |
| and
| |
| | |
| :<math> y_{1'} = {C_1 e^{i k x} - C_1 e^{-i k x} \over 2 i} = C_1 \sin (k x). </math>
| |
| | |
| These last two trigonometric solutions are linearly independent, so they can serve as another basis for the solution space, yielding the following general solution:
| |
| | |
| :<math> y_H = C_0 \cos (k x) + C_1 \sin (k x). </math>
| |
| | |
| ==== Damped harmonic oscillator ====
| |
| Given the equation for the damped [[harmonic oscillator]]:
| |
| | |
| :<math> \left(D^2 + \frac{b}{m} D + \omega_0^2\right) y = 0, </math>
| |
| | |
| the expression in parentheses can be factored out: first obtain the characteristic equation by replacing ''D'' with λ. This equation must be satisfied for all ''y'', thus:
| |
| | |
| :<math> \lambda^2 + \frac{b}{m} \lambda + \omega_0^2 = 0. </math>
| |
| | |
| Solve using the [[quadratic formula]]:
| |
| | |
| :<math> \lambda = \tfrac{1}{2} \left (-\frac{b}{m} \pm \sqrt{\frac{b^2}{m^2} - 4 \omega_0^2} \right ). </math>
| |
| | |
| Use these data to factor out the original differential equation:
| |
| | |
| :<math>\left(D + \frac{b}{2m} - \sqrt{\frac{b^2}{4 m^2} - \omega_0^2} \right) \left(D + \frac{b}{2m} + \sqrt{\frac{b^2}{4 m^2} - \omega_0^2}\right) y = 0. </math>
| |
| | |
| This implies a pair of solutions, one corresponding to
| |
| | |
| :<math> \left(D + \frac{b}{2m} - \sqrt{\frac{b^2}{4 m^2} - \omega_0^2} \right) y = 0 </math>
| |
| :<math> \left(D + \frac{b}{2m} + \sqrt{\frac{b^2}{4 m^2} - \omega_0^2}\right) y = 0 </math>
| |
| | |
| The solutions are, respectively,
| |
| | |
| :<math> y_0 = A_0 e^{-\omega x + \sqrt{\omega^2 - \omega_0^2} x} = A_0 e^{-\omega x} e^{\sqrt{\omega^2 - \omega_0^2} x} </math>
| |
| :<math> y_1 = A_1 e^{-\omega x - \sqrt{\omega^2 - \omega_0^2} x} = A_1 e^{-\omega x} e^{-\sqrt{\omega^2 - \omega_0^2} x} </math>
| |
| | |
| where ω = ''b''/2''m''. From this linearly independent pair of solutions can be constructed another linearly independent pair which thus serve as a basis for the two-dimensional solution space:
| |
| | |
| :<math> y_H (A_0, A_1) (x) = \left(A_0 \sinh \left (\sqrt{\omega^2 - \omega_0^2} x \right ) + A_1 \cosh \left ( \sqrt{\omega^2 - \omega_0^2} x \right ) \right) e^{-\omega x}. </math>
| |
| | |
| However, if |ω| < |ω<sub>0</sub>| then it is preferable to get rid of the consequential imaginaries, expressing the general solution as
| |
| | |
| :<math> y_H (A_0, A_1) (x) = \left(A_0 \sin \left(\sqrt{\omega_0^2 - \omega^2} x \right ) + A_1 \cos \left (\sqrt{\omega_0^2 - \omega^2} x \right ) \right) e^{-\omega x}. </math>
| |
| | |
| This latter solution corresponds to the underdamped case, whereas the former one corresponds to the overdamped case: the solutions for the underdamped case [[oscillation|oscillate]] whereas the solutions for the overdamped case do not.
| |
| | |
| ==Nonhomogeneous equation with constant coefficients==
| |
| To obtain the solution to the '''nonhomogeneous equation''' (sometimes called '''inhomogeneous equation'''), find a particular integral ''y<sub>P</sub>''(''x'') by either the [[method of undetermined coefficients]] or the [[method of variation of parameters]]; the general solution to the linear differential equation is the sum of the general solution of the related homogeneous equation and the particular integral. Or, when the initial conditions are set, use [[Laplace transform]] to obtain the particular solution directly.
| |
| | |
| Suppose we face
| |
| | |
| :<math>\frac {d^{n}y(x)} {dx^{n}} + A_{1}\frac {d^{n-1}y(x)} {dx^{n-1}} + \cdots + A_{n}y(x) = f(x).</math>
| |
| | |
| For later convenience, define the characteristic polynomial
| |
| | |
| :<math>P(v)=v^n+A_1v^{n-1}+\cdots+A_n.</math>
| |
| | |
| We find a solution basis <math>\{y_1(x),y_2(x),\ldots,y_n(x)\}</math> for the homogeneous (''f''(''x'') = 0) case. We now seek a '''particular integral''' ''y<sub>p</sub>''(''x'') by the '''variation of parameters''' method. Let the coefficients of the linear combination be functions of ''x'':
| |
| | |
| :<math>y_p(x) = u_1(x) y_1(x) + u_2(x) y_2(x) + \cdots + u_n(x) y_n(x).</math>
| |
| | |
| For ease of notation we will drop the dependency on ''x'' (i.e. the various ''(x)''). Using the operator notation ''D'' = ''d/dx'', the ODE in question is ''P''(''D'')''y'' = ''f''; so
| |
| | |
| :<math>f=P(D)y_p=P(D)(u_1y_1)+P(D)(u_2y_2)+\cdots+P(D)(u_ny_n).</math>
| |
| | |
| With the constraints
| |
| | |
| :<math>0=u'_1y_1+u'_2y_2+\cdots+u'_ny_n</math>
| |
| :<math>0=u'_1y'_1+u'_2y'_2+\cdots+u'_ny'_n</math>
| |
| :<math> \cdots</math>
| |
| :<math>0=u'_1y^{(n-2)}_1+u'_2y^{(n-2)}_2+\cdots+u'_ny^{(n-2)}_n</math>
| |
| | |
| the parameters commute out,
| |
| | |
| :<math>f=u_1P(D)y_1+u_2P(D)y_2+\cdots+u_nP(D)y_n+u'_1y^{(n-1)}_1+u'_2y^{(n-1)}_2+\cdots+u'_ny^{(n-1)}_n.</math>
| |
| | |
| But ''P''(''D'')''y<sub>j</sub>'' = 0, therefore
| |
| | |
| :<math>f=u'_1y^{(n-1)}_1+u'_2y^{(n-1)}_2+\cdots+u'_ny^{(n-1)}_n.</math>
| |
| | |
| This, with the constraints, gives a linear system in the ''u′<sub>j</sub>''. This much can always be solved; in fact, combining [[Cramer's rule]] with the [[Wronskian]],
| |
| | |
| :<math>u'_j=(-1)^{n+j}\frac{W(y_1,\ldots,y_{j-1},y_{j+1}\ldots,y_n)_{0 \choose f}}{W(y_1,y_2,\ldots,y_n)}.</math> <!-- caution: check my sign -->
| |
| | |
| In the very non-standard notation used above, one should take the i,n-minor of W and multiply it by f. That's why we get a minus-sign. Alternatively, forget about the minus sign and just compute the determinant of the matrix obtained by substituting the j-th W column with (0, 0, ..., f).
| |
| | |
| The rest is a matter of integrating ''u′<sub>j</sub>''.
| |
| | |
| The particular integral is not unique; <math>y_p+c_1y_1+\cdots+c_ny_n</math> also satisfies the ODE for any set of constants ''c<sub>j</sub>''.
| |
| | |
| ===Example===
| |
| Suppose <math>y''-4y'+5y=\sin(kx)</math>. We take the solution basis found above <math>\{e^{(2+i)x}=y_1(x),e^{(2-i)x}=y_2(x)\}</math>.
| |
| | |
| :<math>\begin{align}
| |
| W &= \begin{vmatrix}e^{(2+i)x}&e^{(2-i)x} \\ (2+i)e^{(2+i)x}&(2-i)e^{(2-i)x} \end{vmatrix} = e^{4x}\begin{vmatrix}1&1\\ 2+i&2-i\end{vmatrix} =-2ie^{4x}\\
| |
| u'_1 &=\frac{1}{W}\begin{vmatrix}0&e^{(2-i)x}\\ \sin(kx)&(2-i)e^{(2-i)x}\end{vmatrix} = -\tfrac{i}{2} \sin(kx)e^{(-2-i)x}\\
| |
| u'_2 &=\frac{1}{W}\begin{vmatrix}e^{(2+i)x}&0\\ (2+i)e^{(2+i)x}&\sin(kx)\end{vmatrix} =\tfrac{i}{2} \sin(kx)e^{(-2+i)x}.
| |
| \end{align}</math>
| |
| | |
| Using the [[list of integrals of exponential functions]]
| |
| | |
| :<math>u_1=-\tfrac{i}{2}\int\sin(kx)e^{(-2-i)x}\,dx =\frac{ie^{(-2-i)x}}{2(3+4i+k^2)}\left((2+i)\sin(kx)+k\cos(kx)\right)</math>
| |
| :<math>u_2=\tfrac{i}{2}\int\sin(kx)e^{(-2+i)x}\,dx=\frac{ie^{(i-2)x}}{2(3-4i+k^2)}\left((i-2)\sin(kx)-k\cos(kx)\right).</math>
| |
| | |
| And so
| |
| :<math>\begin{align}
| |
| y_p &= u_1(x) y_1(x) + u_2(x) y_2(x) = \frac{i}{2(3+4i+k^2)}\left((2+i)\sin(kx)+k\cos(kx)\right) +\frac{i}{2(3-4i+k^2)}\left((i-2)\sin(kx)-k\cos(kx)\right) \\
| |
| &=\frac{(5-k^2)\sin(kx)+4k\cos(kx)}{(3+k^2)^2+16}.
| |
| \end{align}</math>
| |
| | |
| (Notice that ''u''<sub>1</sub> and ''u''<sub>2</sub> had factors that canceled ''y''<sub>1</sub> and ''y''<sub>2</sub>; that is typical.)
| |
| | |
| For interest's sake, this ODE has a physical interpretation as a driven damped [[harmonic oscillator]]; ''y<sub>p</sub>'' represents the steady state, and <math>c_1y_1+c_2y_2</math> is the transient.
| |
| | |
| ==Equation with variable coefficients==
| |
| A linear ODE of order ''n'' with variable coefficients has the general form
| |
| :<math>p_{n}(x)y^{(n)}(x) + p_{n-1}(x) y^{(n-1)}(x) + \cdots + p_0(x) y(x) = r(x).</math>
| |
| | |
| ===Examples===
| |
| A simple example is the [[Cauchy–Euler equation]] often used in engineering
| |
| | |
| :<math>x^n y^{(n)}(x) + a_{n-1} x^{n-1} y^{(n-1)}(x) + \cdots + a_0 y(x) = 0.</math>
| |
| | |
| == First order equation ==
| |
| {{ExampleSidebar|35%|2=Solve the equation
| |
| | |
| : <math>y'(x)+3y(x)=2</math>
| |
| | |
| with the initial condition
| |
| | |
| : <math>y(0)=2.</math>
| |
| | |
| Using the general solution method:
| |
| | |
| : <math>y=e^{-3x}\left(\int 2 e^{3x}\, dx + \kappa\right). \,</math>
| |
| | |
| The indefinite integral is solved to give:
| |
| | |
| : <math>y=e^{-3x}\left(2/3 e^{3x} + \kappa\right). \,</math>
| |
| | |
| Then we can reduce to:
| |
| | |
| : <math>y=2/3 + \kappa e^{-3x}. \,</math>
| |
| | |
| where κ = 4/3 from the initial condition.}}
| |
| A linear ODE of order 1 with variable coefficients has the general form
| |
| | |
| :<math>Dy(x) + f(x) y(x) = g(x).</math>
| |
| | |
| Where D is the [[differential operator]]. Equations of this form can be solved by multiplying the [[integrating factor]]
| |
| | |
| :<math>e^{\int f(x)\,dx}</math>
| |
| | |
| throughout to obtain
| |
| | |
| :<math> Dy(x)e^{\int f(x)\,dx}+f(x)y(x)e^{\int f(x)\,dx}=g(x)e^{\int f(x) \, dx},</math>
| |
| | |
| which simplifies due to the [[product rule]] to
| |
| | |
| : <math> D\left (y(x)e^{\int f(x)\,dx} \right )=g(x)e^{\int f(x)\,dx}</math>
| |
| | |
| which, on integrating both sides and solving for ''y''(''x'') gives:
| |
| | |
| : <math> y(x) = \frac{\int g(x)e^{\int f(x)\,dx} \,dx+c}{e^{\int f(x)\,dx}}.</math>
| |
| | |
| In other words: The solution of a first-order linear ODE | |
| | |
| : <math>y'(x) + f(x) y(x) = g(x),</math>
| |
| | |
| with coefficients that may or may not vary with ''x'', is:
| |
| | |
| :<math>y=e^{-a(x)}\left(\int g(x) e^{a(x)}\, dx + \kappa\right)</math>
| |
| | |
| where κ is the constant of integration, and
| |
| | |
| : <math>a(x)=\int{f(x)\,dx}.</math>
| |
| | |
| A compact form of the general solution is (see J. Math. Chem. 48 (2010) 175):
| |
| | |
| : <math> y(x) = \int_a^x \! {[y(a) \delta(t-a)+g(t)] e^{-\int_t^x \!f(u)du}\, dt}\,.</math>
| |
| | |
| where δ(''x'') is the generalized Dirac delta function.
| |
| | |
| === Examples ===
| |
| Consider a first order differential equation with [[constant coefficients]]:
| |
| | |
| :<math>\frac{dy}{dx} + b y = 1.</math>
| |
| | |
| This equation is particularly relevant to first order systems such as [[RC circuit]]s and [[damping|mass-damper]] systems.
| |
| | |
| In this case, ''f''(''x'') = ''b'', ''g''(''x'') = 1.
| |
| | |
| Hence its solution is
| |
| | |
| :<math>y(x) = e^{-bx} \left( \frac{e^{bx}}{b}+ C \right) = \frac{1}{b} + C e^{-bx} .</math>
| |
| | |
| ==Systems of Linear Differential Equations==
| |
| An arbitrary linear ordinary differential equation or even a system of such equations can be converted into a first order system of linear differential equations by adding variables for all but the highest order derivatives. A linear system can be viewed as a single equation with a vector-valued variable. The general treatment is analogous to the treatment above of ordinary first order linear differential equations, but with complications stemming from noncommutativity of matrix multiplication.
| |
| | |
| To solve
| |
| :<math>\left\{\begin{array}{rl}\mathbf{y}'(x) &= A(x)\mathbf{y}(x)+\mathbf{b}(x)\\\mathbf y(x_0)&=\mathbf y_0\end{array}\right.</math>
| |
| (here <math>\mathbf{y} (x)</math> is a vector or matrix, and <math>A( x )</math> is a matrix),
| |
| let <math>U( x )</math> be the solution of <math>\mathbf y'(x) = A(x)\mathbf y(x)</math> with <math>U(x_0) = I</math> (the identity matrix). <math>U</math> is a fundamental matrix for the equation — the columns of <math>U</math> form a complete linearly independent set of solutions for the homogeneous equation. After substituting <math>\mathbf y(x) = U(x)\mathbf z(x)</math>, the equation <math>\mathbf y'(x) = A(x)\mathbf y(x)+\mathbf b(x)</math> simplifies to <math>U(x)\mathbf z'(x) = \mathbf b(x).</math> Thus,
| |
| :<math>\mathbf{y}(x) = U(x)\mathbf{y_0} + U(x)\int_{x_0}^x U^{-1}(t)\mathbf{b}(t)\,dt</math>
| |
| If <math>A(x_1)</math> commutes with <math>A(x_2)</math> for all <math> x_1 </math> and <math> x_2</math>, then
| |
| :<math>U(x) = e^{\int_{x_0}^x A(x)\,dx}</math>
| |
| and thus
| |
| :<math>U^{-1}(x) = e^{-\int_{x_0}^x A(x)\,dx},</math>
| |
| but in the general case there is no closed form solution, and an approximation method such as [[Magnus expansion]] may have to be used. Note that the exponentials are [[matrix exponential]]s.
| |
| | |
| ==See also==
| |
| * [[Matrix differential equation]]
| |
| * [[Partial differential equation]]
| |
| * [[Continuous-repayment_mortgage#Ordinary_time_differential_equation|Continuous-repayment mortgage]]
| |
| * [[Fourier transform]]
| |
| * [[Laplace transform]]
| |
| * [[List of differentiation identities]], Nth Derivatives Section
| |
| | |
| ==External links==
| |
| * http://eqworld.ipmnet.ru/en/solutions/ode.htm
| |
| | |
| == Notes ==
| |
| {{reflist|2}}
| |
| | |
| == References ==
| |
| *{{Citation
| |
| | author = Birkhoff, Garrett and Rota, Gian-Carlo
| |
| | year = 1978
| |
| | title = Ordinary Differential Equations
| |
| | isbn = 0-471-07411-X
| |
| | publisher = John Wiley and Sons, Inc.
| |
| | location = New York
| |
| | oclc =
| |
| }}
| |
| *{{Citation
| |
| | author = Gershenfeld, Neil
| |
| | year = 1999
| |
| | title =The Nature of Mathematical Modeling
| |
| | isbn = 978-0-521-57095-4
| |
| | publisher = Cambridge University Press
| |
| | location = Cambridge, UK.
| |
| | oclc =
| |
| }}
| |
| *{{Citation
| |
| | author = Robinson, James C.
| |
| | year = 2004
| |
| | title = An Introduction to Ordinary Differential Equations
| |
| | isbn = 0-521-82650-0
| |
| | publisher = Cambridge University Press
| |
| | location = Cambridge, UK.
| |
| | oclc =
| |
| }}
| |
| | |
| {{DEFAULTSORT:Linear Differential Equation}}
| |
| [[Category:Differential equations]]
| |