matlab的LMS自适应滤波算法去噪

设定：一个房间中两个麦克风，一个放在远处采集房间噪声，一个放在说话人附近采集带噪语音信号，认为两个音频文件的噪声相似。目标是使用LMS自适应滤波算法来抑制噪声还原语音。 仿真：现给定一录音.mat文件，其中：s是原音频内容；ref_noise是均值为0，方差为1的高斯噪声；mixed是叠加上高斯噪声序列；fs为信号采样率。要求使用LMS自适应滤波法抑制噪声。

clear all

load（‘测试文件.mat‘）;
audiowrite（‘mixed.mp4‘mixedfs）;
audiowrite（‘signal_s.mp4‘sfs）;

num = 547824;

p = 30;

X = zeros（p1）;

u = 1e-2;

W = 0.1*rand（p1）;

e = zeros（p1）;


for k = p:num
   
 for i = 1:p
 
    X（i） = mixed（k-i+1）; 
     
 end

 y（k） = X‘*W;
 
 e（k） =ref_noise（k）-y（k）;
 
 W = W+2*u*e（k）*X;
  
end

y_new = mixed‘-y;

audiowrite（‘signal_new.mp4‘efs）;


% for k = p:num
%  
%     X（k:） = [mixed（k）mixed（k-1）mixed（k-2）mixed（k-3）mixed（k-4）mixed（k-5）mixed（k-6）mixed（k-7）mixed（k-8）mixed（k-p+1）];
% 
% end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        588  2018-05-02 16:36  pj6\pj6.asv

     文件        611  2018-05-02 16:50  pj6\pj6.m

     文件   13691011  2018-04-22 16:56  pj6\测试文件.mat

     目录          0  2019-09-16 16:32  pj6

     文件     204273  2018-05-02 16:50  pj6\mixed.mp4

     文件     207354  2018-04-26 17:27  pj6\signal.mp4

     文件     204816  2018-04-26 17:27  pj6\signalnew.mp4

     文件     204706  2018-05-02 16:50  pj6\signal_new.mp4

     文件     207354  2018-05-02 16:50  pj6\signal_s.mp4

     文件     204273  2018-04-26 17:27  pj6\zao.mp4

----------- ---------  ---------- -----  ----

             14924986                    10