Coupled differential equations

plmokn

1. The problem statement, all variables and given/known data
I need to find the general solution for:
[latex]R \frac{dQ_{1}}{dt} + \frac{1}{C}(Q_{1}-Q_{2})=0[/latex]
[latex]R\frac{dQ_{1}}{dt} + L\frac{d^2Q_{2}}{dt^2}}=0[/latex]
For the case [latex]R=1\times 10^3[/latex], [latex]L=4\times 10^{-3}[/latex], [latex]C=1 \times 10^{-9}[/latex].

3. The attempt at a solution
As an attempt I looked for solutions of the form, [latex]Q_{1}=q_{1}e^{at}, Q_{2}=q_{2}e^{at}[/latex] which then in matrix form gives:
[latex]\begin{pmatrix} Ra+\frac{1}{C} & -\frac{1}{C} \\Ra & La^2\end{pmatrix}\begin{pmatrix}q_{1}\\q_{2}\end{pmatrix}=\begin{pmatri x}0\\0\end{pmatrix}[/latex]

Setting the determinant of the matrix equal to zero, leads to a cubic. There is going to be a constant term corrisponding to the solution a=0, and the remaining quadratic gives a repeated root (using the values of L,C,R given): [latex]a=-\frac{1}{2RC}[/latex]. For this value of a we get a solution:
[latex]\begin{pmatrix}Q_{1}\\Q_{2}\end{pmatrix}=A \begin{pmatrix}2\\1\end{pmatrix} exp(\frac{-t}{2RC})[/latex]
where A is an arbitary constant.

But now I need to find another solution. I think it'll be of the form [latex](Bt+C)exp(\frac{-t}{2RC})[/latex] where B is related to A but I'm not sure how to go about finding it? Most of the results google gives seem to be for a slightly different style of question where there is only one first order derivative in each equation.

Thanks in advance.

HallsofIvy

I would have done it in a slightly different way- not using matrices. First integrate the second equation to get
[tex]RQ_1+ L\frac{dQ_2}{dt}= C_1[/tex]
so
[tex]Q_1= \frac{C_1}{R}-\frac{L}{R} \frac{Q_2}{dt}[/tex]
Also from the second equation,
[tex]R\frac{dQ_1}{dt}= -L\frac{d^2Q_2}{dt^2}[/tex]
so we can replace both instances of Q₁ in the first equation:
[tex]-L\frac{d^2Q_2}{dt^2}+\frac{1}{C}(\frac{C_1}{R}- \frac{L}{R}\frac{dQ_2}{dt}-Q2)= 0[/tex]
or
[tex]-L\frac{d^2Q_2}{dt^2}-\frac{L}{CR}\frac{dQ_2}{dt}- \frac{1}{C}Q_2= -\frac{C_1}{CR}[/itex]
a second order non-homogeneous equation with constant coefficients.
Since that is a second order differential equation, it will have two independent solutions involving two new constants. That's as it should be since you orginally had first and second order equations to solve. After you have solved for Q₂(t) you can immediately get Q₁ from
[tex]Q_1= \frac{C_1}{R}-\frac{L}{R} \frac{Q_2}{dt}[/tex]

plmokn

Thank you.