Научная статья на тему 'ПРИМЕНЕНИЕ МЕТОДА КОНЕЧНЫХ РАЗНОСТЕЙ К РЕШЕНИЮ ЗАДАЧИ ДЛЯ УРАВНЕНИЯ КОЛЕБАНИЯ ПРОДОЛЬНО НАПРЯЖЕННОГО СТЕРЖНЯ С ИНТЕГРАЛЬНЫМИ УСЛОВИЯМИ'

ПРИМЕНЕНИЕ МЕТОДА КОНЕЧНЫХ РАЗНОСТЕЙ К РЕШЕНИЮ ЗАДАЧИ ДЛЯ УРАВНЕНИЯ КОЛЕБАНИЯ ПРОДОЛЬНО НАПРЯЖЕННОГО СТЕРЖНЯ С ИНТЕГРАЛЬНЫМИ УСЛОВИЯМИ Текст научной статьи по специальности «Философия, этика, религиоведение»

CC BY
26
6
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
Разностная задача / интегральные условия / аппроксимация / погрешность аппроксимации / нелокальные условия

Аннотация научной статьи по философии, этике, религиоведению, автор научной работы — Ханкишиев Закир Фарман Оглы

В работе рассмотрена одна задача для уравнения колебания продольно напряженного стержня с интегральными условиями. Известно, что колебание продольно напряженного стержня выражается дифференциальным уравнением в частных производных четвертого порядка. Для таких уравнений ставятся различные задачи, в том числе задачи с интегральными условиями. При наличии интегрального условия в задаче, в ее решении возникают определенные трудности. В настоящей работе после замены интегральных условий нелокальными граничными условиями, к решению полученной новой задачи применяется метод конечных разностей и строится разностная задача, аппроксимирующая эту задачу с первым порядком точности.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «ПРИМЕНЕНИЕ МЕТОДА КОНЕЧНЫХ РАЗНОСТЕЙ К РЕШЕНИЮ ЗАДАЧИ ДЛЯ УРАВНЕНИЯ КОЛЕБАНИЯ ПРОДОЛЬНО НАПРЯЖЕННОГО СТЕРЖНЯ С ИНТЕГРАЛЬНЫМИ УСЛОВИЯМИ»

УДК 519.633

ПРИМЕНЕНИЕ МЕТОДА КОНЕЧНЫХ РАЗНОСТЕЙ К РЕШЕНИЮ ЗАДАЧИ ДЛЯ УРАВНЕНИЯ КОЛЕБАНИЯ ПРОДОЛЬНО НАПРЯЖЕННОГО СТЕРЖНЯ С

ИНТЕГРАЛЬНЫМИ УСЛОВИЯМИ

ХАНКИШИЕВ ЗАКИР ФАРМАН ОГЛЫ

Доцент кафедры Уравнений математической физики Бакинского Государственного

Университета, Баку, Азербайджан

Аннотация. В работе рассмотрена одна задача для уравнения колебания продольно напряженного стержня с интегральными условиями. Известно, что колебание продольно напряженного стержня выражается дифференциальным уравнением в частных производных четвертого порядка. Для таких уравнений ставятся различные задачи, в том числе задачи с интегральными условиями. При наличии интегрального условия в задаче, в ее решении возникают определенные трудности. В настоящей работе после замены интегральных условий нелокальными граничными условиями, к решению полученной новой задачи применяется метод конечных разностей и строится разностная задача, аппроксимирующая эту задачу с первым порядком точности.

Ключевые слова: Разностная задача, интегральные условия, аппроксимация, погрешность аппроксимации, нелокальные условия.

1. Постановка задачи

Задачи для дифференциальных уравнений с интегральными условиями часто встречаются в математической литературе. Это связано с тем, что такие задачи имеют большие прикладные значения [1]. В настоящей работе рассмотрена следующая задача: в замкнутой прямоугольной области О = {0 < х < /, 0 < t < Т} найти функцию и = и(х, t), удовлетворяющую уравнению

а^О^ аЧ^О^ а^О = о < х </,о ^ <т, (1)

аt2 ах4 ах2

граничным условиям

а2ч(о, t) л аV/, t) л

—= о, —V-2 = о, о < t < Т, (2)

ах2 ах2

интегральным условиям " и(

J С (x) u(x, t)dx = jux (t),

0 I

J c2 (x) u( x, t )dx = ju2 (t),

