hc731.m

书籍中matlab第二部分代码

%《数字信号处理教程——MATLAB释义与实现》第七章例7.3.1程序hc731
% 不同长度矩形窗累积符幅特性
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
N0=[10,100];
for k=1:2
    N=N0(k);dw=0.001;               % 设置初值        
    w=[-pi/2:dw:pi/2]+eps;    % 设置自变量数组，偏置微量 
    AR=sin(N*w)./sin(w);         % 求矩形窗符幅函数
    H=cumtrapz(AR)*dw/pi;        % 求累积符幅函数
    subplot(2,2,k),plot(w,AR),grid on  % 画符幅函数曲线
    title('矩形窗符幅特性AR(\omega)'),xlabel('\omega')
    subplot(2,2,k+2),plot(w,H),grid on  % 画累积符幅函数曲线
    title('累积符幅特性H(\omega)'),xlabel('\omega')
    delta1=max(H)-1,    % 求通带最大波动（绝对值和分贝）
    delta2=abs(min(H)),              % 求阻带衰减（绝对值）
    Rp=-20*log10((1-delta1)./(1+delta1))	% 	由delta1及delta2求Rp
    As=-20*log10(delta2./(1+delta1))	    %	由delta1及delta2求As
    deltaw=dw*(find(H==max(H))-find(H==min(H)))    % 求过渡带宽度  
end
set(gcf,'color','w');                % 置图形背景色为白