ГОСТР 54713—2011
-
2 Вычислить комплексный спектр входного сигнала.
Во-первых,
s (0
взвешивается окном Хана. т. е.
sw(i)
=
s(i)
‘
(0.5-0.5 * cos((p/ *(’+0.5)У ibten).
Во-вторых, вычисляется стандартное БПФ
sw(i).
В-третьих, вычисляется полярное представление преобразования, г (iv) и / (w) соответствуют амплитудным
и фазовым частям преобразованного
sw(i)
соответственно.
3 Вычислить предсказание г
(w)
и
( (w).
Предсказанные амплитуда
г_pred
(iv) и фаза
f_pred
(w) вычисляются по двум предыдущим блокам расчета
порога
r(w)
и
f (w):
г_pred (w)
= 2.0 * г (М
)-г
(<-2)
1_pred (w)
=2.0 *
f
(f-1 )-/
(
1 2
),
где
t
— текущий номер блока, 1-1 соответствует данным предыдущего блока и
1-2
соответствует данным
блока перед предыдущим.
4 Вычислить меру неопределенности с
(w).
c(w) = (((г (w) * cos
(f(w ))
-
r_pred
(у/) * cos (
f_pred
(iv))) A2 +
(r
(w )’ sin(/(w)) - г
jx e d (w )
* sin(/jrred(w)) )A2)*0.5 )I (r(w) + abs (г_pred(w)).
Эта формула используется для каждого из коротких блоков с коротким БПФ. для первых 6 строк длинных
блоков мера неопределенности вычисляется по длинному БПФ. для оставшихся строк используется минималь ное
значение неопределенности для коротких блоков БПФ. Если необходимо уменьшить сложность вычислений,
неопределенность для верхней части спектра может быть установлена равной 0.4.
5 Вычислить энергию и неопределенность в разделах вычисления порога.
Энергия каждого раздела,
в
(Ь):
do for each partition b:
e(b)
=
0
do from lower index to upper index
w
ofpartition b e (b)
= e
(b)
+
r (w)
A2
end
do end
do
(e (b) используется в модуле
MIS
(см. Б.6.1): е (6) равно Xengy с X = [R. L. М. S]) и взвешенная неопределен
ность. с (£>):
do for each partition b: c(b)
=
0 do from loner index to upper index
iv
of partition b
c(b) = c
(b)
+
rfw)’’2
■ c
(w) end do end do
Разделы вычисления порога обеспечивают разрешение, примерно равное одной частотной линии БПФ
или 1/3 критической полосы. На нижних частотах одна линия БПФ соответствует одному разделу вычисления. На
верхних частотах несколько строк будут объединены в один раздел вычисления. Набор значений линий
разделов для каждой из трех частот дискретизации задан в таблицах Б.1 — Б.24. Элементы этих таблиц будут
использоваться в процессе вычисления порога маскирования. В каждой таблице есть следующие графы:
1 Индекс раздела вычисления. 6;
2 Самая низкая частотная пиния раздела. w
_
toiv
(by,
3 Самая высокая частотная линия раздела.
w_htgh (b):
4 Среднее значение частоты раздела, барк,
bval (b);
5 Абсолютный порог слышимости
qsthr (b).
6 Максимальное значение
Ь. Ьтах.
равное самому большому индексу, существующему для данной частоты
дискретизации.
6 Вычислить свертку энергии разделов и неопределенности с функцией маскирования.
for each partition b:
ecb(b)
=
0
do for each partition bb:
ecb (b)
=
ecb (b)
+e
(bb)
*
sprdngf (bval (bb) . bval (b))
end do end do do for each partition b:
ct (b) = 0
do for each partition bb:
ct(b)
=
ct(b) +c(bb)’ sprdngf (bval (bb) . bval (b))
end do end do
Поскольку cl (b) взвешивается энергией сигнала, оно должно быть повторно нормализовано к cb (Ь)
cb (b) =cl (b) / ecb (Ь)
Аналогично, из-за ненормализованной природы функции маскирования
ecbb
должно быть повторно нор
мализовано и нормализованная энергия
еп(Ь)
равна:
ел
(Ь)
=
ecb (Ь)
/
толп (Ь)
коэффициент нормализации, rnorni
(b):
do for each partition b tmp (b)
=
0
do for each partition bb
tmp (b)
=
Imp (b)
+
sprdngf (bval (bb) .bval (b)) end do
rnorm (b)
=
1/ tmp (b) end do
104