Всякое уравнение, содержащее более одной переменной, связывает эти переменные. Даже когда мы задаём функцию вида y=f(x), мы пишем уравнение. Это уравнение уже разрешено относительно переменной y, но его можно решить и относительно переменной x.
Дифференциальными называются уравнения, содержащие производные одной переменной по другой.
Если в одном уравнении встречаются производные по разным переменным, например, если уравнение имеет вид F(x,y,z,∂z∂x,∂z∂y,∂2z∂x2,…)=0, то такие дифференциальные уравнения называются уравнениями в частных производных. Если же все производные в уравнении – по одной переменной, то дифференциальное уравнение называется обыкновенным дифференциальным уравнением.
Старший порядок производной, встречающейся в уравнении, называется порядком этого уравнения. Таким образом, обыкновенное дифференциальное уравнение порядка n имеет вид F(x,y,y′,y″,…y(n))=0. Заметим, что в обыкновенном дифференциальном уравнении может встречаться более двух переменных, но все производные, которые в нём есть – по одной.
Мы далее будем рассматривать обыкновенные дифференциальные уравнения первого порядка. Более того, из таких уравнений мы долго будем выбирать уравнения, разрешённые относительно производной – то есть, такие: y′=f(x,y) Общих алгоритмов решения дифференциальных уравнений, как и взятия интегралов, не существует. Свои методы есть только для уравнений отдельно взятых типов. Чем больше таких типов вы знаете – тем лучше вы умеете решать дифференциальные уравнения; хотя никакое сколь угодно глубокое освоение не может гарантированно спасти от фиаско.
Мы начнём с уравнений в частных производных, к которым сводятся многие другие. Выглядят они так: y′=f(x)g(y) Разделив обе части на g(y), проинтегрируем обе части по x. Напоминаю, что если две функции равны, их первообразные равны с точностью до постоянного слагаемого, а их неопределённые интегралы - просто равны. ∫y′g(y)dx=∫f(x)dx В левой части произведём замену y′dx=dy: ∫dyg(y)=∫f(x)dx Осталось взять два интеграла в левой и правой частях, и мы получим формулу вида H(y)=F(x)+C, которую можно решить относительно y и выразить таким образом y через x. Но даже если этого этого не делать, уравнение типа (2) будет уже решением дифференциального уравнения (1). Отметим, что устанавливаемая таким образом связь между y и x неоднозначна: она содержит произвольную константу C.
Пример: решим №53 из Филиппова: (x2−1)y′+2xy2=0, иначе y′=−2xy2(x2−1)=−2x(x2−1)y2. В таком виде уравнение явно имеет вид (1). −y′y2=2xx2−1, −∫dyy2=∫2xx2−1dx, 1y=ln|x2−1|+C. В задаче, однако, есть дополнительное задание: найти в этом семействе решений то, при котором y(0)=1 (т.е. проходящее через точку (0;1)) Тогда 11=ln|02−1|+C,C=1, откуда 1y=ln|x2−1|+1, y=1ln|x2−1|+1. Если бы дополнительного условия y(0)=1 в задаче не было, пришлось бы оставить в ответе всё семейство решений.
Задание: № 54 – 58 (задания, не решённые во время занятия, остаются на дом)
Первым из многочисленных типов уравнений, сводящихся к уравнениям с разделяющимися переменными, будет такой: y′=f(ax+by+c). Он безымянный, не разобран в Филиппове, но уравнения такого типа там попадаются. В таких уравнениях делается замена (и это – первая замена, которую мы делаем в дифференциальных уравнениях) ax+by+c=z. Тогда a+by′=z′, a+by′=a+bf(ax+by+c), z′=a+bf(z), z′a+bf(z)=1, ∫dza+bf(z)=x+C.
Например, №62: y′=cos(y−x) z=y−x, z′=y′−1, y′−1=cos(y−x)−1, z′=cos(z)−1, z′cos(z)−1=1, ∫dzcos(z)−1=−∫dz1−cos(z)=−12∫dzsin2(z2)=ctgz2=x+C. Теперь осталось выполнить обратную замену, и мы получим связь x и y: ctgy−x2=x+C. Задание: № 63, 65
Однородные уравнения имеют вид y′=f(yx) В них выполняется замена y=x⋅u(x), где u – новая искомая функция. При этом с уравнением происходит следующее: xu′+u=f(u), xu′=f(u)−u, u′f(u)−u=1x, ∫duf(u)−u=ln|x|+C, где нужно взять интеграл в левой части и заменить обратно u=yx.
Например, решим №103: (y2−2xy)dx+x2dy=0 Перенося и деля, получим, что y′=dydx=−y2−2xyx2=−y2x2+2yx. Заменим: xu′+u=−u2+2u, xu′=−u2+u, u′u−u2=1x, ln|x|+˜C=∫duu−u2=∫1−u+uu(1−u)du=∫(1u+11−u)du=ln|u|−ln|1−u|, ln|x|+˜C=ln|u1−u|=ln|yx1−yx|=ln|yx−y|, ln|yx−y1x|=˜C. Применим ex к обеим частям: |yx−y1x|=e˜C yx−y1x=±e˜C≡C Решим полученное относительно x: 1xy−1=xC, xy=1xC+1, yx=11xC+1=xC1+xC, и окончательно y=x2C1+xC. Задание: аналогичным способом решить №101, 102, 104-108.