Mètode de variació dels paràmetres

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 comuna en el camp de les matemàtiques pures, però és una eina útil en enginyeria.

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ó $Lu(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 $Lu_{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 intermedis)

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

Lu_{G}=A(x)Lu_{1}(x)+B(x)Lu_{2}(x)+A'(x)u_{1}'(x)+B'(x)u_{2}'(x)\,

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 $Lu_{G}(x)$ perquè $L$ és lineal.