Расчет дифференциального уравнения первого второго и третьего порядка методом Эйлера
Міністерство освіти України
ДАЛПУ
Кафедра
автоматизації
технологічних
процесів і приладобудування
КУРСОВА РОБОТА
з курсу “Математичне моделювання на
ЕОМ”
на тему “Розв’язок диференціального
рівняння
виду апу(п)+ап-1у(п-1)+…+а1у1+а0у=кх при заданих
початкових
умовах з автоматичним вибором кроку
методом Ейлера”
Виконала
студентка групи БА-4-97
Богданова
Ольга Олександрівна
Холоденко
Вероніка Миколаївна
Перевірила Заргун Валентина Василівна
ПОСТАНОВКА ЗАДАЧИ
И МЕТОД РЕШЕНИЯ
Решить
дифференциальное уравнение у/=f(x,y) численным методом -
это значит для заданной последовательности аргументов х0, х1…,
хn и числа у0, не определяя функцию у=F(x), найти такие
значения у1, у2,…, уn, что уi=F(xi)(i=1,2,…,
n) и F(x0)=y0.
Таким образом, численные методы позволяют вместо
нахождения функции
У=F(x) получить таблицу значений этой функции для
заданной последовательности аргументов. Величина h=xk-xk-1
называется шагом интегрирования.
Метод
Эйлера относиться к численным методам, дающим решение в виде таблицы
приближенных значений искомой функции у(х). Он
является сравнительно грубым и применяется в основном для
ориентировочных расчетов. Однако идеи, положенные в основу метода Эйлера,
являются исходными для ряда других методов.
Рассмотрим дифференциальное уравнение первого порядка
y/=f(x,y)
(1)
с
начальным условием
x=x0,
y(x0)=y0 (2)
Требуется найти решение уравнения (1) на отрезке [а,b].
Разобьем отрезок [a,
b] на n равных
частей и получим последовательность х0, х1,
х2,…, хn, где xi=x0+ih (i=0,1,…, n), а h=(b-a)/n-шаг
интегрирования.
В методе Эйлера приближенные значения у(хi)»yi
вычисляются
последовательно по формулам уi+hf(xi,
yi) (i=0,1,2…).
При этом искомая интегральная кривая у=у(х),
проходящая через точку М0(х0, у0), заменяется
ломаной М0М1М2… с вершинами Мi(xi, yi) (i=0,1,2,…); каждое
звено МiMi+1
этой ломаной,
называемой ломаной Эйлера, имеет направление, совпадающее с направлением
той интегральной кривой уравнения (1), которая проходит через точку Мi.
Если правая часть уравнения (1) в
некотором прямоугольнике R£a, удовлетворяет условиям:
|f(x, y1)- f(x, y2)| £ N|y1-y2|
(N=const),
|df/dx|=|df/dx+f(df/dy)|
£ M (M=const),
то имеет место следующая оценка погрешности:
|y(xn)-yn|
£ hM/2N[(1+hN)n-1],
(3)
где у(хn)-значение точного решения
уравнения(1) при х=хn, а уn-
приближенное значение, полученное на n-ом шаге.
Формула (3) имеет в основном теоретическое
применение. На практике иногда оказывается более удобным двойной просчет: сначала расчет ведется с
шагом h, затем
шаг дробят и повторный расчет ведется с шагом h/2.
Погрешность более точного значения уn* оценивается формулой
|yn-y(xn)|»|yn*-yn|.
Метод Эйлера
легко распространяется на системы дифференциальных уравнений и на
дифференциальные уравнения высших порядков. Последние должны быть
предварительно приведены к системе дифференциальных уравнений первого порядка.
Модифицированный метод Эйлера более точен.
Рассмотрим дифференциальное уравнение (1) y/=f(x,y)
с начальным условием y(x0)=y0.
Разобьем наш участок интегрирования на n
равных частей. На малом участке [x0,x0+h]
у
интегральную
кривую заменим прямой
Nk/ y=y(x) линией. Получаем
точку Мк(хк,ук).
Мк Мк/
yk+1
yk
хк хк1/2 xk+h=xk1 х
Через Мк проводим касательную: у=ук=f(xk,yk)(x-xk).
Делим
отрезок (хк,хк1) пополам:
xNk/=xk+h/2=xk+1/2
yNk/=yk+f(xk,yk)h/2=yk+yk+1/2
Получаем
точку Nk/. В этой точке строим
следующую касательную:
y(xk+1/2)=f(xk+1/2,
yk+1/2)=αk
Из точки Мк проводим прямую с угловым
коэффициентом αк и определяем точку пересечения этой прямой с прямой
Хк1. Получаем точку Мк/. В качестве ук+1
принимаем ординату точки Мк/. Тогда:
ук+1=ук+αкh
xk+1=xk+h
(4) αk=f(xk+h/2, yk+f(xk,Yk)h/2)
yk=yk-1+f(xk-1,yk-1)h
(4)-рекурентные формулы метода
Эйлера.
Сначала
вычисляют вспомогательные значения искомой функции ук+1/2
в точках хк+1/2, затем находят значение
правой части уравнения (1) в средней точке y/k+1/2=f(xk+1/2, yk+1/2) и определяют ук+1.
Для оценки погрешности в точке хк
проводят вычисления ук с шагом h, затем с шагом 2h и
берут 1/3 разницы этих значений:
| ук*-у(хк)|=1/3(yk*-yk),
где
у(х)-точное решение дифференциального уравнения.
Таким образом, методом Эйлера можно решать
уравнения любых порядков. Например, чтобы решить уравнение второго порядка y//=f(y/,y,x) c начальными условиями y/(x0)=y/0,
y(x0)=y0, выполняется замена:
y/=z
z/=f(x,y,z)
Тем
самым преобразуются начальные условия: y(x0)=y0, z(x0)=z0, z0=y/0.
РЕШЕНИЕ КОНТРОЛЬНОГО ПРИМЕРА
Приведем расчет дифференциального уравнения первого,
второго и третьего порядка методом
Эйлера
1. Пусть дано дифференциальное уравнение первого порядка:
y/=2x-y
Требуется
найти решение на отрезке [0,1] c шагом h=(1-0)/5=0,2
Начальные
условия: у0=1;
Пользуясь
рекурентными формулами (4), находим:
1). x1=0,2; х1/2=0,1; y(x1)=y(x0)+α0h; y(x1/2)=y(x0)+f(x0,y0)h/2;
3. Чтобы
решить уравнение третьего порядка
y///=2x-y-y/+y//
на
отрезке [0,1], с шагом h=0,2 и начальными условиями
y0//=1
y0/=1
y0=1
необходимо
сделать 3 замены: y/=a y0/=a0=1
y//=a/=b y0//=b0=1
b/=2x-y-a+b
1).x1=0,2; x1/2=0,1
Программа на Turbo Pascal
uses
crt,pram,kurs1_1;
var
yx,xy,l,v,p,ff,ay,by,x:array [0..10] of
real;
y,a,b:array[0..10,0..1] of real;
i,n,o:integer;
c,d,h,k:real;
label
lap1;
begin
screen1;
clrscr;
writeln('введите
наивысший порядок производной не больше трех ');
readln(n);
if
n=0 then begin
writeln('это
прямолинейная зависимость и решается без метода Эйлера');
goto
lap1;end;
writeln('введите
коэффициенты {a0,a1}');
for
i:=0 to n do
readln(l[i]);
if
(n=1) and (l[1]=0) or (n=2) and (l[2]=0) or (n=3) and (l[3]=0) then begin
writeln('деление на ноль');
goto lap1;
end;
writeln('введите
коэффициент при x');
readln(k);
writeln('введите
отрезок ');
readln(c,d);
o:=5;
h:=abs(d-c)/o;
writeln('шаг=',h:1:1);
writeln('задайте
начальные условия y(x)= ');
for
i:=0 to n-1 do
readln(v[i]);
if
n=3 then begin
yx[0]:=v[0];
ay[0]:=v[1];
by[0]:=v[2];
p[0]:=(k*c-l[0]*v[0]-l[1]*v[1]-l[2]*v[2])/l[3];
x[0]:=c;
gotoxy(32,1);
write(' ');
gotoxy(32,2);
write(' x y a b ');
gotoxy(32,3);
write(' ',c:7:7,' ',yx[0]:7:7,'
',ay[0]:7:7,' ',by[0]:7:7,' ');
for i:=0 to o-1 do begin
x[i]:=x[i]+h/2;
y[i,1]:=yx[i]+(h/2)*ay[i];
a[i,1]:=ay[i]+(h/2)*by[i];
b[i,1]:=by[i]+(h/2)*p[i];
p[i+1]:=(k*xy[i]-l[0]*yx[i+1]-l[1]*ay[i+1]-l[2]*by[i+1])/l[3];
end;
for i:=0 to o-1 do begin
gotoxy(32,4+i);
write(' ',xy[i]:7:7,'
',yx[i+1]:7:7,' ',ay[i+1]:7:7,' ',by[i+1]:7:7,' ');
end;
gotoxy(32,4+o);
write(' ');
end;
if
n=2 then begin
x[0]:=c;
yx[0]:=v[0];
ay[0]:=v[1];
p[0]:=(k*c-l[0]*yx[0]-l[1]*v[1])/l[2];
gotoxy(32,1);
write(' ');
gotoxy(32,2);
write(' x y a ');
gotoxy(32,3);
write(' ',c:7:7,'
',yx[0]:7:7,' ',ay[0]:7:7,' ');
for i:=0 to o-1 do begin
x[i]:=x[i]+h/2;
y[i,1]:=yx[i]+(h/2)*ay[i];
a[i,1]:=ay[i]+(h/2)*p[i];
p[i+1]:=(k*xy[i]-l[0]*yx[i+1]-l[1]*ay[i+1])/l[2];
end;
for i:=0 to o-1 do begin
gotoxy(32,4+i);
write(' ',xy[i]:7:7,'
',yx[i+1]:7:7,'
',ay[I+1]:7:7,' ');
end;
gotoxy(32,4+o);
write(' ');
end;
if n=1 then begin
x[0]:=c;
yx[0]:=v[0];
p[0]:=(k*x[0]-l[0]*yx[0])/l[1];
for i:=0 to o-1 do begin
x[i]:=x[i]+h/2;
y[i,1]:=yx[i]+(h/2)*p[i];
xy[i]:=x[i]+h/2;
ff[i]:=(k*x[i]-l[0]*y[i,1])/l[1];
yx[i+1]:=yx[i]+h*ff[i];
x[i+1]:=x[i]+h/2;
p[i+1]:=(k*xy[i]-l[0]*yx[i+1])/l[1];
end;
gotoxy(32,1);
write(' ');
gotoxy(32,2);
write(' x y ');
gotoxy(32,3);
write('
',c:7:7,' ',yx[0]:7:7,' ');
for i:=0 to o-1 do begin
gotoxy(32,4+i);
write(' ',xy[i]:7:7,' ',yx[i+1]:7:7,' ');
end;
gotoxy(32,o+4);
write(' ');
end;
lap1:readln;
pramo;
delay(10000);
clrscr;
end.
ЗАПУСК ПРОГРАММЫ НА ВЫПОЛНЕНИЕ
Программа находится в файле kursova1.pas, и имеет 2 модуля, в которых содержатся
заставки. Модули находятся в файлах pram.tpu и kurs1_1.tpu.
Для запуска файла kursova1.pas
в Turbo Pascal необходимо нажать F9. Появится первая заставка,
далее нажать enter и ввести все необходимые начальные условия: порядок
производной, коэффициенты при членах рада, отрезок и начальные значения у(х0).
На экране выводится шаг вычисления и таблица с ответами. После нажатия enter выводится
вторая заставка, после чего мы возвращаемся к тексту программы.
ОПИСАНИЕ ПРОГРАММЫ
1 – ввод данных, используемых
в программе
2 – использование метки,
очистка экрана, ввод требований, решение
дифференциального уравнения в зависимости от ввода начальных
условий
3 – присвоение начальных
условий для дифференциального уравнения
третьего порядка
4 – вывод таблицы со значениями
5 – ввод формул метода Эйлера
для уравнения третьего порядка
6 – присвоение начальных условий
для решения дифференциального
уравнения второго порядка
7 – вывод таблицы для
уравнения второго порядка
8 – формулы метода Эйлера для
уравнения второго порядка
9 – начальные условия для
дифференциального уравнения первого порядка
10 – формулы метода Эйлера для
решения уравнения первого порядка
11 – вывод таблицы
12 – обращение к метке,
задержка для просмотра результатов, очистка
экрана, конец программы.
Вывод уравнения эйлера из принципа максимума понтрягина для уравнения с запаздыванием. Обыкновенные дифференциальные уравнения в примерах с модификационным методом эйлера. Решение дифференциального уравнения второго порядка метод эйлера алгоритм. Реферат на тему Обыкновенные дифференциальное уравнение первого порядка. Об устойчивости обыкновенных дифференциальных уравнений и метода Эйлера. Модифицированный метод Эйлера блок схема алгоритма Москва Россия Москве. Решение дифференциального уравнения второго порядка методом Эйлера в. Решение дифференциального уравнения третьего порядка методом Эйлера. Пример решения дифференциальных уравнений го порядка методом эйлера. Примеры нахождения значения функции с использованием метода Эйлера. Решение дифференциального уравнения второго порядка методом Эйлера. Решение дифференциальных уравнений второго порядка методом эйлера. Решение дифференциальных уравнений первого порядка методом эйлера. Реферат решение системы дифференциальных уравнений методом Эйлера. Программа дифференциальное уравнение на отрезке методом эйлера.