Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Дифференциальные уравнения.docx
Скачиваний:
2
Добавлен:
01.07.2025
Размер:
122 Кб
Скачать
☆

3. Линейные уравнения.

Линейное уравнение - это уравнение вида http://www.cyberforum.ru/cgi-bin/latex.cgi?\color{red}y'+p(x)y=q(x) . Это уравнение можно решить следующими методами:

1. Метод Бернулли.

Суть метода состоит в разложении искомой функции на произведение двух других - http://www.cyberforum.ru/cgi-bin/latex.cgi?y=u*v , и "подгонкой" их под нужный вид. В этом случае линейное уравнение будет сведено к системе из двух уравнений с разделяющимися переменными.

http://www.cyberforum.ru/cgi-bin/latex.cgi?y'=vu'+uv',\; vu'+uv'+p(x)uv=q(x) \Leftrightarrow v(u'+p(x)u)+uv'=q(x) . Теперь нужно получить функцию http://www.cyberforum.ru/cgi-bin/latex.cgi?u(x) так, чтобы выражение http://www.cyberforum.ru/cgi-bin/latex.cgi?u'+p(x)u равнялось нулю. Для этого решим соответствующее уравнение с разделяющимися переменными http://www.cyberforum.ru/cgi-bin/latex.cgi?u'+p(x)u=0 . Его решением будет http://www.cyberforum.ru/cgi-bin/latex.cgi?u=e^{-\int p(x)\operatorname{d}x} . Также стоит заметить, что функция u не зависит от произвольной постоянной. Теперь функцию http://www.cyberforum.ru/cgi-bin/latex.cgi?v можно найти. http://www.cyberforum.ru/cgi-bin/latex.cgi?v=\int e^{\int p(x)dx}q(x)\operatorname{d}x . Система уравнений, о которой я писал выше имеет вид:

http://www.cyberforum.ru/cgi-bin/latex.cgi?\left\{\begin{matrix}u^'+p(x)u=0\\ v^'=\frac{q(x)}{u}\end{matrix}\right.

Решение исходного уравнения получим просто перемножив http://www.cyberforum.ru/cgi-bin/latex.cgi?u на http://www.cyberforum.ru/cgi-bin/latex.cgi?v:\; y=u(x)v(x,C)

***

2. Метод вариации произвольной постоянной (метод Лагранжа).

Метод состоит в следующем: решим вместо нашего уравнения его вариант http://www.cyberforum.ru/cgi-bin/latex.cgi?y'+p(x)y=0 . Решение записывается в виде: http://www.cyberforum.ru/cgi-bin/latex.cgi?y=Ce^{-\int p(x)\operatorname{d}x} . Положим теперь, что http://www.cyberforum.ru/cgi-bin/latex.cgi?C=C(x) и подставим решение уравнения с нулевой правой частью в исходное. Получим уравнение для C(x):

http://www.cyberforum.ru/cgi-bin/latex.cgi?C'(x)e^{-\int p(x)\operatorname{d}x}-C(x)p(x)e^{-\int p(x)\operatorname{d}x}+C(x)p(x)e^{-\int p(x)\operatorname{d}x}=q(x)

Откуда http://www.cyberforum.ru/cgi-bin/latex.cgi?C'(x)=q(x)e^{\int p(x)\operatorname{d}x}\Leftrightarrow C(x)=\int q(x)e^{\int p(x)\operatorname{d}x}\operatorname{d}x , решение исходного уравнения http://www.cyberforum.ru/cgi-bin/latex.cgi?y=C(x)e^{-\int p(x)\operatorname{d}x} . Сравните его с полученным по методу Бернулли.

***

Пример (метод Бернулли):

http://www.cyberforum.ru/cgi-bin/latex.cgi?y^'+xy=x\;\;y=uv\;\;y^'=vu^'+uv^'\;\;vu^'+uv^'+xuv=x\\v^'+xv=0\;\;\frac{dv}{v}=-xdx\;\;v=e^{-\frac{x^2}{2}}\;\;u^'e^{-\frac{x^2}{2}}=x\Rightarrow u=e^{\frac{x^2}{2}}+C\Rightarrow y=Ce^{-\frac{x^2}{2}}+1

