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

UptoLike

ВУЗ: 

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

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

Рубрика: 

142
Следовательно, предыдущее пробное решение U(x) для 5=n
при t=0; 0.1T; ..., T
имеет вид
UP x s,()V 0 x,()
1
n 1−
k
Vkx,()YP
10s⋅ k,
⋅
∑
=
+:=
Выпишите матрицу U3 получившегося пробного решения, разбив отрезок [a,
b] на 10 частей при t=0, 0.1T,..., T, и график этого решения при t=T
i 010..:=
j 010..:=
U3
ij,
UP a
ba−
10
i⋅+ j,
⎛
⎜
⎝
⎞
⎟
⎠
:=
Матрица предыдущего пробного решения при t=0, 0.1T,..., T
U3
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.229566
0.356331−
0.751021−
0.949740−
0.949631−
0.749740−
0.351021−
0.243669
1.029566
2.000000
1.000000
0.246335
0.334449−
0.729735−
0.930511−
0.931371−
0.730511−
0.329735−
0.265551
1.046335
2.000000
1.000000
0.262140
0.313414−
0.708747−
0.911065−
0.912691−
0.711065−
0.308747−
0.286586
1.062140
2.000000
1.000000
0.277069
0.293155−
0.688046−
0.891445−
0.893656−
0.691445−
0.288046−
0.306845
1.077069
2.000000
1.000000
0.291201
0.273609−
0.667621−
0.871688−
0.874321−
0.671688−
0.267621−
0.326391
1.091201
2.000000
1.000000
0.304609
0.254718−
0.647463−
0.851829−
0.854737−
0.651829−
0.247463−
0.345282
1.104609
2.000000
1.000000
0.317356
0.236430−
0.627562−
0.831897−
0.834952−
0.631897−
0.227562−
0.363570
1.117356
2.000000
1.000000
0.329501
0.218700−
0.607910−
0.811921−
0.815005−
0.611921−
0.207910−
0.381300
1.129501
2.000000
1.000000
0.341097
0.201484−
0.588501−
0.791923−
0.794933−
0.591923−
0.188501−
0.398516
1.141097
2.000000
1.000000
0.352192
0.184744−
0.569326−
0.771927−
0.774770−
0.571927−
0.169326−
0.415256
1.152192
2.000000
⎛
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎝
⎞
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎠
=
Выпишите матрицу сравнения
полученных решений для 5=n и 4=n
Матрица сравнения пробных решений
U23 U2 U3−:=
0
3.025− 10
13−
×
2.54− 10
13−
×
3.104− 10
13−
×
4.396− 10
13−
×
4.032− 10
13−
×
2.043 10
14−
×
8.464 10
13−
×
1.736 10
12−
×
1.892 10
12−
×
0
0
0.003
0.003−
0.003−
0.001
0.003
0.001
0.003−
0.003−
0.003
0
0
0.005
0.004−
0.005−
0.002
0.006
0.002
0.005−
0.004−
0.005
0
0
0.006
0.005−
0.006−
0.002
0.007
0.002
0.006−
0.005−
0.006
0
0
0.006
0.006−
0.006−
0.003
0.008
0.003
0.006−
0.006−
0.006
0
0
0.006
0.007−
0.006−
0.003
0.008
0.003
0.006−
0.007−
0.006
0
0
0.005
0.007−
0.006−
0.003
0.008
0.003
0.006−
0.007−
0.005
0
0
0.005
0.007−
0.006−
0.003
0.008
0.003
0.006−
0.007−
0.005
0
0
0.004
0.007−
0.006−
0.003
0.008
0.003
0.006−
0.007−
0.004
0
0
0.004
0.006−
0.005−
0.003
0.008
0.003
0.005−
0.006−
0.004
0
0
0.003
0.006−
0.005−
0.003
0.007
0.003
0.005−
0.006−
0.003
0
⎛
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎝
⎞
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎠
⎛
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎝
⎞
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎠
=
Максимальное значение |U23
ij
| при
T
t
= равно
K11 max max U23
10
〈〉
(
)
min U23
10
〈〉
(
)
,
(
)
:=

Страницы