Coupled non-homogenous eigenvalue problem help?

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
1 reply · 2K views
de1337ed
Messages
10
Reaction score
0

Homework Statement



Equations:
[itex]\frac{dv_{1}}{dt}[/itex] = [itex]-v_{1}[/itex] - [itex]\frac{2v_{2}}{3}[/itex] + 1 + [itex]\frac{t}{3}[/itex]
[itex]\frac{dv_{2}}{dt}[/itex] = [itex]-2v_{2}[/itex] - 1 - 2t

Initial Conditions:
[itex]v_{1}(0)[/itex] = 6
[itex]v_{2}(0)[/itex] = -6

2. The attempt at a solution

Defined the following:
v(t) =
[ [itex]v_{1}(t)[/itex]
[itex]v_{2}(t)[/itex] ]

[itex]\frac{dv(t)}{dt}[/itex] =
[ [itex]\frac{dv_{1}(t)}{dt}[/itex]
[itex]\frac{dv_{2}(t)}{dt}[/itex] ]

v(0)=
[ 6
-6 ]

s =
[ 1 + [itex]\frac{t}{3}[/itex]
-1 - 2t ]

A =
[ -1 [itex]\frac{-2}{3}[/itex]
0 -2 ]

Particular:

[itex]\frac{dv(t)}{dt}[/itex] = Av + s
Try v = [itex]k_{0}[/itex] + [itex]k_{1}t[/itex]

And now I'm stuck. I know what to do in homogenous cases when we use [itex]cxe^{\lambda t}[/itex] because in that case, you can represent [itex]\frac{dv(t)}{dt}[/itex] as [itex]\lambda v[/itex]
But what do I do in the non-homogeneous case when there is a polynomial driving function? Thank you.
 
Physics news on Phys.org
You do exactly the same thing as in homogeneous case. You want to write your equation as v'=Av+f(t), where v', v, and f(t) are vectors, and A is a 2x2 matrix. You've gotten that far. Now, you need to diagonalize A. So say you found matrix S such that Ad=S-1AS, where Ad is your diagonalized matrix.

Having that, multiply your equation by S-1 from the left on both sides. You get.

S-1v' = S-1Av + S-1 f(t)

Now I'm going to insert I=SS-1 between A and v.

S-1v' = S-1AS S-1v + S-1 f(t)

S-1v' = Ad S-1v + S-1 f(t)

Because Ad is diagonal, you can re-write the above as two completely independent driven first order linear ordinary differential equations. Their driving terms are completely separate, so you should have no trouble solving each independently. Then simply multiply your solution vector by S from the left, and you'll get your answer.