%寫上標(biāo)題 
%設(shè)計(jì)低通濾波器： 
[N,Wc]=buttord() 
%估算得到Butterworth低通濾波器的最小階數(shù)N和3dB截止頻率Wc 
[a,b]=butter(N,Wc); %設(shè)計(jì)Butterworth低通濾波器 
[h,f]=freqz(); %求數(shù)字低通濾波器的頻率響應(yīng) 
figure(2); % 打開窗口2 
subplot(221); %圖形顯示分割窗口 
plot(f,abs(h)); %繪制Butterworth低通濾波器的幅頻響應(yīng)圖 
title(巴氏低通濾波器''); 
grid; %繪制帶網(wǎng)格的圖像 
sf=filter(a,b,s); %疊加函數(shù)S經(jīng)過(guò)低通濾波器以后的新函數(shù) 
subplot(222); 
plot(t,sf); %繪制疊加函數(shù)S經(jīng)過(guò)低通濾波器以后的時(shí)域圖形 
xlabel('時(shí)間 (seconds)'); 
ylabel('時(shí)間按幅度'); 
SF=fft(sf,256); %對(duì)疊加函數(shù)S經(jīng)過(guò)低通濾波器以后的新函數(shù)進(jìn)行256點(diǎn)的基—2快速傅立葉變換 
w= %新信號(hào)角頻率 
subplot(223); 
plot()); %繪制疊加函數(shù)S經(jīng)過(guò)低通濾波器以后的頻譜圖 
title('低通濾波后的頻譜圖');