维纳滤波 MATLAB程序

维纳滤波的程序。有注解，能够实现对语音滤波并输出对比波形

clear all;
clc
load（‘chirp‘‘Fs‘‘y‘）;
p=audioplayer（yFs）;
play（p）;
subplot（211）;
plot（y）;
title（‘原始信号‘）
Py=fft（y）;
subplot（212）;
plot（abs（Py））;
title（‘原始信号频谱‘）     %原始信号时域图、频域图及信号的播放
disp（‘按任意键继续‘）;
pause（）;                  %按任意键程序继续执行
N=length（y）;
% sigma=0.1; 			  %控制噪声强度
% z=zeros（N1）;
% v=randn（N1）*sigma;     %产生噪声
% z=y+v;
z=awgn（y30）;             %将高斯白色噪声添加到信号中.产生高斯白噪声加噪信号z
                          %标量30指定了每一个采样点信号与噪声的比率，单位为dB
q=audioplayer（zFs）;
play（q）;
figure（）;                 %创建新图像窗口
subplot（211）;
plot（z）;
title（‘观测信号‘）
Pz=fftn（z）;
subplot（212）;
plot（abs（Pz））;
title（‘观测信号的频谱‘）     %观测信号时域图、频域图及信号的播放
%加噪后均方误差
err1=mean（（z-y）.^2）;           %对（z-y）列计算（a^2+b^2+...n^2）/n----
disp（[‘加噪后均方误差为‘num2str（err1）]）;      %把数值转换成字符串
disp（‘按任

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

     文件      25276  2004-03-14 10:31  1\chirp.mat

     文件       1688  2018-12-04 19:02  1\winner.m

     目录          0  2018-12-28 15:34  1

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

                26964                    3