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

UptoLike

ВУЗ: 

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

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

Рубрика: 

128
Введем обозначения
Kx() e
a
x
tpt()
⌠
⎮
⌡
d
:=
β
x() Kx()− qx()⋅:=
gx() Kx() fx()⋅:=
Следовательно,
Kx() exp 3− x⋅()→
β
x() 2− exp 3− x⋅()⋅→
gx() exp 3− x⋅()2 x
2
⋅ 6 x⋅− 2+
()
⋅→
Вычислим значения параметров функционала
Ta if a0 0≠
a2
a0
, 0,
⎛
⎜
⎝
⎞
⎟
⎠
:=
Ta 1=
Tb if b0 0≠
b2
b0
, 0,
⎛
⎜
⎝
⎞
⎟
⎠
:=
Tb 4−=
α
aifa00≠ a1 0≠∧
a0
a1
− Ka()⋅, 0,
⎛
⎜
⎝
⎞
⎟
⎠
:=
α
a 1−=
α
bifb00≠ b1 0≠∧
b0
b1
Kb()⋅, 0,
⎛
⎜
⎝
⎞
⎟
⎠
:=
α
b 0.049787=
qa if a0 0 a1 0≠∧
a2
a1
− Ka()⋅, 0,
⎛
⎜
⎝
⎞
⎟
⎠
:=
qa 0=
qb if b0 0 b1 0≠∧
b2
b1
− Kb()⋅, 0,
⎛
⎜
⎝
⎞
⎟
⎠
:=
qb 0=
На йдем коэффициенты системы уравнений AC=B для определения
оптимальных коэффициентов пробных решений C
k
Nkx,()
x
Vk x,()
d
d
:=
i 1 n..:=
()[]
()[]()[]
),(),0(),(),0(
),()(),0()(),(),0()(:
1
aiVqaTaaVabiVqbTbbVb
dxxiVxgxVxxiNxNxKB
b
a
i
⋅−−⋅−⋅+−⋅−
−⋅+⋅+⋅⋅=
∫−
αα
β
B
T
5.784685 5.613800 5.551399 5.522788 5.507581()=
i1n..:=
j1n..:=
[]
),(),(),(),(
),(),()(),(),()(:,
11
ajVaiVabjVbiVb
dxxjVxiVxxjNxiNxKA
b
a
ji
⋅⋅+⋅⋅+
+⋅⋅+⋅⋅=
∫−−
αα
β

Страницы