Решение уравнений i рунге кутта

Метод Рунге — Кутты

Этот онлайн калькулятор реализует классический метод Рунге — Кутты (встречается также название метод Рунге — Кутта) четвертого порядка точности. Метод используется для решения дифференциальных уравнений первой степени с заданным начальным значением

Калькулятор ниже находит численное решение дифференциального уравнения первой степени методом Рунге-Кутты (иногда встречается название метод Рунге-Кутта, а в поисковиках бывает ищут «метод рунге кута», «метод рунги кутта» и даже «метод рунги кута»), который также известен как классический метод Рунге — Кутты (потому что есть на самом деле семейство методов Рунге-Кутты) или метод Рунге — Кутты четвертого порядка.

Для того, чтобы использовать калькулятор, вам надо привести дифференциальное уравнение к форме

и ввести правую часть уравнения f(x,y) в поле y’ калькулятора.

Также вам понадобится ввести начальное значение

и указать точку в которой вы хотите получить численное решение уравнения .

Последнее параметр калькулятора — размер шага с которым вычисляется следующее приближение по графику функции.

Описание метода можно найти под калькулятором.

Метод Рунге-Кутта решения диф. уравнений и их систем.

Метод позволяет решать системы обыкновенных дифференциальных уравнений (ОДУ) первого порядка следующего вида:

которые имеют решение:

где t — независимая переменная (например, время); X, Y и т.д. — искомые функции (зависимые от t переменные). Функции f, g и т.д. — заданы. Также предполагаются заданными и начальные условия, т.е. значения искомых функций в начальный момент.

Одно диф. уравнение — частный случай системы с одним элементом. Поэтому, далее речь пойдет для определенности о системе уравнений.

Метод может быть полезен и для решения диф. уравнений высшего (второго и т.д.) порядка, т.к. они могут быть представлены системой диф. уравнений первого порядка.

Метод Рунге-Кутта заключается в рекурентном применении следующих формул:

Реализация Метода Рунге-Кутта на Delphi может выглядеть так (привожу полностью модуль):

Модуль полностью работоспособен. Возвращаемое функцией Runge_Kutt значение — код ошибки. Вы можете дополнить список ошибок по своему усмотрению. Рассчитанные функции системы помещаются в массив Res. Чтобы не загромождать код, в модуле опущены проверки (типа блоков try). Рекомендую их добавить по своему усмотрению.

Ниже приводится описание функции Runge_Kutt и типов, использующихся в модуле.

• FunArray — вектор функций (правых частей уравнений системы);
• First, Last — начальная и конечная точки расчетного интервала;
• Steps — число шагов по расчетному интервалу;
• InitArray — вектор начальных значений
• var Res — матрица результатов включая независимую переменную.

В модуле описаны типы:

Функция возвращает коды ошибок:

• 0 — нет ошибок;
• 100 — число уравнений не равно числу начальных условий.

Решение содержится в переменной-матрице Res. Первый индекс матрицы относится к переменной (0 — независимая переменная, 1 — первая зависимая и т.д.), второй — к номеру расчетной точки (0 — начальная точка).

Рассмотрим один пример использования модуля. Создадим новое приложение и подключим к нему модуль. На форме приложения разместим кнопку Button1 и область текста Memo1. Поместим в приложение две функции и обработчик нажатия кнопки:

Нажатие кнопки приведет к расчету точек системы, которые будут выведены в текстовую область.

Модуль с примером и справкой можно скачать бесплатно по адресу RK.zip (ZIP, 15,3Kb) (русский вариант). Английский вариант (условно-бесплатный) можно скачать по адресу RK_Eng.zip (ZIP, 23.4Kb)

Ссылки

• http://sadovoya.narod.ru/RK.zip (русский вариант).
• http://sintreseng.narod.ru/RK_Eng.zip (английский, условно-бесплатный вариант)

Оставить комментарий

Комментарии

Скачала по Вашей ссылке русский вариант, изменила для своей системы диф. уравнений, но при запуске выдаёт ошибку :
Project Ex.exe raised exception class EOverflow with message ‘ Floating point overflow ‘
Помогите, пожалуйста .

