ГОСТ Р 59701.1—2022
fu n c tio n у = i s o f i l w k (х ,f s )
%ISOFILWK
F ilte r ISO 8041 Wk, whole body, v e rtic a l direction
%у = i s o f i l w k ( x , f s ) ;
%у output signal vector, acceleration
% x input signal vector, acceleration
% fssam pling frequency Hz
%b ilin e a r transform ation algorithm is used
f1 = 0 . 4;
f2 = 100;
f3 = 12.5;
f4 = 12.5;
Q4 = 0.63;
f5 = 2.37;
Q5 = 0.91;
f6 = 3.35;
Q6 = 0.91;
w3 = 2*pi*f3;
w4 ==2*pi*f4;
w5 - 2*pi*f5;
w6 = 2 *p i*f6 ;
nyq = fs/2;% N yquist frequency
%
%band lim itin g high pass and low pass
%
[b l,a l] =
[b2,a2] =
b u tte r(2,fl/n y q ,’h ig h l);
b u tte r(2,f2/nyq);
% High pass
% Low pass
%
% a-vtra n s itio n
w
%
B3=[l/w 31 ];
A3=[1 /
w
4/
w
4l/Q4/w41 ];
[ЬЗ,аЗ] =b ilin e a r(B 3 ,A3, f s ) ;
%
% upwardstep
%
B4 = [1 /
w
5 / 5 l/Q5/w5 1 ]*w5*w5/w6/w6;
A4 = [1 /
w
6/
w
6 l/Q6/w6 1 ];
[b4,a4] = bilinear(B 4,A 4, f s ) ;
%
%Apply f ilt e r to input signal vector x (output to signal vector y)
%
у = filte r(b 2 ,a 2 ,x );
у = filte r ( b l/a l/y ) ;
у = f ilt e r (ЬЗ,аЗ,у);
у = f ilt e r (Ь4,а4,у);
Рисунок С.1 — Пример программы, реализующей частотную коррекцию Wk
68