Solving a Second-Order ODE with Non-Real Roots

  • Thread starter Thread starter Tom1
  • Start date Start date
  • Tags Tags
    Ode
Click For Summary

Homework Help Overview

The discussion revolves around solving a second-order ordinary differential equation (ODE) with non-real roots, specifically the equation y'' + 2y' + 2y = 0, along with initial conditions y(π/4) = 2 and y'(π/4) = -2.

Discussion Character

  • Mixed

Approaches and Questions Raised

  • Participants explore the characteristic polynomial and its non-real roots, with some questioning the assumptions made in the calculations. There is a focus on transforming complex solutions into real forms using Euler's identity.

Discussion Status

The discussion has progressed through various interpretations of the roots and the general solution. Some participants have offered guidance on expressing complex exponentials in terms of sine and cosine, while others have pointed out potential errors in assumptions and calculations.

Contextual Notes

Participants express confusion regarding the use of complex numbers in ODEs and the application of initial conditions to determine constants in the general solution.

Tom1
Messages
23
Reaction score
0

Homework Statement



Solve the initial value problem:
y''+2y'+2y=0
y([tex]\pi[/tex]/4)=2
y'([tex]\pi[/tex]/4)=-2

Homework Equations



Included in 1

The Attempt at a Solution



I assumed that y=e^rt and came up with the characteristic polynomial r^2+2r+2=0 but when solving it for the two roots, they come back as non-real. Where have I gone wrong? thanks
 
Physics news on Phys.org
You've done nothing wrong, the solutions to ODE's may be complex. However, one may transform them into real solutions using Euler's identity.
 
I have never done one with complex numbers before. Mind walking me through this one?
 
Tom1 said:
I have never done one with complex numbers before. Mind walking me through this one?
No problem, what are your two solutions?
 
Assuming a=1, b=2 c=2 and r1= [tex]\lambda[/tex]+i[tex]\mu[/tex] r2=[tex]\lambda[/tex]-i[tex]\mu[/tex]

r1=(-1+(-4)^1/2)/2
= -(1/2)+i(2)^(1/2)
r2 = -(1/2)-i(2)^(1/2)

Correct?
 
I have different solutions. In your solution, you make the error of assuming that,

[tex]\sqrt{4} = 2\sqrt{2}[/tex]

Which is not the case. Using the quadratic formula,

[tex]r = \frac{-2\pm\sqrt{4-4\cdot1\cdot2}}{2} = \frac{-2\pm i\sqrt{4}}{2}[/tex]

[tex]r = -1\pm i[/tex]

Do you agree?
 
Indeed, that's a mistake I commonly make.

So with those solutions, the general solution of the IVP would be:

y=c1(e^-t)cos(1)+c2(e^-t)sin(-1)

?
 
Tom1 said:
Indeed, that's a mistake I commonly make.

So with those solutions, the general solution of the IVP would be:

y=c1(e^-t)cos(1)+c2(e^-t)sin(-1)

?
I think that you've made a mistake somewhere (perhaps missing a t here and there...?). Let's start from the top,

[tex]y(t) = A\exp\left\{t\left(i-1\right)\right\} + B\exp\left\{-t\left(i+1\right)\right\}[/tex]

Which may be re-written,

[tex]y(t) = e^{-t}\left[Ae^{it}+Be^{-it}\right][/tex]

Now, can you write the two complex exponentials in terms of sines and cosine?
 
Last edited:
Ah, I think I see my mistake. If u=1 then:

y=(Ae^-t)cos(t)+(Be^-t)sin(-t)
 
  • #10
Almost there, note that,

[tex]\sin\left(-\theta\right) = -\sin\theta[/tex]

Furthermore, one must not use the same constants (if you work through the full solution you will see why). Hence, the general solution will be of the form,

[tex]y(t) = e^{-t}C\cos t + e^{-t}D\sin t[/tex]
 
  • #11
Just figured out how to use the tex feature to make it easier on your eyes. I'm a little confused about the variables:

[tex]y(t) = Ae^{-t}C\cos t - Be^{-t}D\sin t[/tex]

Is that the solution?
 
  • #12
I tell you what, since it is clear that you've put the effort in, I'll provide a complete solution in this case. Starting from the second equation in post #8 one may write,

[tex]y(t) = e^{-t}\left[A\left\{\cos t + i\sin t\right\}+B\left\{\cos\left(-t\right) + i\sin\left(-t\right)\right\}\right][/tex]

Noting that cosine and sine are even and odd function respectively, one may write,

[tex]y(t) = e^{-t}\left[A\left\{\cos t + i\sin t\right\}+B\left\{\cos t - i\sin t\right\}\right][/tex]

Collecting the coefficients of sine and cosine,

[tex]y(t) = e^{-t}\left[\left(A+B\right)\cos t + \left(A-B\right)i\sin t\right][/tex]

Let us now define two new constants such that,

[tex]C:= A+B \hspace{2cm} D:=\left(A-B\right)i \hspace{2cm}C,D\in\mathbb{R}[/tex]

Hence, the general solution,

[tex]y(t) = e^{-t}\left[C\cos t + D\sin t\right][/tex]

Do you follow?
 
Last edited:
  • #13
Alright, I see now...but it causes another question to arise: When I go back to find the constants, I apply the initial conditions to solve for C&D, correct?
 
  • #14
Tom1 said:
Alright, I see now...but it causes another question to arise: When I go back to find the constants, I apply the initial conditions to solve for C&D, correct?
Correct indeed :approve:
 
  • #15
Thank you much for all your help.
 
  • #16
Tom1 said:
Thank you much for all your help.
No problem, it was a pleasure :smile:
 

Similar threads

Replies
4
Views
2K
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 16 ·
Replies
16
Views
2K
  • · Replies 33 ·
2
Replies
33
Views
4K
Replies
6
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K
Replies
7
Views
2K
  • · Replies 18 ·
Replies
18
Views
3K
Replies
3
Views
2K