Вот изменённый мною модуль:

unit Unit1;
interface
uses
SysUtils, Forms, StdCtrls, Controls, Classes, Dialogs, Math;
type
TForm1 = class(TForm)
Memo1: TMemo;
rk_But: TButton;
procedure rk_ButClick(Sender: TObject);
private
< Private declarations >
public
< Public declarations >
end;
var
Form1: TForm1;
pn,k,ro,Pzv: Extended;

implementation
uses rk_method, Windows;

<$R *.dfm>
procedure Syst (var t: TFloat; var X: TFloatVector;
var RP: TFloatVector);
const
fdr1=0.503;
fdr2=0.503;
fdr3=0.196;
W1=179.8928;
W2=3773.8568;
W3=2504.1203;
b1=55.9203;
b2=98.6;
b3=98.6;
Ls1=3.78;
Ls2=9;
Ls3=15.3;
Svidj2=1352.438;
Svidj3=1352.438;
my=0.62;
vk=30;
m=1.2;
L1=30.969;
L2=42.131;
delta1=0;

begin
pn:=2.5*Power(10,4);
k:=6*Power(10,-7);
ro:=8.5*Power(10,-7);
Pzv:=3.919*Power(10,7);

RP[0] := (1/(k*W1))*(my*fdr1*sqrt(2/ro)*sqrt(Abs(pn-X[0]))-my*fdr2*sqrt(2/ro)*sqrt(Abs(X[0]-X[1]))-(delta1*delta1*delta1*b1)/(12*ro*vk*Ls1)*X[0]); // dp1/dt
RP[1] := (1/(k*W2))*(my*fdr2*sqrt(2/ro)*sqrt(Abs(X[0]-X[1]))-my*fdr3*sqrt(2/ro)*sqrt(Abs(X[1]-X[2]))-(X[4]*X[4]*X[4]*b2)/(12*ro*vk*Ls2)*X[1]); // dp2/dt
RP[2] := (1/(k*W3))*(my*fdr3*sqrt(2/ro)*sqrt(Abs(X[1]-X[2]))-(X[6]*X[6]*X[6]*b3)/(12*ro*vk*Ls3)*X[2]); // dp3/dt;
RP[3] := (((Svidj2*X[1]*(L1+L2))/L1)-Pzv)*(2/m); // dv2/dt
RP[4] := X[3]; // d delta2/dt
RP[5] := (((Svidj3*X[2]*(L1+L2))/L2)-Pzv)*(2/m); // dv3/dt
RP[6] := X[5]; // d delta3/dt
end;

procedure TForm1.rk_ButClick(Sender: TObject);
var
I, t1, t2: Cardinal;
tOut, InitConds: TFloatVector;
XOuts: TFloatMatrix;
Points: Cardinal;
First, Last: TFloat;
StepsFact: Cardinal;
Count: Word;
begin
Memo1.Clear;
First := 0.0;
Last := 10.0;
Count:= 7;
Points:=10+1; //11 points for output
StepsFact:=1000000; //all steps inside function = 10*StepsFact

try
SetLength(InitConds, Count);
InitConds[0]:=0.0; //x0(0)=0
InitConds[1]:=0.0; //x1(0)=0
InitConds[2]:=0.0; //x2(0)=0
InitConds[3]:=0.0; //x3(0)=0
InitConds[4]:=0.0; //x4(0)=0
InitConds[5]:=0.0; //x5(0)=0
InitConds[6]:=0.0; //x6(0)=0

SetLength(tOut, Points);
SetLength(XOuts, Count, Points);
except
ShowMessage(‘Out of memory. ‘);
exit;
end;

Дифференциальные уравнения первого порядка. Понятие о методе Рунге—Кутта

Содержание:

\

По этой ссылке вы найдёте полный курс лекций по математике:

