«innovative academy» ilmiy tadqiqotlarni


Download 1.67 Mb.
Pdf ko'rish
bet50/70
Sana17.06.2023
Hajmi1.67 Mb.
#1530764
1   ...   46   47   48   49   50   51   52   53   ...   70
Bog'liq
Zamonaviy dunyoda innovatsion tadqiqotlar 1-son

Ключевые слова: DiffSer(), intDiff(), piecewise(), DiffInner(), Функция Тейлора, 
математический маятник. 
Применение рядов для решения ДУ рассмотрим на примере нелинейного 
уравнения второго порядка вида 
)
,
,
(
y
y
x
f
y



удов 
летворяющего начальным условиям 
.
)
(
,
)
(
0
0
0
0
y
x
y
y
x
y




Предположим, что решение поставленной задачи существует и будем его искать в 
форме ряда Тейлора функции y(x). 
...
)
)(
(
2
1
)
)(
(
)
(
)
(
2
0
0
0
0
0








x
x
x
y
x
x
x
y
x
y
x
y
Для построения решения нужно знать производные в точке x0 порядка 1, 2, 3, …, но 
это можно сделать с помощью самого уравнения и его начальных условий. 
Действительно, из начальных условий мы сразу имеем значение функции и ее первой 
производной в точке 
0
x
. Подставляя их в уравнение, сразу находим значение 2-й 
производной 
)
,
,
(
)
(
0
0
0
0
y
y
x
f
x
y



Дифференцируя обе части исходного уравнения по 
аргументу x, получим выражение для 3-й производной 
y
y
y
y
x
f
y
y
y
y
x
f
x
y
y
x
f
y
















)
,
,
(
)
,
,
(
)
,
,
(
Подставляя в правую часть 
0
x
x

и значения предыдущих производных в этой 
точке, находим значение 3-й производной 
)
(
0
x
y

в точке 
0
x
x

. Дифференцируя 
последнее выражение еще раз, и выполняя подстановку известных производных, 
находим значение 4-й производной в точке 
0
x
x

. Для определения производных 
функции y(x) в точке 
0
x
до требуемого порядка n дифференцирование и подстановку 
выполняем необходимое число раз. Затем найденные значения производных 
подставляем в представление решения в форме ряда Тейлора. Для тех значений для 
которых этот ряд сходится, он будет представлять искомое решение задачи Коши.
Для решения задачи Коши 
)
,
,
(
y
y
x
f
y



,
0
0
)
(
y
x
y

,
0
0
)
(
y
x
y



используем
следующую процедуру(см.[1],стр.450).
DiffSer:=proc(f::anything,x::name,y::name,y1::name,\
x0::numeric,y0::numeric,y01::numeric,n::integer)
local p,der,i,j,s:



Download 1.67 Mb.

Do'stlaringiz bilan baham:
1   ...   46   47   48   49   50   51   52   53   ...   70




Ma'lumotlar bazasi mualliflik huquqi bilan himoyalangan ©fayllar.org 2024
ma'muriyatiga murojaat qiling