University of Florida/Egm4313/s12.team5.R2

From testwiki
Jump to navigation Jump to search

Template:Big3

R2.1

Question

Given the two roots and the initial conditions:

λ1=2,λ2=+5

(1.0)

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

(1.1)

Part 1.

Find the non-homogeneous L2-ODE-CC in standard form. Find the solution in terms of the initial conditions and the general excitation r(x).
Now with no excitation and plot the solution:

r(x)=0

(1.2)


Part 2.

Generate 3 non-standard (and non-homogeneous) L2-ODE-CC that admit the two values λ1=2,λ2=+5 as the two roots of the corresponding characteristic equation.

Solution

Part 1.

Characteristic Equation:

(λλ1)(λλ2)=0

(1.3)

(λ(2))(λ5)=0

(1.4)

λ23λ10=0

(1.5)

Non-Homogeneous L2-ODE-CC:

y3y10y=r(x)standardform

(1.6)

Homogenous Solution:

yh(x)=c1e2x+c2e5x

(1.7)

Overall Solution:

y(x)=c1e2x+c2e5x+yp(x)

(1.8)

y(x)=2c1e2x+5c2e5x+y'p(x)

(1.9)

Satisfy Initial Conditions:

y(0)=1=c1+c2+yp(0)

(1.10)

y(0)=0=2c1+5c2+y'p(0)

(1.11)

No excitation:

r(x)=0yp(x)=0y'p(x)=0

(1.12)

From (1.10):

c1+c2=1c1=1c2

(1.13)

From (1.11):

2c1+5c2=0

(1.14)

5c2=2c1

(1.15)

52c2=c1

(1.16)

Plug (1.13) into (1.16):

52c2=1c2

(1.17)

Solve for c2:

c2=27

(1.18)

Plug (1.18) into (1.13) and solve for c1:

c1=57

(1.19)

Therefore, the final solution in terms of the initial conditions and the general excitation r(x) is:

y(x)=57e2x+27e5x

(1.20)


Plot of Solution:

Figure 1File:IEA 2.1 plot.png

Part 2.

Three non-standard and non-homogeneous solutions using the same roots given above:

1.

2(λ(2))(λ5)=0

(1.21)

2λ26λ20=0

(1.22)

2.

3(λ(2))(λ5)=0

(1.23)

3λ29λ30=0

(1.24)

3.

4(λ(2))(λ5)=0

(1.25)

4λ212λ40=0

(1.26)

Author

This problem was solved and uploaded by [Derik Bell]

This problem was proofread by: David Herrick

R2.2

Question

From Section 5 in the notes, pg. 5-6

Solve the L2-ODE-CC from pg. 5-5, equation (4)

 y10y+25y=r(x) 

For initial conditions:

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

Where there is no excitation, i.e.  r(x)=0 

Solution

 λ210λ+25=0 

Factoring:

 (λ5)2=0 

Thus, there is a double root, where  λ=5 

Because there is no excitation, the homogenous solution will be the same as the final solution. Thus, the solution will be:

 yh(x)=y(x)=C1e5x+C2xe5x 

Plugging in the initial condition,

 y(0)=1=C1e0+C2(0)e0 

 C1=1 

In order to solve for the second constant, we need to take a derivative.

 y'h(x)=y(x)=5C1e5x+C2e5x+5C2xe5x 

Plugging in C1 and the other initial condition:

 y'h(0)=y(0)=0=5(1)e0+C2e0+5C2(0)e0 

 5+C2=0C2=5 

So, the final solution is:

 y(x)=e5x5xe5x 

Matlab Code:

x = 0:0.001:1;

f = exp(5*x) - 5*x.*exp(5*x);

plot(x,f)

xlabel ('x')

ylabel ('y(x)')

File:Report2Matlab'.jpg

Author

This problem was solved and uploaded by [David Herrick]
This problem was checked by [William Knapper]

R2.3

Question

Complete problems 3 and 4 from p.59 of K 2011.

Solution

The two problems have the same instructions: "Find a general solution. Check your answer by substitution."

Problem 3

3.y+6y+8.96y=0

We start by using the characteristic equation of this ODE in order to find the roots. We use:

λ2+aλ+b=0

where a = 6 and b = 8.96. Taking the discriminant, we find that

624(8.96)=0.16>0

A positive discriminant implies that there exists two real roots, with a general solution of the form:

y=c1eλ1x+c2eλ2x

