hc742.m

这是数字信号处理——matlab释义与实现一书的源代码

%《数字信号处理教程——MATLAB释义与实现》第七章例7.4.2程序hc742
% 频率样本法过渡带所取样本数的影响
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
N=65;N1=fix(N/6);N2=N-2*N1-1;
A=[ones(1,N1+1),zeros(1,N2),ones(1,N1)];	% 符幅特性样本序列
A1=[ones(1,N1+1),0.5,zeros(1,N2-2),0.5,ones(1,N1)];	% 符幅特性样本序列
T1=0.5925,T2=0.1099,
A2=[ones(1,N1),T1,T2,zeros(1,N2-2),T2,T1,ones(1,N1-1)];	% 符幅特性样本序列
theta = - pi*[0:N-1]*(N-1)/N;			% 相位特性样本序列
H=A.*exp(j*theta);				% 频率特性样本序列
h=real(ifft(H));				% 反变换求出脉冲序列，去掉运算误差造成的虚部
[db,mag,pha,grd,w] = myfreqz(h,[1]);    % 检验设计出的滤波器的频率响应
H1=A1.*exp(j*theta);				% 频率特性样本序列
h1=real(ifft(H1));				% 反变换求出脉冲序列，去掉运算误差造成的虚部
[db1,mag1,pha,grd,w] = myfreqz(h1,[1]);    % 检验设计出的滤波器的频率响应
H2=A2.*exp(j*theta);				% 频率特性样本序列
h2=real(ifft(H2));				% 反变换求出脉冲序列，去掉运算误差造成的虚部
[db2,mag2,pha,grd,w] = myfreqz(h2,[1]);    % 检验设计出的滤波器的频率响应
figure(2),plot(w,mag,'-.',w,mag1,':',w,mag2,'-')
legend('过渡带无样本','过渡带有一样本','过渡带有两样本')
set(gcf,'color','w');                % 置图形背景色为白