%Samp6_5
Wp=2000*2*pi;Ws=3000*2*pi; %濾波器截止頻率
Rp=3;Rs=15;    %通帶波紋和阻帶衰減
Fs=10000;    %采樣頻率
Nn=128;    %調(diào)用freqz所用的頻率點(diǎn)數(shù)
[N,Wn]=buttord(Wp,Ws,Rp,Rs,'s');   %模擬濾波器的最小階數(shù)
[z,p,k]=buttap(N);      %設(shè)計(jì)模擬低通原型Butterworth濾波器
[Bap,Aap]=zp2tf(z,p,k);   %將零點(diǎn)極點(diǎn)增益形式轉(zhuǎn)換為傳遞函數(shù)形式
[b,a]=lp2lp(Bap,Aap,Wn);  %進(jìn)行頻率轉(zhuǎn)換
[bz,az]=impinvar(b,a,Fs);  %運(yùn)用脈沖響應(yīng)不變法得到數(shù)字濾波器的傳遞函數(shù)
figure(1)
[H,f]=freqz(bz,az,Nn,Fs);   %求解數(shù)字濾波器的幅頻特性和相頻特性
subplot(2,1,1),plot(f,20*log10(abs(H)))
xlabel('頻率/Hz');ylabel('振幅/dB');grid on;
subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)))
xlabel('頻率/Hz');ylabel('相位/^o');grid on;
figure(2)
f1=1000;f2=4000;   %輸入信號(hào)的頻率
N=100;    %數(shù)據(jù)長(zhǎng)度
dt=1/Fs;n=0:N-1;t=n*dt;   %采樣間隔和時(shí)間序列
x=sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t); %濾波器輸入信號(hào)
subplot(2,1,1),plot(t,x),title('輸入信號(hào)') %繪制輸入信號(hào)
y=filtfilt(bz,az,x);  %用函數(shù)filtfilt對(duì)輸入信號(hào)進(jìn)行濾波
y1=filter(bz,az,x);   %用filter函數(shù)對(duì)輸入信號(hào)濾波
subplot(2,1,2),plot(t,y,t,y1,':'),title('輸出信號(hào)'),xlabel('時(shí)間/s')
legend( ' filtfilt ', 'filter')  %加圖例