4. Уравнение Бернулли.

Уравнение Бернулли - это уравнение вида http://www.cyberforum.ru/cgi-bin/latex.cgi?\color{red}y^'+p(x)y=q(x)y^s,\;s\in \mathbb{R},\;s\neq 1 .

1. Сведение к линейному уравнению.

Это уравнение сводится к линейному следующим образом:

Домножим исходное уравнение на http://www.cyberforum.ru/cgi-bin/latex.cgi?y^{-s} , получим http://www.cyberforum.ru/cgi-bin/latex.cgi?y^'y^{-s}+p(x)y^{1-s}=q(x) . Сделаем замену http://www.cyberforum.ru/cgi-bin/latex.cgi?u=y^{1-s},\;u^'=(1-s)y^'y^{-s}\Leftrightarrow y^'y^{-sp}=\frac{u^'}{1-s} . В результате получим уравнение вида http://www.cyberforum.ru/cgi-bin/latex.cgi?\frac{u^'}{1-s}+p(x)u=q(x) , которое после домножения на 1-s станет линейным.

Случай, когда s=1 не рассматриваю, т.к. в этом случае уравнение Бернулли - уравнение с разделяющимися переменными: y'+(p(x)-q(x))y=0.

***

2. Метод Бернулли.

Помимо сведения к линейному уравнению, решить уравнение Бернулли можно методом Бернулли.

В этом случае получим http://www.cyberforum.ru/cgi-bin/latex.cgi?u'+p(x)u=0, \;u=e^{-\int p(x)\operatorname{d}x} и http://www.cyberforum.ru/cgi-bin/latex.cgi?v'u = u^sv^sq(x)\Rightarrow \frac{\operatorname{d}v}{v^s}=e^{-(s-1)\int p(x)\operatorname{d}x}q(x)\operatorname{d}x

***

Пример (сведение уравнения Бернулли к уравнению линейному):

http://www.cyberforum.ru/cgi-bin/latex.cgi?y^'+p(x)y=q(x)y^2|:y^2\\\frac{y^'}{y^2}+\frac{p(x)}{y}=q(x)\\u=\frac{1}{y}\;\;u^'=-\frac{y^'}{y^2}\;\Rightarrow \;-u^'+p(x)u=q(x)\\u^'-p(x)u=-q(x)

5. Уравнения в полных дифференциалах.

Уравнение вида P(x,y)dx+Q(x,y)dy=0 называется уравнением в полных дифференциалах, если существует такая функция U(x,y), такая что dU=P(x,y)dx+Q(x,y)dy. Для того, чтобы такая функция существовала необходимо и достаточно, чтобы http://www.cyberforum.ru/cgi-bin/latex.cgi?\frac{\partial P}{\partial y}=\frac{\partial Q}{\partial x} (Этот факт доказывается в курсе математического анализа). Поясню, что это значит:

Если такая функция U существует, то её полный дифференциал имеет вид http://www.cyberforum.ru/cgi-bin/latex.cgi?dU=\frac{\partial U}{\partial x}dx+\frac{\partial U}{\partial y}dy . Кроме того, http://www.cyberforum.ru/cgi-bin/latex.cgi?\frac{{\partial}^{2} U}{\partial x\partial y}=\frac{{\partial}^{2} U}{\partial y\partial x} . http://www.cyberforum.ru/cgi-bin/latex.cgi?P=\frac{\partial U}{\partial x} , http://www.cyberforum.ru/cgi-bin/latex.cgi?Q=\frac{\partial U}{\partial y} . Поэтому нам нужно соблюдение условия http://www.cyberforum.ru/cgi-bin/latex.cgi?\frac{\partial P}{\partial y}=\frac{{\partial}^{2} U}{\partial x\partial y}=\frac{{\partial}^{2} U}{\partial y\partial x}=\frac{\partial Q}{\partial x} .

Решением уравнения в полных дифференциалах будет следующий интеграл: http://www.cyberforum.ru/cgi-bin/latex.cgi?\int_{(x_0;y_0)}^{(x;y)}P(x;y)dx+Q(x;y)dy=0

