czwartek, 14 lutego 2013

Rozwiązywanie równań różniczkowych metodą czynnika całkującego.



Rozważmy równanie różniczkowe w postaci:

\[\frac{{dy}}{{dx}} + p(x)y = q(x)\]
Pomnóżmy obie strony równania przez pewną funkcją J(x):
\[J(x)\frac{{dy}}{{dx}} + J(x)p(x)y = J(x)q(x)\]
Zauważmy, że lewa strona równania przypomina nieco formułę na pochodną iloczynu
\[\frac{d}{{dx}}(J(x)y) = J(x)\frac{{dy}}{{dx}} + y\frac{{dJ(x)}}{{dx}}\]


 
Porównując prawą stronę tego równania z lewą stroną wcześniejszego mamy:
\[J(x)\frac{{dy}}{{dx}} + J(x)p(x)y = J(x)\frac{{dy}}{{dx}} + y\frac{{dJ(x)}}{{dx}}\]
Po prostych przekształceniach mamy:
\[\frac{{dJ(x)}}{{dx}} = J(x)p(x)\]
\[\frac{{dJ}}{J} = p(x)dx\]
Całkując stronami mamy:
\[\int {\frac{{dJ}}{J}}  = \int {p(x)dx} \]
\[\ln J = \int {p(x)dx} \]
\[J(x) = {e^{\int {p(x)dx} }}\]
Funkcję J(x) nazywamy czynnikiem całkującym. Wyjściowe równanie przyjmie postać:
\[\frac{d}{{dx}}(J(x)y) = g(x)J(x)\]
Całkując stronami mamy:
\[J(x)y = \int {g(x)J(x)dx} \]
Ostatecznie
\[y = \frac{1}{{J(x)}}\int {g(x){e^{\int {p(x)dx} }}} dx\]


Przykład:
Rozwiąż równanie
\[\frac{{dy}}{{dx}} + xy = x\]
W naszym przypadku p(x)=x i q(x)=x. Obliczmy czynnik całkujący:
\[J(x) = {e^{\int {xdx} }} = {e^{\frac{1}{2}{x^2}}}\]
Pomnożmy stronami nasze równanie wyjściowe przez czynnik całkujący:
\[{e^{\frac{1}{2}{x^2}}}\frac{{dy}}{{dx}} + xy{e^{\frac{1}{2}{x^2}}} = x{e^{\frac{1}{2}{x^2}}}\]
\[\frac{d}{{dx}}\left( {{e^{\frac{1}{2}{x^2}}}y} \right) = x{e^{\frac{1}{2}{x^2}}}\]
\[{e^{\frac{1}{2}{x^2}}}y = \int {x{e^{\frac{1}{2}{x^2}}}dx} \]
\[\int {x{e^{\frac{1}{2}{x^2}}}dx} = \left\{ {\begin{array}{*{20}{c}}{u = \frac{1}{2}{x^2}}\\{du = xdx}\end{array}} \right\} = \int {{e^u}du = {e^u} + C = {e^{\frac{1}{2}{x^2}}}} + C\]
\[{e^{\frac{1}{2}{x^2}}}y = {e^{\frac{1}{2}{x^2}}} + C\]
\[y = C{e^{ - \frac{1}{2}{x^2}}} + 1\]

Brak komentarzy:

Prześlij komentarz