fgp731.m

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

%《数字信号处理教程——MATLAB释义与实现》第七章演示程序fgp731
% 理想脉冲序列与矩形窗的时域相乘和频域卷积
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
clear
wc=1;N=15;alpha=(N-1)/2;                        % 给出理想频率特性参数
ns=-19;nf=29;n=[ns:nf]+1e-10;                   % 给出序列长度区间
hd=sin(n-alpha)./(n-alpha)/pi;                  % 理想脉冲序列
subplot(3,2,1),stem(n,hd,'.')                   % 画脉冲序列图
win=[zeros(1,0-ns),ones(1,N),zeros(1,nf-N+1)];  % 给出矩形窗函数序列
subplot(3,2,3),stem(n,win,'.')                  % 画矩形窗图        
w=[-1:0.01:1]*pi+1e-10;                         % 给出频率自变量向量
h=seqmult(hd,n,win,n);                   % 理想脉冲序列与矩形窗相乘，得到实际的脉冲响应
Win=sin(0.5*N*w)./(0.5*N*w);                    % 矩形窗函数的频谱
N1=round(1/0.01/pi);
Hd=[zeros(1,100-N1),ones(1,2*N1+1),zeros(1,100-N1)];    % 理想滤波器频谱
H=conv(Hd,Win);                          % 理想滤波器频谱与矩形窗频谱卷积，得到实际频谱
subplot(3,2,5),stem(n,h,'.')                    % 绘实际的脉冲响应图
subplot(3,2,2),plot([-pi,-1,-1,1,1,pi],[0,0,1,1,0,0]) % 绘理想滤波器频谱图
subplot(3,2,4),plot(w,Win);                     % 绘矩形窗函数的频谱图
subplot(3,2,6),plot(w,H(floor(0.5*length(w):1.5*length(w)-1)));  % 绘实际滤波器频谱图
set(gcf,'color','w')                        % 置图形背景色为白