3.3 Quizz 3

  3.3.1 Problem 1
  3.3.2 Problem 2
  3.3.3 Problem 3
  3.3.4 Problem 4

PDF (letter size)

PDF (legal size)

3.3.1 Problem 1

Determine whether or not the given matrix A is diagonalizable. If it is find a diagonalizing matrix P and a diagonal matrix D such that P1AP=DA=[3081] Solution

The first step is to find the eigenvalues. This is found by solvingdet(AλI)=0|3λ081λ|=0(3λ)(1λ)=0

Hence the eigenvalues are λ1=3,λ2=1. Since the eigenvalues are unique, then the matrix is diagonalizable. We need to determine the corresponding eigenvector in order to find P

λ1=3

Solving[3λ1081λ1][v1v2]=[00][330813][v1v2]=[00][0084][v1v2]=[00]

Hence v1=t is free variable and v2 is base variable. Second row gives  8t4v2=0 or v2=2t. Therefore the eigenvector is [v1v2]=[t2t]=t[12] Let t=1, then v1=[12] λ2=1

Solving[3λ2081λ2][v1v2]=[00][3+1081+1][v1v2]=[00][4080][v1v2]=[00]

R2=R22R1 gives[4000][v1v2]=[00] Hence v2=t is free variable and v1 is base variable. First row gives  4v1=0 or v1=0. Therefore the eigenvector is [v1v2]=[0t]=t[01] Let t=1, then v2=[01] Now that both eigenvectors are found, then P=[v1v2]=[1021]

And D is the diagonal matrix of the eigenvalues arranged in same order as the corresponding eigenvectors. (Will verify below). HenceD=[3001] Therefore P1AP=D(1)[1021]1[3081][1021]=[3001]

To verify the above, the LHS of (1) is evaluated directly, to confirm that D is indeed the result and it is diagonal of the eigenvalues. The first step is to find [1021]1. Since this is 2×2 then[1021]1=1det(P)[1021] But det(P)=1. The above becomes[1021]1=[1021] Therefore the LHS of (1) becomes[1021]1[3081]=[1021][3081]=[1×302×3+1×81]=[306+81]=[3021]

And now LHS of (1) becomes[3021][1021]=[302×11×21]=[3001]

HenceD=[3001] Which confirms D is the matrix whose diagonal elements are the eigenvalues of A.

3.3.2 Problem 2

Find the general solution of the homogeneous differential equation y+y10y=0

Solution

This is a linear 3rd order constant coefficient ODE. Hence the method of characteristic equation will be used. Let the solution be y=Aeλx. Substituting this into the ODE givesAλ3eλx+Aλeλx10Aeλx=0Aeλx(λ3+λ10)=0

Which simplifies (for non-trivial y) to the characteristic equation which is a polynomial in λ λ3+λ10=0 By inspection, we see that λ=2 is a root. Therefore a factor of the equation is (λ2). Now doing long division λ3+λ10(λ2) gives λ2+2λ+5.

pict
Figure 3.3:Polynmial long division to find remainder

Hence the above polynomial can be written as(1)(λ2)(λ2+2λ+5)=0 Now the roots for (λ2+2λ+5)=0 are found using the quadratic formula. λ=b2a±12ab24ac=22±1244×5=1±1216=1±12(4i)=1±2i

Hence the roots of the characteristic equation areλ1=2λ2=1+2iλ3=12i

Therefore the basis solution are {e2x,e(1+2i)x,e(12i)x} and the general solution is a linear combination of these basis solutions which givesy=Ae2x+Be(1+2i)x+Ce(12i)x Which can be simplified to (2)y=Ae2x+ex(Be2ix+Ce2ix) By using Euler formula, the above can be simplified further as followsBe2ix+Ce2ix=B(cos(2x)+isin(2x))+C(cos(2x)isin(2x))=(B+C)cos(2x)+sin(2x)(i(BC))

Let (B+C)=B0 a new constant and let i(BC)=C0 a new constant, the above becomesBe2ix+Ce2ix=B0cos(2x)+C0sin(2x) Substituting the above back in (2) gives the general solution asy=Ae2x+ex(B0cos(2x)+C0sin(2x)) The constants A,B0,C0 can be found from initial conditions if given.

3.3.3 Problem 3

Using the method of undermined coefficients, compute the general solution of the given equation y+3y+2y=2sin(x)

Solution

The solution is  y=yh+yp Where yh is the solution to the homogenous ode y+3y+2y=0 and yp is a particular solution to the given ODE. The ode y+3y+2y=0 is linear second order constant coefficient ODE. Hence the method of characteristic equation will be used. Let the solution be yh=Aeλx. Substituting this into y+3y+2y=0Aλ2eλx+Aλ3eλx+2Aeλx=0Aeλx(λ2+3λ+2)=0

And for non trivial solution the above simplifies toλ2+3λ+2=0(λ+1)(λ+2)=0

Hence the roots are λ=1,λ=2. Therefore the basis solutions for yh are {ex,e2x} and yh is linear combination of these basis. Therefore (1)yh=c1ex+c2e2x Now yp is found. Since the RHS is sin(x) then the trial solution is (2)yp=Acos(x)+Bsin(x) This shows that the basis for yp are {sinx,cosx}. There are no duplication between these basis and the basis for yh, so no need to multiply by an extra x. Using (2) gives(3)yp=Asin(x)+Bcos(x)(4)yp=Acos(x)Bsin(x)

Substituting (2,3,4) back into the given ODE givesyp+3yp+2yp=2sin(x)(Acos(x)Bsin(x))+3(Asin(x)+Bcos(x))+2(Acos(x)+Bsin(x))=2sin(x)cos(x)(A+3B+2A)+sin(x)(B3A+2B)=2sin(x)cos(x)(3B+A)+sin(x)(3A+B)=2sin(x)

Comparing coefficients on both sides gives two equations to solve for A,B 3B+A=03A+B=2

Multiplying the second equation by 3 gives 3B+A=09A3B=6

Adding the above two equations gives 10A=6. Hence A=35 and therefore 3B=35 or B=15. Substituting these values of A,B into (2) givesyp=35cos(x)+15sin(x) Hence the solution becomesy=yh+yp=(c1ex+c2e2x)+(35cos(x)+15sin(x))=c1ex+c2e2x35cos(x)+15sin(x)

3.3.4 Problem 4

Show that the given vector functions are linearly independent x1(t)=[et2et]x2(t)=[sintcost] Solution

These functions are defined for all t. Hence domain is t(,). The Wronskian of these vectors isW(t)=|etsint2etcost|=etcost2etsint=et(cost2sint)

We just need find one value t0 where W(t0)0 to show linearly independence. At t=0 the above becomesW(t=0)=1 Therefore the given vector functions are linearly independent. An alternative method is to write c1x1(t)+c2x2(t)=0c1[et2et]+c2[sintcost]=[00]

If the above is true only for c1=0,c2=0 then x1(t),x2(t) are linearly independent. The above can be written as[etsint2etcost][c1c2]=[00] R2=R22R1[etsint0cost2sint][c1c2]=[00] Row two gives c2(cost2sint)=0 For this to be true for any t in the interval t(,), then only solution is c2=0. First row now gives c1et=0 But et is never zero which means c1=0.

Since the only solution to c1x1+c2x2=0 is c1=c2=0, then this shows that x1,x2 are linearly independent.