University of Florida/Egm4313/s12.team11.R2

From testwiki
Jump to navigation Jump to search

Problem R2.1


Part 1


Problem Statement


Given the two roots and the initial conditions:

λ1=2,λ2=5
y(0)=1,y(0)=0

Find the non-homogeneous L2-ODE-CC in standard form and the solution in terms of the initial conditions and the general excitation r(x).
Consider no excitation:
r(x)=0
Plot the solution

Solution


Characteristic Equation:


(λλ1)(λλ2)=0
(λ+2)(λ5)=λ2+2λ5λ10=0

 λ23λ10=0


Non-Homogeneous L2-ODE-CC


 y3y10=r(x)
Homogeneous Solution:


yh(x)=c1e2x+c2e5x
y(x)=c1e2x+c2e5x+yp(x)
Since there is no excitation,
yp(x)=0

 y(x)=c1e2x+c2e5x
Substituting the given initial conditions:


y(0)=1

 1=c1+c2

y(0)=0

 0=2c1+5c2

Solving these two equations for c1 and c2 yields:

 c1=5/4,c2=1/4
Final Solution


 y(x)=(5/4)e2x(1/4)e5x

File:2.1fig1.jpg

Part 2


Problem Statement


Generate 3 non-standard (and non-homogeneous) L2-ODE-CC that admit the 2 values in (3a) p.3-7 as the 2 roots of the corresponding characteristic equation.

Solutions


2(λ+2)(λ5)=2λ2+4λ10λ20=0

 2λ26λ20=0


3(λ+2)(λ5)=3λ2+6λ15λ30=0

 3λ29λ30=0


4(λ+2)(λ5)=4λ2+8λ20λ40=0

 4λ212λ40=0


--Egm4313.s12.team11.gooding 02:01, 7 February 2012 (UTC)

Report 2, Problem 2

Problem Statement

Find and plot the solution for the homogeneous L2-ODE-CC

y(x)10y(x)+25y(x)=0

with initial conditions y(0)=1 ,and y(0)=0

Characteristic Equation

λ210λ+25=0

(λ5)(λ5)=0

λ=5

Homogeneous Solution

The solution to a L2-ODE-CC with real double root is given by

y(x)=c1eλx+c2xeλx

First initial condition

y(0)=1

y(0)=c1e5*0+c2*0*e5*0=1

c1=1

Second initial condition

y(0)=0

ddxy(x)=y(x)=5e5x+c2e5x(5x+1)

y(0)=5e5*0+c2e5*0(5*0+1)=0

5+c2=0

c2=5

The solution to our L2-ODE-CC is

                       y(x)=e5x(15x)

Plot

y(x)=e5x(15x)

File:PlotR2 2.jpg

Egm4313.s12.team11.imponenti 00:30, 8 February 2012 (UTC)

Problem R2.3

Problem Statement (K 2011 p.59 pb. 3)

Find a general solution. Check your answer by substitution.

Given

y+6y+8.96y=0

Solution

We can write the above differential equation in the following form:

d2ydx2+6dydx+8.96y=0

Let ddx=λ.

The characteristic equation of the given DE is

λ2+6λ+8.96=0

Now, in order to solve for λ, we can use the quadratic formula:

λ=(6)±(6)24(1)(8.96)2(1)

λ=6±3635.842

λ=6±0.162

λ=6±0.42

Therefore, we have:

λ1=6+0.42=5.82=2.8

and

λ1=60.42=6.42=3.2

Thus, we have found that the general solution of the DE is actually:

y=c1e2.8x+c2e3.2x

Check:

To check if y is indeed the solution of the given DE, we can differentiate the

what we found to be the general solution.

dydx=y=2.8xc1e2.8x+3.2xc2e3.2x

d2ydx2=y=7.84xc1e2.8x+10.24c2e3.2x

Substituting the values of y,y and y in the given equation, we get:

7.84c1e2.8x+10.24c2e3.2x+6(2.8c1e2.8x3.2c2e3.2x)+8.96(c1e2.8x+c2e3.2x)=0

7.84c1e2.8x+10.24c2e3.2x16.8c1e2.8x19.2c2e3.2x+8.96c1e2.8x+8.96c2e3.2x=0

and thus:

00.

Therefore, the solution of the given DE is in fact:


y=c1e2.8x+c2e3.2x

Problem Statement (K 2011 p.59 pb. 4)

Find a general solution to the given ODE. Check your answer by substituting into the original equation.

Given

y'+4y+(π2+4)y=0

Solution

The characteristic equation of this ODE is therefore:

λ2+aλ+b=λ2+4λ+(π2+4)=0

