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

UptoLike

ВУЗ: 

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

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

Рубрика: 

154
AAP
0
0
0
0
2.837−
4.549
2.937−
1.027 10
14−
×
0
0
0
0
7.582
26.749−
36.102
18.479−
0
0
0
0
7.048
28.254−
47.418
28.628−
0
0
0
0
9.878
48.398−
88.611
53.498−
1
0
0
0
0
0
0
0
0
1
0
0
0
0
0
0
0
0
1
0
0
0
0
0
0
0
0
1
0
0
0
0
⎛
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎝
⎞
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎠
=
i1n1−..:=
D2P
ni+ 2−
N2P
i 1−
:=
D2P
T
3.194 2.104 10
12−
× 1.156− 10
12−
× 1.926 10
13−
× 0000
()
=
i1n1−..:=
BBP
i 1−
0:=
BBP
ni+ 2−
B1P
i 1−
:=
BBP
T
00000000()=
HP D2P:=
Dt HP,()AAP HP⋅
B
BP+:=
YP rkfixed HP 0, T, 100, D,():=
Следовательно, предыдущее пробное решение 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.
i 010..:=
j 010..:=
U3
ij,
UP a
ba−
10
i⋅+ j,
⎛
⎜
⎝
⎞
⎟
⎠
:=
Матрица предыдущего пробного решения при t=0, 0.1T,..., T

Страницы