% ch3example3A.m
f_N=8000;                            % 采樣率
f_p=1000; f_s=700; R_p=3; R_s=20;    % 設計要求指標
Ws=f_s/(f_N/2); Wp=f_p/(f_N/2);      % 計算歸一化頻率
[n, Wn]=cheb1ord(Wp,Ws,R_p,R_s);     % 計算階數和截止頻率
[b,a]=cheby1(n, R_p, Wn, 'high');    % 計算H(z)
freqz(b,a, 1000, 8000) % 作出H(z)的幅頻相頻圖, freqz(b,a, 計算點數, 采樣率)
subplot(2,1,1); axis([0 4000 -30 3])