Обыкновенным дифференциальным уравнением называется уравнение вида связывающее независимую переменную х, искомую функцию у = у(х) и ее производные у'(х),у»(х)>. у(п)(ж) (наличие хотя бы одной производной обязательно). Здесь — заданная функция своих аргументов. Замечание. Обозначения зависимой и независимой переменных через х и у, используемые в приведенном определении, не являются жесткими; часто в качестве независимой удобно брать переменную t, иными буквами обозначают и зависимую переменную (см. ниже пример 2).

В обыкновенном дифференциальном уравнении искомая функция у = у(х) есть функция одной независимой переменной х. Если искомая функция есть функция двух (и более) независимых переменных, то имеем дифференциальное уравнение с частными производными. В этой и двух следующих главах мы будем рассматривать только обыкновенные дифференциальные уравнения.

Простейшим дифференциальным уравнением является уравнение вида ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ ПЕРВОГО ПОРЯДКА Эквивалентные дифференциальные уравнения Задача Коши еорема существования и единственности решения задачи Коши для уравнения Приближенные методы интегрирования уравнения Метод последовательных приближений Численные методы решения задачи Коши Метод Эйлера Понятие о методе Рунге—Кутта где f(x) — известная непрерывная на некотором интервале (а, 6) функция, а у = у(х) — искомая функция.

С таким уравнением мы уже встречались в интегральном исчислении, когда поданной функции f(x) требовалось найти ее первообразную F(x). Как известно, всякая функция, удовлетворяющая уравнению (2), имеет вид где F(x) — какая-нибудь первообразная для функции /(ж) на интервале (а, 6), а С — произвольная постоянная. Таким образом, искомая функция у = у(х) определяется из уравнения (2) неоднозначно. Порядком дифференциального уравнения называется порядок наивысшей производной, входящей в уравнение.

Например, — дифференциальное уравнение 1-го порядка; дифференциальное уравнение 2-го порядка; — дифференциальное уравнение пятого порядка. Решением дифференциального уравнения n-го порядка на интервале (а, Ь) называется всякая функция у = (р(х), имеющая на этом интервале производные до n-го порядка включительно и такая, что подстановка функции у = tp(x) и ее производных в дифференциальное уравнение обращает последнее в тождество по х на интервале (а, Ь).

Например, функция у = sin х является решением дифференциального уравнения второго порядка на интервале (-оо, +оо). В самом деле, . Подставив в данное уравнение найденные значения у и у», получим — Задача. Найти совпадающие решения двух дифференциальных уравнений (не решая самих уравнений): График решения дифференциального уравнения называется интегральной кривой этого уравнения. Процесс нахождения решения дифференциального уравнения называется интегрированием дифференциального уравнения.

К составлению и интегрированию дифференциальных уравнений приводят многочисленные задачи как самой математики, так и других наук (физики, химии, биологии и т. п.). Пример 1. Найти такую кривую, чтобы тангенс угла наклона касательной в каждой ее точке численно равнялся ординате точки касания. Пусть — уравнение искомой кривой.

Как известно, , значит, определяющее свойство кривой есть — дифференциальное уравнение первого порядка.

Нетрудно видеть, что функция _ть решение этого уравнения. Оно также имеет очевидное решение у = 0. Кроме того, решениями оудут функции где С — произвольная постоянная, так что уравнение имеет бесконечное множество решений. ^ Пример 2. Найти закон прямолинейного движения материальной точки, движущейся с постоянным ускорением а.

А Требуется найти формулу з = s(t). выражающую пройденный путь как функцию времени. По условию имеем — дифференциальное уравнение второго порядка. Последовательно находим: Произвольные постоянные можно определить, если положить В самом деле, полагая первом из соотношений (*). получаем , откуда — Из второго соотношения (*) при t = tq имеем откуда Подставляя найденные значения С\ и С? в выражение для функции *(f). приходим к известному закону движения материальной точки с постоянным ускорением: § 2. Эквивалентные дифференциальные уравнения.

Задача Коши Пусть имеем дифференциальное уравнение первого порядка Если в этом уравнении удается выразить производную у’ через х и у, то получаем уравнение разрешенное относительно производной. Здесь / — заданная функция своих аргументов. Наряду с уравнением (1) рассматривают эквивалентное ему дифференциальное уравнение или уравнение более общего вида получаемое из (Г) путем умножения на некоторую функцию известные функции своих аргументов).

