Решение линейных задач математической физики на основе методов взвешенных невязок. Анкилов А.В - 129 стр.

UptoLike

129
AAP
i 1 nj 2
if ij1 0()
i1n1
D2P
ni 2
N2P
i 1

i1n1
BBP
i 1
0

BBP
ni 2
B1P
i 1

HP D2
P

DtHP()AAP HP
B
B
P

YP rkfixed HP 0 T 100 D()
Следовательно, предыдущее пробное решение U(x) для
n 5
при t=T имеет
вид
UP x() V 0 x()
1
n 1
k
Vkx()YP
100
k

Cравним полученные решения для
n 5
и 4
n
при t=T
02
0
0.002
0.004
Ux()UP x()
Замените старое значение меры точности
2
3
наибольшим значением
Ux() UP x()
на отрезке [a,b]
23 3.183 10
3

Найдем невязки полученного пробного решения.
При t=T получим невязку



)(),0()(),0()(2),0()(),()(
),()(),()(),()(
),(:)(
,100
111
,1001,1
11
1,1001,1,1001,1
xgxVxxV1xKxV2xK1YxkVx
xkV1xK2xkV2xK1YM1xkVx
B1YC1YM1xkVxR1
k
n
k
n
k
n
z
znzk
n
k
n
z
kzzkznzk



