fig734.m

高西全2009年最新版<数字信号处理>源码程序

T=input('T= ')  %键入过渡采样值T
datB=pi/16;wc=pi/3;   %过渡带宽度pi/16,通带截止频率为pi/3；
m=1;
N=(m+1)*2*pi/datB+1     %按（7.3.15）式估算采样点数N
%N=N+mod(N+1,2);          %确保h(n)长度N为奇数
F=[0,wc/pi,wc/pi+2/N,wc/pi+4/N,1];
A=[1,1,T,0,0];
subplot(3,2,1);plot(F,A);
axis([0,1,0,1.2]);grid on
xlabel('ω/π');ylabel('H_d_g(ω)');