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

UptoLike

117
AA
i 1 j 1
0

AA
ni 1 nj 1
M1
i 1 j 1

AA
ni 1 j 1
C1
i 1 j 1

AA
i 1 nj 1
if ij1 0()
i1
n

BB
ni 1
B1
i 1

BB
i 1
0

Найдем решение системы дифференциальных уравнений dH/dt=AA*H+BB.
HD
2

DtH()AA H
BB

Y rkfixed H 0 T 100 D()
Следовательно, при t=T получим следующие коэффициенты
k 1
n

Y
100 k
1.4891
2.0906
9.0874
11.0421
4.1935
Подставив коэффициенты Y
100,k
, наберите в файле отчета получившееся
пробное решение.
Для примера решение имеет вид U(x,1)=U
0
(x)+1.4891U
1
(x)–
–2.0906U
2
(x)+9.0874U
3
(x)–11.0421U
4
(x)+ 4.1935U
5
(x).
Пробное решение U(x) для
n 5
при t= T имеет вид
Ux() V 0 x()
1
n
k
Vkx()Y
100

График пробного решения
02
0
1
2
Ux()
x
Сравним решения, полученные методом Галеркина и с помощью метода Фурье
при t=T