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

UptoLike

ВУЗ: 

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

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

Рубрика: 

94
Подставив коэффициенты Y
100,k
, наберите в файле отчета получившееся
пробное решение.
Для примера решение имеет вид U(x,1)=U
0
(x)–2.8878U
1
(x)–
–0.0481U
2
(x)–0.0021U
3
(x)–0.000069U
4
(x)–0.000001U
5
(x).
Пробное решение U(x) для
n 5
при t= T имеет вид
Ux() V 0 x()
1
n
k
Vkx()Y
100
k

График пробного решения
02
2
0
2
Ux()
x
Сравним решения, полученные методом Галеркина и с помощью метода Фурье
при t=T
02
0
5
10
9
1
10
8
1.5
10
8
UT x T()Ux()
x
Замените старое значение меры точности
1
3
наибольшим значением
UT x() Ux()
на отрезке [a,b]
13 1.068 10
8

Получим матрицу предыдущего (для 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
C
P

HP AP
1
D1

DtHP()A1P H
P

YP rkfixed HP 0 T 100 D()
Следовательно, предыдущее пробное решение U(x) для
n 5
имеет вид

Страницы