Evaluating the discriminant:


a24b=424(π2+4)=4π2<0


Therefore the equation has two complex conjugate roots and a general homogenous solution of the form:



y=eax/2(Acos(ωx)+Bsin(ωx))


Where: ω=b14a2=π2+414(4)2=π2=π

And finally we find the general homogenous solution:


                                                  y=e2x(Acos(πx)+Bsin(πx))


Check:

We found that:


y=e2x(Acos(πx)+Bsin(πx))

Differentiating y to obtain y' and y' respectively:

y'=2e2x(Acos(πx)+Bsin(πx))+e2x(πAcos(πx)πBsin(πx))

                                   y'=2e2x(Acos(πx)+Bsin(πx))+πe2x(Acos(πx)Bsin(πx))


y'=4e2x(Acos(πx)+Bsin(πx))2πe2x(Acos(πx)Bsin(πx))2πe2x(Acos(πx)Bsin(πx))π2e2x(Acos(πx)+Bsin(πx))


                 y'=4e2x(Acos(πx)+Bsin(πx))4πe2x(Acos(πx)Bsin(πx))π2e2x(Acos(πx)+Bsin(πx))


Substituting these equations into the original ODE yields:

4e2x(Acos(πx)+Bsin(πx))4πe2x(Acos(πx)Bsin(πx))π2e2x(Acos(πx)+Bsin(πx))

+4(2e2x(Acos(πx)+Bsin(πx))+πe2x(Acos(πx)Bsin(πx)))+(π2+4)(e2x(Acos(πx)+Bsin(πx)))=0

4e2x(Acos(πx)+Bsin(πx))4πe2x(Acos(πx)Bsin(πx))π2e2x(Acos(πx)+Bsin(πx))

8e2x(Acos(πx)+Bsin(πx))+4πe2x(Acos(πx)Bsin(πx))+π2e2x(Acos(πx)+Bsin(πx))+4e2x(Acos(πx)+Bsin(πx))=0

(48+4)e2x(Acos(πx)+Bsin(πx))+(4+4)πe2x(Acos(πx)Bsin(πx))+(π2π2)e2x(Acos(πx)+Bsin(πx))=0

00

Therefore, the solution is correct.

--Gonzalo Perez

Problem R2.4

K 2011 p.59 pb. 5

Problem Statement

Find a general solution to the given ODE. Check your answer by substituting into the original equation.

Given

y'+2πy+π2y=0

Solution

The characteristic equation of this ODE is therefore:

λ2+aλ+b=λ2+2πλ+π2=0

Evaluating the discriminant:

a24b=(2π)24(π2)=4π24π2=0

Therefore the equation has a real double root and a general homogenous solution of the form:

y=(c1+c2x)eax/2 [1]

And finally we find the general homogenous solution:

                                                       y=(c1+c2x)eπx


Checking:

We found that:

y=(c1+c2x)eπx

Differentiating y to obtain y' and y' respectively:

                                                  y'=c2eπxπ(c1+c2x)eπx


And,

y'=πc2eπxπc2eπx+π2(c1+c2x)eπx

                                                 y'=2πc2eπx+π2(c1+c2x)eπx


Substituting these equations into the original ODE yields:

2πc2eπx+π2(c1+c2x)eπx+2π(c2eπxπ(c1+c2x)eπx)+π2((c1+c2x)eπx)=0

2πc2eπx+π2(c1+c2x)eπx+2πc2eπx2π2(c1+c2x)eπx+π2(c1+c2x)eπx=0

(2π+2π)c2eπx+(π22π2+π2)(c1+c2x)eπx=0

00

Therefore this solution is correct.

  1. Kreyszig 2011, p.54-57.

K 2011 p.59 pb. 6

Problem Statement

Find a general solution to the given ODE. Check your answer by substituting into the original equation.

Given

10y32y+25.6y=0

Solution

Initially we modify the original ODE to put it in the form of a second-order homogenous linear ODE with constant coefficients:

10y32y+25.6y=0

Dividing both sides by 10:

y'3.2y+2.56y=0

The characteristic equation of this ODE is now therefore:

λ2+aλ+b=λ23.2λ+2.56=0

Evaluating the discriminant:

a24b=(3.2)24(2.56)=10.2410.24=0

Therefore the equation has a real double root and a general homogenous solution of the form:

y=(c1+c2x)eax/2 [1]

And finally we find the general homogenous solution:

                                                       y=(c1+c2x)e1.6x


Checking:

We found that:

y=(c1+c2x)e1.6x

Differentiating y to obtain y' and y' respectively:

                                                  y'=1.6c1e1.6x+c2e1.6x+1.6c2xe1.6x


