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

UptoLike

66
;
7
1
1)(;
6
1
1)(;
5
1
1)(
;
4
1
1)(;
3
1
1)(;56)(
6
5
5
4
4
3
3
2
2
10
xxxuxxxuxxxu
xxxuxxxuxxu
причем

0
0
uR .
В результате расчета по программе ODE.mcd при 5
n получим вектор
коэффициентов
)130778,1133009,0466128,2564241,2140938,1(
С
.
Подставив коэффициенты
k
C , набираем в файле отчета получившееся пробное
решение:
).(130778,1)(133009,0
)(466128,2)(564241,2)(140938,1)()(
54
32105
xuxu
xuxuxuxuxy
Определяем значения мер точности:

,0004801,0)()(max
45
,
13
xyxy
ba

0,004606,,,...,max
51
,
23
xCCR
ba

,00003062.0)()(max
5
,
33
xyxY
ba
4.2. В качестве пробных функций используем функции, полученные в
примере 4 раздела 2.5:
);4cos(4)4sin()();3cos(3)3sin()();2cos(2
)2sin()();cos()sin()(;)(;56)(
54
3210
xxxuxxxux
xxuxxxuexuxxu
x
причем

0
0
uR .
Замечание. Процедуру получения функций вида (2.34) – (2.36)
необходимо описать в файле отчета.
В результате расчета по программе ODE.mcd при 5
n получим вектор
коэффициентов
)000202,0001314,0008526,0209285,0539001,4(
С .
Подставив коэффициенты
k
C
, набираем в файле отчета получившееся пробное
решение:
).(000202,0)(001314,0
)(008526,0)(209285,0)(539001,4)()(
54
32105
xuxu
xuxuxuxuxy
Определяем значения мер точности:

,006344,0)()(max
45
,
14
xyxy
ba

0,229,,,...,max
51
,
24
xCCR
ba

,006182.0)()(max
5
,
34
xyxY
ba
5. Найдем приближенное решение краевой задачи (2.41) интегральным
методом наименьших квадратов.