JUST THE MATHS SLIDES NUMBER ORDINARY DIFFERENTIAL EQUATIONS 4 (Second order equations (A)) A.J.Hobson

Transcription

1 JUST THE MATHS SLIDES NUMBER 15.4 ORDINARY DIFFERENTIAL EQUATIONS 4 (Second order equations (A)) by A.J.Hobson Introduction Second order homogeneous equations Special cases of the auxiliary equation

2 UNIT ORDINARY DIFFERENTIAL EQUATIONS 4 SECOND ORDER EQUATIONS (A) INTRODUCTION A second order ordinary linear differential equation (with constant coefficients) has the general form a d2 y dx 2 + bdy dx + cy = f(x), where a, b and c are the constant coefficients The various cases of solution which arise depend on the values of the coefficients, together with the type of function, f(x), on the right hand side SECOND ORDER HOMOGENEOUS EQUATIONS The term homogeneous, in the present context means that f(x) 0. That is, we consider a d2 y dx 2 + bdy dx + cy = 0. 1

3 Note: A very simple case of this equation is d 2 y dx 2 = 0, which, on integration twice, gives the general solution y = Ax + B, where A and B are arbitrary constants. We should therefore expect two arbitrary constants in the solution of any second order linear differential equation with constant coefficients. THE STANDARD GENERAL SOLUTION The equivalent of a d2 y dx 2 + bdy dx + cy = 0 in the discussion of first order differential equations would have been That is, b dy dx + cy = 0. 2

4 dy dx + c b y = 0. This could have been solved using an integrating factor of giving the general solution e c b dx = e c b x, y = Ae c b x, where A is an arbitrary constant. We therefore make a trial solution of the form y = Ae mx, where A 0, in the second order case. We shall need dy dx = Amemx and d2 y dx 2 = Am2 e mx. Hence, aam 2 e mx + bame mx + cae mx = 0. 3

5 In other words, am 2 + bm + c = 0. This quadratic equation is called the auxiliary equation, having the same (constant) coefficients as the orginal differential equation. In general, the auxiliary equation will have two solutions, say m = m 1 and m = m 2, giving corresponding solutions, of the differential equation. y = Ae m 1x and y = Be m 2x, The linearity of the differential equation implies that the sum of any two solutions is also a solution. Thus, is another solution. y = Ae m 1x + Be m 2x Since this contains two arbitrary constants, we take it to be the general solution. 4

6 Notes: (i) It may be shown that there are no solutions other than those of the above form, though special cases are considered later. (ii) It will be possible to determine particular values of A and B if an appropriate number of boundary conditions for the differential equation are specified. (iii) Boundary conditions usually consist of a set of given values for y and dy dx at a certain value of x. EXAMPLE Determine the general solution of the differential equation d 2 y dx 2 + 5dy dx + 6y = 0 and also the particular solution for which y = 2 and = 5 when x = 0. dy dx 5

7 Solution The auxiliary equation is which can be factorised as m 2 + 5m + 6 = 0, (m + 2)(m + 3) = 0. Its solutions are therefore m = 2 amd m = 3. Hence, the differential equation has general solution y = Ae 2x + Be 3x, where A and B are arbitrary constants. Applying the boundary conditions, we shall also need dy dx = 2Ae 2x 3Be 3x. Hence, 2 = A + B and 5 = 2A 3B, giving A = 1, B = 1. The required particular solution is y = e 2x + e 3x. 6

8 SPECIAL CASES OF THE AUXILIARY EQUATION (a) The auxilary equation has coincident solutions Suppose that both solutions of the auxiliary equation are the same number, m 1. In other words, am 2 + bm + c is a perfect square and hence a(m m 1 ) 2 Apparently, the general solution of the differential equation is y = Ae m 1x + Be m 1x. But this does not genuinely contain two arbitrary constants since it can be rewritten as y = Ce m 1x where C = A + B. It will not, therefore, count as the general solution. The fault seems to lie with the constants A and B rather than with m 1. 7

9 We now examine a new trial solution of the form y = ze m 1x, where z denotes a function of x. We shall also need and dy dx = zm 1e m1x + e m 1x dz dx d 2 y dx = 2 zm2 1e m1x + 2m 1 e m 1x dz dx + em 1x d2 z dx 2. Substituting these into the differential equation, we obtain e m 1x or a zm 2 dz 1 + 2m 1 dx + d2 z dx 2 + b zm 1 + dz + cz = 0 dx z(am bm 1 + c) + dz dx (2am 1 + b) + a d2 z dx 2 = 0. The first term on the left hand side is zero since m 1 is already a solution of the auxiliary equation. 8

10 The second term on the left hand side is also zero since the auxiliary equation is equivalent to That is, Thus b = 2am 1. We conclude that with the result that a(m m 1 ) 2 = 0. am 2 2am 1 m + am 2 1 = 0. d 2 z dx 2 = 0, z = Ax + B, where A and B are arbitrary constants. Summary The general solution of the differential equation in the case of coincident solutions to the auxiliary equation y = (Ax + B)e m 1x. 9

11 EXAMPLE Determine the general solution of the differential equation Solution The auxilary equation is 4 d2 y dx 2 + 4dy dx + y = 0. 4m 2 + 4m + 1 = 0 or (2m + 1) 2 = 0. It has coincident solutions at m = 1 2. The general solution is therefore y = (Ax + B)e 1 2 x. 10

12 (b) The auxiliary equation has complex solutions If the auxiliary equation has complex solutions, they will appear as a pair of complex conjugates, say m = α ± jβ. Using m = α ± jβ instead of m = m 1 and m = m 2, the general solution of the differential equation will be y = P e (α+jβ)x + Qe (α jβ)x, where P and Q are arbitrary constants. By properties of complex numbers, a neater form of this result is obtainable as follows: y = e αx [P (cos βx + j sin βx) + Q(cos βx j sin βx)] or y = e αx [(P + Q) cos βx + j(p Q) sin βx]. Replacing P + Q and j(p Q) (which are just arbitrary quantities) by A and B, we obtain the standard general solution for the case in which the auxiliary equation has complex solutions. The general solution is 11

13 EXAMPLE y = e αx [A cos βx + B sin βx]. Determine the general solution of the differential equation Solution d 2 y dx 2 6dy dx The auxiliary equation is + 13y = 0. m 2 6m + 13 = 0, which has solutions given by m = ( 6) ± ( 6) = 6 ± j4 2 = 3 ± j2. The general solution is therefore y = e 3x [A cos 2x + B sin 2x]. 12