Coupled ODE with missing connecting derivatives

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:

w_1^{''}=aw_2^{''}+bw_1+cw_2

w_2^{''}=dw_1^{''}+ew_1+fw_2

Now move everything over:

w_1^{''}-bw_1-aw_2^{''}-cw_2=0

w_2^{''}-fw_2-dw_1^{''}-ew_1=0

Now convert to operators and line-up everything:

(D^2-b)w_1-(aD^2+c)w_2=0

(dD^2+e)w_1-(D^2-f)w_2=0

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

Your answers helped me a lot.
 
Thread 'Direction Fields and Isoclines'
I sketched the isoclines for $$ m=-1,0,1,2 $$. Since both $$ \frac{dy}{dx} $$ and $$ D_{y} \frac{dy}{dx} $$ are continuous on the square region R defined by $$ -4\leq x \leq 4, -4 \leq y \leq 4 $$ the existence and uniqueness theorem guarantees that if we pick a point in the interior that lies on an isocline there will be a unique differentiable function (solution) passing through that point. I understand that a solution exists but I unsure how to actually sketch it. For example, consider a...
Back
Top