Составители:
5. О программе сжатия-восстановления числовых потоков
Программа содержит ряд процедур, которые описаны ниже.
1. FUN(t:extended): — аппроксимируемая функция
extended; t — аргумент (extended);
2. DifFUN(t:extended): — производная аппроксимируемой
extended; функции;
3. APPR_B_spl — линейная комбинация
(var cm2,cm1,c0, трех B-сплайнов ω
−2
, ω
−1
, ω
0
,
t,xm1,x0,x1,x2: eu(t) = c
−2
ω
−2
(t) + c
−1
ω
−1
(t)+
t,xm1,x0,x1,x2: +c
0
ω
0
(t) при t ∈ (x
0
, x
1
);
extended):extended; здесь участвуют узлы
x
−1
, x
0
, x
1
, x
2
. Вычисления
проводятся по формуле
APPR_B_spl= c
−2
(t−x
1
)
2
(x
1
−x
0
)(x
1
−x
−1
)
+c
−1
((x
−1
− x
1
− x
2
+ x
0
)t
2
−
−2(x
0
x
−1
− x
1
x
2
)t + x
−1
x
0
x
2
−
−x
−1
x
1
x
2
+ x
−1
x
0
x
1
− x
0
x
1
x
2
)
×((x
1
− x
−1
)(x
1
− x
0
)(x
2
− x
0
))
−1
+ c
0
(t−x
0
)
2
(x
1
−x
0
)(x
2
−x
0
)
;
4. fg(var x1,x2:extended): — вычисление функционала:
extended hg, ui = u(x
1
)+
+(x
2
− x
1
)u
0
(x
1
)/2;
5. pr_DECOMP_partly(var — процедура отыскания a
−1
и b
0
:
xm1,x0,x1,x2,x3, a
−1
= −c
−2
x
2
−x
1
x
1
−x
−1
+ c
−1
x
2
−x
−1
x
1
−x
−1
;
cm2,cm1,c0,c1,am1,b0 b
0
= (c
−2
(x
3
− x
1
)(x
2
− x
1
)−
:extended) −c
−1
(x
2
− x
−1
))+(x
1
− x
−1
)×
×(c
0
(x
3
− x
0
) − c
1
(x
1
− x
0
))×
×((x
3
− x
0
)(x
1
− x
−1
))
−1
;
6. PR_cm1c0(var — процедура отыскания c
−1
и c
0
:
xm1,x0,x1,x2,x3, c
−1
= (a
−2
(x
2
− x
1
)+
46