基于matlab各滤波器源代码,基于matlab各滤波器源代码

汉明窗的FIR低通滤波：

Fs=22050;

[x,FS,bits]=wavread('G:\Users\DP\Desktop\SoundTest.wav');%G:\Users\DP\Desktop \SoundTest.wav G:\Users\DP\Desktop\好铃网-湖人掘金宣传片.wav

x=x(:,1);

figure(1);

subplot(2,1,1);

plot(x);

%sound(x,FS,bits); %回放语音

title('语音信号时域波形图')

y=fft(x,3260);

f=(FS/3260)*[1:1630];

subplot(2,1,2);

plot(f(1:1630),abs(y(1:1630)));

title('语音信号频谱图');

%产生噪声信号并加到语音信号

t=0:length(x)-1;

zs0=0.05*cos(2*pi*100*t/22050);

figure(2);

subplot(2,1,1)

plot(zs0)

title('噪声信号波形');

zs1=fft(zs0,1200);

%sound(zs0,FS,bits);%回放噪音

f=(FS/1200)*[1:600];

subplot(2,1,2)

plot(f(1:600),abs(zs1(1:600)));

title('噪声信号频谱');