程序exn914(2)-画波德图 %下面的语句用来绘制波德图,如果用bode函数,只要一句 figure(2),bode(prod(abs(po)),poly(po));hold on bh=by;ah=poly(p0); %求H(S)的分子分母系数 w=logspace(2,4);%给出频率范围和分度 H polyval(bh,j*w)./polyval(ah,j*w);%H(jw) aH=unwrap(angle(H))*180/pi; %求出相角 fH=20*log10(abs(H)); %求出振幅 figure(2), subplot(2,1,1),semilogx(w,fH),grid on,hold on%绘幅频图 subplot(2,1,2),semilogx(w,aH),grid on,hold on%相频图 end
程序exn914(2)-画波德图 %下面的语句用来绘制波德图,如果用bode函数,只要一句 % figure(2),bode(prod(abs(p0)), poly(p0));hold on bh=by;ah= poly(p0); % 求H(s)的分子分母系数 w=logspace(2,4); % 给出频率范围和分度 H = polyval(bh,j*w)./polyval(ah,j*w); % 求H(jw) aH=unwrap(angle(H))*180/pi; % 求出相角 fH=20*log10(abs(H)); % 求出振幅 figure(2), subplot(2,1,1),semilogx(w,fH),grid on,hold on % 绘幅频图 subplot(2,1,2),semilogx(w,aH),grid on,hold on% 相频图 end
程序en914运行结果(时域) ·运行此程序,设N=4, 可得过渡过程如右图, 0.8 Ln= 3 从中看出输出信号达 到0.6处所需的时间约 0.6 为单级时常数乘以级 0.4 数。由于极点互相接 近,此程序在N>4时又 0.2 会出现很大误差。 0 0.005 0.01
程序exn914运行结果(时域) • 运行此程序,设N=4, 可得过渡过程如右图, 从中看出输出信号达 到0.6处所需的时间约 为单级时常数乘以级 数。由于极点互相接 近,此程序在N>4时又 会出现很大误差。 0 0.005 0.01 0 0.2 0.4 0.6 0.8 1 n=1 2 3 4