简单的回波对消程序（LMS）

Agoni

1. 条件：不考虑近端干扰，不设双端检测器，自适应算法（LMS）
   
2. %% 回波对消  简单的仿真实例
   
3. % 仿真语音对消
   
4. % 不考虑近端干扰
   
5. % 不设双端检测器
   
6. % 采用LMS算法
   
7. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   
8. clc
   
9. clear all
   
10. close all
    
11. 
    
12. Fs = 1000;
    
13. t = 0:1/Fs:2;
    
14. fc = 10;
    
15. x = sin(2*pi*fc*t);                 % 以sin信号作为远端信号
    
16. h = [0.3 0.5 0.3];                  % 回波干扰冲激响应
    
17. y = filter(h,1,x);                  % 回波
    
18. y = awgn(y,30,'measured');          % 加噪声
    
19. 
    
20. M = 21;
    
21. mu = 0.04;
    
22. len = length(x);             % 信号的长度
    
23. output = zeros(1, len);      % LMS滤波器的输出
    
24. 
    
25. wn = zeros(M, len);          % 初始化权向量wn
    
26. ref = zeros(M, 1);
    
27. 
    
28. for n = 1:len-1
    
29.   
    
30.     if n < M
    
31.         ref = [y(n:-1:1), zeros(1, M - n)].';
    
32.     else
    
33.         ref = y(n:-1:n - M + 1).';
    
34.     end
    
35.     output(n) = x(n) - wn(:, n)' * ref;
    
36.     wn(:,n+1) = wn(:,n) + mu*ref*conj(output(n));
    
37. end
    
38. ref = y(len:-1:len-M+1).';
    
39. output(len) = x(len) - wn(:, len)' * ref;
    
40. plot(abs(output))

复制代码

转自：http://blog.sina.com.cn/s/blog_5def5a660100u2hn.html