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

UptoLike

ВУЗ: 

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

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

Рубрика: 

119
v2
0
1:=
load1 a v1,()
if a0
0≠
a2
a0
a1
a0
v1
0
⋅−, v1
0
,
⎛
⎜
⎝
⎞
⎟
⎠
if a0 0≠ v1
0
,
a2
a1
a0
a1
v1
0
⋅−,
⎛
⎜
⎝
⎞
⎟
⎠
⎛
⎜
⎜
⎜
⎝
⎞
⎟
⎟
⎟
⎠
:=
load2 b v2,()
if b0
0≠ v2
0
,
b2
b1
b0
b1
v2
0
⋅−,
⎛
⎜
⎝
⎞
⎟
⎠
if b0 0≠
b2
b0
b1
b0
v2
0
⋅−, v2
0
,
⎛
⎜
⎝
⎞
⎟
⎠
⎛
⎜
⎜
⎜
⎝
⎞
⎟
⎟
⎟
⎠
:=
score xf y,()y:=
S bvalfit v1 v2, a, b,
ab+
2
, D, load1, load2, score,
⎛
⎜
⎝
⎞
⎟
⎠
:=
S 0.153224 0.062349()=
Первый коэффициент вектора S – точ ное значение y(a), если a0=0 и y'(a) в
противном случае, второй коэффициент – то чное значение y(b), если b0=0 и
y'(b) в противном случае. Т. е. свели краевую задачу к задаче Коши.
На йдем решение этой задачи Коши, разбив отрезок на N частей
y
0
if a0 0≠
a2
a0
a1
a0
S
0
⋅−, S
0
,
⎛
⎜
⎝
⎞
⎟
⎠
:=
y
1
if a0 0≠ S
0
,
a2
a1
a0
a1
S
0
⋅−,
⎛
⎜
⎝
⎞
⎟
⎠
:=
Введите число точек разбиения отрезка [a, b]
N 100:=
Получ им матрицу решений
Y rkfixed y a, b, N, D,():=
Вычислим значения компьютерного решения в 10 промежуточных то чках
n 010..:=
T
n 0,
Y
N
10
nba−()⋅ 0,
:=
T
n 1,
Y
N
10
nba−()⋅ 1,
:=
Матрица компьютерного решения, где первая строка значения переменной
x
,
вторая –
y
, имеет вид
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
=
062349.0402849.0632877.0779347.0863987.0
0.19.08.07.06.0
904302.0914371.0905491.0886703.0865220.0846776.0
5.04.03.02.01.00.0
T
T
График точ ного компьютерного решения

Страницы