hc733.m

matlab数字信号处理代码

%《数字信号处理教程——MATLAB释义与实现》第七章例7.3.3程序hc733
% 三角窗累积符幅特性
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
clear
N=51; n=1:N; wd=triang(N);		% 给出三角窗参数，产生窗函数
[A,omega]=amplres(wd');			% 调用amplres，求符幅函数，变wd为行向量
Ad=[fliplr(A),A(2:end)];		% 将符幅函数对称扩展至负频率区
omegad=[-fliplr(omega),omega(2:end)]; % 将频率对称扩展至负频率区
Gd=cumtrapz(Ad)*(omega(2)-omega(1))/pi; % 积分求累积符幅函数Gd
AdB=20*log10(abs(Ad));			% 将符幅函数转换为分贝
GdB=20*log10(Gd); 				% 将累积符幅函数转换为分贝
subplot(2,2,1),plot(omegad,Ad)  % 绘图
subplot(2,2,2),plot(omegad,Gd)
subplot(2,2,3),plot(omegad,AdB)
axis([-pi,pi,-50,50]),grid on
subplot(2,2,4),plot(omegad,GdB)       % 绘分贝幅特性
axis([-pi,pi,-50,10]),grid on
set(gcf,'color','w');                % 置图形背景色为白