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

UptoLike

Рубрика: 

130
U3
i 0,
a
ba
10
i+:=
U3
i 1,
UP a
ba
10
i+
:=
Матрица предыдущего пробного решения
=
061840.0402690.0633110.0779613.0864049.0
0.19.08.07.06.0
904175.0914222.0905463.0886786.0865260.0846705.0
5.04.03.02.01.00.0
T
U3
Выпишите матрицу сравнения полученных решений для 5=n и 4=n
U23
0〈〉
U2
0〈〉
:=
U23
1〈〉
U2
1〈〉
U3
1〈〉
:=
Матрица сравнения приближенных решений U
n
(x) и U
n-1
(x)
=
00048.000016.000021.000026.000008.0
0.19.08.07.06.0
00012.000015.000004.000008.000005.000008.0
5.04.03.02.01.00.0
T
U23
Максимальное значение |U23
ij
| равно
K31 max max U23
1〈〉
(
)
min U23
1〈〉
(
)
,
(
)
:=
K31 0.000479=
Выпишите это значение для n=1,, n=5 и сделайте вывод.
На йдем невязку полученного решения
Rx() L 0 x, V,()fx()
1
n
i
C
i 1
Lix, V,()
=
+:=
i 010..:=
UN
i 0,
a
ba
10
i+:=
UN
i 1,
Ra
ba
10
i+
:=
Матрица невязок
=
0176.00025.00039.00002.00021.0
0.19.08.07.06.0
0013.00008.00016.00003.00021.00055.0
5.04.03.02.01.00.0
T
UN