Два дифференциальных уравнения называются эквивалентными в некоторой области D изменения величин х, у, у’, если всякое решение у(х) 6 D одного из этих уравнений является решением другого уравнения и наоборот. При преобразовании дифференциальных уравнений надо следить за тем, чтобы преобразованное уравнение было эквивалентным исходному. Если дифференциальное уравнение имеет решение, то, как правило, множество его решений оказывается бесконечным.

Впрочем, дифференциальное уравнение имеет только одно решение а уравнение Чтобы выделить определенное решение уравнения (1), надо задать начальное условие, которое заключается в том, что при некотором значении независимой переменной х заранее дано значение уо искомой функции у(х): Геометрически это означает, что задается точка Мо(х0, уо), через которую должна проходить искомая интегральная кривая.

Возможно вам будут полезны данные страницы:

Задачу отыскания решения у(х) уравнения (1), удовлетворяющего начальному условию (2), называют задачей Коши (начальной задачей) для уравнения (1). §3. Теорема существования и единственности решения задачи Коши для уравнения Теорема 1 (существования и единственности решения). Пусть имеем дифференциальное уравнение и пусть функция f(x, у) определена в некоторой области D на плоскости хОу. Выберем произвольную точку Мо(хо,уо) € D.

точки, в которой функция 1) непрерывна по совокупности аргументов; 2) имеет ограниченную частную производную ^, то найдется интервал оси Ох, на котором существует, и притом единственная, функция у = , являющаяся решением уравнения(1) и принимающая при х = жо значение уо (рис. Геометрически это означает, что через точку ) проходит одна и только одна интегральная кривая уравнения (1). Теорема 1 имеет локальный характер: она гарантирует существование единственного решения у = tp(x) уравнения (1) лишь в достаточно малой окрестности точки х0-

Из теоремы 1 вытекает, что уравнение (1) имеет бесконечное множество различных решений (например, одно решение, график которого проходит через точку (жо, г/о); другое решение, когда график проходит через точку Пример 1. В уравнении функция определена и непрерывна во всех точках плоскости хОу и имеет всюду щ = 1. В силу теоремы 1 через каждую точку (xq, уо) плоскости хОу проходит единственная интегральная кривая этого уравнения. Пример 2.

В уравнении функция определена и непрерывна на всей плоскости хОу. Здесь так что второе условие теоремы 1 нарушается в точках оси Ох. Нетрудно проверить, что функция где С — любая постоянная, является решением данного уравнения. Кроме того, уравнение имеет очевидное решение Если искать решения этого уравнения, соответствующие условию у(0) = 0, то таких решений найдется бесчисленное множество, в частности, следующие (рис. 2): Таким образом, через каждую точку оси Ох проходят по крайней мере две интегральные кривые и, следовательно, в точках Этой оси нарушается единственность.

Если взять точку Л/|(1, 1), то в достаточно малой ее окрестности выполнены все условия теоремы 1. Следовательно, через данную точку в малом квадрате П проходит единственная интегральная кривая ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ ПЕРВОГО ПОРЯДКА Эквивалентные дифференциальные уравнения Задача Коши еорема существования и единственности решения задачи Коши для уравнения Приближенные методы интегрирования уравнения Метод последовательных приближений Численные методы решения задачи Коши Метод Эйлера Понятие о методе Рунге—Кутта уравнения у’ = 3у >.

Если квадрат Q взять достаточно большим (подумайте, каким), то в нем единственность решения уже не будет иметь места. Это подтверждает локальный характер теоремы 1. Теорема 1 дает достаточные условия существования единственного решения уравнения у’ = . Это означает, что может существовать единственное решение у = у(х) уравнения у’ = f(x> у)* удовлетворяющее условию у\х=го = у0, хотя в точке (жо> уо) не выполняются условия 1) или 2) теоремы или оба вместе.

Пример 3:

