Un sistema d'equacions diferencials lineals és una EDO (equació diferencial ordinària) del tipus: $$$x'(t)=A(t)\cdot x+ b(t)$$$ on, $$A(t)$$ és una matriu, $$n \times n$$, de funcions en la variable $$t$$, $$b (t)$$ és un vector de dimensió $$n$$ de funcions en la variable $$t$$, i $$x$$ és un vector de mida $$n$$ que és la funció que volem trobar.
Un exemple de sistema d'EDO's lineal seria: $$$\begin{pmatrix}x \\ y \end{pmatrix}'= \begin{pmatrix} 0 & \ln t \\ -e^t & 3 \cos t \end{pmatrix} \cdot \begin{pmatrix} x \\ y \end{pmatrix}+\begin{pmatrix} e^t \\ 3e^t \end{pmatrix}$$$
Un sistema lineal de dimensió $$n$$ té $$n$$ solucions linealment independents i resoldre el sistema significa trobar-les totes. Tota solució és una combinació lineal d'aquestes $$n$$ solucions. Així, donat un vector de condicions inicials ($$n$$), determinarem les $$n$$ constants trobant una única solució.
Quan resolem un sistema lineal posarem els $$n$$ vectors solució (linealment independents) a les columnes d'una matriu, l'anomenada matriu fonamental del sistema (de mida $$n \times n$$). Per tant, entendrem per resoldre el sistema trobar una matriu fonamental. Multiplicant aquesta matriu per un vector de constants arbitràries tindrem la solució general.
Una propietat important de les matrius fonamentals és que, si multipliquem una matriu fonamental per una matriu constant amb determinant diferent de zero, el resultat és una altra matriu fonamental (és important que la matriu constant es multipliqui per la dreta, si no, no és cert).
Per resoldre aquest tipus d'equacions no hi ha mètodes explícits (només en dimensió $$1$$). Tot i això, hi ha alguns casos particulars que sí que sabrem resoldre: Sistemes d'EDO's homogenis a coeficients constants, sistemes d'EDO's a coeficients constants no homogenis, i sistemes triangulars d'equacions diferencials.
En aquest tema ens centrarem en Sistemes d'EDO's a coeficients constants no homogenis.
Suposem que tenim el següent sistema a resoldre:$$$x'=A \cdot x+b(x)$$$
Un exemple seria: $$$\begin{pmatrix} x \\ y \end{pmatrix}'=\begin{pmatrix} 0 & 1 \\ -1 & 2 \end{pmatrix} \cdot \begin{pmatrix} x \\ y \end{pmatrix} +\begin{pmatrix} e^t \\ 3e^t\end{pmatrix}$$$
Per resoldre-ho, seguirem els mateixos passos que en una equació lineal, és a dir, trobarem la solució general de la part homogènia i després trobarem una solució particular del sistema no homogeni.
Considerem el problema homogeni: $$$x'=A \cdot x$$$ Seguint els passos descrits anteriorment, trobem la matriu principal del sistema $$\phi_h(t)$$.
Sabem que, aleshores, la solució general de la part homogènia és: $$$x_h(t)=\phi_h(t) \cdot C$$$ on $$C$$ és un vector de constants.
Ara busquem una solució particular de la forma $$$x_p(t)=\phi_h(t)\cdot u(t)$$$ Imposem que sigui solució: $$$\left .\begin{array}{r} x_p'=A\cdot x_p+b(t)=A\cdot \phi_h(t)\cdot u(t)+b(t) \\ x_p'= \phi_h(t) \cdot u(t)+\phi_h(t) \cdot u'(t)=A \cdot \phi_h(t)\cdot u(t)+\phi_h(t)\cdot u'(t) \end{array} \right\} $$$ $$$\Rightarrow \phi_h(t) \cdot u'(t)=b(t) \Rightarrow u(t)=\int \Big(\phi_h(t)\Big) ^{-1} \cdot b(t) \cdot dt $$$
Així, la solució general és: $$$x(t)=x_h(t)+x_p(t)=\phi_h(t) \cdot C+ \phi_h(t) \cdot \int \Big(\phi_h(t)\Big) ^{-1} \cdot b(t) \cdot dt$$$
Reprenem l'exemple inicial.
Tenim, segons la nostra notació que: $$$A=\begin{pmatrix} 0 & 1 \\ -1 & 2 \end{pmatrix} \ b(t)=\begin{pmatrix} e^t \\ 3e^t \end{pmatrix}$$$ Sabem com resoldre la part homogènia d'aquest problema; una matriu fonamental era: $$$\phi_h(t)=e^t\begin{pmatrix} t & 1 \\ t+1 & 1 \end{pmatrix}$$$ Ara busquem una solució particular.
La teoria ens diu que busquem una solució del tipus $$$\begin{pmatrix} x_p(t) \\ y_p(t) \end{pmatrix}= \phi_h (t) \cdot \begin{pmatrix} u_1(t) \\ u_2(t) \end{pmatrix}$$$ sabent que $$$\begin{pmatrix} u_1(t) \\ u_2(t) \end{pmatrix}=\int \phi_h^{-1} (t) \cdot b(t) \cdot dt$$$ Calculant primer $$$\displaystyle \phi_h^{-1}(t)=\frac{1}{-e^{2t}}\begin{pmatrix} e^t & -e^t \\ -e^t-t\cdot e^t & t \cdot e^t\end{pmatrix}=\begin{pmatrix} -e^{-t} & e^{-t} \\ e^{-t}-t\cdot e^t & -te^{-t}\end{pmatrix} = e^{-t} \begin{pmatrix} -1 & 1 \\ 1+t & -t \end{pmatrix}$$$ tenim: $$$\displaystyle \phi_h^{-1}(t)\cdot b(t)= e^{-t} \begin{pmatrix} -1 & 1 \\ 1+t & -t \end{pmatrix}\cdot e^t \begin{pmatrix}1 & 3\end{pmatrix}=\begin{pmatrix}2 \\ 1+t-3t\end{pmatrix}=\begin{pmatrix}2 \\ 1-2t\end{pmatrix}$$$ Per tant: $$$\displaystyle \begin{pmatrix}u_1(t) \\ u_2(t) \end{pmatrix}= \begin{pmatrix}\int 2 \cdot dt \\ \int (1-2t)\cdot dt \end{pmatrix}=\begin{pmatrix}2t \\ t-t^2\end{pmatrix}$$$ i, finalment: $$$\begin{pmatrix}x_p(t) \\ y_p(t) \end{pmatrix}=\phi_h(t) \cdot \begin{pmatrix}u_1(t) \\ u_2(t) \end{pmatrix}=e^t\begin{pmatrix} t & 1 \\ t+1 & 1 \end{pmatrix} \cdot \begin{pmatrix}2t \\ t-t^2\end{pmatrix}=$$$ $$$= e^t \cdot \begin{pmatrix}2t^2+t-t^2 \\ 2t^2+2t+t-t^2 \end{pmatrix}=e^t \begin{pmatrix}t^2+t \\ t^2+3t \end{pmatrix} $$$
D'aquesta manera, la solució general del sistema serà la suma de la solució de la homogènia i de la solució particular. $$$ \begin{pmatrix}y(t) \\ x(t) \end{pmatrix}= e^t\begin{pmatrix}t & 1 \\ t+1 & 1 \end{pmatrix}\cdot \begin{pmatrix} C_1 \\ C_2 \end{pmatrix}+e^t\begin{pmatrix}t^2+t \\ t^2+3t\end{pmatrix}$$$