hc735.m

书籍中matlab第二部分代码

%《数字信号处理教程——MATLAB释义与实现》第七章例7.3.5程序hc735
% 用哈明窗设计FIR低通滤波器
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
wp = 0.2*pi; ws = 0.3*pi;  deltaw= ws - wp;		% 过渡带宽Δω的计算
N0 = ceil(6.6*pi/ deltaw) + 1;	% 按哈明窗(表7.3.2)计算所需的滤波器阶数N0
N=N0+mod(N0+1,2);		% 为了实现第一类偶对称滤波器，应使其长度N为奇数
wdham = (hamming(N))';	% 求窗函数
wc = (ws+wp)/2	;		% 截止频率取为两边界频率的平均值
hd = ideallp(wc,N);		% 求理想脉冲响应	
h = hd .* wdham;		% 设计的脉冲响应(即系数)为理想脉冲响应与窗函数乘积
[db,mag,pha,grd,w] = myfreqz(h,[1]);	% 对设计结果进行检验
dw = 2*pi/1000;
Rp = -(min(db(1:wp/dw+1)))        	% 检验通带波动
As = -round(max(db(ws/dw+1:501))) 	% 检验最小阻带衰减
%绘图
n=0:N-1;
subplot(2,2,1); stem(n,hd,'.'); title('理想脉冲响应')
axis([0 N-1 -0.1 0.3]); ylabel('hd(n)');text(N+1,-0.1,'n')
subplot(2,2,2); stem(n,wdham,'.');title('哈明窗')
axis([0 N-1 0 1.1]); ylabel('w(n)');text(N+1,0,'n')
subplot(2,2,3); stem(n,h,'.');title('实际脉冲响应')
axis([0 N-1 -0.1 0.3]); xlabel('n'); ylabel('h(n)')
subplot(2,2,4); plot(w/pi,db);title('幅度响应（单位： dB）');grid
axis([0 1 -100 10]); xlabel('频率（单位：pi）'); ylabel('分贝')
set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1])
set(gca,'YTickMode','manual','YTick',[-50,0])
set(gca,'YTickLabelMode','manual','YTickLabels',['50';' 0'])
set(gcf,'color','w');                % 置图形背景色为白