Лабораторная работа №5

Численное решение дифференциальных уравнений

Цель работы. Изучение численных методов решения дифференциальных уравнений.

Задание. 1. Решить дифференциальное уравнение аналитически и численно указанными методами для двух значений шага интегрирования h=0.01; 0.001 . Результаты расчета вывести на экран и распечатать в виде таблицы (табл.8).

2. Построить графики функций y(x) (5 графиков).

Варианты уравнений. Варианты уравнений и методов их решения приведены в табл.9.

Таблица 8

D	Решения уравнения, у(x)
	аналит.	численное
х		метод 1			метод 2
		h=0.01	h=0.001	h=0.01		h=0.001

Таблица 9

Вар.	Вид уравнения	Метод	Вар.	Вид уравнения	Метод
1	у'=(xy2+x)/(y-x2y)	1,4	14	у'=cos(t)-y	3,5
2	у'=(1-2x)/y2	2,4	15	у'=exp(bx)-ay	1,4
3	у'=(1-x2)/xy	3,4	16	у'=-2y/(y2-6x)	2,4
4	у'=(y2-y)/x	1,5	17	у'=1/(2x-y2)	3,4
5	у'=(1+y)/(tg(x)	2,5	18	у'=sec(x)- y tg(x)	1,5
6	у'=exp(x)-1	3,5	19	у'=(exp(x)-y)/x	2,5
7	у'=y ln(y)/sin(x)	1,4	20	у'=1+y/(x(x+1))	3,5
8	у'=(1+y2)/(1+x2)	2,4	21	у'=(y+yx2-x2)/(x(1+x2))	1,4
9	у'=4x-2y	3,4	22	у'=cos(x-y)	2,4
10	у'=x exp(-x2)-2xy	1,5	23	у'=3x-2y+5	3,4
11	у'=2x-y	2,5	24	у'=sin(x)-y	1,5
12	у'=exp(-x)-2y	3,5	25	у'=exp(x)-y	2,5
13	у'=exp(-x)-2x	1,4	26	у'=exp(2x)-1	3,5

Примечание. Значение параметров a,b и начальные условия y|_x=xo=y_o выбрать cамостоятельно.

Математическое описание. Дифференциальным уравнением первого порядка называется уравнение вида F(x,y,у')=0 или у'=f(x,y). Функция y(x), при подстановке которой уравнение обращается в тождество, называется решением дифференциального уравнения.

Рассмотрим несколько численных методов решения дифференциальных уравнений первого порядка. Описание численных методов приводится для уравнения в виде у'=f(x,y). Расчетные зависимости для одного шага интегрирования имеют следующий вид.

1.Метод Эйлера.

у_i+1=у_i+hf(x_i,y_i),

x_i+1=x_i+h.

2.Модифицированный метод Эйлера (вариант 1).

у_i+1=у_i+hf(x_i+h/2,y_i+hf(x_i,y_i)/2),

x_i+1=x_i+h.

3.Модифицированный метод Эйлера (вариант 2).

у_i+1=у_i+(h/2)[f(x_i,y_i)+f(x_i,+h,y_i+hf(x_i,y_i))],

x_i+1=x_i+h.

4.Метод Рунге-Кутта третьего порядка.

у_i+1=у_i+(k₁+4k₂+k₃)/6,

k₁=hf(x_i,y_i),

k₂=hf(x_i+h/2,y_i+k₁/2),

k₃=hf(x_i+h,y_i+2k₂-k₁),

x_i+1=x_i+h.

5.Метод Рунге-Кутта четвертого порядка.

у_i+1=у_i+(k₁+2k₂+2k₃+k₄)/6,

k₁=hf(x_i,y_i),

k₂=hf(x_i+h/2,y_i+k₁/2),

k₃=hf(x_i+h/2,y_i+k₂/2),

k₄=hf(x_i+h,y_i+k₃),

x_i+1=x_i+h,

где у_i+1,у_i - значения искомой функции в точках x_i+1,x_i соответственно, индекс i показывает номер шага интегрирования, h - шаг интегрирования. Начальные условия при численном интегрировании учитываются на нулевом шаге: i=0, x=x₀, y=y₀ .

Содержание отчета:

1. Название, цель работы и задание.

2. Математическое описание, алгоритм (структограмма) и текст программы.

3. Результаты расчета, пять графиков зависимости y(Dx) и выводы по работе .