ВЕСТНИК ОШСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА МАТЕМАТИКА. ФИЗИКА. ТЕХНИКА. 2024, № 1(4)
УДК: 517.928
DOI: https://doi.org/10.52754/16948645 2024 1(4) 6
О МЕТОДЕ ГАЛЕРКИНА ПОСТРОЕНИЯ ПЕРИОДИЧЕСКИХ РЕШЕНИЙ КВАЗИЛИНЕЙНОЙ ИНТЕГРО-ДИФФЕРЕНЦИАЛЬНОЙ УРАВНЕНИИ ВТОРОГО
ПОРЯДКА
Алымбаев Асангул Темиркулович доктор ф.-м.н., профессор
asangul1952@gmail. сот Государственный университет им. И. Арабаева
Бишкек, Кыргызстан Бапа кызы Айнура, магистр математики, ст.преподаватель
abapakyzy@gmail. сот
Иссык-Кульский государственный университет им.К.Тыныстанова
Каракол, Кыргызстан
Аннотация. В статье рассматривается задача построения и установления существования, приближенного 2п- периодического решения квазилинейного интегро-дифференциального уравнения второго порядка с периодической правой частью. Решение ищется в виде ряда Фурье. Построена система квазилинейных алгебраических уравнений относительно ряда Фурье. Доказано однозначной разрешимости системы алгебраического уравнения, что равносильно существованию 2п- периодического решения квазилинейного интегро-дифференциального уравнения второго порядка. Оценена погрешность разности между точным и приближенными решениями.
Ключевые слова: Метод Галеркина, 2п- периодическое решение, квазилинейное интегро-дифференциальное уравнение второго порядка, точное и приближенное решение.
ЭКИНЧИ ТАРТИПТЕГИ ИНТЕГРО-ДИФФЕРЕНЦИАЛДЫК ТЕЦДЕМЕНИН МЕЗГИЛДИК ЧЫГАРЫЛЫШЫН ГАЛЕРКИНДИН МЕТОДУ МЕНЕН ЧЫГАРУУ
Алымбаев Асангул Темиркулович, ф.-м.и.доктору, профессор
asangul1952@gmail. com И.Арабаев атындагы КМУ Бишкек, Кыргызстан Бапа кызы Айнура, математиканын магистри, ага окутуучу
ahapakyzyagmail. com К. Тыныстанов атындагы ЫМУ Каракол, Кыргызстан
Аннотация. Макалада экинчи тартиптеги квазисызыктуу интегро-дифференциалдык тецдеменин 2п — мезгилдYY жакындаштырылган чыгарылышын табуу жана анын жашашын далилдвв маселеси каралат. Чыгарылыш Фурьенин катары тYPYндв изделет. Катардын коэффициенттерине карата квазисызыктуу алгебралык тецдеме тYЗYЛYп, анын бир маанилYY чечилиши далилденет жана тыянак катары квазисызыктуу интегро-дифференциалдык тецдеменин Галеркиндин ыкмасы боюнча тургузулган, 2 п — мезгилдик чыгарылышынын жашашы далилденет. Так жана жакындаштырылган чыгарылыштардын ортосундагы айырманын чени аныкталат.
Негизги свздвр: Галеркиндин методу, 2 п —мезгилдYY чыгарылыш, экинчи тартиптеги квазисызыктуу интегро-дифференциалдык тецдеме, так жана жакындаштырылган чыгарылыштар.
ON GALERKIN'S METHOD FOR CONSTRUCTING PERIODIC SOLUTIONS OF A
SECOND-ORDER QUASILINEAR INTEGRO-DIFFERENTIAL EQUATION
29
Alymbaev Asangul Temirkulovich, Doctor of Ph. and Math. Sc., Professor
asangul1952@gmail. com State University named after I. Arabaev Bishkek, Kyrgyzstan
Bapa kyzy Ainura, Master of Physics and Mathematics Education, Senior Lecturer
abapakyzy@gmail. com Issyk-Kul State University named after K. Tynystanov
Karakol, Kyrgyzstan
Annotation. The paper considers the problem of constructing and establishing the existence of an approximate 2n-periodic solution of a second-order quasilinear integro-differential equation with a periodic right-hand side. The solution is sought in the form of a Fourier series. A system of quasilinear algebraic equations with respect to the Fourier series is constructed. The unique solvability of a system of an algebraic equation is proved, which is equivalent to the existence of a 2n-periodic solution of a second-order quasilinear integro-differential equation. The error of the difference between the exact and approximate solutions is estimated.
Keywords: Galerkin method, 2n-periodic solution, second-order quasilinear integro-differential equation, exact and approximate solution.
Изучение периодических решений сильно нелинейных систем (систем, не содержащих малого параметра) классическими асимптотическими методами не всегда возможно. Для таких систем разработаны ряд аналитических и численно-аналитических методов. Среди существующих методов, имеются методы наряду с доказательством теорем существования периодических решений позволяет построить этих решений. Одним из подобных методов является метод Бубнова-Галеркина, получивший достаточно полное обоснование для нелинейных систем обыкновенных дифференциальных уравнений в работе М.Урабе [4] и Л.Чезари [7]. В дальнейшем идеи работы М.Урабе распространены для изучения периодических решений системы дифференциальных уравнений с запаздыванием, интегро-дифференциальных уравнений с конечным и бесконечным последействием, а также системы автономных дифференциальных и интегро-дифференциальных уравнений с малым параметром [1,2,3,5,6]. Настоящая работа посвящена задаче построения и доказательству теоремы существования -периодического решения квазилинейного интегро-дифференциального уравнения второго порядка с конечным последействием.
Рассмотрим 2 ^-периодическое по t и дифференцируемое по x,u функцию f ( t,x,u) в области ( t,x,u)eTxD1xD2, где Т = (—оо,+оо ) ,D1,D2 - ограниченные области в R = (—о, + оо ) . Введем нормы:
lf(t,x,u)l0= max \\f(t,x,u)\\,
J XD;^ X.D2
1
2n
\\f (t,x,u)\\o =
2i\f\:
dt
Пусть - д ( t) непрерывная 2n -периодическая и разлагающая в ряд Фурье функция
со
д (t) = с0 + V2^ (ск с os к t + d ks ink t). ( 1 )
k=1
Введем на множестве периодических функций операторы , такие, что
^тд ( ° = с0 + ^ J( CkC0Skt + dks ink 0, (2 )
I
k=1
l(m + l)4 + (m + 2 )4 +
, , V2 V2 V2
ffi(m) = 7——tj , 7—
(m + 1 у (m +1)^ m^
Рассмотрим интегро-дифференциальную уравнению второго порядка
(5)
(0 = V2 ^ ( cfccо skt + d fcsink t). ( 3 )
k=m+l
Заметим, что
Rmg(t) = g(t) — Smg(t). Рассмотрим дифференциальное уравнение второго порядка вида
d2x(t)
= (4)
где б ( t) — 2 п- периодическая функция, представимое в ряд Фурье вида (2). С учетом (1) уравнение (4) записываем в виде
СО
d2x(t) /-W ч
—¡т;— = с о + V 2 > ( cfcc о s к t + dfcs ink t) .
k=l
Теорема 1. Пусть x = x ( t) — 2 п периодическое решение (5). Если
со со
со = 0, x(0) = — V^. ^ = ^
к=1 fc=l тогда периодическое решение представимо в виде
со
x ( t) = V2 ^ — (—cfc с os k t + d fcs ink t) . fc=l
CO
йшх ( t) = x ( t) — 5шх ( t) = — V2 ^ — (cfcc о s k t + d fcs ink t) .
k=m+1
Теорема 2. Для разности x ( t) — x?n ( t) - справедлива оценка
|x( 0 — S,„x( 01о < k( 0 — Xш ( OU < о" (жНЯо ,
||x( t)— x„,( t)y о <ffi(™)11Л1о
где
2 2 i1^
<j(m) =
^ = Л x ( t ) + / ( t, x ( t) , T ( ( t, s, x (s ) ds ) , (6)
где А-вещественное число, /(t, x, и( t) ), ( ( t, s, x) -периодическое по t, s с периодом 2 п функции,
Предполагается, непрерывно-дифференцируемость ( ) - по своим
переменным. Периодическое решение уравнения (6) ищем в виде
m
xjn ( t) = ао + V2 afcc os k t + bfcs ink t). ( 7 )
fc=i
Поставляя (7) в (6) получим равенство
771 771
V2 —k2afcc о s k t — k2 bfcs ink t ) = Л ао + V2 ^(Л afcc о s k t + Л bfcs ink t) + fc=i fc=i
lit
+A (m) + VlSiA (т cos kt + в(т s ink t), ( 8)
k=1
где
271
A0 ) = 2^12 j f (t,xт( ^ ,um( ^) d^ 0
2n
A^ =—1=~ I f(t,xm(t),um(t)) сosktdt, 2 V 2 n J
о
2n
B(m->=—-=- I f(t,xm(t),um(t))sinktdt, ,M=k.B
t-T
Приравнивая коэффициенты при одинаковых гармониках функции получим из (8) систему алгебраических уравнений относительно неизвестных а0, ак, Ьк разложении (7).
Аа0 + 0 ■ ак + 0 ■ bk + А¡,т) = О, О ■ а0 + (А + к2)а0 + О ■ Ьк + А¡¿т) = О, ,0 ■ а0 + О ■ ак + (А + к2)Ък + = 0.
или в матричной форме
'А О 0 \ /aoN
/4т)\
О A + k2 0 )( ак ) + vO О А + к2J \bk
A^) = 0 ,
Отсюда получим
где
(А О О
D m = (0 A + k2 0 ), F (m)( а) = vO 0 А + к2,
\<7
D(m)a + F(m)(а) = 0, (9)
/4т)\
A(m ), а = (а0^. k = 1,т. ув(т) J \bkJ
Обозначим через х = х (t), 2n — периодическое решение уравнения (6) и представим его в виде ряда Фурье
со
х (t) = а0 + ак с os k t + bks ink t). (1 0 )
k=1
Поставляя (10) в уравнение (6) и с учётом свойства оператора Sm, получим следующее равенство
d2x (t) d2 хт (t) if 1
Sm dt2 d12 AXm( ^ + Sm f I t,xxm( 0 , I t-,s,xm(s^) ds I +
t-T
/
\
/(г,Х (0, |<р(г,5,Х (5) )(5) -5т/( ^Хт(г), |<р(г,5,Хш (5) )(
Г—Т \ Г—Т
Отсюда имеем
О + (ш)( й) =
/ '
=
V
/( г,Х ( 0 , г, 5 ,Хт (5 ) ) (¿5 )-/( ( 0 , ^5,Хт (5 ) ) (5 ) ) = -р (т) ,
С-т
С-Т
где а = ( хо , х^, Ь^), к = 1, т . Представим разность
/ ( г, X ( О, /_т <р(г, 5 , X ( 5 ) ) ( 5 ) - / ( г, Хш ( О, /_т <р(г, 5 , Хш (5)) ( 5)
в виде
/( г,х ( 0 ,/_ т<р( г,5 ,х(5)) (5 ) - /( гДш ( 0,/_<?(г,5,Хш (5 )) (5) = = /( г,х ( г),/_ г,5,х (5 )) (5 ) - /(г,йп ( г),/_ г,5,х (5 )) (5) +
+/(г,Хш (г), |<р(г,5,х (5) )(5^-/(г,Жш(г), |<р(г,5,Хш (5) =
а/ ( г, Хш + 0!(х - «ш) , /_ т <Р ( г, 5, /_ т <?( г, 5, х (5 ) ) (5))
дх
( а а )
а/( г,хш( г) ,й + б2(и- йп)) г а <р( г,5,х(5 ) + б2(х (5)- хш(5) ) )
ди
Оценим разность: Оценим разность:
/
С-Т
5х
х х
О < вг < 1, 0 < 02 < 1, 0 < 03 < 1.
|5ш (/( г,* ( г) ,й ( г) ) -/( г,Хш ( г) ,йш ( г) ))|
<
д/( + 6>1(:х ^тХм ( О)
9х
|х(г)
х |
| д /( С, хт ( С) , и + 02 (м _ % ) )
ди
С-т
д <¡0 ( С, 5, :й (я) + 03 (:й (я) _ хт (я)})
дх
|Х (5) - Хш(5)|о(5 . ( 1 1 )
Так как, согласно теореме 2, для разности й( г) - хш ( г) имеют место оценки: |х(г) х ( г)|о<а (т)|/)|о, ||х ( г) - Хш( г)|| <0! (т)||/||о,
С учетом этих оценок, из (2.4.6) получим
||р(т)|| < а (т)[|/|о1/1 1 + |/111111/1о= а (т)|/|о1/к( 1 + ^кт) (12)
или
||р(т)| <а1(т)|/|о1/11 ( 1 + |<ркт). Представим уравнение (9) в виде
а + (О (ш))_ V (ш)( а) = 0, и положив а = <х, отсюда имеем
х ( ) х ( )
из которого следует оценка
(13) (14)
а + (Б(т))-V(т)(а) < К\\р(т)\\ < а(т)К^^Щ^ 1 + 1р^т) (15)
или
а + (Б(т))~V(т(а) < К\\р(т)\\ < аг(т)К\П\01Пг( 1 + 1р^т), (16) (БЫ)-.
Решаем уравнение (14) методом последовательных приближений
^+1 = - (Б(т)) ~ V(т) (ак), к — 0,1, 2..........(17)
За начальное приближение берем число а0 — а т.е коэффициенты частной суммы
где К —
ряда
а
(1) = а0+^2у (аксозкt + Ькзтк^.
к=1
Докажем равномерную сходимость итерационного процесса (17). Представим разность а к+1 — а к в виде
ак^ 1 — ак — —(Б(т))'1 (р(т) (ак) — Р(т) (а^)) —
= —(Б (т)) - 1дд Р(т)(ак - ^ (ак — ак . О )(ак—ак_ О.
Предположим, что существует достаточно большое число т0, так ое, чт о п р и т > т0 выполняется условие
(18)
дР(т\а)
да
X
< —, при 0 < х < 1-
(19)
Тогда из (18) следует оценка
\\ак+1 — ак\\< (Б(т)) - 1
<К — 1ак — ак - 11 <х№к — ак-11.
К1
дР (т)(ак - 1 + в( ак — ак -1) )
да
\\ак — ак -1\\ <
Отсюда получим
\\ак+1 — ак\\ < х\\ак — ак - 1У < Х2\\ак -1 — ак-2\\ < • • • < Хк\\а1 — ао\\. (20) Оценим разность а 1 — а0 , положив а0 — а. Из (17) при к — 0 имеем
а1 — а0 — — *а0 + (Б(т)) -1Р (т)( а0 ) + — — [а + (Б(т)) - 1Р (т)(а) + — —(Б(т)) - 1р (т).
Отсюда с учетом неравенства (15) или (16) получим
\\а1 — а0\\ < ||(Б(т)-\\р(т)\\<а(т)КШП1(1 + 1р^т) (21)
или
\\а1 — а0\\<а1 ^К^У/к(1+1Р^т), при т - оо. (22)
С учётом неравенств (21), (22), получим из (20) оценку \\ак+1 — ак\\ < Хка1 (т)К1/^1/11 (1 + 1р^т), для к — 0,1,2,... Нетрудно показать справедливости выполнения неравенства
X
\\ак+Р — ак\\ < 1 ОтН/иЯ^ 1 + Ы1Т).
Отсюда, при получим
\\а — а к\\ <
X*
1-Х
а 1( т) К\\/\\01/11( 1+1р^т),
(23)
которая обеспечивает равномерную сходимость последовательности (17) к решению уравнения (14) т.е. а = — Так как — = ( -0 , -bfc), к = 1, m, то приближенное решение интегро-дифференциального уравнения (6) согласно по методу Галеркина записывается в виде
хт ( t) = -0 + V2 Z( -kc os к t + -ks ink t) .
Из (23) при к = О
fc=i
,, Л ^ ffi(m)Wllo^li ( 1 + lyliT)
На - «У <-:-. ( 2 4)
1
Оценим разность ||xm ( t) — хт ( t)||0 :
2.71
dt =
ухт ( 0 — хт ( t)yo = -1 J ухт ( 0 — хт ( t)||2d t =
о
27Т , т ч 2
= — J (a 0 + VlJ ( ( âk — a k) с o s к t + ( bk — -k)s i n к t) ) о V k=l /
7TL
= »■Z«! +1J (V2 (âk — a-k) 2 + V2(ik — âk)2} = +
fc=l
m
+ 2 Z (y^2âk — Vl-klf + yvibk — V2-k|2) = 1 llaâk — -J2.
k=1
Следовательно
1
Ухш ( г) - Хш ( г)|о = Уй - аУ,
то с учетом неравенства (24) получим оценку
Ух ( г)-х- ( г)| ^|Л|о|/|1 (
И^ш ( г ) Хш ( г )|о < /77-. . .
V 2 ( 1 -х)
Так как, согласно теореме (2)
|х ( г) - Хш ( г)|о < а ОтНЯо
а также
|х ( г) - йш( г)|о < а1 (т)|/|о, тогда для разности | х х | справедлива оценка:
|х ( г) - Хш( г)|о = |х ( г) - 4 (г) + 4( г) - Хш( г)|о < |х ( г) - 4( г)|о +
+ |Хш ( г) - Хш ( г)|о < < а (тпЛЯо +--. (25)
\[2 л/2
Поскольку а (т) < —, а-^т) < — и ||/||о < |/|о то из (25) следует неравенства
|х ( г) - хш( г)и < т | Л о ( 1 + ^ ( 1-х) ) =
= 1 - х) + К-|/| 1 ( 1 + М^) |/|о[У2 + /ТЩ 1 + |«р|1т)]
т 2|/ |о VI ( 1-х) " т 2( 1-х) ,
при
Таким образом доказано утверждение.
Теорема 3. Пусть интегро-дифференциальное уравнение (6) имеет 2п — периодическое решение а(Ь) и удовлетворяет следующим требованиям: а) выполняется требование теоремы 2;
б) Il â + (D (т)) 1 F (т)(â)
Л 5F(m)(a)
в)
<ai(m) K\\f\\0lfli( 1+Whj) ;
да
К
ai(m) < Щ.
Tïl
о < X < 1.
Тогда, алгебраическое уравнение (6) имеет единственное решение а — (а0,а1 ,Ь1 ,• ■ • ,ат,Ьт) , такое, что между точным а( Ь) и приближенным решением хт (Ь) найденного методом Галеркина, справедлива оценка
|х (t)-Xm(<
|f|0[V2 + K|f| i ( l+Wh-c) ] rn2(l-/)
при m -> 00.
Литература
1. Алымбаев А.Т. Метод Бубнова-Галеркина построения периодических решений автономных систем интегро-дифференциальных уравнений //А.Т.Алымбаев// Математическая физика. - Киев, 1993. - вып.22. -С.3-7.
2. Алымбаев А.Т. Численные, численно-аналитические и асимптотические методы исследования краевых задач: Монография // А.Т.Алымбаев. - Бишкек, 2015. - 205 с.
3. Алымбаев А.Т. О методе гармонического баланса построения периодического решения системы автономных интегро-дифференциальных уравнений с бесконечным последействием // А.Т. Алымбаев, Бапа кызы Айнура // Алатоо academic studies. -Бишкек, 2022. - №2. - С.459-464.
4. Урабе М. Метод Галеркина для нелинейных периодических систем. - Москва, 1966. - 97, №3. - С.3-34.
5. Самойленко А.М. Метод Бубнова-Галеркина построения периодических решений интегро-дифференциальных уравнений типа Вольтерра // М.А.Самойленко., О.Д.Нуржанов// Дифференц.уравнения. - 1979. -15, №8. - С.1503-1517.
6. Мартынюк Д.Н. Метод Бубнова-Галеркина отыскания периодических и квазипериодических решений систем с запаздыванием. //Д.Н.Мартынюк// В кн.: Тр. Второй конф. «Дифференциальные уравнения и применения». - 1982. - ч.2. - С.445-448.
7. Cesari L. Functional analysis and periodic solutions of nonlinear differential equations// L. Cesari//In: Contributions to differential equations. - New-York, 1963,vol.1. - Р.145-197.