基于RLS算法的语音噪声回声消除处理MATLAB代码

基于RLS算法的语音噪声回声消除处理MATLAB代码，可以对有回声或者有噪音的语音信号进行自适应滤波处理重建原始语音信号。

clc
clear all
close all
[s fs] = audioread（‘1.wav‘）; 
N=length（s）;
time=（0:N-1）/fs;
figure（1）;   
plot（times）;
xlabel（‘时间‘）;ylabel（‘幅度‘）;    % x轴、y轴标题
title（‘原始信号波形‘）;
hold on 

a=mean（s）;
s=s-a;
b=max（abs（s））;
s=s/b;
% s=awgn（s10）;
% s1=s.*b;
% s1=s1+a;
% audiowrite（‘2.wav‘s1fs）;

% 
% clean=s‘;
% ref_noise=.1*randn（1length（s））;
% mixed = clean+ref_noise;
[echo f1s] = audioread（‘1echo.wav‘）; 


% s1=mixed‘*b;
% s1=s1+a;
figure（2）; 
plot（timeecho）;
xlabel（‘时间‘）;ylabel（‘幅度‘）;    % x轴、y轴标题
title（‘有回声信号波形‘）;
hold on 
audiowrite（‘2.wav‘echof1s）;

a=mean（echo）;
echo=echo-a;
b=max（abs（echo））;
echo=echo/b;

mu=0.05;M=2;espon=1e-4;
% [enwnyn]=lmsFunc（muMref_noisemixed）;
% [enwnyn]=nlmsFunc（muMref_noisemixedespon）;
delta = 1e-7;
lambda = 1;
[enwy]=rls（lambdaMsechodelta）;


s2=y‘*b;
s2=s2+a;
figure（3）; 
plot（times2）;
xlabel（‘时间‘）;ylabel（‘幅度‘）;    % x轴、y轴标题
title（‘去回声后信号波形‘）;
hold on 
audiowrite（‘3.wav‘s2fs）;

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

     文件    1097778  2018-11-21 21:34  回声处理\1.mp3

     文件     558196  2018-11-22 20:08  回声处理\1.wav

     文件     558196  2018-11-22 20:13  回声处理\1echo.wav

     文件     558196  2018-12-01 14:26  回声处理\2.wav

     文件     558196  2018-12-01 14:26  回声处理\3.wav

     文件       1048  2018-12-01 14:31  回声处理\audio.m

     文件        141  2018-11-22 19:58  回声处理\getaudio.m

     文件        162  2018-11-22 20:13  回声处理\getecho.m

     文件        858  2018-11-22 16:56  回声处理\rls.m

     文件       3000  2018-11-22 20:57  回声处理\Untitled.m

     文件        632  2018-11-22 16:43  回声处理\wavread.asv

     文件      29492  2018-12-01 14:27  回声处理\原始信号波形.png

     文件      30219  2018-12-01 15:48  回声处理\去回声后信号波形.png

     文件      27832  2018-12-01 14:28  回声处理\有回声信号波形.png

     目录          0  2018-12-01 15:48  回声处理

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

              3423946                    15