jidiantu.m

经典《信号与系统》教程的matlab例程,对深入理解信号与系统相关概念有很大帮助

%=========================================================
%   程序名：jidiantu
%   绘制系统函数的拉普拉斯变换曲面图
%=========================================================
clear;close all;
a=-6:0.2:6;
b=-6:0.2:6;
[a1,b1]=meshgrid(a,b);			
s=a1+i*b1;					  %确定绘制曲面图的复平面区域
c=ones(length(b),length(a));
n=input('请输入系统函数表达式中的分子多项式：');
m=input('请输入系统函数表达式中的分母多项式：');
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);                 %确定三维视角
%==========================================================