University of Florida/Egm4313/s12.team11.imponenti/R5.5

From testwiki
Jump to navigation Jump to search

R5.6

solved by Luca Imponenti

Problem Statement

Complete the solution to the following problem

y+4y+13y=2e2xcos(3x)

where

yh=e2x[Acos(3x)+Bsin(3x)]

and

yp=xe2x[Mcos(3x)+Nsin(3x)]

Find the overall solution y(x) corresponds to the initial condition:

y(0)=1 , y(0)=0
Plot the solution over 3 periods.

Taking the derivatives of the particular solution yp(x)

Particular Solution

yp=xe2x[Mcos(3x)+Nsin(3x)]

y'p=e2x[sin(3x)(N2Nx3Mx)+cos(3x)(3Nx+M2Mx)]

y'p=e2x[sin(3x)(12Mx6M5Nx4N)+cos(3x)(6N5Mx4M12Nx)]

Plugging these into the ODE yields

sin(3x)(6M12Mx+5Nx+4N)+cos(3x)(5Mx+4M+12Nx6N)+ 4[sin(3x)(N2Nx3Mx)+cos(3x)(3Nx+M2Mx)]+ 13x[Mcos(3x)+Nsin(3x)]=2cos(3x)

Equating like terms allows us to solve for M and N

sin(3x)[(12Mx6M5Nx4N)+4(N2Nx3Mx)+13Nx]=0

cos(3x)[(6N5Mx4M12Nx)+4(3Nx+M2Mx)+13Mx]=2cos(3x)

6M=0

6N=2

M=0 , N=13

So the particular solution is

yp=13xe2xsin(3x)

Overall Solution

The overall solution in the sum of the homogeneous and particular solutions

y(x)=yh(x)+yp(x)

y(x)=e2x[Acos(3x)+Bsin(3x)]+13xe2xsin(3x)

To find A and B we apply the initial conditions

y(0)=1 , y(0)=0

y(0)=1=A

Taking the derivative

y(x)=ddx[e2x[cos(3x)+Bsin(3x)]+13xe2xsin(3x)]

y(x)=e2x[(3B+x2)cos(3x)(2B+23x+83)sin(3x)]

y(0)=0=3B2

B=23

Giving us the overall solution

   y(x)=e2x[cos(3x)+23sin(3x)+13xsin(3x)]

Plot

The period for cos(3x) , sin(3x) is 2π3

Plotting the solution y(x) over 3 periods yields

File:Plot56.jpg

Egm4313.s12.team11.imponenti (talk) 01:56, 30 March 2012 (UTC)

Template:CourseCat