Для уравнения имеем В точках оси Ох функции / и разрывны, причем Но через каждую точку (хо, 0) оси Ох проходит единственная интегральная кривая Замечание. Если отказаться от ограниченности , то получается следующая теорема существования решения. Теорема 2. Если функция f(x, у) непрерывна в некоторой окрестности точки (х0, уо), то уравнение имеет в этой окрестности по крайней мере одно решение , принимающее при х = х0 значение уо.

Теорема существом*и« и единственности решение задачи Кошм для уравнения Задача. Найти интегральную кривую уравнения проходящую через точку . Задача. Найти решение задачи Коши Определение 1. Общим решением дифференциального уравнения в некоторой области Q существования и единственности решения задачи Коши называется однопараметрическое семейство S функций у = ip(x, С), зависящих от переменной ж и одной произвольной постоянной С (параметра), такое, что 1) при любом допустимом значении постоянной С функция является решением уравнения (1): 2) каково бы ни было начальное условие .

Можно подобрать такое значение Со постоянной С, что решение у = Со) будет удовлетворять начальному условию При этом предполагается, что точка (жо, уо) принадлежит области П существования и единственности решения задачи Коши. Пример 4. Показать, что общим решением дифференциального уравнения является функция где С — произвольная постоянная. 4 В данном случае /, и условия теоремы 1 выполняются всюду. Следовательно, через каждую точку (zo, уо) плоскости хОу проходит единственная интегральная кривая данного уравнения. Проверим, что функция удовлетворяет условиям 1) и 2), содержащимся в определении общего решения.

Действительно, при любом С имеем так что у = х + С есть решение данного уравнения. Потребовав, чтобы при х = хо решение принимало значение уо, приходим к соотношению уо = х0 + Со, откуда Решение . или удовлетворяет поставленному начальному условию. Частным решением дифференциального уравнения (1) называется решение, получаемое из общего при каком-либо конкретном значении произвольной постоянной С (включая ±оо).

Таким образом, общее решение этого дифференциального уравнения можно определить как множество всех частных решений уравнения. В процессе интегрирования дифференциального уравнения мы часто приходим к уравнению неявно задающему общее решение уравнения. Уравнение (2) называют общим интегралом дифференциального уравнения (\). Уравнение где Со — некоторое конкретное значение постоянной С, называется частным интегралом. Замечание.

Название происходит от того, что для простейшего дифференциального уравнения вида его общее решение действительно записывается при помощи обычного неопределенного интеграла Пример. Общий интеграл уравнения имеет следующий вид или В дальнейшем для краткости мы будем иногда говорить, что решение уравнения проходит через некоторую точку Mq(xq, уо), если точка М0 лежит на графике этого решения.

Определение 2:

Решение дифференциального уравнения (1) называется особым, если в каждой его точке нарушается свойство единственности, т. е. если через каждую его точку (so, уо) кроме этого решения проходит и другое решение уравнения (1), не совпадающее в сколь угодно малой окрестности точки График особого решения называют особой интегральной кривой уравнения. Геометрически это — огибающая семейства интегральных кривых дифференциального уравнения, определяемых его общим интегралом1*.

Если для дифференциального уравнения (1) в некоторой области D на плоскости хОу выполнены условия теоремы 1, то через каждую точку (хо, уо) € D проходит единственная интегральная кривая уравнения. Эта кривая входит в однопа-раметрическое семейство кривых образующих общий интеграл уравнения (1), и получается из этого семейства при конкретном значении параметра С, т.е. является частным интегралом уравнения (1). Никаких других решений, проходящих через точку (жо> 2/о)> здесь быть не может.

Следовательно, для существования особого решения у уравнения (1) необходимо, чтобы не выполнялись условия теоремы 1. В частности, если правая часть уравнения (1) непрерывна в рассматриваемой области D, то особые решения могут проходить только через те точки, где производная Ц становится бесконечной. Напомним, что огибающей семейства кривых Ф(х, у, С) = 0 называется такая кривая, которая в каждой своей точке касается некоторой кривой семейства и каждого отрезка которой касается бесконечное множество кривых из этого семейства.

