% 第七章: 例 7.16
% 頻率采樣技術(shù): 低通, 最優(yōu)法T1 & T2
%            wp = 0.2pi, ws=0.3pi, Rp=0.25dB, As=50dB
%            T1 = 0.5925, T2=0.1099
M = 60; alpha = (M-1)/2; l = 0:M-1; wl = (2*pi/M)*l;
Hrs = [ones(1,7),0.5925,0.11,zeros(1,43),0.11,0.5925,ones(1,6)];
Hdr = [1,1,0,0]; wdl = [0,0.2,0.3,1];
k1 = 0:floor((M-1)/2); k2 = floor((M-1)/2)+1:M-1;
angH = [-alpha*(2*pi)/M*k1, alpha*(2*pi)/M*(M-k2)];
H = Hrs.*exp(j*angH);
h = real(ifft(H,M));
[db,mag,pha,grd,w] = freqz_m(h,1);
[Hr,ww,a,L] = Hr_Type2(h);
subplot(1,1,1)
subplot(2,2,1);plot(wl(1:31)/pi,Hrs(1:31),'o',wdl,Hdr); 
axis([0,1,-0.1,1.1]); title('低通: M=60,T1=0.59, T2=0.109')
xlabel(''); ylabel('Hr(k)')
set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1])
set(gca,'YTickMode','manual','YTick',[0,0.59,0.109,1]); grid
subplot(2,2,2); stem(l,h); axis([-1,M,-0.1,0.3])
title('脈沖響應(yīng)'); ylabel('h(n)');text(M+1,-0.1,'n')
subplot(2,2,3); plot(ww/pi,Hr,wl(1:31)/pi,Hrs(1:31),'o');
axis([0,1,-0.1,1.1]); title('振幅響應(yīng)')
xlabel('頻率（單位：pi）'); ylabel('Hr(w)')
set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1])
set(gca,'YTickMode','manual','YTick',[0,0.59,0.109,1]); grid
subplot(2,2,4);plot(w/pi,db); axis([0,1,-100,10]); grid
title('幅度響應(yīng)'); xlabel('頻率（單位：pi）');
ylabel('分貝數(shù)');
set(gca,'XTickMode','Manual','XTick',[0;0.2;0.3;1]);
set(gca,'YTickMode','Manual','YTick',[-63;0]);
set(gca,'YTickLabelMode','manual','YTickLabels',['63';' 0'])
%