0

0 < t < T, (3)

и начальным условиям

и(х,о) = ^(х), аи(хо) = ^(х), о < х < /. (4)

аг

Здесь /(х, t), ¡их (t), ¡лг (^), (рх (х), (рг (х) - известные непрерывные функции, а > о, Ь > о -действительные числа. Предполагается, что функции с (х) и с2 (х) удовлетворяют системе дифференциальных уравнений

|с;(х) = а11с1 (х)+а12 с2 (хХ [с;(х) = а21с1 (х) + а22с2 (x),

где ак, /, к = 1,2 - известные действительные числа.

2. Замена интегральных условий (3) нелокальными граничными условиями

Рассмотрим первое интегральное условие в (3) и продифференцируем его по t два раза:

i

J Ci(x)

d2u(x, t)

dt2

dx = jul(t).

С учетом уравнения (1), это равенство можно переписать в виде

J Ci( x)

- a

д 4u(x, t) 2 д 2u(x, t)

dx4

- + b

dx2

+ f (x, t)

dx = jul(t),

или в виде

- a

J Ci(x) d U(t t) dx + b2 J Ci(x)

d2 u(x, t)

dx + J c (x)f (x, t )dx = ju"(t).

(6)

-х4 -х2

0 дх 0 дх 0 Применяя формулу интегрирования по частям, первым двум интегралам в этом равенстве, с учетом равенства (5) и граничных условий (2), после элементарных преобразований, получим:

i

J Ci( x)

d4u(x, t) , . 4d3u(x, t) dx = c (x)

0

dx4

dx3

+ (aiici(x) + ai2c2( x))

du(x, t)

dx

(a„c;(x) + ai2c2(x)Hx,t)lx=0 +(ai2i + ai2a2i )/"i(t) + (aiiai2 + ai2a22 (t).

-(a,c

l ^2

J Ci( x)

d u(x, t)

dx2

dx = c (x)

du(x, t)

dx

- c[( x)u(x, t } x=0 + an^i (t) + ai2^2 (t).

x=0

Подставив найденные выражения для этих интегралов в левую часть равенства (6), после простых преобразований, получим:

d3u(l,t) d3u(0,t) , . du(l,t) J/V^du(0,t) Ci (l)-^ - ci (0)-+ di (l) —^^ - di (0) v 7 7 - d; (l)u(l, t) +

dx3 1 dx3 + d[ (0)u(0, t) = Д (t),

dx

dx

(7)

где

d ( x ) =

f

b

aii - 2 v a J

2 ^ - i if

- ci (x) + ai2C2 (x), Mi (t) =--2 ^r(t) - J Ci (x)f (X, t)dx

a i

2

2

aii ^ ai2a2i aii ' 2

v a J

Mi(t) -

f

0

2

aiiai2 ^ ai2a22 ai2 ' 2 v a J

M2(t).

Используя второе интегральное условие, аналогичным образом, справедливость равенства

/;ч-3и(/,-) ,Л.д3и(0,Г) , ._-и(/,Г) , ^-и(0,Г) с2 (/) -иг^ - ^ (0) д') + а 2 (/) —^ - а2 (0) —^ _ а2 (/)и(/, г) +

получим

dx3 2 4 ' dx3 + d2 (0)u(0, t) = M (t),

dx

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

dx

(8)

где

Г

d2 (x) = a2iCi(x) +

f

a22 ' 2 a2

b2 ^ - i Г f '

— C2( x), M2(t) =--2 M2O) - I C2( x)f (X, t )dx

™ a J

,2 Л

aiia2i ^ a2ia22 a2i ' 2 v a J

f

Mi(t) -

2

0

2

a2iai2 ^ a22 a22 ' 2 v 2i i2 22 22 a2 J

M2(t).

Таким образом, интегральные условия (3), при выполнении условий (5), заменили нелокальными граничными условиями (7) и (8).

- 3и (/, -) - 3и(0, -)

Если исключить из условий (7)-(8), сначала ——;—, затем ——-—, то вместо этих

dx3

dx3

условии получим условия

0

0

x=0

x=0

0

д3u(0,t) ди(0,t) SU(l,t) . .

Л + а11 \ + «12 + ^13U(0, t) + «!4U(l, t) = ft (t), (9)

д. д. д.

^ + A ^ + A + AU0 t) + £„u(l, t) = д,(,), (10)

где

a,, =

д. д. д.

с2 (l)di (0) - Ci (l)d2 (0) _ Ci (l)d2 (l) - C2 (l)di (l)

Ci(0)c2(l) - Ci(l)C2(0) '

Ci(l )d 2 (0) - C2(l )d1(0)

Ci(0)c2(l) - Ci(l)C2(0)'

C2(0)di(0) - Ci(0)d 2(0)

Ci(0)c2(l) - Ci(l)C2(0) '

ci(0)d; (0) - -C2(0)d1(0)

Ci(0)c2(l) - -Ci(l)C2(0) '

_ C1(l)^2(t) - C2(l)Ai(t)

11 ------- ------ i2 Ci(0)c2(l) - cl(l)c2(0Y

= ^ ^ ...........= C2(l)d1(l)-Ci(l)d2(l)

13 14 Ci(0)c2(l) - Ci(l)c2(0)'

= C2(Q)di(Q) clVy;u2^ ff = C1(0)d2 (l) - C2 (0)d1 (l)

Ai1 ------„ ------ Ai2 Ci(0)c2(l)-Ci(l)c2(0)'

_C2(0)d1(l) - Ci(0)d; (l)

Pi3 = /74 /ГЧЧ , A14 =

Ci(0)c2(l) - Ci(l )C2(0)'

> J~2 (t) =

q(0)^(t) - C2(0)^i(t)

Ci (0)C2 (l) - Ci (l)C2 (0) 2 Ci (0)C2 (l) - Ci (l)C2 (0)

3. Построение разностной задачи

В дальнейшем рассмотрим задачу нахождения решения уравнения (1), удовлетворяющего граничным условиям (2), (9)-(10) и начальным условиям (4). К решению этой задачи будем применять метод конечных разностей. С этой целью определим в замкнутой области Б ={0 < х < 1,0 < t < Т} сеточную область

= {(хп , tJ X хп = пк> П = 0,1,-) Я, tJ = - = 0,1,-) jo, ЯИ = ^ ЛТ = Т }.

Рассмотрим сначала уравнение (1) в узлах (хи, tj), сеточной области аИт при п = 2,3,..., N — 2, j = 1,2,...,- — 1, и заменим в нем, значения частных производных соответствующими разностными выражениями:

и(хп , tj+1 ) — 2и(хп , tj ) + и(хп , t j—1 )

г2

■ +

2 и(Хп—2 , tj ) — 4и(Хп—1 , tj ) + 6и(Хп , tj ) — 4и(Хп+1 , tj ) + и(Хп+2 , tj )

+ а ---------—

И2

— Ъ 2 и( ХП—1, ^ ) — 2и( ^ , ^ ) + ^ ( ХП+1, ^ ) = / ( ^ , ^ ) + ОИ 2 Л

п = 2,3,...,N — 2, - = 1,2,...,-о -1. (11)

Очевидно, что эти равенства справедливы тогда, когда решение уравнения (1) в области Б = {0 < х < I, 0 < t < Т} имеет ограниченные частные производные по х до шестого и по t до четвертого порядков.

Если отбросим в равенствах (11) слагаемые порядка 0(И2 + т2), и при этом обозначим приближенное значение и (хп, tj) через у-, то относительно функций у- получим разностные уравнения

уГ1— 2уп + у-1, 2у-—2— 4у-—1 + 6у п — 4уп+1 + уп+2 ,2уп—1— 2у п + у^ _

-п--+ а ----ъ ---=

т2 И4 И2

= / -, п = 2,3,...,N — 2, - = 1,2,...,-0 — 1, (12)

где / п = / (х п, tJ).

Impact Factor: SJIF 2020 - 5.497 2021 - 5.81

В этих разностных уравнениях при каждом фиксированном значении j, участвует N +1 значение yj, n = 0,i,..., N, а число уравнений равно N - 3. Недостающих четырех уравнений мы получим с использованием условий (2) и (9)-(10).

Пусть x = ah, i < а < 2.Рассмотрим уравнение (1) в точке (x,t.):

d2u(x, tj) 2 d4u(x, tj) 2 d2u(x, tj)

1 " — b

dt2

■ + a

dx4

dx2

= f (x, t,).

(13)

Разложим функцию u(x, t ) по формуле Тейлора в окрестности точки (x, t ) :

u(x, tj) = u(x,t.) + (x - x)

du(x, tj) (x - x)2 d u(x, tj) (x - x)3 d3u(x, t,)

dx

- +

2!

dx2

+ -

3!

dx3

- +

+

(x - x)4 d u(x,tj) (x - x)5 d u(x,tj) (x - x)6 d6u(x,t})

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

4!

dx4

+

5!

dx5

+ -

6!

dx6

(14)

где х — точка, лежащая между точками х и х.

- 2и( х, -.)

Аналогичным образом, для функции

dx2

имеем:

d2u(x,tj) d2u(x,tj) d3u(x,tj) (x - x)2 d4u(x,tj) (x - x)3 d5u(x,tj)

dx2

dx2

- + (x - x)

dx3

+ -

2!

dx4

3!

dx5

■ +

(x - x)4 d6u(~, tj)

+

(15)

4! -х6 ' где ~ — точка, лежащая между точками х и х.

Взяв в равенстве (14), по очереди, х = х0, хх, х2, х3 а в равенстве (15) х = х0, получим пять равенств. Пусть а, Ь, с, а, е — произвольные числа. Умножим, обе части полученных первых четырех равенств, последовательно, на числа а, Ь, с, а, а последнего равенства на

е^2. После сложения этих равенств, получится равенство, которое, легко можно проверить,

24 , 60 48 , 12 12 12

что при а =--, Ь = —, с =--, а = —, е = —, а = —, принимает вид

11 11 11 11 11 11

Я4 Гi2h

d u

ii

i2

(- 2u( x0, t j) + 5u( xj, t j) - 4u( x2, t j) + u( x3, tj))+ o(h2).

dx4 iih4 v 'j> ^ 3> j

Аналогичным образом можно получить справедливость равенства

(16)

d 4ul l -

i2h

dx4

(- 2u(xn , tJ ) + 5u(xN-i , tJ ) - 4u(xN-2 , tj ) + u(xN-3 , tj )) +

iih J J J J

+ o(h2).

(17)

Равенство (13), с учетом равенства (16) можно переписать в виде

л2 Г i2h ,

d ul-, t,.

V ii 1

dt2 i2h

+ (- 2u(x0, tj) + 5u(xi, tj) - 4u( x2, tj) + u(x3, tj))- b'

db^

dx2

= f 1 1ТГ''1 )+o(h!)

так как x = — h.

ii

(18)

Легко можно показать справедливость равенства

ОФ "Международный научно-исследовательский центр "Endless Light in Science"

1 v ^ tj) + v [^, t, ) = v Г t,) + o(h 2), (19)

если функция v(x, t) имеет ограниченные частные производные по x до второго порядка. При

. . d 2u( x, t) v (x, t) = ——— это равенство принимает вид

дUl2k t I

1 д2u(0,tj) | д U 1,j J_ 12 д2u(xt,tj) 1 1 dt2 dt2 1 1 dt2

o(h2 ). (20)

„ , N д u(x, t)

Если v ( x, t ) =--—, то из равенства (19) получим справедливость равенства

дх

12h

1д 2u(0, tj )+^ ¿u^^u*

д 2u| ^^, t,

0(h2 )

11 дх2 дх2 11 дх2

или

52мf12h , |

Ч 11 , j )_ 12 дtj)+Ыи2

