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

UptoLike

ВУЗ: 

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

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

Рубрика: 

151
Скопируйте в файл отчета этот вектор. Подставив коэффициенты C
k
,
наберите в файле отчета получившееся пробное решение.
Следовательно, пробное решение U(x,y) для
n1 3
имеет вид
Uxy()V 00 x y()
1
n1
k
1
n1
m
C
k 1 n1 m 1()
Vkm x y()

Построим таблицу U2 получившегося пробного решения, разбив область D на
100 частей, и график этого решения
i 01
0

j 01
0

U2
ij
Ua
i
10
b
j
10

Сравним точное и приближенное (при
n1 3
) решения, для этого найдем
разность матриц этих решений U1 и U2
Таблица сравнения точного и приближенного решения
U12 U1 U2
0
0
0
0
0
0
0
0
0
0
0
0
0.003
0.006
0.009
0.01
0.011
0.01
0.009
0.006
0.003
0
0
0.011
0.019
0.025
0.028
0.029
0.028
0.025
0.019
0.011
0
0
0.015
0.026
0.032
0.035
0.036
0.035
0.032
0.026
0.015
0
0
0.013
0.021
0.024
0.025
0.025
0.025
0.024
0.021
0.013
0
0
0.004
0.006
0.004
0.001
0
0.001
0.004
0.006
0.004
0
0
0.005
0.011
0.018
0.025
0.027
0.025
0.018
0.011
0.005
0
0
0.012
0.022
0.032
0.041
0.045
0.041
0.032
0.022
0.012
0
0
0.011
0.019
0.029
0.038
0.041
0.038
0.029
0.019
0.011
0
0
0.003
0.006
0.011
0.017
0.019
0.017
0.011
0.006
0.003
0
0
0
0
0
0
0
0
0
0
0
0
Максимальное значение |U12
ij
| равно
12 max max U12()min U12()()
12 0.045
Найдем предыдущее пробное решение
i 1 n1 1
j 1 n1 1
B1
i 1 n1 1()j 1()
0
a
0
b
xfxy()L 00 x y V()()Wi j x y()
d
d
i1 1 n1 1
j1 1 n1 1
i2 1 n1 1
j2 1 n1 1

Страницы