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

UptoLike

ВУЗ: 

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

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

Рубрика: 

60
0 0.5 1
0
1
10
5
2
10
5
3
10
5
Yk x()Ux()
x
Замените старое значение меры точности
3
5
наибольшим значением
Yx() Ux()
на отрезке [a,b]
35 2.334 10
5

2. В качестве пробных функций возьмем систему функций (2.29).
Введите систему пробных функций:
12
32
)1(,
4
1,2,
3
1,1,56,0:),(
k
xx
x
xkif
x
xkifxkifxkV2
Найдем коэффициенты системы уравнений AC=B для определения
коэффициентов пробных решений C
k
i1
n

j1
n

B
i 1
a
b
x
fx() L 0 x V2()()Li x V2()
d
A
i 1 j 1
a
b
x
Li x V2()Ljx V2()
d
Решая систему уравнений AC=B матричным методом, получим вектор
коэффициентов C
k
CA
1
B

C
T
4.207857 9.361081 1.023365 0.337991 0.171763()
Скопируйте в файл отчета этот вектор. Подставив коэффициенты C
k
,
наберите в файле отчета получившееся пробное решение.
Пробное решение U(x) для
n 5
имеет вид
Ux() V2 0 x()
1
n
k
C
k 1
V2 k x()

Найдем вектор коэффициентов C
k
для предыдущего пробного решения.
Для этого решим систему уравнений A
1
C=B
1
, где A
1
– угловая матрица ( 1
n)-
го порядка матрицы A, а B
1
– вектор-столбец, содержащий первые ( 1
n)
элементы столбца B.

Страницы