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

UptoLike

ВУЗ: 

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

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

Рубрика: 

141
)(665031.0
)
(
75110
9
.1
)
(
19784
2
.2
)
(
50384
3
.1
)
(
45401
3
.2
)
(
)
1,
(
xU5
x
U
4
x
U
3
x
U
2
x
U1
x
U
0
x
U
−
−+−++=
Выпишите аналогичным образом решение задачи своего варианта при 5=n .
Пробное решение U(x) для 5
=n при t=0; 0.1T; 0.2T;..., T имеет вид
Uxs,()V 0 x,()
1
n
k
Vkx,()Y
10 s⋅ k,
⋅
∑
=
+:=
Выпишите матрицу U2 получившегося пробного решения, разбив отрезок [a,
b] на 10 частей при t=0; 0.1T; 0.2T;..., T.
i 010..:=
j 010..:=
U2
ij,
Ua
ba−
10
i⋅+ j,
⎛
⎜
⎝
⎞
⎟
⎠
:=
Матрица пробного решения при t=0; 0.1T; 0.2T;..., T
U2
1.000000
0.211736
0.379137−
0.772617−
0.968705−
0.967401−
0.768705−
0.372617−
0.220863
1.011736
2.000000
1.000000
0.232938
0.358927−
0.753854−
0.948639−
0.946292−
0.748639−
0.353854−
0.241073
1.032938
2.000000
1.000000
0.251526
0.338787−
0.734303−
0.928632−
0.925849−
0.728632−
0.334303−
0.261213
1.051526
2.000000
1.000000
0.268133
0.318883−
0.714301−
0.908646−
0.905799−
0.708646−
0.314301−
0.281117
1.068133
2.000000
1.000000
0.283214
0.299317−
0.694078−
0.888664−
0.885965−
0.688664−
0.294078−
0.300683
1.083214
2.000000
1.000000
0.297098
0.280153−
0.673791−
0.868678−
0.866232−
0.668678−
0.273791−
0.319847
1.097098
2.000000
1.000000
0.310024
0.261420−
0.653549−
0.848689−
0.846531−
0.648689−
0.253549−
0.338580
1.110024
2.000000
1.000000
0.322170
0.243131−
0.633423−
0.828703−
0.826821−
0.628703−
0.233423−
0.356869
1.122170
2.000000
1.000000
0.333666
0.225287−
0.613460−
0.808729−
0.807082−
0.608729−
0.213460−
0.374713
1.133666
2.000000
1.000000
0.344613
0.207878−
0.593691−
0.788776−
0.787307−
0.588776−
0.193691−
0.392122
1.144613
2.000000
1.000000
0.355086
0.190891−
0.574135−
0.768857−
0.767498−
0.568857−
0.174135−
0.409109
1.155086
2.000000
⎛
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎝
⎞
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎠
=
Получ им матрицу предыдущего (для 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−
CP⋅:=
A1P
0.283699−
0.454931
0.293657−
1.082121− 10
14−
×
0.758218
2.674879−
3.610220
1.847938−
0.704776
2.825390−
4.741831
2.862814−
0.987765
4.839839−
8.861091
5.349758−
⎛
⎜
⎜
⎜
⎜
⎝
⎞
⎟
⎟
⎟
⎟
⎠
=
HP AP
1−
D1P⋅:=
Dt HP,()A1P HP⋅:=
YP rkfixed HP 0, T, 100, D,():=

Страницы