Например, для уравнения ) функция / = Зу2/3 непрерывна всюду, но производная Ц обращается в бесконечность при у = 0. т.е. на оси Ох плоскости хОу. Уравнение (3) имеет общее решение — семейство кубических парабол — и очевидное решение проходящее через те точки, где производная ^ не ограничена. Решение у = 0 — особое, так как через каждую его точку проходит и кубическая парабола, и сама эта прямая у = 0 (см. рис. 2).

Таким образом, в каждой точке решения у = 0 нарушается свойство единственности. Особое решение у = 0 не получается из решения у = (х +С)3 ни при каком числовом значении параметра С (включая ±оо). Из теоремы 1 можно вывести только необходимые условия для особого решения. Множество тех точек, где производная не ограничена, если оно является кривой, может и не быть особым решением уже потому, что эта кривая, вообще говоря, не является интегральной кривой уравнения (1).

Если, например, вместо уравнения (3) взять уравнение то в точках прямой у = 0 по-прежнему нарушается условие ограниченности производной ^, но эта прямая, очевидно, не является интегральной кривой уравнения (4). Итак, чтобы найти особые решения уравнения (1), надо 1) найти множество точек, где производная ^ обращается в бесконечность; 2) если это множество точек образует одну или несколько кривых, проверить, являются ли они интегральными кривыми уравнения (1);

3) если это интегральные кривые, проверить, нарушается ли в каждой их точке свойство единственности. При выполнении всех этих условий найденная кривая представляет собой особое решение уравнения (1). Задача. Найти особые решения уравнения Сделать рисунок. §4. Приближенные методы интегрирования уравнения 4.1. Метод изоклин Пусть имеем дифференциальное уравнение где функция /(ж, у) в некоторой области D на плоскости хОу удовлетворяет условиям теоремы I.

Это уравнение определяет в каждой точке (ж, у) области D значение у’, т. е. угловой коэффициент касательной к интегральной кривой в этой точке. Говорят, что уравнение (1) определяет в области D поле направлений. Чтобы его построить, надо в каждой точке (ж0, уо) € D представить с помощью некоторого отрезка направление касательной к интегральной кривой в этой точке, определяемое значением /(жо, уо). Совокупность этих отрезков дает геометрическую картину поля направлений.

Задача интегрирования дифференциального уравнения (1) может быть теперь сформулирована так: найти такую кривую, чтобы касательная к ней в каждой точке имела направление, совпадающее с направлением поля в этой точке. Такое истолкование дифференциального уравнения и его интегрирования дает графический способ решения уравнения. Для построения интегральных кривых пользуются изоклинами. Изоклиной называется множество всех точек плоскости хОу, в которых касательные к искомым интегральным кривым имеют одно и то же направление (у’ = const).

Из этого определения следует, что семейство изоклин дифференциального уравнения (1) задается уравнением где к — числовой параметр. Если придать параметру к близкие числовые значения, можно найти достаточно густую сеть изоклин и приближенно построить интегральные кривые дифференциального уравнения. Пример 1. Проинтегрировать уравнение по способу изоклин. Семейство изоклин данного уравнения определяется уравнением ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ ПЕРВОГО ПОРЯДКА Эквивалентные дифференциальные уравнения Задача Коши еорема существования и единственности решения задачи Коши для уравнения Приближенные методы интегрирования уравнения.

Метод последовательных приближений Численные методы решения задачи Коши Метод Эйлера Понятие о методе Рунге—Кутта Полагая , получаем изоклины , по которым строим интегральные кривые уравнения (рис. 4). Нулевая изоклина определяет множество возможных точек экстремума интегральных кривых (прямая х = О в примере 1).

Для большей точности построения интегральных кривых определяют направление вогнутости и точки перегиба этих кривых (если такие точки существуют). Для этого находят у» в силу уравнения (1): Знак правой части определяет знак у», т.е. направление вогнутости интегральных кривых. Линия, заданная уравнением есть множество всех возможных точек перегиба интегральных кривых. В примере 1 имеем поэтому все интегральные кривые обращены вогнутостью вверх, и точек перегиба интегральных кривых нет.

