feifei-weinalvbo.m

信号从发送者传送到接收者往往受到几种形式的变形而消弱

close all
clc;
load('chirp','y','Fs');
sound(y);
M=length(y);
z=awgn(y,10);   %加噪
figure(1);
plot(y);
title('语音信号y的时域图');
figure(2);
plot(z);
title('加噪后的信号');
rZ=xcorr(z,M-1,'unbiased');         %观测数据的自相关估计
rY=xcorr(y,M-1,'unbiased');         %原始数据的自相关估计
Rz=zeros(M,N);
for k=0:N-1
    Rz(:,k+1)=rZ(M-k:M+M-1-k);
end
H_opt=pinv(Rz)*rY(M:2*M-1);         %用广义逆的方法解上面的方程 
Z=filter(H_opt,1,z);        %加噪信号通过滤波器后的输出
figure(3);
plot(Z);title('经滤波后的信号');
pause(3);
soundsc(Z);