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

UptoLike

90
Найдем решение получившейся системы дифференциальных уравнений:
HD
2

DtH()A1 H
B
1
Y rkfixed H 0 T 100 D()
Следовательно, при t=T получим следующие коэффициенты
Y
100 k
2.4678
1.4224
2.0088
1.5563
0.591
Подставив коэффициенты Y
100,k
, наберите в файле отчета получившееся
пробное решение.
Для примера решение имеет вид U(x,1)=U
0
(x)+2.4678U
1
(x)+
+1.4224U
2
(x)–2.0088U
3
(x)+1.5563U
4
(x)–0.591U
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
4
0.001
0.0015
UT x T()Ux()
x
Замените старое значение меры точности
1
наибольшим значением
UT x() Ux()
на отрезке [a,b]
12 1.072 10
3
