lpp.m

来自「我的数字信号处理课程设计」· M 代码 · 共 23 行

[y,fs,bits]=wavread('b.wav');
%信号的频谱
m=fft(y,512);  %取y的512点fft变换。
p=(0:255)/256*(fs/2);  
%低通滤波器
wp=0.2;ws=0.4;rp=0.5;rs=40;
[N,Wc]=buttord(wp,ws,rp,rs);
[num,den]=butter(N,Wc);
[h,f]=freqz(num,den,256,1000);

lp=filter(num,den,y);
u=fft(lp,512);  %取y的512点fft变换。
lpp=(0:255)/256*(fs/2);  

%subplot(211);plot(y);%信号时域波形
%subplot(212);plot(p,abs([m(1:256)]));%信号的频谱图
%subplot(313);plot(f,abs(h));%低通滤波器的频谱图
%figure(2);
%subplot(211);plot(lp);%滤波后的信号时域波形
%subplot(212);
plot(lpp,abs([u(1:256)]));%滤波后的信号的频谱图

%wavplay(y,fs);