Метод последовательных приближений

Пусть имеем дифференциальное уравнение где функция f(x, у) в некоторой области D изменения х, у удовлетворяет условиям теоремы 1, и пусть точка (хо, уо) € D. Решение задачи Коши равносильно решению некоторого интегрального уравнения, т. е. уравнения, в которое неизвестная функция входит под знаком интеграла. В самом деле, пусть — решение уравнения (2), заданное в некоторой окрестности точки х0 и удовлетворяющее начальному условию (3).

Тогда при имеет место тождество Проинтегрируем это тождество по х Отсюда учитывая (3), получаем так что решение y(x) задачи Коши удовлетворяет интегральному уравнению (4) Обратно: если непрерывная функция , удовлегворяет интегральному уравнению (4), то, как легко проверить, у(х) является решением задачи Коши (2)-(3). Решение интегрального уравнения (4) для всех х, достаточно близких к хо» может быть построено методом последовательных приближений по формуле причем в качестве (po(t) можно взять любую непрерывную на отрезке функцию, в частности, ( Прммер 2.

Методом последовательных приближений решить задачу Коши Сводим данную задачу к интегральному уравнению Выбирая за нулевое приближение функцию последовательно находим: Легко видеть, что функция у = ех есть решение задачи. 4.3. Численные методы решения задачи Коши Метод Эйлера Пусть требуется найти приближенное решение дифференциального уравнения удовлетворяющее начальному условию.

Будем предполагать, что в некотором прямоугольнике функция f(x, у) непрерывна и имеет непрерывные частные производные достаточно высокого порядка по всем аргументам, так что решение задачи Коши (1)-(2) существует, единственно и является функцией, дифференцируемой достаточное число раз. Численное решение задачи состоит в построении таблицы приближенных значений , уп решения задачи в точках Чаще всего выбирают . Точки хк называют узлами сетки, а величину — шагом сетки.

Геометрически искомая интегральная кривая у = у(х), проходящая через точку заменяется ломаной Эйлера MqM\Mi . с вершинами в точках Метод Эйлера относится к группе одно-шаговых методов, в которых для вычисления точки требуется знание только предыдущей вычисленной точки . Для оценки погрешности метода на одном шаге сетки разложим точное решение в окрестности узла х = Хк по формуле Тейлора Сравнение формул () показывает, что они совпадают до членов первого порядка по h включительно, а погрешность формулы (4) равна 0(h2). Поэтому говорят, что метод Эйлера имеет первый порядок.

Пример. Методом Эйлера решить задачу Коши на отрезке с шагом В данном случае . Пользуясь формулой последовательно находим Результаты вычислений сведем в таблицу Точное решение Замечание. Если рассмотреть задачу Коши на любом отрезке с любым шагом , то получим и т.д., так что в этом случае ломаная Эйлера «распрямляется» и совпадает с прямой точным решением поставленной задачи Коши. Понятие о методе Рунге—Купа Метод Эйлера весьма прост, но имеет низкую точность.

Точность решения можно повысить путем усложнения разностной схемы. Весьма распространенными на практике являются схемы Рунге—Кутта. Пусть опять требуется решитьзадачу Коши . Будем строить таблицу приближенных значенийу решения уравнения вточках(узлах сетки). Рассмотрим схему равноотстоящих узлов — шаг сетки. В методе Рунге—Кутта величины yi+\ вычисляются по следующей схеме

Присылайте задания в любое время дня и ночи в ➔

Официальный сайт Брильёновой Натальи Валерьевны преподавателя кафедры информатики и электроники Екатеринбургского государственного института.

Все авторские права на размещённые материалы сохранены за правообладателями этих материалов. Любое коммерческое и/или иное использование кроме предварительного ознакомления материалов сайта natalibrilenova.ru запрещено. Публикация и распространение размещённых материалов не преследует за собой коммерческой и/или любой другой выгоды.

Сайт предназначен для облегчения образовательного путешествия студентам очникам и заочникам по вопросам обучения . Наталья Брильёнова не предлагает и не оказывает товары и услуги.