0(h2 ) (21)

дх2 11 дх2

д 2u (0, tt ) дх2

Учитывая равенства (20) и (21) в равенстве (18), приходим к равенству

так как в силу граничного условия в (2) -^ ' 1 = 0.

d2u(0,t,) d u(xl,t,) 12a2 / \ --—r~ +12-~2 + "^T" !" 2u(xo, t, ) + 5u( Xj, tj ) - 4u( X2, t, ) + u( x3, t, )j-

-12*'^ = 1 / f, tj V Oh 2 )■ (22)

Пусть ~ = ph, 0<p< 1.

d2u(x, t )

Разложим функцию u(x, t ) и --- по формуле Тейлора в окрестности точки

j дх

(~, t-):

, , ~ du(~, t, ) (x - ~)2 d2u(~, t, ) (x - ~)3 d3u(~, t, )

u(x, t.) = u(x, t .) + (x - x)-— + -----;Г— + -----T-O— +

^ ^ dx 2! dx2 3! dx3

(х -~)4 д u(x,tj) (х -~)5 д u(х,tj) (х -~)6 дХ^,tj) + 4! дх4 + 5! дх5 + 6! дх6 '

— и(х,г ) -2и(~,г ) ~д3и(~,г ) (х — ~)2 —4и(~,г ) (х — ~)3 —5и(~,г )

