wsyn.m

来自「ITU-T G.723.1 Speech Coder: Matlab imple」· M 代码 · 共 33 行

function [y, WSynMem] = WSyn (x, WSynCof, WSynMem)
% Filter with three stages of filtering
% - All-pole filter (synthesis filter)
% - Pole-zero filter (formant weighting filter)
% - Recursive pitch filter (harmonic noise weighting filter)

% $Id: WSyn.m 1.1 2003/11/21 G.723.1-v2r1a $

[yAP, WSynMem.APMem] = PZFilter (1, WSynCof.aQ, x, WSynMem.APMem);

[yWF, WSynMem.WFMem] = PZFilter (WSynCof.bW, WSynCof.aW, yAP, WSynMem.WFMem);

[y, WSynMem.HNWMem] = HNWFilter (yWF, WSynCof.G, WSynCof.L, WSynMem.HNWMem);

return

%--------------------
function [y, HNWMem] = HNWFilter (x, G, L, HNWMem)
% FIR filter, lag L, gain -G

LMem = length (HNWMem.xp);

% Form the extended input vector
xe = [HNWMem.xp; x];

% Filter with the harmonic noise weighting filter
Nx = length (x);
y = x - G * xe(LMem-L+1:LMem-L+Nx);

HNWMem.xp = xe(end-LMem+1:end);

return