УДК 517.928.4
ВЛИЯНИЕ ВОЗМУЩЕНИЯ ПОДВИЖНОЙ ОСОБОЙ ТОЧКИ НА ПРИБЛИЖЕННОЕ РЕШЕНИЕ ОДНОГО КЛАССА НЕЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ ВТОРОГО ПОРЯДКА
INFLUENCE OF PERTURBATIONS OF THE MOVABLE SINGULAR POINT ON THE APPROXIMATE SOLUTION OF A CLASS OF NONLINEAR DIFFERENTIAL
EQUATIONS OF SECOND ORDER
В.Н. Орлов 1, Т.Ю. Леонтьева 2 V.N. Orlov 1, T.Yu. Leont'eva 2
1) Гуманитарно-педагогическая Академия (филиал) «КФУ им. В.И. Вернадского»,
298635, г. Ялта, ул. Севастопольская, д. 2а, тел/факс (0654У32-30-13. Humanity and Pedagogical Academy (Branch) «V. Vernadsky named CFU» in Yalta, 298635, city of Yalta, st. Sevastopol, 2a, tel/fax: (0654) -32-30-13,
2) Чувашский государственный педагогический университет им. И.Я. Яковлева,
428000, г. Чебоксары, ул. К. Маркса, д. 38, тел (8352^62-03-12 I. Yakovlev named Chuvash State Pedagogical University, 428000, Cheboksary, st. K. Marx, 38, tel (8352)^2-03-12
E-mail: [email protected] [email protected]
Аннотация. В современном мире с нелинейными дифференциальными уравнениями можно встретиться почти во всех областях науки и техники. Но решение этих уравнений связано с большими трудностями в связи с наличием подвижных особых точек. В данной работе представлена апробация одного из шагов метода приближенного решения нелинейных дифференциальных уравнений для одного класса дифференциальных уравнений. Данный метод включает решение шести задач. Решение первых двух задач: доказательство теоремы существования и единственности решения нелинейного дифференциального уравнения; построение приближенного решения и исследование влияния возмущения начальных условий на приближенное решение, опубликованы ранее. В статье рассмотрено влияние возмущения подвижной особой точки на приближенное решение.
Resume. In today's world with nonlinear differential equations can be found in almost all fields of science and technology. But the solution to these equations is connected with great difficulties due to the presence of moving singular points. This paper presents the testing of one of the steps of the method of approximate solution of nonlinear differential equations for a class of differential equations. This method involves the solution of six problems. The decision of the first two tasks: the proof of the existence and uniqueness of solutions of nonlinear differential equations; construction of an approximate solution and investigation of the influence of the perturbation of the initial conditions on an approximate solution, published earlier. The article considers the influence of perturbations of movable singular point on the approximate solution of this class of nonlinear differential equations of second order. These results are accompanied by estimates.
Ключевые слова: подвижная особая точка, нелинейное дифференциальное уравнение второго порядка, приближенное решение, метод мажорант, окрестность подвижной особой точки, возмущение подвижной особой точки, апостериорная погрешность
Key words: movable singular point, non-linear second-order differential equation, approximate solution, majorant method, a neighborhood of the movable singular point, the perturbations of the movable singular point, a posteriori error.
Введение
Нелинейные дифференциальные уравнения представляют большой интерес в связи с их приложением во многих областях науки и техники [Kalman, 1961; Axford, 1970; Hill, 1977; Ockendon, 1978; Shi, 2005]. Решения нелинейных дифференциальных уравнений связаны с большими трудностями, вызванными наличием подвижных особых точек у интегралов этих уравнений, которые и являются препятствием к использованию известных на данный момент приближенных
численных и аналитических методов решения [Березин, Жидков, 1960; Бахвалов, 1970; Фильчаков, 1970]. Упомянем о работах, направленных на разрешение в квадратурах нелинейных дифференциальных уравнений с подвижными особыми точками, но это удается сделать лишь в частных случаях [Яблонский, 1964; Еругин, 1967; Самодуров, 1983; Кондратеня и др., 1988; Мататов, Сабынич, 1991; Лукашевич, 1995; Чичурин, Швычкина, 2014]. В данной статье дано исследование влияния возмущения подвижной особой точки на приближенное решение в окрестности подвижной особой точки.
Объекты и методы исследования. Объектом исследования является один класс нелинейных дифференциальных уравнений второго порядка с полиномиальной правой частью пятой степени. Существующие методы решения дифференциальных уравнений можно классифицировать: 1) точные; 2) приближенные; 3) асимптотические. Большинство работ посвященных решению дифференциальных уравнений можно отнести в основном к точным и частично к асимптотическим методам, и только незначительная их часть использует формальный аппарат приближенных методов, без строгого доказательства соответствующих пунктов этих методов. В данной работе приведен очередной этап приближенного аналитического метода решения нелинейных дифференциальных уравнений с подвижными особыми точками, позволяющий получить решение рассматриваемого класса нелинейных дифференциальных уравнений с заданной точностью в окрестности подвижной особой точки. Решение данной задачи для других классов нелинейных дифференциальных уравнений представлены в работах [Орлов, 2006; Орлов, 2008а, б; Орлов, 2009]. Предыдущие этапы аналитического приближенного метода решения рассматриваемого класса нелинейного дифференциального уравнения были опубликованы ранее [Орлов, Леонтьева, 2013а, б].
Результаты и их обсуждение. Для задачи Коши
/(*)=у5 М+Ах), (1)
У(хо) = Уо> У'(хо) = У (2)
в случае точного значения подвижной особой точки в работе [Орлов, Леонтьева, 2014] было
*
получено приближенное решение в окрестности подвижной особой точки X в виде
у(х) = (х* - х)-2-± Ся (х* - х)2, Со * 0. (3)
п=0
Возмущение подвижной особой точки ~ * оказывает влияние на структуру аналитического приближенного решения (3), которое принимает следующий вид:
N ~ / ~ \(п-1)/2
У»(x)=Ncn-*г)/2, Co *0, (4)
n=0
где Cn - возмущенные значения коэффициентов, ~ - возмущенное значение подвижной особой точки.
Теорема. Пусть выполняются следующие условия: r (x )е C * в области
i ) |r (n )(~У *)
K = \x :|Г-X <P0| P0 = const > 0; 3M0;I-< M0, M 0 = const, n = 0,1,2,....; x*< x*;
1 1 n!
известны оценки погрешности У* и с : — х* <А~*, с — Сс < Ас~
где
Ау* < у ^4 • 5{М + I)2 ^ . Тогда для аналитического приближенного решения (4) задачи (1)-(2) в областях
х — р3< х < х — Ах
справедлива оценка погрешности
где
АуК {х )<Л0 +Л! +Л2 +А3,
{5) {6)
Ах
А =--4
А 0 ~ * 4
х — х
2кМ{М +1), 5 ] У* — х 2 8 2'^(М +1),5]• У* — х 2
Л = "" Х {ж +1 + 2){ж +1 — б)
+1) У *— х|52 '=0
1 — 25 {М +1) ~
4
А2 = 26М{М+ У \ Г¿22{М +1)7 р1 +2А1/2¿22' {М +1)7 р'
1 — 210 {М +1)2 р
Л =
26 (АЛ~ + ])//Р
1=0 2
1=0
1 — 2>р5
4 , 4
X 221 !ЛП р1 +2р1/2 X 221 ¡л72 р1
Р0 = тЦл, Р2},
а =■
1
4^ {М +1)2
(из [Орлов, Леонтьева, 2014]), р2 =
8{м + ам +1)
Р
~ * — х|, х е {5) А~ *, х е{б)
1 =
|1 +1,1 = 0,1,2,3,4, 9 — 1, г = 5,6,7,8
! = М + АМ +1,
71 =■
(0,1 = 0,1,2 1,1 = 3,4 !
72 =
0, Г = 0,1 1,1 = 2,3,4!
М = тах <
с, вир
гу 7 х
и!
М 2 !> АМ
вир
г
Л
и!
А~ , с -
параметр,
зависящий от условий (3), п = 0,1,2,.... Доказательство. Используя классический подход, имеем
Аук{х)=|Ях)—{х) < |Ях )—~{х)+{х)—Ук{х).
Оценим |у{х )— ~ {х):
|Ях) — У{х)
<
сп {х — х) сп {~ — х)
<
<
V с {х* х){п—1У2 V с х Лп—1У2 с {х* х){п—1)/2 V с {~* Лп—1У2
X сп {х — х) сп {х — х) сп {х — х) сп {х — х)
<
~ — А~ < х < ~
N—1
1=0
1=0
1
*
*
п
п
и=0
п=0
и=0
п=0
п=0
и=0
<
К- ( - Уп-1^2 -А ~ // - \(п-1У2 /-- \(п-1У2^
АСп (х - х) + Х Сп ((х - Х) -(-Х - х) )
п=0 п=0
<
<£ дС. (г - х +дт !п-1)"|+]Г|^|.|(х- - хГ2 -(г - х)(п-1
п=0
п=0
- \(п-1)/ 2 X
Рассмотрим ряд п
п=0
VI- (- )(п-1У2 /-- V X Сп '(х - х) -(- - х)
- |(п-1 у 2 X
Так как х < г * < х* и Со | =
С
= 2 — , то при п = 0 V 4
~ I - \-12 уу2
С • (х - х) -(х - х)
Дх
<--4
- х V
х - х
Учитывая |С| = С = 0, |С2| = С2| = 0, |С3| = |сз| = 0, |С4| = |С4| = 0, получаем:
|у(х)-- (х )< 3 +]Г| I ](х' - хГ» "(Г - х)
х - х
- \(п-1)/2
х
п=5
+
п=5
■А ~ ( „ Лп-1)/ 2
+ ХДС„ (- - х + Д- )
Далее, в случае п = 5,6,7,...
\(п-1)/ 2
/ - \(п-1>2 /-- \( (х - х) - (х - х)
,(п-1)/ 2
<
•---- - « ----- -
х - х + дх
)(п-0/2 _ (-- _ Л)(п-1У2
<
<
Д- (- - х + Д-
Следовательно, для оценки приближенного решения (4) имеем:
|у(х)- Ум(х )
ДТ
. • 4
х - х
3+ Х С-.
п=N+1
х - х
(п-1)/ 2
да /
+Е 1ч-Д- *(- -
п=5
х
+ Дх'
\(п-2У 2
+
п=5
- / „ Л(п-1)/2
+X дс •(- -х + Д- )
= Д0 + Д + Д2 + Д,
где С - С = дС .
п п
Таким образом,
Дх-
Д =--4
х - х
Выражение оценки Д следует из теоремы 2 работы [Орлов, Леонтьева, 2014].
Перейдем к оценке Д2. Проведем суммирование отдельно по целым и дробным степеням,
учитывая, что Д- < - - х :
(п-2)/2
да
*
*
А 2 =Х\Сп\ * {У х + АУ 2 2 2 = Е |(У2п—^ ■АУ * (•
п=5 п=3
+ Х\с2п •ах*{х* — х + ау*) =а21
./* I ~ *
{х — х
+ А~'
\{2и—3)/ 2
+
п=3
:А2,1 +А2,2 .
С учетом закономерности получения оценок для Си , следует:
А2Д = X
п=3
с.
|{2п—3)/2
к=1 ТО
к=1
с
10к—3
Ах ~ { у * — х| + Ах * )
10к—5)/ 2
=х
к=1
с
10к —5
Ау * {у * — х| + Ау * )
+ X с10к—1
к=1
АУМ|У* | ЛУ*^10к—7 У2
Ах {х — х\ + Ах )
10к—3У 2
+
Ау * { у * — х| + ах *)
+
с
10к+1
У-/!У» | х» ¥10к—1 Ах {.х — х + Ах )
-1)/ 2
+ Е
к=1
с
10к+3
Ах * {у * — х| + Ах *)
10к+1У2
4 то
= ЕЕ
г=0 к=1 >6 ЛУ*
с
10к—5+2г
•|Ах * {у * — х| + Ах
10к—7+2г У 2
<
<
26 АУ *М {М +1) У * —
х
3/2
\2 ~ * • х — х
4
•¿22г {М +1)7
х — х
г=0
1 — 210 •{М +1)
I К и-\Т> [0,1 = 0,1,2
при условии |х *— х| < 1 ^4 • у {М +1) J, где 7\=< ^ ^ .
В случае х * — х < Ах * получим:
А 2,1 < ■
26 М {М + 1)АУ
>52
1 — 210•{М +1)2 •АУ*5
•¿22г{М +1)71 АУ*г.
г=0
Аналогичным способом получим оценку для А
2,2
А2,2 < 27М{М + ^• ¿¿22г {М +1)72 Ау*г,
4
1 — 210 •{М +1)2 •АУ
г=0
0,1 = 0,1
где 72 = ■
[1, г = 2,3,4
Найдем оценку для А3 . Для выражения ас предполагаем оценки
ЛХ5И <25,1 (АЛЛ + &+АМ + 1)и , лс5„+1 <
25и+1 {АМ + ИМ + АМ + 1)п {5и + 3){5и — 5) '
___25и+3 {ам + Им + АМ + 1)п
5и+2 < {5и + 4){5и — 4) ' 5и+3 < {5и + 5){5и — 3)
25и+4 {ам + им + АМХ + 1)и {5и + 6){5и — 2) '
А^о <
{5и + 2){5и — 6)
25и+2 {алл+Им+ам+1)и
Ас5п+4 <
ТО
то
То'
ТО
где М = Бир
г(п
(
п!
, ДМ =
Бир
п!
Дх . Докажем оценку для Д(-5п в случае
N +1 = 5(2п +1)
ДС
10п+5
С
10п+5 - С10п+5
210п+5 М (М + 1)2п+1 210п+5 М (М +1)
2п+1
>10и+5
(10п + 7)(10п -1)
(10п + 7)(10п -1) (10п + 7)(10п -1)
(М + дм )М + ДМ + 1)2п+1 - М (М +1)2
\2и+1
210п+5 (м + ДМ + 1)2п+1
(10п + 7)(10п -1)
(м + дм)-м|
М +1
Лп+1
М + ДМ +1,
210п+5 (м + ДМ +1)
(10п + 7)(10п -1)
2п+1
(
М + ДМ - М
ДМ
,2п+1
1 -
ч М + ДМ + 1,
<
<
210п+5 (дм- + 1)(м + ДМ +1)
2п+1
(10п + 7)(10п -1)
Аналогичные выражения имеем и в случаях N +1 = 5п +1, N +1 = 5п + 2 , N +1 = 5п + 3 и N +1 = 5п + 4 . Таким образом получаем оценку
ДСп+1 <
2п+1 (дм-+цм+Дм+1|(п+1); 5]
(п + 3)п - 5)
Разделяя целые и дробные степени в выражении Д , имеем:
д з =хК
п-Ъ
да
+ХД(-
* А '—'*
х - х +
Дх ) =ХДС2п-1
п=3
х - х + Дх
п—1
+
х - х + Дх
(2и-1У 2 4 ш —
= X X ДС10£-5+2г
t=0 к=1
+ХХДС
10к -4+2t
t=0 к=1
— - . х - х + Дх
х - х + Дх
(10к -5+2t У 2
5к-3+t
+
4 да
< X X
t=0 к=1
210к-5+2t (дМ- + 1)(М + ДМ— +1)[(]
10к - 5+2t )/ 5]
(10к - 3 + 2t )(10к -11 + 2t)
<
х - х
+ Д—
5к-3+t
+
4 да 210к -4+2t (дм- + 1)(М + дМ + 1)[(10к -4+2t^5]
+ X X
t=0 к=1
4
(10к - 2 + 2Г )(10к -10 + 2Г)
х - х
+ Д—
(10к - 5+2t У 2
.¿2* (дл- + 210к-5 (М + ДЛ— + Г''-5+2>5] ^о 1 ¿X (10к - 3 + 2t )(10к -11 + 2t)
-
х - х
+ Д-
5к-3+t
+
*
г
п
п
п=3
*
| - 4 M + ЛМ + - 4+2t V5] + ¿1 (l°k - 2 + 2t )(l0k -10 + 2t)
— *
x - x
+ Л—
(l°k-5+2t )/ 2
<
< +2P12¿2*^p<\,
1 - 2 ¡Л p v-=° ¿=0
- 0.x x\, x e (5) ~ , [0,i = 0,1,2 [ 0,i = 0,1
ГДе ? = Лх',,^e(6) , * = m + ш +1 • * = i 1,- = 3,4 , * =-[1,- = 2,3,4'
Рассматриваемые в работе ряды являются сходящимися в областях
х -Р0 < x < — -Л— и — -Л~ < x < — ,
где
—* * x - x
< Л— * и р0 = min
1
1
4^{М +1)2 ' 8{м + аЛХ +1)2 ] Замечание 1. Теорема 3 справедлива в области
У* А У* ^ У*
х +Ах < х < х + Р0
— * < х < — * + Ау *.
* .—' *
При этом пункт 3 теоремы примет вид х < х .
(7)
(8)
Пример. Найдем приближенное решение задачи (1)-(4) в случае г{х) = 0 при начальных данных
у{0,5) = 1, у'{0,5) = ^= и с = 0,0001. Величина возмущения не превышает е = 0,8 -10 3. -V 3
Данная задача имеет точное решение y =
1 + л/э - 2x
. Найдем радиус окрестности подвижной
„ „ „ * 1+V3
особой точки Р0 « 0,0312497 . Точное значение подвижной особой точки х =——— .
Произведем расчет в случае — = 1,366025, л— = 0,0000004, значение аргумента х = 1,34 и
структуры приближенного решения —12 . Результаты представлены в таблице. Таблица.
Сравнительный анализ погрешностей приближенного решения Comparative analysis of the approximate solution error
у
x y —12 ¿У Л—12 Ay
1,34 5,768549 5,768594 0,000045 0,059 0,0008
где —12 - приближенное решение (4); у - значение точного решения; Л—12 - оценка погрешности приближенного решения, полученная по теореме; Ау - абсолютная погрешность
приближенного решения _—12; Д1 у - апостериорная оценка погрешности. В нашем случае для N = 25 априорная оценка будет удовлетворять требуемой точности £ = 0,8 -10 . Добавки в структуре приближенного решения для N = 13,...,25 не превышают требуемой точности. Поэтому в структуре приближенного решения можем ограничиться значением N = 12, при котором приближенное решение будет иметь погрешность £ = 0,8 • 10 3.
Заключение
В статье сформулирована и доказана теорема, отражающая влияние возмущения подвижной особой точки на приближенное решение одного класса нелинейных дифференциальных уравнений второго порядка в вещественной области.
Список литературы References
1. Axford R. A. 1970. Differential equations invariant urber two-parameter Lie groups with applications to nonlinear diffusion. Los Alamos Report. (LA-4517, UC-34).
2. Hill J. M. 1977. Radial deflections of thin precompressed cylindrical rubber bush mountings. International Journal of Solids and Structures, 13: 93-104.
3. Kalman R. 1961. New results in linear filtering and predication theory. Journal of Basic Engineering (American Society of Mechanical Engineers Trans.), 83D: 95-108.
4. Ockendon J. R. 1978. Numerical and analytical solutions of moving boundary problems. Proceedings of Symposium on Moving Boundary Problems edited by D. G. Wilson, A. D. Solomon and P. T. Boggs, New York: 129-145.
5. Shi M. 2005. On the solution of a one-dimensional Riccati equation related to risk-sencitive portfolio optimization problem. Repts Fac. Sci. and Eng. Soga Univ. Math. 34 (1): 17-24.
6. Bakhvalov N. S. 1970. Numerical methods. M., Nauka, 632.
7. Berezin I. S., Zhidkov N. P. 1960. Methods of calculation: in 2 v. M., Fizmatgiz.
8. Fil'chakov P. F. 1970. Numerical and graphical methods of applied mathematics. Kiev, Naukova dumka, 800.
9. Erugin N. P. 1967. Analytical theory and the real problems of the theory of differential equations associated with the first method and the analytic theory. Differentsial'nye uravneniya, 3 (11): 1821-1864.
10. Kondratenya S. G., Prolisko E. G., Shilo T. I. 1988. On the question of the existence of polar solutions of first order differential equations. Differentsial'nye uravneniya, 24 (10): 1824-1826.
11. Lukashevich N. A. 1995. The simplest differential equation of the third order of P-type. Differentsial'nye uravneniya, 31 (6): 955-961.
12. Matatov V. I., Sabynich L. V. 1991. On mobile features of autonomous Hamiltonian systems. Vestnik Belorusskogo universiteta. Seriya 1, Minsk, 8. Dep. v VINITI 09.04.91, № 1532-V91.
13. Samodurov A. A. 1983. On the integrability of Abel differential equation in parametric form. Vestnik Belorusskogo universiteta. Seriya 1, 2: 57-59.
14. Chichurin A. V., Shvychkina E. N. 2014. On the integrability of a third-order system, the equivalent equation of Chazy six fixed poles. Vestnik Chuvashskogo gosudarstvennogo pedagogicheskogo universiteta im. I. Ya. Yakovleva. Seriya: Mekhanika predel'nogo sostoyaniya, 4 (22): 176-187.
15. Yablonskiy A. I. 1964. The asymptotic expansion of the right solutions of certain classes of differential equations. Doklad AN BSSR, 8 (2): 77 80.
16. Orlov V.N. 2006. Criteria for the existence of mobile singular points of solutions of differential equations of Riccati. Vestnik Samarskogo gosudarstvennogo universiteta. Estestvenno nauchnaya seriya, 6/1 (46): S. 64-69.
17. Orlov V.N. 2008. The approximate solution of the first Painleve equation. Vestnik KGTU im. AN. Tupoleva, 2:
42-46.
18. Orlov V.N. 2008. A method for the approximate solution of matrix differential Riccati equations. Vestnik MAI, 15 (5): 128-135.
19. Orlov V.N. 2009. A study of the approximate solution of Abel differential equation in the neighborhood of a movable singular point. Vestnik MGTU im. N. E. Baumana. Seriya: Estestvennye nauki, 4 (35): 102-108.
20. Orlov V. N., Leont'eva T. Yu. 2013. The influence of perturbation of the initial data on the approximate solution of a nonlinear differential equation of second order in the analyticity region. Vestnik Chuvashskogo gosudarstvennogo pedagogicheskogo universiteta im. I. Ya. Yakovleva. Seriya Mekhanika predel'nogo sostoyaniya, 3 (17): 103-109.
21. Orlov V. N., Leont'eva T. Yu. 2013. Construction of the approximate solution of a nonlinear second-order differential equation in the region of holomorphy. Vestnik Chuvashskogo gosudarstvennogo pedagogicheskogo universiteta im. I. Ya. Yakovleva. Seriya: Estestvennye i tekhnicheskie nauki, 4 (80): 156-162.
22. Orlov V. N., Leont'eva T. Yu. 2014. Construction of the approximate solution of a nonlinear second-order differential equations in the neighborhood of a movable singular point in a complex region. Vestnik Chuvashskogo gosudarstvennogo pedagogicheskogo universiteta im. I. Ya. Yakovleva. Seriya: Mekhanika predel'nogo sostoyaniya, 4 (22): 157-166.