ГОСТ Р 56047— 2014
% выделение памяти
XwN2 = zeros (1, NF/2+-1);
% фильтрация на основе модели внешнего и среднего уха
Xw2(1,:) = W2 .* X2(1.1:NF/2+1);
Xw2(2,:) = W2 .* X2(2.1:NF/2+1);
for (k = 0:NF/2)
XwN2(k+1)= (Xw2(1.k+1) - 2 • sqrt (X\v2(1 .k+1) * Xw2(2.k+1))...
+ Xw2(2.k+1)):
end
Eb(1.:) = PQgroupCB (Xw2(1.:), Version):
Eb<2.:) = PQgroupCB (Xw2(2.:), Version);
EbN= PQgroupCB (XwN2. Version);
E(1.:) = Eb<1.:) + EIN:
E(2,:) = Eb(2.:) + EIN:
Es(1,:)= PQspreadCB (E(1,:), Version);
Es{2.:) = PQspreadCB (E(2.:>, Version):
%---------------------
function (Ehs. Ef] = PQ_timeSpread (Es. Ef)
persistent Nc a b
if (isempty (Nc))
[Nc. fc] = PQCB (’Basic’):
Fs = 48000;
NF = 2048;
Nadv = NF / 2:
Fss = Fs/Nadv;
1100 = 0.030:
tmin = 0.008;
[a. b) = PQtConst (1100. tmin. fc. Fss);
end
Ehs = zeros (1. Nc);
for (m = 0:Nc-1)
Ef(m+1) = a(m+1) * Ef(m+1) + b(m+1) ’ Es(m+1);
Ehs(m+1) = max(Ef(m+1). Es(m+1)):
end
function X2 = PQDFTFrame (x)
persistent h\v
NF = 2048:
if (isempty (hw))
Amax = 32768;
fc = 1019.5:
Fs = 48000;
Lp = 92;
GL = PQ_GL (NF. Amax. ftVFs. Lp);
hw = GL ’ PQHannWin (NF):
end
xw = hw .* x;
X = PQRFFT (xw. NF. 1);
X2 = PQRFFTMSq (X. NF);
%-------------------------------------------
function GL = PQ_GL (NF. Amax. fcN. Lp)
W = NF -1;
gp = PQ_gp (fcN. NF. W);
GL = 10*(Lp / 20) / (gp * Amax/4 * W):
33