% ch3example4A.m
f_N=10000;                                        % 采樣率
f_p=[1000, 1500]; f_s=[600, 1900]; R_p=3; R_s=20; % 設計要求指標
Ws=f_s/(f_N/2); Wp=f_p/(f_N/2);                   % 計算歸一化頻率
[n, Wn]=ellipord(Wp,Ws,R_p,R_s);                  % 計算階數和截止頻率
[b,a]=ellip(n, R_p, R_s, Wn);                     % 計算H(z)
freqz(b,a, 1000, 10000) % 作出H(z)的幅頻相頻圖, freqz(b,a, 計算點數, 采樣率)
subplot(2,1,1); axis([0 5000 -30 3])