To find the unknowns, we find the roots of the above equation:

λ2+6λ+8.96=0

λ1=3.2,λ2=2.8

Therefore, the correct general solution is:

y=c1e3.2x+c2e2.8x 

To check that this solution is correct, we take the first and second derivatives:

y'=3.2c1e3.2x2.8c2e2.8x

y'=10.24c1e3.2x+7.84c2e2.8x

Now, we plug these values into the original equation:

10.24c1e3.2x+7.84c2e2.8x19.2c1e3.2x16.8c2e2.8x+8.96c1e3.2x+8.96c2e2.8x=0

Upon further inspection, all the terms on the left side of the equation cancel out and equal zero.

Problem 4

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

We start by using the characteristic equation of this ODE in order to find the roots. We use:

λ2+aλ+b=0

where a = 4 and b = (π^2 + 4). Taking the discriminant, we find that

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

A negative discriminant implies that there exists complex conjugate roots to this equation. The general solution for this type of equation is

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

Where

ω=b14a2=(π2+4)(14)42=π

This yields

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

To verify that this is the correct solution, we first find the first and second derivatives of the solution:

y=e2x(Acos(πx)+Bsin(πx))
y'=e2x(Aπsin(πx)+Bπcos(πx))2e2x(Acos(πx)+Bsin(πx))
y'=e2x[Aπ2cos(πx)Bπ2sin(πx)+2Aπsin(πx)2Bπsin(πx)]2e2x[Aπsin(πx)+Bπcos(πx)2Acos(πx)2Bsin(πx)]

We then plug these values into the original problem to get this somewhat lengthy equation:

e2xAπ2cos(πx)e2xBπ2sin(πx)+2e2xAπsin(πx)2e2xBπsin(πx)+2e2xAπsin(πx)2e2xBπcos(πx)+
4e2xAcos(πx)+4e2xBsin(πx)4e2xAπsin(πx)+4e2xBπcos(πx)
8e2xAcos(πx)8e2xBsin(πx)e2xAπ2cos(πx)+e2xBπ2sin(πx)+4e2xAcos(πx)+4e2xBsin(πx)=0

Upon further inspection, all the terms on the left side of the equation cancel out to 0.

Author

This problem was solved and uploaded by: (Will Knapper)

This problem was proofread by: David Herrick

R2.4

Question

K 2011 p.59 pbs.5,6. Find a general solution. Check your answer by substitution.

Problem 5

y+2πy+π2y=0 

Solution

Determine the characteristic equation and the nature of the determinate:

y=eλx 
y=λeλx 
y=λ2eλx 

Characteristic Equation:

λ22πλ+π2λ=0 
a24b=(2π)24(π2)=4π24π2=0 

0 indicates a real double root, therefore:

yg=ea2x(C1+C2x) 

a=2π  from the original equation. Therefore:

yg=eπx(C1+C2x) 
y'g=eπx(C2πC1πC2x) 
y'g=eπx(C1π2+C2π2x2C2π) 

Substituting into the original equation gives:

eπx[(C1π2+C2π2x2C2π)+2π(C2πC1πC2x)+π2(C1+C2x)]=0 
=eπx[C1π2+C2π2x2C2π+2C2π2C1π22C2π2x+C1π2+C2π2x] 

Canceling gives 0 = 0, confirming the general solution found above.

Author

This problem was solved and uploaded by: (John North)

This problem was proofread by: David Herrick

Question 2

Problem 5

10y32y+25.6y=0 

Solution

Reduce original equation so that the first coefficient is 1:

y3.2y+2.56y=0 

Determine the characteristic equation and the nature of the determinate:

y=eλx 
y=λeλx 
y=λ2eλx 

Characteristic Equation:

λ23.2λ+2.56=0 
a24b=3.224(2.56)=10.2410.24=0 

0 indicates a real double root, therefore:

yg=ea2x(C1+C2x) 

a=(3.2)2=1.6 from the original equation. Therefore:

yg=e1.6x(C1+C2x) 
y'g=e1.6x(1.6C1+1.6C2x+C2) 
y'g=e1.6x(2.56C1+2.56C2x+3.2C2) 

Substituting into the original equation gives:

