ГОСТ ИСО 8041—2006
2x1с
2 - 1
s = -
l
-
g
--
(
-
*
--
/
-
0
=
r
-
J
---
)
-
2
----
г
---
1
. где 1. — частота среза фильтра.
function у = isofilwk(x.fs)
<*lSOFILW K
Filter ISO 8041 Wk. whole body, vertical direction
%у = isofilwklx. fs) ;
%у output signal vector, acceleration
%x input signal vector, acceleration
%fs sampling frequency Hz
°tbilinear transformation algorithm is used
П = 0.4.
f2 = 100;
f3 = 12.5;
f4 = 12.5;
Q4 = 0.63;
f5 = 2.37;
05 = 0.91;
f6 = 3.35;
06 = 0.91;
w3 = 2*pi*f3,
w4 = 2*pi*f4.
w5 = 2*pi*f5;
w6 = 2*pi*fl>.
% Nyquist frequency
nyq = fs/2;
%
band lim iting high pass and low pass
%
|bl.al| = butter(2. П/nyq. ’highI):
% High pass
% Low pass
|b2,a2] = butlcrl2, f2/nyq);
*
% a-v transition
%
B3 = |l/w3 11;
A3 = [l/w 4/w 4 1/Q4/W4 l|;
|ЬЗ.аЗ| = bilinear! B3,A3.fs),
%
°t upward step
%
B4 =
|1 /
w
5 /
w
5
!/Q 5/w 5 I |*w5*w5/w6/w6;
A4 = |l/w 6/w 6 l/Q 6/w6 I );
|b4.a4| = bilinear) B4.A4.fs):
%
% Apply filter to input signal vector x (output to signal vector y)
%
у = filtcr(b2.a2.x);
у = filtcr(bl/al/y);
у = filtcr!b3,a3.y);
у = filtcr(b4,a4,y):
Рисунок C.1 — Пример программы, реализующей фильтр
понизить, если изменить программу соответствующим образом, например использовать преобразование
61