Coupled ODE with missing connecting derivatives

Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
4 replies · 2K views
FrankST
Messages
23
Reaction score
0
Hi,

I have a coupled system of ODE like:

w1'' = A w2'' + B w1 + C w2
w2'' = D w1'' + E w1 + F w2

I need to solve it analytically but it seems it cannot be solved using eigenvalue method. My concern is first that if this system have sufficient equations and if so how it can be solved analytically. If there is no analytical solution, is there any numerical solution for that?


Thanks for your assistance,

Frank
 
Physics news on Phys.org
Welcome to PF!

Hi Frank! Welcome to PF! :smile:

Isn't that just PW'' = QW, with P and Q matrices? :wink:
 
Solving is possible by substitution.
w1'' = A w2'' + B w1 + C w2
w2'' = D w1'' + E w1 + F w2
w2''= D(A w2'' + B w1 + C w2) + E w1 + F w2
w1 = (w2'' - D A w2''- D C w2 - F w2 ) / ( D A + E )
w1'' = (w2'''' - D A w2''''- D C w2'' - F w2'' ) / ( D A + E )
Then, bringing back w1 and w1'' into w2'' = D w1'' + E w1 + F w2 leads to a fourth order linear ODE with only one unknown function w2
 
I personally prefer the more elegant operator approach. First, I'll use lower-case constants and subscripts to make it easier to read. The operator approach is more methodical and easy to apply:

[tex]w_1^{''}=aw_2^{''}+bw_1+cw_2[/tex]

[tex]w_2^{''}=dw_1^{''}+ew_1+fw_2[/tex]

Now move everything over:

[tex]w_1^{''}-bw_1-aw_2^{''}-cw_2=0[/tex]

[tex]w_2^{''}-fw_2-dw_1^{''}-ew_1=0[/tex]

Now convert to operators and line-up everything:

[tex](D^2-b)w_1-(aD^2+c)w_2=0[/tex]

[tex](dD^2+e)w_1-(D^2-f)w_2=0[/tex]

Now, to eliminate [itex]w_2[/tex], operate on the first by [itex](D^2-f)[/itex] and on the second by [itex]-(aD^2+c)[/itex] and add. I get then:<br /> <br /> [tex](D^2-f)(D^2-b)w_1-(aD^2+c)(dD^2+e)w_1=0[/tex]<br /> <br /> Do something similar to get the equation in [itex]w_2[/itex][/itex]
 
Thank you guys for you valuable advices.

Your answers helped me a lot.