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

UptoLike

ВУЗ: 

УлГТУ | Ульяновск

Составители: 

Рубрика: 

123
порядка матрицы A, а B
1
вектор-столбец, содержащий первые (1−n ) элементы
столбца B.
C1 if n 1≠ submatrix A 0, n 2−, 0, n 2−,()()
1−
submatrix B 0, n 2−, 0, 0,()⋅, 10,
⎡
⎣
⎤
⎦
:=
C1
T
1.283640 3.626226− 0.295870 3.106352−()=
Получ им матрицу предыдущего (для 4
=n ) пробного решения, разбив
отрезок [a, b] на 10 частей
UP x() if n 1≠ V 0 x,()
1
n 1−
k
C1
k 1−
Vk x,()⋅
∑
=
+, V
0 x,(),
⎛
⎜
⎜
⎝
⎞
⎟
⎟
⎠
:=
i 010..:=
U3
i 0,
a
ba−
10
i⋅+:=
U3
i 1,
UP a
ba−
10
i⋅+
⎛
⎜
⎝
⎞
⎟
⎠
:=
Матрица предыдущего пробного решения
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
=
062772.0403509.0633662.0779985.0864404.0
0.19.08.07.06.0
904636.0914814.0906106.0887337.0865612.0846932.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)
⎟
⎟
⎠
⎞
−−−−−
⎜
⎜
⎝
⎛
−−−−−−
=
0005.00007.00008.00007.00005.0
0.19.08.07.06.0
0004.00005.00006.00007.00004.00002.0
5.04.03.02.01.00.0
T
U23
Максимальное значение |U23
ij
| равно
K11 max max U23
1〈〉
(
)
min U23
1〈〉
(
)
,
(
)
:=
K11 0.000806=
Выпишите это значение для n=1, …, n=5 и сделайте вывод.
На йдем невязку полученного решения
Rx() L 0 x, V,()fx()−
1
n
i
C
i 1−
Lix, V,()⋅
∑
=
+:=
i 010..:=

Страницы