拉普拉斯变换及其曲面图

>> 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)')
>>。