% 第八章: 例 8.30
%            切比雪夫-2帶阻濾波器設計
%                 CHEBY2函數的應用
%
% 數字低通濾波器指標:            % 切比雪夫-2型帶阻:
ws = [0.4*pi 0.7*pi];                   % 數字阻帶邊緣頻率
wp = [0.25*pi 0.8*pi];                  % 數字通帶邊緣頻率
Rp = 1;	                            % 通帶波動(dB) 
As = 40;                             % 阻帶衰減(dB) 
Ripple = 10 ^ (-Rp/20);                 % 通帶波動
Attn = 10 ^ (-As/20);                   % 阻帶衰減

% 計算切比雪夫-2型濾波器參數
 [N,wn] = cheb2ord(wp/pi,ws/pi,Rp,As);

% 數字切比雪夫-2型帶阻濾波器設計:
[b,a] = cheby2(N,As,ws/pi,'stop');

% 級聯型式實現:
[b0,B,A] = dir2cas(b,a)
%%b0 = 0.1558
%%B = 1.0000    1.1456    1.0000
%%    1.0000    0.8879    1.0000
%%    1.0000    0.3511    1.0000
%%    1.0000   -0.2434    1.0000
%%    1.0000   -0.5768    1.0000
%%A = 1.0000    1.3041    0.8031
%%    1.0000    0.8901    0.4614
%%    1.0000    0.2132    0.2145
%%    1.0000   -0.4713    0.3916
%%    1.0000   -0.8936    0.7602

% 繪圖:
[db,mag,pha,grd,w] = freqz_m(b,a);
subplot(2,2,1);plot(w/pi,mag);grid;title('頻率響應')
xlabel(''); axis([0 1 0 1])
set(gca,'XTickMode','manual','XTick',[0;0.25;0.4;0.7;0.8;1])
set(gca,'YTickMode','manual','YTick',[0;Ripple;1])
subplot(2,2,3);plot(w/pi,db);grid;title('幅度(dB)')
xlabel('數字頻率（單位：pi）'); axis([0 1 -50 0]);
set(gca,'XTickMode','manual','XTick',[0;0.25;0.4;0.7;0.8;1])
set(gca,'YTickMode','manual','YTick',[-40;0])
subplot(2,2,2);plot(w/pi,pha/pi);grid;title('相位響應')
xlabel('');ylabel('相位（單位：pi）') 
set(gca,'XTickMode','manual','XTick',[0;0.25;0.4;0.7;0.8;1])
subplot(2,2,4);plot(w/pi,grd);grid;title('群延遲')
xlabel('數字頻率（單位：pi）'); 
set(gca,'XTickMode','manual','XTick',[0;0.25;0.4;0.7;0.8;1])
subplot;