Приложение Б
(обязательное)
ГОСТ Р 56066—2014
Требования к программному обеспечению
Б.1 Цифровая фильтрация зарегистрированных данных может быть выполнена во время обработки
данных и должна выполняться только один раз. Частота среза *
р
применяемого цифрового фильтра должна
соответствовать 5 Гц. В качестве цифровых фильтров могут быть применены преобразование Фурье или
четырехполюсный бесфазный цифровой фильтр Багтерворта.
Б.2 В качестве цифрового фильтра рекомендуется применять алгоритм прямого и обратного
преобразования Фурье с отсечением шумовой составляющей верхних частот спектра ускорения:
Прямое преобразование:
_12я/
к = 0
........
N -
1.
л-О
Обратное преобразование:
1
N
- 12х/
кп
п = 0
........
N -
1,
/v к -О
где
N
- число значений сигнала, измеренных за период, а также число компонентов разложения;
х„. (п
= 0.......
N -
1) - измеренные значения сигнала (в дискретных временных точках с номерами
п
= 0
........
N
- 1. которые являются входными данными для прямого преобразования и выходными для обратного:
х», (к
= 0
.......
N -
1) -
N
комплексных амплитуд синусоидальных сигналов, слагающих исходный сигнал;
являются выходными данными для прямого преобразования и входными для обратного;
■
-
к I
-----
- обычная (вещественная) амплитуда Л-ro синусоидального сигнала;
N
агд(х*) - фаза
к-
го синусоидального сигнала (аргумент комплексного числа);
к -
индекс частоты. Частота
к-
го сигнала равна Л/Г. где Г - период времени, в течение которого брались
входные данные.
Б.З Алгоритм цифрового бесфазового фильтра Баттерворта 4 порядка реализуется путем повторного
прохода через фильтр Баттерворта 2 порядка отфильтрованных данных. Алгоритм фильтра Баттерворта 2
порядка определяется разностным уравнением
т =оаХ[1-\]+а2X[t - 2]+ ЬхY[t - 1]+Ь2Y[t - 2].(Б.
1
)
где YW - поток отфильтрованных данных;
X I/]- поток исходных данных;
/-текущ ее время;
а;. а
1
. too.
b,. Wa,W. -
константы, зависящие от класса частот канала и определяемые следующими
выражениями:
IV, = 2п (КЧК) 2.0775.
W.
=
sin-(
W , •
772)
/
cos
(W< ■TI2).
a0 = W ;/(\,0 + p W a +W a2 ),
3i=2fti,
aj = Эй.
bt = -2 (W ;-1)/(1 + ^21A/i + W ;).
b: = (-! + f t Wa- W -)/(1 + J2Wa+ W ;):
T -
период выборки:
КЧК = 10 - для класса частоты канала 10.
7