拉普拉斯变换及其曲面图.docVIP

x1=-0.2:0.01:0.2; y1=-0.2:0.01:0.2; [x,y]=meshgrid(x1,y1); s=x+i*y; fs=abs(1./s); Warning: Divide by zero. (Type warning off MATLAB:divideByZero to suppress this warning.) mesh(x,y,fs); 拉普拉斯变换及其曲面图： x1=-0.2:0.01:0.2; y1=-0.2:0.01:0.2; [x,y]=meshgrid(x1,y1); s=x+i*y; fs=abs(1./s); mesh(x,y,fs); title(单位阶跃信号拉普拉斯变换曲面图); colormap(hsv); axis([-0.2,0.2,-0.2,0.2,0,60]); 绘制连续系统零、极点图： A=[1 3 4]; p=roots(A) p = -1.5000 + 1.3229i -1.5000 - 1.3229i B=[1,2,5]; q=roots(B) q = -1.0000 + 2.0000i -1.0000 - 2.0000i 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,实轴) 连续系统零、极点分析： p=p; q=q; f=1:2:3; 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,lenght(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))