Mètode de variació dels paràmetres

De Viquipèdia

En matemàtiques, la variació dels paràmetres és una tècnica usada per resoldre certes equacions diferencials ordinàries de segon ordre no homogènies. La variació de paràmetres no és una tècnica molt comú en el camp de les matemàtiques pures, però és una eina útil en enginyeria.

[edita] Tècnica

Donada una equació diferencial de la forma

$u''+p(x)u'+q(x)u=f(x)\,$

es defineix l'operador lineal

$L=D^2+p(x)D+q(x)\,$

on D representa l'operador diferencial. S'ha de resoldre, doncs, l'equació $L u (x) = f (x)$ per $u (x)$ , on $L$ i $f (x)$ són conegudes.

Suposant que es tenen dues solucions linealment independents per l'equació diferencial donada, u₁ i u₂. Sigui W el Wronskià d'aquestes dues funcions, i W sigui diferent de zero (les solucions són linealment independents).

Es busca la solució general a l'equació diferencial $u G (x)$ que serà de la forma

$u_G(x)=A(x)u_1(x)+B(x)u_2(x).\,$

Aquí, $A (x)$ i $B (x)$ són desconegudes, i $u 1 (x)$ i $u 2 (x)$ són les solucions de l'equació homogènia. Es pot observar que si $A (x)$ i $B (x)$ són constants, llavors $L u G (x) = 0$ . És desitjable que A=A(x) i B=B(x) siguin de la forma

$A'(x)u_1(x)+B'(x)u_2(x)=0.\,$

Ara,

$u_G'(x)=(A(x)u_1(x)+B(x)u_2(x))'=(A(x)u_1(x))'+(B(x)u_2(x))'\,$

$=A'(x)u_1(x)+A(x)u_1'(x)+B'(x)u_2(x)+B(x)u_2'(x)\,$

$=A'(x)u_1(x)+B'(x)u_2(x)+A(x)u_1'(x)+B(x)u_2'(x)\,$

i com que es requereix la condició de sobre, llavors es té que

$u_G'(x)=A(x)u_1'(x)+B(x)u_2'(x)\,$

Derivant un altre cop (i ometent passos intermitjos)

$u_G''(x)=A(x)u_1''(x)+B(x)u_2''(x)+A'(x)u_1'(x)+B'(x)u_2'(x)\,$

Ara es pot escriure l'acció de L sobre u_{G _{com a}}

_{Com que u₁ i u₂ són solucions, llavors}

$Lu_G=A'(x)u_1'(x)+B'(x)u_2'(x)\,$

Es té el sistema d'equacions

$\begin{pmatrix} u_1(x) & u_2(x) \\ u_1'(x) & u_2'(x) \end{pmatrix} \begin{pmatrix} A'(x) \\ B'(x)\end{pmatrix} = \begin{pmatrix} 0\\ f\end{pmatrix}$

Desenvolupant,

$\begin{pmatrix} A'(x)u_1(x)+B'(x)u_2(x)\\ A'(x)u_1'(x)+B'(x)u_2'(x)\end{pmatrix} = \begin{pmatrix} 0\\f\end{pmatrix}$

Per tant, el sistema de sobre determina les condicions

$A'(x)u_1(x)+B'(x)u_2(x)=0\,$

$A'(x)u_1'(x)+B'(x)u_2'(x)=Lu_G=f\,$

Es troben A(x) i B(x) d'aquestes condicions, per tant, donades

$\begin{pmatrix} u_1(x) & u_2(x) \\ u_1'(x) & u_2'(x) \end{pmatrix} \begin{pmatrix} A'(x) \\ B'(x)\end{pmatrix} = \begin{pmatrix} 0\\ f\end{pmatrix}$

es pot resoldre per (A′(x), B′(x))^T, i per tant

$\begin{pmatrix} A'(x) \\ B'(x)\end{pmatrix}= \begin{pmatrix} u_1(x) & u_2(x) \\ u_1'(x) & u_2'(x) \end{pmatrix}^{-1} \begin{pmatrix} 0\\ f\end{pmatrix}$

$={1\over W} \begin{pmatrix} u_2'(x) & -u_2(x) \\ -u_1'(x) & u_1'(x) \end{pmatrix} \begin{pmatrix} 0\\ f\end{pmatrix}$

Finalment,

$A'(x) = - {1\over W} u_2(x) f(x),\; B'(x) = {1 \over W} u_1(x)f(x)$

$A(x) = - \int {1\over W} u_2(x) f(x)\,dx,\; B(x) = \int {1 \over W} u_1(x)f(x)\,dx$

Mentre les equacions homogènies són relativament fàcils de resoldre, aquest mètode permet el càlcul dels coeficients de la solució general de l'equació particular, i per tant es pot determinar la solució general completa.

Cal tenir en compte que $A (x)$ i $B (x)$ es determinen per només una constant arbitràries addicional (la constant d'integració); es podrien esperar dues constants d'integració perquè l'equació original era de segon ordre. Afegir una constant a $A (x)$ o a $B (x)$ no canvia el valor de $L u G (x)$ perquè $L$ és lineal.