问题一:以下代码后续没有用到,但从注释来看,变量u和u1又是很重要,后边必须要使用的参量。这个矛盾怎么解释?
u = zeros(1,ITER); % Initialize input of adaptive filter to zerou1 = zeros(1,ITER); % Initialize input of feedforward path to zero
if n > (delay+1)u(n) = S.ff(end)*en(n-delay-1);elseu(n) = 0;end
问题二:自适应滤波器的输入信号是un,也就是系统纯净的语音信号,这与作者给出的系统逻辑图、问题一中定义的参量是不相符合的。
for k = 1:1:M-1; % Updates of tapped-delay line of adaptive filterfftap(M-k+1) = fftap(M-k);endfftap(1) = un(n);yn(n) = w'*fftap'; % Compute output by inner producten(n) = un(n)+yfb(n)-yn(n); % Compute error signalw = w + (mu*en(n))*fftap'; % LMS algorithm
问题一:以下代码后续没有用到,但从注释来看,变量u和u1又是很重要,后边必须要使用的参量。这个矛盾怎么解释?
u = zeros(1,ITER); % Initialize input of adaptive filter to zero
u1 = zeros(1,ITER); % Initialize input of feedforward path to zero
if n > (delay+1)
u(n) = S.ff(end)*en(n-delay-1);
else
u(n) = 0;
end
问题二:自适应滤波器的输入信号是un,也就是系统纯净的语音信号,这与作者给出的系统逻辑图、问题一中定义的参量是不相符合的。
for k = 1:1:M-1; % Updates of tapped-delay line of adaptive filter
fftap(M-k+1) = fftap(M-k);
end
fftap(1) = un(n);
yn(n) = w'*fftap'; % Compute output by inner product
en(n) = un(n)+yfb(n)-yn(n); % Compute error signal
w = w + (mu*en(n))*fftap'; % LMS algorithm
最后希望作者能够贴出书中原始的代码>~<。