ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 07.12.2021
Просмотров: 2281
Скачиваний: 1
91
Получим
матрицу
предыдущего
(
для
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
CP
HP
AP
1
D1P
D t HP
(
)
A1P HP
YP
rkfixed HP
0
T
100
D
(
)
Следовательно
,
предыдущее
пробное
решение
U(x)
для
n
5
имеет
вид
UP x
( )
V
0
x
(
)
1
n
1
k
V k x
(
)
YP
100 k
C
равним
полученные
решения
для
n
5
и
4
n
при
t=T
0
2
0
0.01
0.02
U x
( )
UP x
( )
x
Замените
старое
значение
меры
точности
22
наибольшим
значением
U x
( )
UP x
( )
на
отрезке
[a,b]
22
0.019
Найдем
невязки
полученного
пробного
решения
.
При
t=T
получим
невязку
)
(
)
,
0
(
)
(
)
,
0
(
1
)
(
)
,
0
(
2
)
(
)
,
(
)
(
)
,
(
1
)
(
)
,
(
2
)
(
1
)
,
(
:
)
(
1
,
100
1
1
1
,
100
1
,
1
x
g
x
V
x
x
V
x
K
dx
d
x
V
x
K
Y
x
k
V
x
x
k
V
x
K
dx
d
x
k
V
x
K
Y
A
x
k
V
x
R
k
n
k
n
k
n
z
z
z
k
0
2
0
0.005
0.01
0.015
R1 x
( )
x
92
Замените
старое
значение
меры
точности
32
наибольшим
значением
R1 x
( )
на
отрезке
[a,b]
32
0.012
При
t=
0
получим
невязку
R2 x
( )
V
0
x
(
)
f x
( )
1
n
k
D2
k
1
V k x
(
)
0
2
0
5
10
15
R2 x
( )
x
Замените
старое
значение
меры
точности
42
наибольшим
значением
R2 x
( )
на
отрезке
[a,b]
42
7.55 10
15
3.
Введите
систему
пробных
и
поверочных
функций
(
для
примера
в
качестве
пробных
и
поверочных
функций
возьмем
систему
тригонометрических
функций
):
V k x
(
)
if k
0
sin
2
k
1
(
)
b a
x
2
b0 a2
b
b2 a0
a
a0 b0
b a
(
)
b2 a0
b0 a2
(
)
x
a0 b0
b a
(
)
k
1
n
V1 k x
(
)
if k
0
2
k
1
(
)
b
a
cos
2
k
1
(
)
b
a
x
2
b2 a0
b0 a2
a0 b0
b
a
(
)
V2 k x
(
)
if k
0
2
k
1
(
)
2
b
a
2
sin
2
k
1
(
)
b
a
x
2
0
W k x
(
)
V k x
(
)
93
Найдем
коэффициенты
системы
дифференциальных
уравнений
A
t
H
d
d
C H
B
для
отыскания
функций
H
k
(t)
с
начальными
условиями
A H
0
( )
D1
i
1
n
B
i
1
a
b
x
K x
( )
V2
0
x
(
)
x
K x
( )
V1
0
x
(
)
d
d
x
( )
V
0
x
(
)
g x
( )
W i x
(
)
d
i
1 n
j
1 n
A
i
1
j
1
a
b
x
V j x
(
)
W i x
(
)
d
C
i
1
j
1
a
b
x
K x
( )
V2 j x
(
)
x
K x
( )
d
d
V1 j x
(
)
x
( )
V j x
(
)
W i x
(
)
d
i
1 n
D1
i
1
a
b
x
f x
( )
V
0
x
(
)
(
)
W i x
(
)
d
Приведем
систему
к
виду
t
H
d
d
A1 H
B1
с
начальными
условиями
H
0
( )
D2
A1
A
1
C
B1
A
1
B
D2
A
1
D1
Найдем
решение
получившейся
системы
дифференциальных
уравнений
:
H
D2
D t H
(
)
A1 H
B1
Y
rkfixed H
0
T
100
D
(
)
Следовательно
,
при
t=T
получим
следующие
коэффициенты
Y
100
k
2.8878
0.0481
2.0958
10
3
6.9289
10
5
1.3289
10
6
94
Подставив
коэффициенты
Y
100,k
,
наберите
в
файле
отчета
получившееся
пробное
решение
.
Для
примера
решение
имеет
вид
U
(
x
,1)=
U
0
(
x
)–2.8878
U
1
(
x
)–
–0.0481
U
2
(
x
)–0.0021
U
3
(
x
)–0.000069
U
4
(
x
)–0.000001
U
5
(
x
).
Пробное
решение
U(x)
для
n
5
при
t= T
имеет
вид
U x
( )
V
0
x
(
)
1
n
k
V k x
(
)
Y
100 k
График
пробного
решения
0
2
2
0
2
U x
( )
x
Сравним
решения
,
полученные
методом
Галеркина
и
с
помощью
метода
Фурье
при
t=T
0
2
0
5
10
9
1
10
8
1.5
10
8
UT x T
(
)
U x
( )
x
Замените
старое
значение
меры
точности
13
наибольшим
значением
UT x
( )
U x
( )
на
отрезке
[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
CP
HP
AP
1
D1P
D t HP
(
)
A1P HP
YP
rkfixed HP
0
T
100
D
(
)
Следовательно
,
предыдущее
пробное
решение
U(x)
для
n
5
имеет
вид
95
UP x
( )
V
0
x
(
)
1
n
1
k
V k x
(
)
YP
100 k
C
равним
полученные
решения
для
n
5
и
4
n
при
t=T
0
2
0
5
10
7
1
10
6
1.5
10
6
U x
( )
UP x
( )
x
Замените
старое
значение
меры
точности
23
наибольшим
значением
U x
( )
UP x
( )
на
отрезке
[a,b]
23
1.06 10
6
Найдем
невязки
полученного
пробного
решения
.
При
t=T
получим
невязку
)
(
)
,
0
(
)
(
)
,
0
(
1
)
(
)
,
0
(
2
)
(
)
,
(
)
(
)
,
(
1
)
(
)
,
(
2
)
(
1
)
,
(
:
)
(
1
,
100
1
1
1
,
100
1
,
1
x
g
x
V
x
x
V
x
K
dx
d
x
V
x
K
Y
x
k
V
x
x
k
V
x
K
dx
d
x
k
V
x
K
Y
A
x
k
V
x
R
k
n
k
n
k
n
z
z
z
k
0
1
2
3
0
5
10
17
R1 x
( )
x
Замените
старое
значение
меры
точности
33
наибольшим
значением
R1 x
( )
на
отрезке
[a,b]
33
5 10
17
При
t=
0
получим
невязку
R2 x
( )
V
0
x
(
)
f x
( )
1
n
k
D2
k
1
V k x
(
)