-1— =-;Г— + (х — х )--О— + ------О— + -----+

—х2 —х2 —х3 2! —х4 3! —х5

, (х — ~)4 —6и(х2, ^ )

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

4! —х6 ' где х1, .х2 — точки, лежащие между точками х и ~.

—3и(х, г )

Одновременно, функцию --- также разложим по формуле Тейлора в окрестности

—х

точки (~, г.):

д 3u( x, t,) д 3u(~, t,)

д4u(x, tj) (x - X)2 д5u(x, tj) (x - X)3 д6u(x.3, tj)

-=-+ (х — х)-+ --Ч- + --^^, (25)

сх3 сх3 дх4 2! дх5 3! дх6

где хз, — точка, лежащая между точками х и ~.

Взяв в равенстве (23), х = х0, хх, х2, а в равенствах (24) и (25) х = х0, получим снова пять равенств. Пусть а, Ъ, с, ё, е — произвольные числа. Умножим, обе части полученных первых трех равенств, последовательно, на числа а, Ъ, с, четвертого равенства на ёИ2 а последнего

равенства на еИ3 . После сложения этих равенств, получится равенство, которое, легко можно

12 , 24 12 , 12 12 „ 3

проверить, что при а = —, Ъ =--, с = —, ё =--, е =--, В = —,

