ep731b.m

来自「高西全2009年最新版<数字信号处理>源码程序」· M 代码 · 共 21 行

datB=pi/16;wc=pi/3;   %过渡带宽度pi/16,通带截止频率为pi/3；
m=1;
N=(m+1)*2*pi/datB+1     %按（7.3.15）式估算采样点数N
F=[0,wc/pi,wc/pi+2/N,wc/pi+4/N,1];
A=[1,1,0.38,0,0];
hn=fir2(N-1,F,A,boxcar(N));
Hw=fft(hn,1024);    %计算频率响应函数:DFT[h(n)]
wk=2*pi*[0:1023]/1024;
Hgw=Hw.*exp(j*wk*(N-1)/2);  %计算幅度响应函数Hg(ω)
%=====================================================
%计算通带最大衰减Rp和阻带最小衰减Rs
Rp=max(20*log10(abs(Hgw)))
hgmin=min(real(Hgw));Rs=20*log10(abs(hgmin))
%==============================================================================================
%以下绘图7.3.4
subplot(3,2,1);plot(F,A);axis([0,1,0,1.3]);grid on
xlabel('ω/π');ylabel('H_d_g(ω)');
n=0:N-1;
subplot(3,2,3);stem(n,hn,'.');axis([0,65,-0.2,0.5]);xlabel('n');ylabel('h(n)');
subplot(3,2,5);plot(wk/pi,20*log10(abs(Hgw)));axis([0,1,-60,3]);grid on
xlabel('ω/π');ylabel('20lg|Hg(ω)|');