Действующий
Таблица В.9 - Функция частотной коррекции 
для общей вибрации в зданиях во всех направлениях (на основе ИСО 2631-2)
для общей вибрации в зданиях во всех направлениях (на основе ИСО 2631-2)n | Среднегеометрическая частота, Гц | Характеристика полосового фильтра | Функция частотной коррекции W_m | Допуск | |||||||||||
Номинальная | Истинная | Модуль | дБ | Фаза, ...° | Модуль | дБ | Фаза, ...° | Модуль, % | дБ | Дельта фи_0,...° | |||||
-10 | 0,1 | 0,1000 | 0,01585 | -36,00 | 169,7 | 0,01584 | -36,00 | 168,7 | + 26/-100 | + 2/ - бескон. | + бескон. / - бескон. | ||||
-9 | 0,125 | 0,1259 | 0,02511 | -32,00 | 166,9 | 0,02510 | -32,00 | 165,7 | + 26/-100 | + 2/ - бескон. | + бескон. / - бескон. | ||||
-8 | 0,16 | 0,1585 | 0,03978 | -28,01 | 163,5 | 0,03976 | -28,01 | 161,9 | + 26/-100 | + 2/ - бескон. | + бескон. / - бескон. | ||||
-7 | 0,2 | 0,1995 | 0,06297 | -24,02 | 159,1 | 0,06293 | -24,02 | 157,1 | + 26/-100 | + 2/ - бескон. | + бескон. / - бескон. | ||||
-6 | 0,25 | 0,2512 | 0,0995 | -20,04 | 153,4 | 0,09941 | -20,05 | 150,8 | + 26/-100 | + 2/ - бескон. | + бескон. / - бескон. | ||||
-5 | 0,315 | 0,3162 | 0,1565 | -16,11 | 146,0 | 0,1563 | -16,12 | 142,8 | + 26/-100 | + 2/ - бескон. | + бескон. / - бескон. | ||||
-4 | 0,4 | 0,3981 | 0,2436 | -12,27 | 136,3 | 0,2430 | -12,29 | 132,2 | + 26/-100 | + 2/ - бескон. | + бескон. / - бескон. | ||||
-3 | 0,5 | 0,5012 | 0,3699 | -8,64 | 123,6 | 0,3684 | -8,67 | 118,6 | + 26/-100 | + 2/ - бескон. | + бескон. / - бескон. | ||||
-2 | 0,63 | 0,6310 | 0,5336 | -5,45 | 107,7 | 0,5304 | -5,51 | 101,3 | + 26/-21 | + 2/-2 | +12/-12 | ||||
-1 | 0,8 | 0,7943 | 0,7071 | -3,01 | 89,36 | 0,7003 | -3,09 | 81,40 | + 26/-21 | + 2/-2 | +12/-12 | ||||
0 | 1 | 1,000 | 0,8457 | -1,46 | 71,00 | 0,8329 | -1,59 | 42,49 | + 26/-21 | + 2/-2 | + 12/-12 | ||||
1 | 1,25 | 1,259 | 0,9291 | -0,64 | 54,98 | 0,9071 | -0,85 | 26,56 | + 12/-11 | + 1/-1 | + 6/-6 | ||||
2 | 1,6 | 1,585 | 0,9699 | -0,27 | 42,14 | 0,9342 | -0,59 | 12,83 | + 12/-11 | + 1/-1 | + 6/-6 | ||||
3 | 2 | 1,995 | 0,9877 | -0,11 | 32,17 | 0,9319 | -0,61 | 0,5459 | + 12/-11 | + 1/-1 | + 6/-6 | ||||
4 | 2,5 | 2,512 | 0,9950 | -0,04 | 24,39 | 0,9101 | -0,82 | -10,89 | + 12/-11 | + 1/-1 | + 6/-6 | ||||
5 | 3,15 | 3,162 | 0,9980 | -0,02 | 18,20 | 0,8721 | -1,19 | -21,86 | + 12/-11 | + 1/-1 | + 6/-6 | ||||
6 | 4 | 3,981 | 0,9992 | -0,01 | 13,15 | 0,8184 | -1,74 | -32,52 | + 12/-11 | + 1/-1 | + 6/-6 | ||||
7 | 5 | 5,012 | 0,9997 | 0,00 | 8,884 | 0,7498 | -2,50 | -42,85 | + 12/-11 | + 1/-1 | + 6/-6 | ||||
8 | 6,3 | 6,310 | 0,9999 | 0,00 | 5,135 | 0,6692 | -3,49 | -52,73 | + 12/-11 | + 1/-1 | + 6/-6 | ||||
9 | 8 | 7,943 | 0,9999 | 0,00 | 1,680 | 0,5819 | -4,70 | -62,07 | + 12/-11 | + 1/-1 | + 6/-6 | ||||
10 | 10 | 10,00 | 0,9999 | 0,00 | -1,680 | 0,4941 | -6,12 | -70,84 | + 12/-11 | + 1/-1 | + 6/-6 | ||||
11 | 12,5 | 12,59 | 0,9999 | 0,00 | -5,135 | 0,4114 | -7,71 | -79,15 | + 12/-11 | + 1/-1 | + 6/-6 | ||||
12 | 16 | 15,85 | 0,9997 | 0,00 | -8,884 | 0,3375 | -9,44 | -87,25 | + 12/-11 | + 1/-1 | + 6/-6 | ||||
13 | 20 | 19,95 | 0,9992 | -0,01 | -13,15 | 0,2738 | -11,25 | -95,45 | + 12/-11 | + 1/-1 | + 6/-6 | ||||
14 | 25 | 25,12 | 0,9980 | -0,02 | -18,20 | 0,2203 | -13,14 | -104,2 | + 12/-11 | + 1/-1 | + 6/-6 | ||||
15 | 31,5 | 31,62 | 0,9950 | -0,04 | -24,39 | 0,1760 | -15,09 | -114,0 | + 12/-11 | + 1/-1 | + 6/-6 | ||||
16 | 40 | 39,81 | 0,9877 | -0,11 | -32,17 | 0,1396 | -17,10 | -125,7 | + 12/-11 | + 1/-1 | + 6/-6 | ||||
17 | 50 | 50,12 | 0,9699 | -0,27 | -42,14 | 0,1093 | -19,23 | -139,8 | + 12/-11 | + 1/-1 | + 6/-6 | ||||
18 | 63 | 63,10 | 0,9291 | -0,64 | -54,98 | 0,08336 | -21,58 | -156,9 | + 12/-11 | + 1/-1 | + 6/-6 | ||||
19 | 80 | 79,43 | 0,8457 | -1,46 | -71,00 | 0,06036 | -24,38 | -176,1 | + 26/-21 | + 2/-2 | +12/-12 | ||||
20 | 100 | 100,0 | 0,7071 | -3,01 | -89,36 | 0,04013 | -27,93 | -195,1 | + 26/-21 | + 2/-2 | +12/-12 | ||||
21 | 125 | 125,9 | 0,5336 | -5,46 | -107,7 | 0,02407 | -32,37 | -211,5 | + 26/-21 | + 2/-2 | +12/-12 | ||||
22 | 160 | 158,5 | 0,3699 | -8,64 | -123,6 | 0,01326 | -37,55 | -224,6 | + 26/-100 | + 2/ - бескон. | + бескон. / - бескон. | ||||
23 | 200 | 199,5 | 0,2436 | -12,27 | -136,3 | 0,006937 | -43,18 | -234,7 | + 26/-100 | + 2/ - бескон. | + бескон. / - бескон. | ||||
24 | 250 | 251,2 | 0,1565 | -16,11 | -146,0 | 0,003541 | -49,02 | -242,3 | + 26/-100 | + 2/ - бескон. | + бескон. / - бескон. | ||||
25 | 315 | 316,2 | 0,09950 | -20,04 | -153,4 | 0,001788 | -54,95 | -248,3 | + 26/-100 | + 2/ - бескон. | + бескон. / - бескон. | ||||
26 | 400 | 398,1 | 0,06297 | -24,02 | -159,1 | 0,000899 | -60,92 | -252,8 | + 26/-100 | + 2/ - бескон. | + бескон. / - бескон. | ||||
Для получения среднеквадратичного значения корректированного ускорения 
могут быть использованы любые методы частотного анализа [аналоговые или цифровые, в реальном масштабе времени, в третьоктавных полосах или с использованием быстрого преобразования Фурье (БПФ)], позволяющие получить среднеквадратичные значения спектральных составляющих ускорения 
для их последующего возведения в квадрат и суммирования с соответствующими весовыми коэффициентами:
могут быть использованы любые методы частотного анализа [аналоговые или цифровые, в реальном масштабе времени, в третьоктавных полосах или с использованием быстрого преобразования Фурье (БПФ)], позволяющие получить среднеквадратичные значения спектральных составляющих ускорения для их последующего возведения в квадрат и суммирования с соответствующими весовыми коэффициентами:
, (C.1)
- значение функции частотной коррекции на среднегеометрической частоте i-й полосы частот.
Примечание - Анализ вибрации в частотной области нельзя применять для получения значения дозы вибрации VDV. Его нельзя использовать также для получения текущего среднеквадратичного значения корректированного ускорения по данному стандарту, поскольку время интегрирования или постоянная времени интегрирования (1 с) сопоставимы с постоянной времени фильтра (величиной, обратно пропорциональной ширине полосы частот).
Значения среднегеометрических частот при третьоктавном анализе - в соответствии с таблицами В.1-В.9 (приложение В). В диапазон анализа должны входить по крайней мере по одной третьоктавной полосе выше и ниже номинального диапазона, определяемого границами 
и 
таблицы 3.
и таблицы 3.
Перед проведением операций возведения в квадрат и суммирования по формуле (С.1) значения ускорения следует умножить на соответствующее значение функции частотной коррекции (см. 5.6 и таблицы В.1-В.9 приложения В).
Среднеквадратичное значение корректированного ускорения может быть получено на основе среднеквадратичных значений составляющих БПФ сигнала вибрации в соответствии с формулой (С.1) или на основе составляющих спектральной плотности мощности PSD в соответствии с формулой (С.2). В каждом из этих случаев весовые коэффициенты (значения функции частотной коррекции) следует рассчитывать по формулам 5.6, а не брать из таблиц приложения В:
, (C.2)
При суммировании спектральных составляющих следует принимать во внимание эффект просачивания, связанный с использованным для формирования выборки временным окном. При анализе в широкой полосе частот рассчитанное в соответствии с формулой (С.2) среднеквадратичное значение корректированного ускорения 
следует разделить на коэффициент, соответствующий ширине полосы эквивалентного идеального фильтра, пропускающего белый шум той же мощности (см. таблицу С.1).
следует разделить на коэффициент, соответствующий ширине полосы эквивалентного идеального фильтра, пропускающего белый шум той же мощности (см. таблицу С.1).Форма окна* | Коэффициент | Типичное применение |
| Хэннинг | 1,50 | Общего применения, нестационарный случайный процесс |
| Плосковершинное | 3,77 | Периодический или синусоидальный сигнал (например, при калибровке) |
| * В других типичных ситуациях возможно применение других временных окон. | ||
Разрешение по частоте для БПФ должно быть менее 40% (предпочтительно 20%) нижней границы номинального диапазона частот. Частота выборки должна по крайней мере в 5 раз превышать верхнюю границу номинального диапазона частот.
Оценка сигналов ускорения с точки зрения их воздействия на человека включает в себя процедуру частотной коррекции с использованием одного из фильтров, описанных в 5.6. При определении среднеквадратичного значения корректированного ускорения процедура умножения на функцию частотной коррекции может быть осуществлена как до интегрирования, так и после вычисления среднеквадратичных значений спектральных составляющих - результат будет одним и тем же. Но для определения таких параметров, как максимальное кратковременное среднеквадратичное значение MTTV [см. формулы (3) и (4)], необходимо измерять максимальное значение текущего среднеквадратичного значения ускорения, что требует выполнения процедуры умножения на функцию частотной коррекции до интегрирования по времени.
Применение при анализе во временной области цифровых фильтров позволяет ограничить использование дорогостоящих и громоздких (особенно в многоканальных системах) аналоговых устройств.
Подобно тому, как для построения аналоговых фильтров в частотной области используется преобразование Лапласа, для цифровой фильтрации часто применяют z-преобразование. Передаточная функция цифрового фильтра может быть представлена в виде его z-преобразования H(z). В z-области преобразование Y(z) выходного сигнала цифрового фильтра связано с z-преобразованием входного сигнала Х(z) формулой
. (C.3)
, (C.4)
и 
- постоянные коэффициенты;
, (С.5)
и 
- выборочные значения входного и выходного сигналов соответственно в момент времени 
.
Коэффициенты фильтра 
и 
могут быть получены методом билинейного преобразования или методом импульсных инвариантов [5]. Метод билинейного преобразования наилучшим образом подходит для фильтров Баттерворта, в частности фильтров верхних и нижних частот, описанных в 5.6; z - преобразование этих двухполюсных фильтров может быть получено из преобразования Лапласа передаточной функции в 5.6 заменой переменной Лапласа s:
и могут быть получены методом билинейного преобразования или методом импульсных инвариантов [5]. Метод билинейного преобразования наилучшим образом подходит для фильтров Баттерворта, в частности фильтров верхних и нижних частот, описанных в 5.6; z - преобразование этих двухполюсных фильтров может быть получено из преобразования Лапласа передаточной функции в 5.6 заменой переменной Лапласа s:
, (C.6)
- период выборки.
Аналогичный расчет или альтернативный метод импульсных инвариантов может быть использован для переходных и ступенчатых фильтров.
Фильтры с бесконечной импульсной характеристикой по очереди применяют к последовательности оцифрованных данных в порядке их поступления в соответствии с уравнением (С.5).
В качестве примера ниже приведена программа, моделирующая процедуру фильтрации для функции частотной коррекции 
в кодах MATLABR, где используют встроенную функцию 'fillter.m', а также стандартные функции анализа сигналов 'butter.m' и 'bilinear.m'*.
в кодах MATLABR, где используют встроенную функцию 'fillter.m', а также стандартные функции анализа сигналов 'butter.m' и 'bilinear.m'*.
Примечание - Приведенная на рисунке С.1 программа требует, чтобы частота выборки по крайней мере в 9 раз превышала значение верхней частоты диапазона измерений f_2 (по таблице 2), чтобы удовлетворить требованиям настоящего стандарта по допуску на функцию частотной коррекции. Требование к частоте выборки можно понизить, если изменить программу соответствующим образом, например использовать преобразование
,
- частота среза фильтра.| function у - isofilwk(x,fs) | |||
| %ISOFILWK | |||
| Filter ISO 8041 Wk, whole body, vertical direction | |||
| % | у = isofilwk(x, fs); | ||
| % | у output signal vector, acceleration | ||
| % | x input signal vector, acceleration | ||
| % | fs sampling frequency Hz | ||
| % | bilinear transformation 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*pi*f6; | |||
| nyq = fs/2; | % Nyquist frequency | ||
| % | |||
| % band limiting high pass and low pass | |||
| % | |||
| [b1,a1] = butter(2, f1/nyq, 'high1); | % High pass | ||
| [b2,а2] =butter(2, f2/nyq); | % Low pass | ||
| % | |||
| % a-v transition | |||
| % | |||
| B3 = [1/w3 1]; | |||
| A3 = [1/w4/w4 1/Q4/w4 1]; | |||
| [b3,а3] =biIinear(B3,A3,fs); | |||
| % | |||
| % upward step | |||
| % | |||
| B4 = [1/w5/w5 1/Q5/w5 1]*w5*w5/w6/w6; | |||
| A4 = [1/w6/w6 1/Q6/w6 1]; | |||
| [b4,а4] = bilinear(B4,A4,fs); | |||
| % | |||
| % Apply filter to input signal vector x (output to signal vector y) | |||
| % | |||
| y = filter(b2,a2,x); | |||
| y = filter(b1/a1/y); | |||
| y = filter(b3,a3,y); | |||
| y = filter(b4,a4,y); | |||