7 7 7 17 7 7

принимает вид

д 4u

3h

12

дх4

7h4

u(х0, t .) - 2u(Xj, t,) + u(x2, t.) - h

д 3u( x0, tj)

ax^

+

o(h2).

(26)

Аналогичным образом можно получить справедливость равенства

д 4u

, 3h

l--, tj

7 1

12

ax4

7h4

u(xw, tj) - 2u(xN-i, tj) + u(xw-2, tj) - h

3h

ди( xn , tj)

Л

ax3

+

o(h2). (27)

Если рассмотреть уравнение (1) в точке (~, tj) = 1 , ^ I, то это уравнение, с учетом

равенства (26) можем переписать в виде

12а2

f 3h

д 2ul —, t

V 7 -

дt2 ' 7h4

V

= f i f, tj )+ Oh)

u(x0, t .) - 2u(xx, t ) + u(x2, t.) - h

д 3u( x0, t j)

дx3

Л д 2u| - b2

3h

дx2

На этот раз нам понадобится равенство

4 - (^ t])+3 ^ и tj)=V (3И, ^ )+о(и 2),

которое справедливо, если функция V (х, t) имеет ограниченные частные производные по х

(28) (29)

до второго порядка. При V (х, t) =

4 д2u(0, tj) 3 д 2u(x, tj) 7 dt2 + 7 д2 д 2u( x, t)

д2u( x, t)

^ f3h

д 2ul —, t

a2

это равенство принимает вид

■ +

o(h2).

(30)

Если v (x, t) =

дx2

то из равенства (29) получим справедливость равенства

4 д 2u(0, t ) 3 д 2u(x, t )

7 дx2 7 дx

13h

д 2ul —, t

дx2

■ +

o(h2)

или

2 f 3h ^ IT-^ J_ 3 д2u(xi,tj)

o(h2)

(31)

дx2 7 дx2

ОФ "Международный научно-исследовательский центр "Endless Light in Science"

t

7

3

3

t

7

3

7

7

так как в силу граничного условия в (2)

д 2u (0, tj) дx2

Учитывая равенства (30) и (31) в равенстве (28), после несложных вычислений получим справедливость равенства

4

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

д2 u(0, tj) д 2u(xi, tj) 12а ■ + 3-—— + ■

дt2

дt2

h4

u(x0, t,) - 2u(xx, t.) + u(x2, t.) - h

д 3u( x0, tj)

Л

дx3

- 3b:

д 2u (xx, tj)

aX2

= 7fГf, tj ) + 0(h2).

(32)

Если исключить из равенств (22) и (32), сначала

этих равенств получим равенства

д2 u( x, tj)

дt2

затем

д 2u(0, tj)

дt2

, то вместо

д 2u( x0, tj) 12а:

д'2 11

17 h4

6u(x0, t.) - 13u(x:, t.) + 8u(x2, t.) - u(x3, t.) - 4h

