MATLAB参考于程序.doc

附录 MATLAB参考程序 1. 绘制连续系统零极点图程序 function [p,q]=sjdt(A,B) %A:系统函数分母多项式系数向量 %B:系统函数分子多项式系数向量 %p:函数返回的系统函数极点位置行向量 %q:函数返回的系统函数零点位置行向量 p=roots(A); 　%求系统极点 q=roots(B);; 　%求系统零点 p=p; %将极点列向量转置为行向量 q=q; ％将零点列向量转置为行向量 x=max(abs([p q])); %确定纵坐标范围 x=x+0.1; y=x; %确定横坐标范围 clf hold on axis([-x x -y y]); %确定坐标轴显示范围 axis(square) plot([-x x],[0 0]) %画横坐标轴 plot([0 0],[-y y]) %画纵坐标轴 plot(real(p),imag(p),x) %画极点 plot(real(q),imag(q),o) %画零点 title(连续系统零极点图) %标注标题 text(0.2,x-0.2,虚轴) text(y-0.2,0.2,实轴) 2.根据连续系统零极点分布绘制系统频率响应曲线程序 function splxy(f1,f2,k,p,q) %f1、f2：绘制频率响应曲线的频率范围（即频率起始和终止点，单位为赫兹） %p、q：系统函数极点和零点位置行向量 %k：绘制频率响应曲线的频率取样间隔 p=p; q=q; f=f1:k:f2; %定义绘制系统频率响应曲线的频率范围 w=f*(2*pi); y=i*w; n=length(p); m=length(q); if n==0 %如果系统无极点 yq=ones(m,1)*y; vq=yq-q*ones(1,length(w)); bj=abs(vq); ai=1; elseif m==0 %如果系统无零点 yp=ones(n,1)*y; vp=yp-p*ones(1,length(w)); ai=abs(vp); bj=1; else yp=ones(n,1)*y; yq=ones(m,1)*y; vp=yp-p*ones(1,length(w)); vq=yq-q*ones(1,length(w)); ai=abs(vp); bj=abs(vq); end Hw=prod(bj,1)./prod(ai,1); plot(f,Hw); title(连续系统幅频响应曲线) xlabel(频率w（单位：赫兹）) ylabel(F(jw)) 3. 绘制离散系统零极点图程序 function ljdt(A,B) p=roots(A)； %求系统极点 q=roots(B)； %求系统零点 p=p； %将极点列向量转置为行向量 q=q； %将零点列向量转置为行向量 x=max(abs([p q 1]))； %确定纵坐标范围 x=x+0.1； y=x； %确定横坐标范围 clf hold on axis([-x x -y y]) %确定坐标轴显示范围 w=0:pi/300:2*pi； t=exp(i*w)； plot(t) %画单位园 axis(square) plot([-x x],[0 0]) %画横坐标轴 plot([0 0],[-y y]) %画纵坐标轴 text(0.1,x,jIm[z]) text(y,1/10,Re[z]) plot(real(p),imag(p),x) %画极点 plot(real(q),imag(q),o) %画零点 title(pole-zero diagram for discrete system) %标注标题 hold off 4.根据离散系统零极点分布绘制系统频率响应曲线程序 fun