And,

y'=(1.6)2c1e1.6x+1.6c2e1.6x+1.6c2e1.6x+(1.6)2c2xe1.6x

                                                 y'=(1.6)2c1e1.6x+3.2c2e1.6x+(1.6)2c2xe1.6x


Substituting these equations into the original ODE yields:

10[(1.6)2c1e1.6x+3.2c2e1.6x+(1.6)2c2xe1.6x]32[1.6c1e1.6x+c2e1.6x+1.6c2xe1.6x]+25.6[(c1+c2x)e1.6x]=0

25.6c1e1.6x+32c2e1.6x+25.6c2xe1.6x51.2c1e1.6x32c2e1.6x51.2c2xe1.6x+25.6c1e1.6x+25.6c2xe1.6x=0

(25.651.2+25.6)c1e1.6x+(3232)c2e1.6x+(25.651.2+25.6)c2xe1.6x=0

00

Therefore this solution is correct.

  1. Kreyszig 2011, p.54-57.

Report 2, Problem 2.5

Problem Statement

Problem 2.5 Find an Ordinary Differential Equation. Use y+ay+by=0 for the given basis.

ODE Solutions

Two Real Roots: y(x)=c1eλ1x+c2eλ2x
Double Real Roots: y(x)=c1eλ1x+c2xeλ2x

Reverse Engineering

eλx(λ2+aλ+b)=0

Case 1: a24b>0
λ2(λ1+λ2)λ+λ1λ2
y(λ1+λ2)y+λ1λ2y=0

Case 2: a24b=0
λ2(λ1+λ2)λ+λ1λ2
y(λ1+λ2)y+λ1λ2y=0

Problem Solutions

R2.5 K2011 p.59 pbs.16

e2.6x,e4.3x
(λ2.6)(λ+4.3)
λ2+1.7λ11.18


ODE Form:

                                                            y+1.7y11.18=0

R2.5 K2011 p.59 pbs.17

e5x,xe5x
(λ+5)2
λ2+25λ+5


ODE Form:

                                                            y+25y+5y=0


Created by [Daniel Suh] 20:57, 7 February 2012 (UTC)

Problem 2.6

Solved by: Andrea Vargas

Problem Statement

For the following spring-dashpot-mass system (in series) find the values for the parameters k,c,m knowing that the system has the double real root λ=3

Figure

File:R2.6.jpg

Solution

Previously, we have derived the following equation for such a system:
(From Sec 1 (d), (3) p.1-5)
m(yk+kcyk)+kyk=f(t)
We can write this equation in standard form by diving through by m:
yk+kcyk+kmyk=f(t)m
Here, we can take the coefficients of yk and yk as a and b:

yk+kcayk+kmbyk=f(t)m


Next,considering the double real root:
λ=3
We can find the characteristic equation to be:
(λ+3)2=λ2+6λ+9=0
Which is in the form:
λ2+aλ+b=0

Then, we know that a=6 and b=9:

Setting a and b from the first equation equal to these, we obtain:

                                                             kc=6andkm=9



Clearly, there is an infinite amount of solutions to this problem because we have 2 equations but 3 unknowns. This can be solved by fixing one of the values and finding the other two.

Example of Solution

An example of fixing one of the constants to find the other two is provided here. By solving the simple equations above, we can illustrate how to find k,c,m . We had:
kc=6andkm=9

If we fix the mass to m=10kg. We find:
k10=9

k=90

Then,
90c=6

c=15

Finally, we obtain:

                                                             m=10,k=90,andc=15


--Andrea Vargas 21:44, 7 February 2012 (UTC)


Problem 2.7: McLaurin Series

Problem Statement

Develop the McLaurin Series (Taylor Series at t=0) for et,cost,sint

Solution

et=n=0tnn!
=t00!+t11!+t22!+t33!+...tnn!
                                      et=1+t+t24+t36+...tnn!


cost=k=0(1)kt2k(2k)!
=(1)0t2(0)(2(0))!+(1)1t2(1)(2(1))!+(1)2t2(2)(2(2))!+(1)3t2(3)(2(3))!+...(1)kt2k(2k)!
                                      cost=1t22+t424+t6720...(1)t2k(2k)!


sint=k=0(1)kt2k+1(2k+1)!
=(1)0t2(0)+1(2(0)+1)!+(1)1t2(1)+1(2(1)+1)!+(1)2t2(2)+1(2(2)+1)!+(1)3t2(3)+1(2(3)+1)!+...(1)kt2k+1(2k+1)!
                                      sint=tt36+t5120t75040+...(1)kt2k+1(2k+1)!

