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

UptoLike

91
Получим матрицу предыдущего (для 4
n) пробного решения
AP submatrix A 0 n 2 0 n 2
()
CP submatrix C 0 n 2 0 n 2
()
D1P submatrix D1 0 n 2 0
0
()
A1P AP
1
C
P

HP AP
1
D1
P

DtHP()A1P H
P

YP rkfixed HP 0 T 100 D()
Следовательно, предыдущее пробное решение U(x) для
n 5
имеет вид
UP x() V 0 x()
1
n 1
k
Vkx()YP
100

Cравним полученные решения для
n 5
и 4
n
при t=T
02
0
0.01
0.02
Ux()UP x()
x
Замените старое значение меры точности
2
2
наибольшим значением
Ux() UP x()
на отрезке [a,b]
22 0.019
Найдем невязки полученного пробного решения.
При t=T получим невязку


)(),0()(),0(1)(),0(2)(),()(
),(1)(),(2)(1),(:)(1
,100
111
,1001,1
xgxVxxVxK
dx
d
xVxKYxkVx
xkVxK
dx
d
xkVxKYAxkVxR
k
n
k
n
k
n
z
zzk
02
0
0.005
0.01
0.015
R1 x()
x