Differential Equations Systems problem (3X3 matrix) Please help

AI Thread Summary
The discussion revolves around solving a system of differential equations represented by a 3x3 matrix. The characteristic polynomial yields eigenvalues r1 = 1 and r2 = r3 = 2, with corresponding eigenvectors derived for each eigenvalue. The user is unsure how to proceed after finding the eigenvectors, particularly for the repeated eigenvalue. The general solution is proposed as a combination of the eigenvalues and eigenvectors, but there is confusion regarding the specific eigenvectors used. Ultimately, the correct general solution incorporates the eigenvectors associated with each eigenvalue, leading to a complete representation of the system's behavior.
VinnyCee
Messages
486
Reaction score
0
Here is the problem:

X'\,=\,\left(\begin{array}{ccc}3 & -1 & -1\\1 & 1 & -1\\1 & -1 & 1\end{array}\right)\,X

Here is what i have so far:

det(A\,-\,rI)\,=\,0

det\left[\left(\begin{array}{ccc}3 & -1 & -1\\1 & 1 & -1\\1 & -1 & 1\end{array}\right)\,-\,\left(\begin{array}{ccc}r & 0 & 0\\0 & r & 0\\0 & 0 & r\end{array}\right)\right]\,=\,\left(\begin{array}{ccc}3\,-\,r & -1 & -1\\1 & 1\,-\,r & -1\\1 & -1 & 1\,-\,r\end{array}\right)

-(r\,-\,1)\,(r^2\,-\,4r\,+4)\,=\,0

r_1\,=\,1,\,\,\,\,r_2\,=\,2,\,\,\,\,r_3\,=\,2

r\,=\,2 is repeated once.

(A\,-\,r_1I)\,\xi\,=\,0

\left[\left(\begin{array}{ccc}3 & -1 & -1\\1 & 1 & -1\\1 & -1 & 1\end{array}\right)\,-\,\left(\begin{array}{ccc}1 & 0 & 0\\0 & 1 & 0\\0 & 0 & 1\end{array}\right)\right]\,\left(\begin{array}{ccc} \xi_1 \\ \xi_2 \\ \xi_3 \end{array}\right)=\,\left(\begin{array}{ccc}2 & -1 & -1\\1 & 0 & -1\\1 & -1 & 0\end{array}\right)\,\left(\begin{array}{ccc} \xi_1 \\ \xi_2 \\ \xi_3 \end{array}\right)\,=\,0

Now, multiplying out the right hand side of the equation above, and rref'ing, I get:

\left(\begin{array}{cccc} 1 & 0 & -1 & 0 \\ 0 & 1 & -1 & 0 \\ 0 & 0 & 0 & 0 \end{array}\right)

Which means that,

\xi_1\,-\,\xi_3\,=\,0
\xi_2\,-\,\xi_3\,=\,0

So then,

\xi_1\,=\,\xi_2\,=\,\xi_3?

Here is where I am stuck! What am I supposed to do now to solve this system?
 
Last edited:
Physics news on Phys.org
Just pick a number, like \xi_1 = \xi_2 = \xi_3 = 1 or something. That's just saying that any vector with all three components the same is an eigenvector corresponding to eigenvalue r = 1. Is that what you are asking?
 
Ok, \xi_1=\xi_2=\xi_3=1.

What do I do with this information?

Do I just solve for the r_2=r_3=2 part and that is the answer?

I am still not understnding exactly what it is that I am supposed to be doing with the \xi_n values.
 
(note: unfortunately, there's a little difference between VinnyCee's notation and what follows. VinnyCee uses ξ1, ξ2, and ξ3 to represent eigenvector components. In what follows, ξ1, ξ2, and ξ3 are the actual eigenvectors themselves. sorry.)

from your results (& standard textbook): for single unique eigenvalue r1=1 and its eigenvector ξ1=<1, 1, 1>, a specific solution is:

X_{1}(t) \ = \ C_{1} \cdot \exp(r_1 t) \cdot \xi_{1} \ = \ C_1 \cdot \exp(t) \cdot \left ( \begin{array}{c} 1 \\ 1 \\ 1 \end{array} \right )

then for repeated eigenvalue r23=2 and its eigenvectors ξ2=<2, 1, 1> and ξ3=<1.33, -1, 2.33>, a specific solution is:

X_{23}(t) \ = \ C_{2} \cdot \exp(r_{23} t) \cdot \xi_{2} \ + \ C_{3} \cdot \exp(r_{23} t) \cdot ( \xi_{2} t \ + \ \xi_{3} ) \ = \ C_{2} \cdot \exp(2 t) \cdot \left ( \begin{array}{c} 2 \\ 1 \\ 1 \end{array} \right ) \ + \ C_3 \cdot \exp(2t) \cdot \left ( \left ( \begin{array}{c} 2 \\ 1 \\ 1 \end{array} \right ) \cdot t \ + \ \left ( \begin{array}{c} 1.33 \\ -1 \\ 2.33 \end{array} \right ) \right )

\mbox{General Solution} \ = \ \textbf{X(t)} \ = \ X_{1}(t) \ + \ X_{23}(t)

general solution is (X1(t) + X23(t)), and constants C1, C2, and C3 obtained from initial cond.
 
Last edited:
The general solution is...

Is this the correct general solution?:

X\,=\,\ C_1 \cdot e^t \cdot \left ( \begin{array}{c} 1 \\ 1 \\ 1 \end{array} \right ) \,+\,\ C_2 \cdot e^{2t} \cdot \left ( \begin{array}{c} 1 \\ 0 \\ 1 \end{array} \right ) \,+\,\ C_3 \cdot e^{2t} \cdot \left [ \left( \begin{array}{c} 1 \\ 0 \\ 1 \end{array} \right )\cdot t\,+\,\left ( \begin{array}{c} 1 \\ 1 \\ 0 \end{array} \right ) \right]

This is using different eigenvectors than you calculated. I got:

\xi^{(2)}\,=\,\left ( \begin{array}{c} 1 \\ 0 \\ 1 \end{array} \right ) and \xi^{(2)}\,=\,\left ( \begin{array}{c} 1 \\ 1 \\ 0 \end{array} \right ).

I think that they are just multipules of your solutions.
 
Thread 'Variable mass system : water sprayed into a moving container'
Starting with the mass considerations #m(t)# is mass of water #M_{c}# mass of container and #M(t)# mass of total system $$M(t) = M_{C} + m(t)$$ $$\Rightarrow \frac{dM(t)}{dt} = \frac{dm(t)}{dt}$$ $$P_i = Mv + u \, dm$$ $$P_f = (M + dm)(v + dv)$$ $$\Delta P = M \, dv + (v - u) \, dm$$ $$F = \frac{dP}{dt} = M \frac{dv}{dt} + (v - u) \frac{dm}{dt}$$ $$F = u \frac{dm}{dt} = \rho A u^2$$ from conservation of momentum , the cannon recoils with the same force which it applies. $$\quad \frac{dm}{dt}...

Similar threads

Replies
5
Views
915
Replies
2
Views
1K
Replies
5
Views
2K
Replies
7
Views
2K
Replies
4
Views
1K
Replies
3
Views
3K
Replies
1
Views
1K
Back
Top