q724.m

matlab在电子信息中的应用和实例中的源代码

%《MATLAB及在电子信息课程中的应用》第七章例7.24程序q724
% 用各种窗函数设计FIR数字滤波器q724
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
clear;close all
N=21; wc=pi/4;						    % 理想低通滤波器参数
n=0:N-1; r=(N-1)/2;
hdn=sin(wc*(n-r))/pi./(n-r); 	        % 计算理想低通单位脉冲响应hd(n)
if rem(N,2)~=0   hdn(r+1)=wc/pi; end	% N为奇数时，处理n=r点的0/0型
wn1=boxcar(N); 		                    % 矩形窗
hn1=hdn.*wn1'; 		                    % 加窗
% 以上两条语句可代以fir函数：hn1=fir1(N-1,wc/pi,boxcar(N));
wn2=hamming(N);		                    % hamming 窗
hn2=hdn.*wn2'; 		                    % 加窗
% 以上两条语句可代以fir函数：hn2=fir1(N-1,wc/pi,hamming(N));k=3;
k=1;                                    % 绘图函数p751f的位置参数
hnwplot(hn1,k);title('矩形窗设计的h(n)')
set(gcf,'color','w')                    % 设置图形背景色为白色
%hamming 窗
wn2=hamming(N);
hn2=hdn.*wn2';
k=3;
hnwplot(hn2,k);title('hamming设计的h(n)')