CC BY
9
УДК 517.928.4
ПОСТРОЕНИЕ ПРИБЛИЖЕННОГО РЕШЕНИЯ ОДНОГО НЕЛИНЕЙНОГО ДИФФЕРЕНЦИАЛЬНОГО УРАВНЕНИЯ ВТОРОГО ПОРЯДКА
В ОБЛАСТИ ГОЛОМОРФНОСТИ
APPROXIMATE SOLUTION FOR A NONLINEAR DIFFERENTIAL EQUATION OF THE SECOND ORDER IN THE REGION OF HOLOMORPHY
В. Н. Орлов, Т. Ю. Леонтьева V. N. Orlov, T. Y. Leontyeva
ФГБОУВПО «Чувашский государственный педагогический университет им. И. Я. Яковлева», г. Чебоксары
Аннотация. В работе рассматривается нелинейное дифференциальное уравнение второго порядка с подвижными особыми точками. Дано доказательство теоремы существования и единственности решения этого уравнения в области голоморфности, построено приближенное решение в случае точного значения начальных условий.
Abstract. The article considers the nonlinear differential equation of the second order with movable special points. The proof of the theorem of existence and uniqueness of the solution for this equation in the region of holomorphy is given and the approximate solution in case of exact value of entry conditions is made.
Ключевые слова: подвижная особая точка, нелинейное дифференциальное уравнение второго порядка, приближенное решение, область голоморфности.
Keywords: movable special point, nonlinear differential equation of the second order, approximate solution, region of holomorphy.
Актуальность исследуемой проблемы. Уравнения, содержащие неизвестные функции и их производные в степени выше первой, называются нелинейными. В последние годы они привлекают все большее внимание. Например, уравнения физических явлений обычно линейны лишь в первом приближении; дальнейшее и более точное исследование, как правило, требует использования нелинейных уравнений. Решения нелинейных уравнений связаны с большими трудностями, вызванными наличием подвижных особых точек у интегралов этих уравнений, которые и являются препятствием к использованию известных на данный момент приближенных методов решения.
Материал и методика исследований. В данной работе используется метод решения нелинейных дифференциальных уравнений, предложенный в работах [2] и [3], который в целом состоит из шести этапов. В статье [4] был рассмотрен первый шаг вышеуказанного метода, который упоминается и в работах [1], [5].
Результаты исследований и их обсуждение. Рассматривается нелинейное дифференциальное уравнение второго порядка с полиномиальной правой частью пятой степени в комплексной области:
У "(2) = Ь0 (2)У 5 (7) + Ъ, (2)У 4 (2) + Ъ2 (2)У 3 (7) + 63 (2)У 2 (2) + Ъ4 (2)у (2) + Ь5 (7) , (1)
где Ь, I = 0,1,...,5 - голоморфные функции в рассматриваемой области. С помощью замены переменной
У(2) = 4^ -А*2! (2)
У( ' 5Ъ0(2)
при выполнении условии
2
Мz) b2(z) Ьз(Z) _
п / п ^2^
bo (z) ч bo(z) ,
b4(z) + b0(z) 5
4bo( z) 16
5bo(z) 2bi(z) b2(z) 63 (z) (3)
уравнение (1) приводится к нормальной форме [5]:
u"( z) = u 5( z) + r( z) (4)
где
, ч , . . ft , ч ._,__ ft „ ч . . . f
'5/
()=_ bi5 (zХ/6Д7) 3bo (zЩТ) bi (z2bo (z)bi (z
^ = 55bo4(z) 20bo2(z) 5bo(z) 5bo2(z)
^0
» N 2
'0 (z) I b1'
+ b + 2l bo(z )J ^ ^ 1*0 (z)J b- + bo"(z) u ,(z) (5)
+b5 (zX,b° (z)+-5Щ---lib3«--u(z) (5)
в каждой области, в которой b0 (z) Ф 0 . Рассмотрим задачу Коши:
У"(z) = y5(z) + r(z), (6)
У( z0) = У0> У'( z0) = У1. (7) Используя идею, предложенную в работах [2], [3], докажем теорему 1. Теорема 1. Пусть выполняются следующие условия: 1. r(z) е C1 в области
|z "z0\<Р0, (8)
где р0 = const;
\r(n}( z0)\
3M0 : J-1 <M0, где M0 = const, n = 0,1,2,....
n!
Тогда решение задачи (6)-(7) является голоморфной функцией
ад
У(2) = £С„ (2 - 2о)" (9)
п=0
в области
¡2 - 2^ , (10)
где
р =
1
М = тах
Ы^Р"
г (п)( г,
( 2 о)
п!
1 22(М +1)4 Доказательство. По условию теоремы имеем
ад
г(2) = Х Ап (2 - 2о)п .
п=0
Подставив выражения (9) и (11) в уравнение (6), получим
X Сп (2 - 2о) п = X СП (2 - 2о) п Ап (2 - 2о ) П .
V п=о у V п=о у п=о
После выполнения соответствующих операций получаем
ад ад
ХСпп(п - 1)( 2 - 2о) п-2 =Х(^п" + Ап )(2 - 2о)п ,
, п = о,1,2,....
(11)
(12)
п п п
где В** = X В* С. , В* = X В В , В = XС С , п = о,1,2,....
^ п / > п-г г ' п / , п-г г ' п / > п-г г ' ? ? ?
г=о
г=о
г=о
Равенство (12) обратится в тождество при условии, что
п(п - 1)Сп = В*-2 + Ап-2 , п = 2,3,... . (13)
Данное соотношение позволяет однозначно определить все коэффициенты Си . Таким образом, получено формальное представление решения задачи (6)-(7) в некоторой окрестности точки ^ в виде (9), единственность которого следует из однозначности определения С .
Покажем сходимость полученного ряда в области (10). Положим
М = тах
Ы^ир
\г (п)( 2о)|
п!
>, п = о,1,2,...,
(14)
что возможно в силу п. 2 данной теоремы. Тогда из выражения (13) с учетом (14) получаем оценку для коэффициентов С :
1
-2М(М + 1)4п, п = 2,3,....
(15)
п(п -1)
С помощью метода математической индукции докажем эту оценку в случае п +1 = 2k + 1:
2к-1
2к-1
2к(2к + = ВЦ-! + Аи-! = x в^-д. + А^ = x x в*-1-,-а с, + Аи-! =
2к-1-г
Л
2к-1^2к-1-г (2к-1-г-1 ]
= x xI x С2к-1-г-]-тСт II x С1 -1С1
г=о V 1=о V т=о /V 1=о
откуда с учетом оценки (15) следует, что
о V 1=о
Сг + А2к-1
5
со
со
со
п=2
п=о
п
о
f
2k+11 —
Y 2k-1Í2k-1-i í2k-1-i-j 22(2k-1-i-j-m) 2 /д^ j\4(2k-1-i-j-m) Л
-У У У -(-^-x
2k(2k +1) У | y (2k - i - j - m - 1)(2k - i - j - m - 2) J
f j 22^j-г2
y 2 M 2(M + 1)4(j-г) ^
1у (j - i)*((1 - /)* -1)
JJ
22iM (M +1)4 i *(¿ * -1)
■ +
+
M
1 2k-1Í 2k-1-i, =-1-УI У ^22(2k-1-i-j)M2 (M + 1)4(2k-1-i-j)
2k(2k +1) 2k(2k +1) У |
2k -1-1 - j ( x у
m=0
1
1
(2k - i - j - m - 1)(2k - i - j - m - 2) m * (m * -1)
v
JJ
Г
22 }M 2(M +1)41У
1
1
2 M (M + 1)4i
x-^-— +
i (i -1)
2k-1Í 2k-1-i^ 2k-1-i- j
1=0 V(j - /)*((/ - O* -1) l (l -1)
M
JJJ
2k (2k +1)
<22(2k-1) M 5(M + 1)4(2k-1)
2k (2k +1)
i=0
У У I У
v i=0
1
1
Vv
m=0
(2k - i - j - m - 1)(2k - i - j - m - 2) m* (m* -1)
1
1
. l =0
d-l) (Ü-l) -1) l (l -1)
1
■ + ■
1
JJJ
i (i -1) 2k(2k +1)
-M —
—
1
2k (2k +1)
где (2k -1 - i - j - m)* =
22(2k+1 M (M + 1)4(2k+1,
1, (2k -1 - i - j - m) = 0,1
(2k -1 - i - j - m),(2k -1 - i - j - m) = 2,3,...'
! . V I -1 -* - -Я = 0,1 * I 1,т = 0,1
(2к -1 - г - я) = { , т = { ,
[(2к -1 - г - Я),(2к -1 - г - Я) = 2,3,../ [т, т = 2,3,...
(2к -1 - г)* = { ',(2к -1-0 = , /• ={ = 01 , Я* = { = ,
[(2к -1 - г),(2к -1 -1) = 2,3,... {/,/ = 2,3,... [я,я = 2,3,...
( я -1)* ={ 1(Я -1) =04 , ={11 =0Д.
(Я ) {(Я -1 ),(Я-I) = 2,3,..., {1,1 = 2,3,...
Аналогичное соотношение имеем и в случае п +1 = 2к. Таким образом, подтверждаем оценку для коэффициентов Сп ряда (9). Рассмотрим ряд
1
n=2
n(n - 1)
22n • M(M + 1)4n|z - z0|n
(16)
x
x
x
x
x
x
x
x
oo
который является мажорирующим для ряда (9). На основании признака Даламбера заключаем сходимость ряда (16) в области |2 - 20| < —1-— . Тогда, положив
р = ттр
1
22(М +1)4
22(М +1)4
', получаем сходимость ряда (9) в области (10), что доказывает
нашу теорему.
Полученные в данной теореме оценки позволяют построить приближенное решение задачи (6)-(7):
N
УЫ (2 ) = XCn (2 - 2о) п .
(17)
Теорема 2. Пусть выполняются пп. 1 и 2 теоремы 1, тогда для приближенного решения (17) задачи (3)-(4) в области (10) справедлива оценка погрешности
2Ж+2 • М(М +1)4Л+4|2-2,
4 N+41
aУN (2) <
М+1
1
N (N +1)
1 - 22 • (М +1)4 • 2 - 2о
(18)
где р = тт<!ро,-
1
М = тах
1 22(М +1)4 j Доказательство. По определению
ДУN (2) = |у(2) - УN (2) =
|УоИиР-
\г (п)( 2,
( 2о)|
п!
п = о,1,2,
ад N
XCn (2 - 2о)п -XCn (2 - 2о)
п=о
п=о
XCn (2 - 2о)п < X|Cn
п=^1 n=N+1
Учитывая оценки для Сп из теоремы 1, имеем
• 2 - 2п
ад -1
AУN(2) < ^ ——-22п • М(М +1)
n=N+\п(п 1)
+ 1)4п ¡2 - 2ПГ <■
1 -. у-^-. у- .-,4N+41 N+1
N (N +1)
•22"+2 • М(М +1)4"+4 2 - 20
1
1 - 22 • (М + 1)4 • \2 - 20| Пример. Найдем приближенное решение задачи (6)-(7) в случае г(2) = о для
2-\/з 2л/э
начальных условий у(о,5 + о,5г) = 1 +1 и у'(о,5 + о,5г) =--1--г .
Данная задача имеет точное решение у =
л/3
1 + Г,+* I 2 • г - 2 2
Найдем радиус
аналитичности р = о,оо7359312. Выберем значение 2 = о,5о 1 + о,5о1г из области голоморфности. Применяя (17), N = 3, вычислим приближенное значение функции. Произведенные расчеты приведены в таблице 1,
п=о
п
п
п
X
X
где у3 - приближенное решение (17); у - значение точного решения;
Ду3 - оценка погрешности приближенного решения (18);
Ду - истинная величина погрешности приближенного решения у3;
А У - апостериорная оценка погрешности, которая определяется путем решения обратной задачи теории погрешности. На основании апостериорной оценки убеждаемся в том, что в структуре приближенного решения (17) значение N = 5. Добавки в структуре приближенного решения для N = 4 и N = 5 в общей сумме не превышают требуемой точности. Поэтому в структуре приближенного решения можем ограничиться значением N = 3, при котором приближенное решение будет иметь погрешность ^ = 0,00001.
Таблица 1
1 Уз У ¿Уз Л.У Л1У
0,501+0.5011 1,000004015+1,0023054011 0,999995985+1,0023134011 0,00019894 0,000011335 0,00001
Резюме. В статье доказана теорема существования и единственности решения рассматриваемого нелинейного дифференциального уравнения второго порядка, что позволяет построить приближенное решение для данного уравнения в области голоморфности.
ЛИТЕРАТУРА
1. Орлов, В. Н. Аналитическое приближенное решение одного нелинейного дифференциального уравнения в комплексной области / В. Н. Орлов, М. П. Гузь // Вестник Чувашского государственного педагогического университета им. И. Я. Яковлева. Серия: Механика предельного состояния. - 2012. - № 2 (12). - С. 75-82.
2. Орлов, В. Н. Метод приближенного решения дифференциального уравнения Риккати / В. Н. Орлов // Научно-технические ведомости Санкт-Петербургского государственного политехнического университета. -2008. - № 4. - С. 102-108.
3. Орлов, В. Н. Об одном методе приближенного решения матричных дифференциальных уравнений Риккати / В. Н. Орлов // Вестник Московского авиационного института. - 2008. - Т. 15. - № 5. - С. 128-135.
4. Орлов, В. Н. Построение приближенного решения нелинейного дифференциального уравнения в области аналитичности / В. Н. Орлов, Т. Ю. Леонтьева // Фундаментальные и прикладные проблемы механики деформируемого твердого тела, математического моделирования и информационных технологий : сб. ст. по мат-лам международ. науч.-практ. конф., Чебоксары, 12-15 августа 2013 г. Ч. 2. Математическое моделирование и информационные технологии. - Чебоксары : Чуваш. гос. пед. ун-т, 2013. - С. 47-52.
5. Орлов, В. Н. Приближенное решение одного нелинейного дифференциального уравнения в области голоморфности / В. Н. Орлов, А. З. Пчелова // Вестник Чувашского государственного педагогического университета имени И. Я. Яковлева. - 2012. - № 4 (76). - С. 133-139.
