Текст подпрограммы и версий de26r_p.zip , de26e_p.zip |
Тексты тестовых примеров tde26r_p.zip , tde26e_p.zip |
Выполнение одного шага численного интегрирования системы обыкновенных дифференциальных уравнений первого порядка методом Адамса без контроля точности.
Выполняется один шаг численного интегрирования системы М обыкновенных дифференциальных уравнений первого порядка
Y ' = F (X, Y) , Y = ( y1, ..., yM ) , F = ( f1 (X, y1,..., yM),..., fM (X, y1,..., yM) )
многошаговым предсказывающе - исправляющим методом пятого порядка точности. Суть метода состоит в том, что сначала по явной формуле Адамса строятся предсказанные значения приближенного решения, которые затем исправляются по неявной формуле Адамса.
По заданному значению YX решения в узле Xn вычисляется значение этого решения в узле Xn + H. Значение H шага интегрирования задается пользователем при обращении к подпрограмме.
Березин И.С., Жидков Н.П. Методы вычислений, т.2. Физматгиз, М., 1960.
Беленький В.З. Стандартная программа для интегрирования системы дифференциальных уравнений методом Адамса. Сб. "Вычислительные методы и программирование", вып 3. Изд-во МГУ, 1965
procedure DE26R(F :Proc_F70; M :Integer; var JSTART :Integer; var YX :Array of Real; var X :Real; H :Real; var DELTY :Array of Real; var DF :Array of Real; var RF :Array of Real; var YP :Array of Real; var RY :Array of Real; var RFN :Array of Real; var IERR :Integer);
Параметры
F - |
имя подпрограммы вычисления значений правой
части дифференциального уравнения. Первый
оператоp подпрограммы должен иметь вид: procedure F (X :Real; var Y :Array of Real; var DY :Array of Real; M :Integer); Здесь: X, Y - значения независимой и зависимой переменных, соответственно. Вычисленное значение правой части должно быть помещено в DY. B случае системы уравнений, т.е. когда M ≠ 1, параметры Y и DY представляют массивы длины M (тип параметров X, Y и DY: вещественный); |
M - | количество уравнений в системе (тип: целый); |
JSTART - | целый указатель режима использования подпрограммы, имеющий следующие значения: |
0 - | первое обращение к подпрограмме должно быть выполнено с нулевым значением JSTART; |
+ 1 - | выполнить один шаг интегрирования системы дифференциальных уравнений для значений независимой и зависимой переменных и шага интегрирования, заданных параметрами X, YX и H, соответственно; |
- 1 - | повторить последний шаг интегрирования с новым значением параметра H; |
на выходе из подпрограммы JSTART полагается равным + 1; | |
X, YX - | начальные вещественные значения аргумента и решения уравнения в нем; в результате работы подпрограммы в X получается новое значение аргумента, а в YX - соответствующее значение решения; в случае системы уpавнений, т.е. когда M ≠ 1, YX задается одномерным массивом длины M; |
H - | вещественное значение шага интегрирования; |
DELTY - RF, YP RFN | одномерные вещественные рабочие массивы длины M; |
DF, RY - | двумерные вещественные рабочие массивы размеpа M*5; |
IERR - | целая переменная, значение которой в pезультате работы подпрограммы полагается равным 65, если при заданном значении шага интегрирования не может быть достигнута сходимость численного решения к точному значению решения; в этом случае интегрирование системы можно повторить обращением к подпрограмме с новым значением параметра H и со значением JSTART = 1 . |
Версии
DE26E - | выполнение одного шага численного интегрирования системы обыкновенных дифференциальных уpавнений первого порядка методом Адамса без контроля точности, при этом все вычисления над действительными числами ведутся с удвоенным числом значащих цифр. В этом случае параметры YX, X, H, DELTY, DF, RF,YP, RY, RFN и параметры X, Y и DY в подпрограмме F должны иметь тип Extended. |
Вызываемые подпрограммы
DE33R - | построение начальных значений при интегрировании системы обыкновенных дифференциальных уpавнений первого порядка методом Адамса без контроля точности; вызывается при работе подпрограммы DE26R. |
DE33E - | построение начальных значений при интегрировании системы обыкновенных дифференциальных уpавнений первого порядка методом Адамса без контроля точности, при этом все вычисления над действительными числами ведутся с удвоенным числом значащих цифр; вызывается при работе подпрограммы DE26E. |
UTDE12 - | подпрограмма выдачи диагностических сообщений при работе подпрограммы DE26R. |
UTDE13 - | подпрограмма выдачи диагностических сообщений при работе подпрограммы DE26E. |
Kpоме того, DE26R и DE26E используют рабочие подпрограммы DE28RP, DE28RS и DE28EP, DE28ES, соответственно. |
Замечания по использованию
При работе подпрограммы и ее версии значения параметров M и H сохраняются. При многократном использовании подпрограммы и ее версии для вычисления решения на отрезке значения параметров M, YX, X, DF, YP не должны изменяться в вызывающей программе между последовательными обращениями к ним. |
y1' = y2 y2' = -y1 y1 (3/4 π) = √2 /2 , y2 (3/4 π) = -√2 /2
Приводятся подпрограмма вычисления значений правой части системы и фрагмент вызывающей программы, выполняющей несколько шагов из одной и той же точки, а также результаты счета.
Unit TDE26R_p; interface uses SysUtils, Math, { Delphi } Lstruct, Lfunc, UtRes_p, FDE26R_p, DE26R_p; function TDE26R: String; implementation function TDE26R: String; var M,JSTART,IH,IERR :Integer; XN,H,X,Y1,Y2 :Real; RFN :Array [0..1] of Real; YN :Array [0..1] of Real; Y :Array [0..1] of Real; DELTY :Array [0..1] of Real; DF :Array [0..9] of Real; RF :Array [0..1] of Real; YP :Array [0..1] of Real; RY :Array [0..9] of Real; label _6,_7,_8,_20; begin Result := ''; { результат функции } M := 2; XN := 0.75*3.14159265359; Н := 0.01; H := -H; Y[0] := Sqrt(2.0)/2.0; Y[1] := -Y[0]; JSTART := 0; IH := 0; Х := XN; _6: DE26R(FDE26R,M,JSTART,Y,X,H,DELTY,DF,RF,YP,RY,RFN,IERR); IH := IH+1; Y1 := Sin(X); Y2 := Cos(X); Result := Result + Format(' %20.16f %20.16f %20.16f %20.16f %20.16f %20.16f ', [X,Y[0],Y[1],H,Y1,Y2]) + #$0D#$0A; if ( IH = 1 ) then goto _6; if ( IH = 5 ) then goto _20; JSTART := -1; if ( IH = 4 ) then goto _8; if ( IH = 3 ) then goto _7; H := 0.005; H := -H; goto _6; _7: H := 0.02; goto _6; _8: H := 0.01; goto _6; _20: UtRes('TDE26R',Result); { вывод результатов в файл TDE26R.res } exit; end; end. Unit fde26r_p; interface uses SysUtils, Math, { Delphi } Lstruct, Lfunc; procedure fde26r(X :Real; var Y :Array of Real; var DY :Array of Real; M :Integer); implementation procedure fde26r(X :Real; var Y :Array of Real; var DY :Array of Real; M :Integer); begin DY[0] := Y[1]; DY[1] := -Y[0]; end; end. Результаты: после первого обращения к подпрограмме - X Y(1) Y(2) 2.346194490+00 7.141423761-01 -7.000004762-01 H Y1 Y2 1.000000000-02 7.141423761-01 -7.000004762-01 после второго обращения к подпрограмме - X Y(1) Y(2) 2.336194490+00 7.211065574-01 -6.928241717-01 H Y1 Y2 1.000000000-02 7.211065574-01 -6.928241717-01 после третьего обращения к подпрограмме - X Y(1) Y(2) 2.341194490+00 7.176334371-01 -6.964210292-01 H Y1 Y2 5.000000000-03 7.176334371-01 -6.964210292-01 после четвертого обращения к подпрограмме - X Y(1) Y(2) 2.366194490+00 7.000004763-01 -7.141423759-01 H Y1 Y2 2.000000000-02 7.000004762-01 -7.141423761-01 после пятого обращения к подпрограмме - X Y(1) Y(2) 2.356194490+00 7.071067812-01 -7.071067812-01 H Y1 Y2 1.000000000-02 7.071067812-01 -7.071067812-01