Ввиду того, что dU=0, U=const. По аналогии с определёнными интегралами: http://www.cyberforum.ru/cgi-bin/latex.cgi?\int dC=C,\;\int_{a}^{b}dC=\int_{a}^{b}0dx=0 .

Интеграл от U берётся по любой кривой, соединяющей точки http://www.cyberforum.ru/cgi-bin/latex.cgi?(x_0;y_0) и http://www.cyberforum.ru/cgi-bin/latex.cgi?(x;y) . По этому удобнее взять ломаную ABC, где AB - отрезок, соединяющий http://www.cyberforum.ru/cgi-bin/latex.cgi?(x_0;y_0)\;c\;(x;y_0) , BC - отрезок, соединяющий http://www.cyberforum.ru/cgi-bin/latex.cgi?(x;y_0)\;c\;(x;y) . В этом случае интеграл будет выглядеть наиболее просто:http://www.cyberforum.ru/cgi-bin/latex.cgi?\int_{x_0}^{x}P(\xi ;y_0)d\xi +\int_{y_0}^{y}Q(x;\eta )d\eta

***

Пример:

http://www.cyberforum.ru/cgi-bin/latex.cgi?(ye^{xy}-sin(x))dx+xe^{xy}dy=0\\P(x;y)=ye^{xy}-sin(x)\;\;Q(x;y)=xe^{xy}\\\frac{\partial P}{\partial y}=e^{xy}+xye^{xy}\;\;\frac{\partial Q}{\partial x}=e^{xy}+xye^{xy}\\U=\int_{x_0}^{x}(y_0e^{\varepsilon y_0}-sin(\varepsilon ))d\varepsilon+\int_{y_0}^{y}xe^{x\eta }d\eta \\U=e^{xy_0}+cos(x)-e^{x_0y_0}-cos(x_0)+e^{xy}-e^{xy_0}\\U=e^{xy}+cos(x)-e^{x_0y_0}-cos(x_0)\\U=0\;\;e^{xy}+cos(x)+C=0\;\;C=-e^{x_0y_0}-cos(x_0)

6. Уравнения, неразрешенные относительно производной.

Формально, в этот класс можно отнести все уравнения первого порядка, если перенести все слагаемые из левой (правой) части в правую (левую) и получить уравнение вида F(x,y,y')=0. Но в этот класс я буду относить только уравнения, в которых нельзя производную выделить алгебраически, т.е. записать y'=f(x,y). Например: y=x*cos(y')+y'.

Для решения уравнений этого класса часто используется следующий метод:

Делаем замену y'=p и дифференцируем уравнение по х. Во многих случаях (хотя и далеко не всех)

эти действия позволяют свести решение нашего уравнения к уравнениям вышеописанных классов.

Так, например, можно решить уравнение, которое я дал в качестве примера выше:

p=cos(p)-xp'sin(p)+p', которое является линейным относительно x(p): x'(p-cos(p))=1-x*sin(p).

Естественно, данная замена (y'=p) не единственная. Можно применить и такую: cos(y')=p, но решать получающееся в этом случае уравнение сложнее (но в каком-нибудь другом уравнении именно эта замена может быть лучше).

***

Пример:

Рассмотрим уравнение вида y=x*f(y')+g(y'), которое называется уравнением Лагранжа. Его можно решить методом, который я кратко изложил в I.7). Сделаем замену y'=p и продифференцируем по х:

y'=f(p)+x*p'f'(p)+p'g'(p) <=> p=f(p)+x*p'f'(p)+p'g'(p) <=> p-f(p)=p'(xf'(p)+g'(p))

http://www.cyberforum.ru/cgi-bin/latex.cgi?p-f(p)=p^'(xf^'(p)+g^'(p))|*\frac{x^'}{p-f(p)}\\x^'=\frac{xf^'(p)+g^'(p)}{p-f(p)}\;\;\;\;\;x^'-x\frac{f^'(p)}{p-f(p)}=\frac{g^'(p)}{p-f(p)}

Последнее является линейным относительно x(p).