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

UptoLike

ВУЗ: 

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

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

Рубрика: 

153
Таблица невязки
U4
0
0.182
0.301
0.373
0.411
0.423
0.411
0.373
0.301
0.182
0
0.023
0.042
0.053
0.057
0.059
0.059
0.059
0.057
0.053
0.042
0.023
0.01
0.063
0.107
0.131
0.143
0.146
0.143
0.131
0.107
0.063
0.01
1.323 10
3
0.115
0.171
0.195
0.203
0.204
0.203
0.195
0.171
0.115
1.323 10
3
0.011
0.108
0.151
0.155
0.147
0.142
0.147
0.155
0.151
0.108
0.011
0.054
0.052
0.067
0.042
0.012
1.039 10
13
0.012
0.042
0.067
0.052
0.054
0.122
0.028
0.04
0.095
0.147
0.167
0.147
0.095
0.04
0.028
0.122
0.193
0.094
0.121
0.195
0.262
0.288
0.262
0.195
0.121
0.094
0.193
0.232
0.093
0.107
0.182
0.253
0.281
0.253
0.182
0.107
0.093
0.232
0.189
0.042
0.083
0.037
0.021
0.045
0.021
0.037
0.083
0.042
0.189
0
0.395
0.543
0.565
0.544
0.531
0.544
0.565
0.543
0.395
0
Максимальное значение |U4
ij
| равно
32 max max U4() min U4()()
32 0.56477
3. Введите пробные функции
V1 k m x y()sin
k x
a
sin
m y
b

Нормируем их. Для этого вычислим нормировочные коэффициенты
i 1 n1
j 1 n1
VV
i 1 j 1
0
a
y
0
b
xV1 i j x y()()
2
d
d

Получили нормированные пробные функции
Vkm x y()if k m 0
V1 k m
x
y
()
VV
k 1 m 1
d

Введите поверочные функции (для примера в качестве поверочных возьмем
пробные функции)
Wkm x
y
()Vkm
x
y
()
Найдем коэффициенты системы уравнений AC=B для определения
коэффициентов пробных решений C
k
i 1 n1
j 1 n1
B
i 1 n1 j 1()
0
a
y
0
b
xfxy()L 00 x y V()()Wi j x y()
d
d
i1 1 n1
j1 1 n1
i2 1 n1
j2 1 n1

Страницы