e1.6x[(2.56C1+2.56C2x+3.2C2)3.2(1.6C1+1.6C2x+C2)+2.56(C1+C2x)]=0 
=e1.6x[(2.56C1+2.56C2x+3.2C25.12C15.12C2x3.2C2+2.56C1+2.56C2x]=0 

Canceling gives 0 = 0, confirming that general solution found above.

Author

This problem was solved and uploaded by: (John North)

This problem was proofread by: David Herrick

R2.5

Question

Problems 16 and 17 from pg. 59 of the textbook. The two questions have the same instructions: Find an ODE of the form
y'+ay+by=0
for the given basis.

Solution

Problem 16

16.e2.6x,e4.3x
This basis implies that there are two distinct, real roots: 2.6 and -4.3. We put these roots into the characteristic equation to get: (λ2.6)(λ+4.3)=0λ2+1.7λ11.18=0
Then, you convert the characteristic equation to a regular ODE:

y'+1.7y11.18y=0 

Problem 17

17.e5x,xe5x
A basis of this form implies that the solution is of the form:
y=(c1+c2x)eax/2
The real double root can be found by this equation:
λ=a/2=5
That means that the characteristic equation can be found by solving:
(λ+5)2=λ2+25λ+5
Converting the characteristic equation to the ODE, we get:

y'+25y+5y=0 

Author

This problem was solved and uploaded by: (Will Knapper)
This problem was proofread by (Josh House)

R2.6

Question

Realize spring-dashpot-mass system in series as shown in Fig. p.1-4 with the similar characteristic as in (3)p.5-5, but with the double real root λ=3 , i.e., find values for parameters k,c,m.

Solution

First we put the double root into the characteristic equation:
(λ+3)2=0
λ2+6λ+9=0

The corresponding L2-ODE is:
y'+6y+9y=0

The equation of motion of the spring-dashpot-mass system is:
m(yk+kcyk)+kyk=0 [1]
or:
myk+mkcyk+kyk=0

Solving for the parameters k,c,m:
m=1

mkc=6

k=9

We get:

m=1 

k=9

c=32

Author

This problem was solved and uploaded by: (Radina Dikova)

This problem was proofread by: Mike Wallace

R2.7

Question

Develop the McLauren series (Taylor series at t=0) for e^t, cos(t), and sin(t).

Solution

The general form of the Taylor series expansion at point a is

f(a)+f(a)1!(xa)+f(a)2!(xa)2+f(a)3!(xa)3+...

for et at t=0:

e0+e01!(x)+e02!(x)2+e03!(x)3+...


=1+x+x22!+x33!+...


for cos(t) at t=0:

cos(0)+sin(0)1!(x)+cos(0)2!(x)2+sin(0)3!(x)3+cos(0)4!(x)4+...


=1x22!+x44!+...


for sin(t) at t=0:

sin(0)+cos(0)1!(x)+sin(0)2!(x)2+cos(0)3!(x)3+sin(x)4!(x)4+cos(0)5!(x)5+...


=xx33!+x55!+...



Author

This problem was solved and uploaded by: (John North)
This problem was checked by: William Knapper

R2.8

Question

Find a general solution. Check your answer by substitution.

8. y'+y'+3.25y=0


15. y'+0.54y'+(0.0729+π)y=0

Solution

8. y'+y'+3.25y=0

Writing the characteristic equation:

λ2+λ+3.25=0


Which is now in the form:

λ2+aλ+b=0


Solving for the two roots:

λ1=12(a+a24b),λ2=12(aa24b)


We will find the discriminant to be less than 0, leading to complex conjugate roots:

a24b=124(3.25)=12


Which leads to a solution of the form:

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


Where ω2=b14a2, therefore:

ω2=3.2514(1)2=3ω=3


Giving us a solution of:

y=ex/2(Acos(3x)+Bsin(3x))


Checking our solution:

y=ex/2(Acos(3x)+Bsin(3x))
y'=12ex/2(Acos(3x)+Bsin(3x))+ex/2(3sin(3x)+3Bcos(3x))
y'=14ex/2(Acos(3x)+Bsin(3x))12ex/2(3Asin(3x)+3Bcos(3x))12ex/2(3Asin(3x)+3Bcos(3x))+ex/2(3Acos(3x)3Bsin(3x))


Plugging all of this into the original differential equation and combining like terms gives us:

[14312+3.25]ex/2Acos(3x)
[32+323]ex/2Asin(3x)
[3232+3]ex/2Bcos(3x)
[14312+3.25]ex/2Bsin(3x)


The terms in the brackets all add up to zero, thus verifying we have a correct solution to our problem.



15. y'+0.54y'+(0.0729+π)y=0


Writing the characteristic equation:

λ2+0.54λ+(0.0729+π)=0


Which is now in the form:

λ2+aλ+b=0


Solving for the two roots:

λ1=12(a+a24b),λ2=12(aa24b)


We will find the discriminant to be less than 0, leading to complex conjugate roots:

a24b=0.5424(0.0729+π)=4π


Which leads to a solution of the form:

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


Where ω2=b14a2, therefore:

ω2=(0.0729π)14(0.54)2=πω=π


Giving us a solution of:

y=e0.27x(Acos(πx)+Bsin(πx))


Checking our solution:

y=e0.27x(Acos(πx)+Bsin(πx))
y'=0.27e0.27x(Acos(πx)+Bsin(πx))+e0.27x(πAsin(πx)+πBcos(πx))
y'=0.0729e0.27x(Acos(πx)+Bsin(πx))0.27e0.27x(πAsin(πx)+πBcos(πx))0.27e0.27x(πAsin(πx)+πBcos(πx))+e0.27x(πAcos(πx)πBsin(πx))


Plugging all of this into the original differential equation and combining like terms gives us:

[0.0729π0.1458+0.0729+π]e0.27xAcosπx
[0.27π+0.27π0.54π]e0.27xAsinπx
[0.27π0.27π+0.54π]e0.27xBcosπx
[0.0729π0.1458+0.0729+π]e0.27xBsinπx


The terms in the brackets all add up to zero, thus verifying we have a correct solution to our problem.

Author

This problem was solved and uploaded by: (Josh House)
This problem was proofread by: (Radina Dikova)

R2.9

Question

For the given ODE, find and plot the solution for the L2-ODE-CC corresponding to

 λ2+4λ+13=0 

In another figure superpose this figure, Fig. from R2.6 p.5-6, and the Fig. from R2-1 p.3-7.

This corresponds to the L2-ODE-CC in standard form

 y+4y+13y=r(x) 

Given

Initial conditions

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

Where there is no excitation

 r(x)=0 

Solution

Solving the quadratic equation for  λ  determines the roots of the ODE


 λ1=b+b2(4ac)2a=4+(4)24(1)(13)2(1)


 λ2=bb2(4ac)2a=4(4)24(1)(13)2(1)

The roots of the quadratic equation come out to be:

 λ=(2±3i) 


The general homogeneous solution of the ODE will resemble the following form.

 y=C1ea1xcosb1x+C2ea2xsinb2x 

Plugging in the roots into the equation we get

 y=C1e2xcos3x+C2e2xsin3x 


Implementing initial conditions where:  y(0)=1,y(0)=0 

 y(0)=C1e2(0)cos3(0)+C2e2(0)sin3(0)=1 

 y(0)=2C1e2xcos3x3C1e2xsin3x2C2e2xsin3x+3C2e2xcos3x=0 

 y(0)=2(1)e2(0)cos3(0)+3C2e2(0)cos3(0)=2+3C2=0; 


So from our initial conditions the constants are determined to be

 C1=1,C2=2/3 


The final solution becomes

 y(x)=e2xcos3x+23e2xsin3x 

Figure 1 for R2.9 Superposed Figs. for R2.9

Author

This problem was solved and uploaded by: Mike Wallace
This problem was proofread by John North

Contribution Summary

Problem 2 was solved and Problems 1, 3, and 4 were proofread by David Herrick 15:56, 6 February 2012 (UTC)

Problem 8 was solved and Problem 5 was proofread by Josh House 23:41, 7 February 2012 (UTC)

Problem 6 was solved and Problem 8 was proofread by Radina Dikova 16:26, 7 February 2012 (UTC)

Problem 9 was solved and Problem 6 was proofread by Mike Wallace 23:35, 7 February 2012 (UTC)

Problems 3 and 5 were solved and Problems 2 and 7 were proofread by William Knapper 05:04, 8 February 2012 (UTC)

Problem 1 was solved by Derik Bell 15:44, 8 February 2012 (UTC)

Problems 4 and 7 were solved and Problem 9 was proofread by John North 15:47, 8 February 2012 (UTC)

References

  1. Vu-Quoc , Loc. 2012. IEA S12 Lecture notes, audios, videos. Section 1(d) https://docs.google.com/Doc?id=dc82nfgb_858gmw855hj

Template:CourseCat