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

UptoLike

ВУЗ: 

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

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

Рубрика: 

55
Ux() V 0 x()
1
n
k
C
k 1
Vkx()

Найдем вектор коэффициентов C
k
для предыдущего пробного решения.
Для этого решим систему уравнений A
1
C=B
1
, где A
1
– угловая матрица ( 1
n)-
го порядка матрицы A, а B
1
– вектор-столбец, содержащий первые ( 1
n)
элементы столбца B.
C1 if n 1 submatrix A 0 n 2 0 n 2()()
1
submatrix B 0 n 2 0 0() 10

UP x() if n 1 V 0 x()
1
n 1
k
C1
k 1
Vkx()
V 0 x()

Cравним полученные решения для
n 5
и 4
n
0 0.5 1
0
2
10
4
4
10
4
6
10
4
Ux()UP x()
x
Замените старое значение меры точности
1
наибольшим значением
Ux() UP x()
на отрезке [a,b]
13 4.801 10
4

Найдем невязку полученного решения
Rx() L1 0 x V()gx()
1
n
i
C
i 1
L1 i x V()

0 0.5 1
0
0.002
0.004
0.006
Rx()
x
Замените старое значение меры точности 23
наибольшим значением
Rx()
на отрезке [a,b]
23 4.606 10
3


Страницы