--Egm4313.s12.team11.arrieta 17:07, 6 February 2012 (UTC)


Problem R2.8


Problem Statement


Find a general solution. Check your answer by substitution.

Problem 8


y+y+3.25y=0

Let:
λ=d/dx

Characteristic Equation


λ2+λ+3.25=0
Using the quadratic equation to find roots we get:
λ1=1+i(12)2
λ2=1i(12)2
Therefore:

yh(x)=e12x(c1cos(x3)+c2sin(x3))
Check By Substitution


y(x)=12e12x(c1cos(x3)+c2sin(x3)+e12x(3c1sin3x+3c1cos3x)
y(x)=14e12x(c1cos(x3)+c2sin(x3)12e12x(3c1sin3x+3c1cos3x)
12e12x(3c1sin3x+3c1cos3x)e12x(3c1cos(x3)3c2sin(x3)
Substituting y,y,y into the original equation, the result is

 y+y+3.25y=0

Problem 15


y+0.54y+(0.0729+π)y=0
Let: ddx=λ

Characteristic Equation


λ2+0.54λ+(0.0729+π)=0
Using the quadratic equation to find roots we get:
λ1=0.27+i(π)2
λ2=0.27i(π)2
Therefore:

yh(x)=e0.27x(c1cos(xπ)+c2sin(xπ)
Check By Substitution


y(x)=0.27e0.27x(c1cos(xπ)+c2sin(xπ)+e0.27x(πc1sinπx+πc1cosπx)

y(x)=0.0729e0.27x(c1cos(xπ)+c2sin(xπ)0.27e0.27x(πc1sinπx+πc1cosπx)

0.27e0.27x(πc1sinπx+πc1cosπx)+e0.27x(π(c1cos(xπ))π(c2sin(xπ)))

Substituting y,y,y into the original equation, the result is

 y+0.54y+(0.0729+π)y=0

Egm4313.s12.team11.gooding 03:41, 7 February 2012 (UTC)

Report 2, Problem 9

Problem Statement

Find and plot the solution for the L2-ODE-CC corresponding to

λ2+4λ+13

with r(x)=0

and initial conditions y(0)=1, y(0)=0

In another figure, superimpose 3 figs.:(a)this fig. (b) the fig. in R2.6 p.5-6, and (c) the fig. in R2.1 p.3-7

Quadratic Equation

λ=b±(b24ac)2a with a=1,b=4,c=13

λ=4±(424*1*13)2*1=2±3i

λ=2±3i

Homogeneous Solution

The solution to a L2-ODE-CC with two complex roots is given by

y(x)=ea2x[Acos(ωx)+Bsin(ωx)]

where λ=a2±ωi=2±3i

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

Solving for A and B

first initial condition y(0)=1

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

y(0)=e2*0[Acos(3*0)+Bsin(3*0)]=1

A=1

second initial condition y(0)=0

y(x)=ddxy(x)=ddxe2x[cos(3x)+Bsin(3x)]

y(x)=e2x[(2B3)sin(3x)+(3B2)cos(3x)]

y(0)=e2*0[(2B3)sin(3*0)+(3B2)cos(3*0)]

0=3B2

B=23

so the solution to our L2-ODE-CC is

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

Solution to R2.6

After solving for the constants kc and km we have the following homogeneous equation

y(x)+6y(x)+9y=0

Characteristic Equation and Roots

λ2+6λ+9=0

(λ+3)(λ+3)=0

We have a real double root λ=3

Homogeneous Solution

We know the homogeneous solution to a L2-ODE-CC with a double real root to be

y(x)=c1eλx+c2xeλx

Assuming object starts from rest

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

Plugging in λ and applying our first initial condition

y(0)=c1e3*0+c2*0*e3*0=1

c1=1

Taking the derivative and applying our second condition

y(x)=ddxy(x)=ddxe3x+c2xe3x

y(x)=3e3x+c2e3x3c2xe3x

y(0)=3e3*0+c2e3*03c2*0*e3*0=0

c2=3

Giving us the final solution

                 y(x)=e3x+3xe3x

Plots

Solution to this Equation

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

File:Plotr2 9.jpg

Superimposed Graph

Our solution: y(x)=e2x[cos(3x)+23sin(3x)] shown in blue

Equation for fig. in R2.1 p.3-7: y(x)=54e2x14e5x shown in red

Equation for fig. in R2.6 p.5-6:y(x)=e3x+3xe3x shown in green

File:R2superposed.jpg

Egm4313.s12.team11.imponenti 03:38, 8 February 2012 (UTC)

Template:CourseCat