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

UptoLike

Рубрика: 

161
U1
min U1()10=
max U1()12.043=
Получение приближенного решения
Введите порядок пробного решения
=
+=
n
k
kn
yxkVCyxVU
1
),,(),,0(.
n 3:=
Введите пробные функции
V1 k x, y,()x
k
ax() y
k
by():=
Нормируем их. Для этого вычислим нормировочные коэффициенты и
выпишем их аналитические выражения
i 1 n..:=
VV
i 1
0
a
y
0
b
xV1 i x, y,()()
2
d
d:=
VV
T
1
30
π
5
1
105
π
7
1
252
π
9
Получили нормированные пробные функции
Vkx, y,()if k 0
V1 k x, y,()
VV
k 1
, d,
:=
Выпишем первые четыре функции
V 0 x, y,()10
V 1 x, y,()30 x
π
x
()
y
π
y
()
π
5
V 2 x, y,()105 x
2
π
x
()
y
2
π
y
()
π
7
V 3 x, y,()252 x
3
π
x
()
y
3
π
y
()
π
9
Введем оператор, равный левой части уравнения