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

UptoLike

ВУЗ: 

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

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

Рубрика: 

57
Найдем вектор коэффициентов 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 V1 0 x()
1
n 1
k
C1
k 1
V1 k x()
V1 0 x()

Cравним полученные решения для
n 5
и
n 1
()4
0 0.5 1
0
0.005
0.01
Ux()UP x()
x
Замените старое значение меры точности
1
наибольшим значением
Ux() UP x()
на отрезке [a,b]
14 6.344 10
3

Найдем невязку полученного решения
Rx() L1 0 x V1()gx()
1
n
i
C
i 1
L1 i x V1()

0 0.5 1
0
0.1
0.2
0.3
Rx()
x
Замените старое значение меры точности 24
наибольшим значением
Rx()
на отрезке [a,b]
24 0.229
Сравним решения, полученные методом Ритца и с помощью стандартной
функции системы MathCAD

Страницы