%=========================================================
%   程序名：jidiantu
%   繪制系統(tǒng)函數(shù)的拉普拉斯變換曲面圖
%=========================================================
clear;close all;
a=-6:0.2:6;
b=-6:0.2:6;
[a1,b1]=meshgrid(a,b);			
s=a1+i*b1;					  %確定繪制曲面圖的復(fù)平面區(qū)域
c=ones(length(b),length(a));
n=input('請輸入系統(tǒng)函數(shù)表達(dá)式中的分子多項式：');
m=input('請輸入系統(tǒng)函數(shù)表達(dá)式中的分母多項式：');
l=length(n);r=length(m);
p=n(l)*c;q=m(r)*c;
if l==1;
    p=n(l)*c;
else for i=1:l-1;
       p=p+n(i)*s.^(l-1);
     end
end
if r==1;
    q=m(r)*c;
else for i=1:r-1;
         q=q+m(i)*s.^(r-1);
     end
end
H=p./(q+eps);
H=abs(H);				      %計算拉普拉斯變換的樣值
%subplot(221)
mesh(a,b,H);                  %繪制曲面圖
surf(a,b,H);
axis([-6,6,-6,6,-0,3]);
title('拉普拉斯曲線變換圖');
colormap(hsv);
%view(-20,30);                 %確定三維視角
%==========================================================