3 д 3u(tj)

дx3

17'

12h

11 , i2h 1 28 ,

= -- f —, t, 1 + 28 f

11

17

3h 1

T, ^ ) +

0(h2)

(33)

д 2u( xi, tj) 4а2

дt2

17h4

- 7u(x0, t.) + 18u(x:, t.) - 15u(x2, t ) + 4u(x3, t ) - h

3 д

дx3

д2u(xi, tj) 44 / 12h

- ьj

дx2

=41 f f , 'j)+51 f f 3h,'')+Oh2)

(34)

Аналогичным образом можно получить справедливость следующих равенств: 5 2и(хЯ—1, tj) 4а 2 /

--2—— + —4 (— 7и(хя, tj) +18и(хя—1, tj) — 15и(хя—2, tj) + 4и(хя—3, tj)

дt 17И J J J J

- h

д 3u( xn , tj)

)

д2u(xn , t,) 12а

дx3

- b2

7

' 2 ^u^ = 44ffl - 12h,j 1+Lffl - ,j )+Oh2) (35)

дx 51 Г 11 1) 51 Г 7 1) V Л

дt2 11

17 h4 12h

6u(xn ,tj) - 13u(xn-i, tj) + 8u(xn-2, tj) - u(xn-3, tj) - 4h

3 д u(xN , tj )

Л

ax3

=-17 f f l - Tt , 'i)+28- f (l - f, 'I )+o(h •).

(36)

В этих четрех равенствах участвуют значения частной производной функции и( х, t) по х третьего порядка в точках х = 0 и х = I. Используя равенства (9)-(10), их можно заменить следующими разностными выражениями:

д 3u( x0, tj) u( xi, tj) - u( x0, tj) u( xn , tj) - u( xn-i, tj)

дx

3 = -aii

h

-a

12

h

ai3u(x0 ,

tf) -

