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

UptoLike

ВУЗ: 

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

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

Рубрика: 

153
Выпишите матрицу U2 получившегося пробного решения, разбив отрезок
[a,b] на 10 частей при t=0, 0.1T,..., T.
i 010..:=
j 010..:=
U2
ij,
Ua
ba−
10
i⋅+ j,
⎛
⎜
⎝
⎞
⎟
⎠
:=
Таблица пробного решения при t=0, 0.1T,..., T
U2
1
0.211736
0.379137−
0.772617−
0.968705−
0.967401−
0.768705−
0.372617−
0.220863
1.011736
2
1
0.222982
0.369046−
0.763454−
0.958653−
0.956663−
0.758653−
0.363454−
0.230954
1.022982
2
1
0.255069
0.33871−
0.735383−
0.92855−
0.924937−
0.72855−
0.335383−
0.26129
1.055069
2
1
0.303519
0.288078−
0.686925−
0.87852−
0.873467−
0.67852−
0.286925−
0.311922
1.103519
2
1
0.362292
0.217478−
0.616388−
0.80867−
0.80368−
0.60867−
0.216388−
0.382522
1.162292
2
1
0.425439
0.128049−
0.522747−
0.718993−
0.716472−
0.518993−
0.122747−
0.471951
1.225439
2
1
0.488581
0.022063−
0.406384−
0.609314−
0.611655−
0.409314−
6.384336− 10
3−
×
0.577937
1.288581
2
1
0.549782
0.097031
0.269432−
0.479332−
0.487772−
0.279332−
0.130568
0.697031
1.349782
2
1
0.609527
0.224919
0.115568−
0.328764−
0.34243−
0.128764−
0.284432
0.824919
1.409527
2
1
0.669893
0.357026
0.050668
0.157581−
0.173056−
0.042419
0.450668
0.957026
1.469893
2
1
0.733225
0.489275
0.22484
0.033727
0.022119
0.233727
0.62484
1.089275
1.533225
2
⎛
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎝
⎞
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎠
=
Получ им матрицу предыдущего (для 4
=n ) пробного решения
AP submatrix A 0, n 2−, 0, n 2−,():=
MP submatrix M 0, n 2−, 0, n 2−,():=
CP submatrix C 0, n 2−, 0, n 2−,():=
BP submatrix B 0, n 2−, 0, 0,():=
D1P submatrix D1 0, n 2−, 0, 0,():=
N1P submatrix N1 0, n 2−, 0, 0,():=
M1P AP
1−
MP⋅:=
C1P AP
1−
CP⋅:=
B1P AP
1−
BP⋅:=
D2P AP
1−
D1P⋅:=
N2P AP
1−
N1P⋅:=
i1n1−..:=
j1n1−..:=
AAP
i 1− j 1−,
0:=
AAP
ni+ 2− nj+ 2−,
M1P
i 1− j 1−,
:=
AAP
ni+ 2− j 1−,
C1P
i 1− j 1−,
:=
AAP
i 1− nj+ 2−,
if ij1, 0,():=

Страницы