3.6. Примеры

№1. Найти общий интеграл уравнения 6e^x cos²y dx + (1 – 2e^x) ctg y dy=0.

Решение. Разделим переменные в данном уравнении, поделив обе его части на выражение cos²y (1 – 2e^x):

.

Интегрируя обе части уравнения, имеем

–3 ln|1–2e^x|+ln|tg y| = ln |C|, C

(поскольку C – произвольная постоянная, то для удобства дальнейших преобразований мы заменили С на ln|C|). Отсюда

или .

Получили общий интеграл данного уравнения. При делении на

cos²y (1–2e^x) мы могли потерять решения , k – целое число, x= –ln2, но они содержатся в общем интеграле, если подставить значение С = 0.

№2. Найти частное решение уравнения если у=–1 при х=1.

Решение. Перепишем уравнение в виде

х²dy = (xy + y²)dx (*)

и воспользуемся подстановкой у=uх. Тогда dy= udх+ хdu. Подставив выражение у и dy в уравнение (*), имеем

х²(udх + хdu)=(х^.uх+u²х²)dх;

x²(udx+xdu)=x²(u+u²)dx;

udx+xdu= udx+ u²dx; т.е.

xdu= u²dx.

Получили уравнение с разделяющимися переменными. Разделив переменные и интегрируя, получим

.

Так как u=, то .

Используя начальные условия х=1, у= –1, имеем 1=ln1+C, откуда С=1. Следовательно,

,

Отсюда получаем искомое частное решение данного уравнения .

№3. Решить уравнение .

Решение. Здесь P(x)=–ctgx, Q(x)=sinx. Решим уравнение двумя методами.

I. Метод Лагранжа

Найдем сначала общее решение соответствующего ЛОДУ , т.е. . Предположим, что (у=0 – решение данного уравнения), разделяя переменные и интегрируя, получим

, ,

Отсюда .

Общее решение ЛНДУ ищем в виде .

Найдем .

Подставим у и в исходное уравнение:

или .

Получили или . Тогда .

Следовательно, общее решение исходного уравнения есть y=(x+C) sin x.

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

Пусть. Тогда и уравнение принимает вид

,

или

.

Подберем функцию u(x) так, чтобы выражение в скобках равнялось нулю, т.е. решим первое дифференциальное уравнение с разделяющимися переменными ,

.

Откуда u=С₁sin x.

Пусть С₁=1, u=sin x.

, отсюда , т.е. .

Итак, y=(x+C)·sin x, есть общее решение данного ЛНДУ.

№4. Решить методом Рунге – Кутта дифференциальное уравнение при начальном условии у(0) = 1 на отрезке [0; 0,5] с шагом 0,1.

Решение.

Для i = 0 вычислим коэффициенты k_i:

Последующие вычисления приводить не будем, а результаты представим в виде таблицы.

i	xi	k		yi	yi
0	0	1	0,1000	0,1104	1
		2	0,1100
		3	0,1105
		4	0,1211
1	0,1	1	0,1210	0,1325	1,1104
		2	0,1321
		3	0,1326
		4	0,1443
2	0,2	1	0,1443	0,1569	1,2429
		2	0,1565
		3	0,1571
		4	0,1700
3	0,3	1	0,1700	0,1840	1,3998
		2	0,1835
		3	0,1842
		4	0,1984
4	0,4	1	0,1984	0,2138	1,5838
		2	0,2133
		3	0,2140
		4	0,2298
5	0,5				1,7976

№5. Решим предыдущий пример методом Эйлера.

Решение.

Применяем формулу

Производя аналогичные вычисления далее, получаем таблицу значений:

i	0	1	2	3	4	5
xi	0,0	0,1	0,2	0,3	0,4	0,5
yi	1	1,1	1,22	1,362	1,528	1,721

Для сравнения точности приведенных методов численного решение данного уравнения решим его аналитически и найдем точные значения функции у на заданном отрезке.

Уравнение является линейным неоднородным дифференциальным уравнением первого порядка. Решим соответствующее ему однородное уравнение.

Решение неоднородного уравнения имеет вид .

Общее решение: .

C учетом начального условия:

Частное решение: .

Для сравнения полученных результатов составим таблицу.

i	xi	yi
		Метод Эйлера	Метод Рунге-Кутта	Точное значение
0	0	1	1	1
1	0,1	1,1	1,1104	1,1103
2	0,2	1,22	1,2429	1,2428
3	0,3	1,362	1,3998	1,3997
4	0,4	1,528	1,5838	1,5837
5	0,5	1,721	1,7976	1,7975

Как видно из полученных результатов метод Рунге – Кутта дает наиболее точный ответ. Точность достигает 0,0001. Кроме того, следует обратить внимание на то, ошибка (расхождение между точным и приближенным значениями) увеличивается с каждым шагом вычислений. Это обусловлено тем что, во-первых, полученное приближенное значение округляется на каждом шаге, а во-вторых – тем что, в качестве основы вычисления принимается значение, полученное на предыдущем шаге, т.е. приближенное значение. Таким образом, происходит накопление ошибки.

Это хорошо видно из таблицы. С каждым новым шагом приближенное значение все более отличается от точного.