ai4u(xN , ' j

( xn , tj) + ~i(tj) + o(h2),

д3u(xw, t)

дx3

= -Ai

u(X1, ',, ) - u(x 0 , tj ) „ u(xN , 'j ) (xN-1, 'j )

h

h

-A

12

^^ X0, 'j ) -

-Ai4"(xn , 'j ) + /X2('; ) + O(h 2 ).

3

d2u(x;, t ) d 2u(xn-1, t )

и

, участвующих в

С другой стороны значение частных производных и

dx dx

этих равенствах, можно заменить обычными разностными выражениями с точностью o(h2). В результате таких замен, равенства (33)-(36) примут следующий вид:

d2u(x0,tf) i2a2 / \ -0— +--1 (6u( x0, tj) - 13u( x, tj) + 8u( x2, tj) - u( x3, tj))+

+ -

dt2 17hj 48a2 [ u(x;,tj) -u(x0,tj) u(Xn,tj) -u(x^i,tj)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

17 h

a

h

+ a

h

+ auu( x0, tj) + a;4u( xn , tj)

11 / 12h 1 28 / 3h 1 48a2 ~ , .

=--fl-, t, l + — f l —, t, 1 +-цМ,) + o(h),

17 l 11 j J 17 | 7 j) 17 h 1 j V '

d2u(x;,tj) 4a2 / ч -+ Г7/Й (- 7u( x0, tj) + 18u(x;, tj) - 15u( x2, tj) + 4u(x3, tj))+

+ ■

4a2 [ u(x;,tj) -u(x0,tj) u(Xn,tj) -U(XN-;,tj)

17 h

a

h

+ a

h

+ a;3u( x0, tj) + a;4u( xn , tj)

J

- b 2 »(x0, 'j)- 2u<x;,,,) + u( x2,,,) = 44 112Г, , ) + 2. f [ 3Г, ,)+ 4ai , + ^

h2 51 | 11 jJ 51 | 7 jJ 17h j w

d2u(xN-1, tj) 4a2 / \

-172-+ ^TT (- 7u(x N, t j) + i8u(xN-1, t j) - i5u( XN-2, t j) + 4u(xN- 3, t j))+

dt 17h

4a2 [_ u(x;,tj) -u(x0,tj) u(Xn,tj) -U(XN-;,tj)

+ ■

17 h

- b

Pi;

h

+ P;

h

■ + P13u(X0, tj ) + P14u(xN , tj )

2 u(xN, tj) - 2u(xN-i, tj) + u(xn-2, tj) = 44 f[r_ 12Г ^

= ' ii ,tj У 5У Г 7 j

h2

= 44 f 11 - ;2T, tj 1+7 f[ l - ?, tj 1+

+ — M2(t;) + o(h),

17 h j W d2 u( xn , tj) 12a2

+ ;27T (6u(xn , tj) - i3u(xN-i > tj) + 8u(xn-2, tj) - u(xN-3 > tj )) +

17h4 j j J J

dt2 17hv N-1' j

48a2 [rt u(x;,tj)-u(x0,tj) u(Xn,tj)-U(XN-;,tj)

+ ■

17 h

Pi,

h

+ P;

h

+ P13u(^ tj ) + P14u(xN , tj )

11 / 12Л 1 28 / 3й 1 48а2 ~ , .

= — л—, г, 1 + —Л1—, г, 1 +-/~2(г;) + 0Ш.

17 ^ 11 1) 17 ^ 7 1) 17 И 1 У ' Заменив в этих равенствах частные производные по г, соответствующими разностными выражениями, отбросив слагаемые порядка 0(И), обозначив, как в начале, приближенное значение и (хп, tj), через у-, получим еще четыре разностных уравнения, которых если

присоединим к разностным уравнениям (12), то придем к следующей разностной задаче, аппроксимирующей исходную задачу (1)-(4):

У0 ~ 2yj + У0- + ^ (6У0 - 13yj + 8yi - yj )+

+ ■

48a

17 h

2 a y;j , y° + a yN yN-; + a,3y0 + a,4yN

v

h

h

11 / 12h 1 28 13h 1 48a2 ~ . .

= -iif 1 ;11h• 'j )+18f(3T• 'J )+Mi(,j',

^-j11 + t40t (- 7yj +18yi -15 У2 + 4y3)+

+ ■

г2 17h ^

4а2f W _ W W _ W ^

aii ZL-y°L + ai2 yN /N-i + ai3yj + ai4yj

V h h )

17 h

-b2 y° 2yi ' y2 = 44fl —,t, 1 + -fl-,t, 1 + —£(',), h2 51 V 11 1) 51 V 7 1) 17h 1

JiMlZNL^ yJ-2-^-ynf'-i + 6УП-4yj+i + y;f+2 ,2 ynf-i-2УП + yj+i _

-^--+ а ---b ---=

r2 h4 h2

= fj, n = 2,3,...,N - 2, j = 1,2,..., j -1, (37)

yN+-1 - 2yj-1 +1 + ^T (- 7yN +18yj-i -15yNj-2 + 4yN-3)+

+

2 117 4 V s N s N-1 У N-2 У N-3 ,

г 17h

4а 2 f~ yi-yj^R yN - y^-i , A J

17 h

Ai —;-+ P12-;-+ A3 Уь + Ai4 yj

V h h )

,2 yjJ -2yN-1 + yjJ-2 44 12h 1 7 3A 1 4а2 x , .

-b ---= — f l l--,t, I + — fl l--,t, I + +-fi2(t,),

h2 51 г 11 1) 51 г 7 1) 17h 1

,J+1 -1 i"»-2

yj - 2yN + yj- + 12а (6yj 13yj + 8yj yj )+

-2-+ ^Г7Т (6-yN - 13yN-1 + 8yN-2 - yN-3 ) +

г 1/h

48а

+

2 Ay—^ + A12 + A13 У0 + A yN ^ h h

)

17 И

11 / 12И ^ 28 /3И ^ 48а2 ~ , .

= — п / [ "57,'-)+17/ [ Т,]+Т7И ^).

у,0 =«(х,), уп =«(х„) + тр2(х,), п = 0,1,2,...,N. (38)

Отметим, что первые два и последние два уравнения в (37) имеют погрешность О(И), а погрешность аппроксимации остальных уравнений равна о(и 2).

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Разностная задача (37)-(38) есть явная разностная схема. На основе этой схемы можно написать однопараметрическую неявную схему, которую можно решить известными методами.

ЛИТЕРАТУРА

1. Ханкишиев З.Ф. Решение одной задачи для линейного нагруженного диффе-ренциального уравнения параболического типа с интегральными условиями // Вестник Бакинского Университета, серия физико-математических наук.- 2021.- № 4.- с. 25-38.

i Надоели баннеры? Вы всегда можете отключить рекламу.