%Samp6_9
Fs=1000;     %采樣頻率
wp=[100 200]*2/Fs;   %通帶邊界頻率（歸一化頻率）（6-20式）
ws=[50 250]*2/Fs;   %阻帶邊界頻率（歸一化頻率）（6-20式）
Rp=1;Rs=30;Nn=128;  %通帶波紋和阻帶衰減以及繪制頻率特性的數據點數
[N,Wn]=cheb1ord(wp,ws,Rp,Rs);%求得數字濾波器的最小階數和歸一化截止頻率
[b,a]=cheby1(N,Rp,Wn);  %按最小階數、通帶波紋和截止頻率設計數字濾波器
figure(1)
[H,f]=freqz(b,a,Nn,Fs);   %求得濾波器的頻率特性
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=30;f2=100;f3=270;   %輸入信號的三種頻率成分
N=100;   %輸入信號的數據點數
dt=1/Fs;n=0:N-1;t=n*dt; %時間序列
x=sin(2*pi*f1*t)+0.3*cos(2*pi*f2*t)+0.1*sin(2*pi*f3*t); %輸入信號
subplot(2,1,1),plot(t,x),title('輸入信號')  %繪制輸入信號
y=filtfilt(b,a,x);   %對輸入信號進行濾波
subplot(2,1,2),plot(t,y)  %繪制輸出信號
ylim([-0.2 0.3])
title('輸出信號'),xlabel('時間/s')