[a,b]=meshgrid(a,b); s=a+i*b; s=1./((s+2).*(s+2)+1); Fs=abs(s); mesh(a,b,Fs); suf(a,b,Fs); title('单边正弦信号f (t) = e-2t sin( t)u(t)拉氏变换曲面图'); colormap(hsv); 2520151050210-1信号拉氏变换曲面图-1-1.5-2-2.5-2-3 2、已知连续时间信号f(t)?e?3tu(t)，试求出该信号的拉普拉斯变换F(s)和傅立叶变换F(j?)，用MATLAB绘出拉普拉斯变换曲面图F(s)及幅频曲线F(j?，观察曲面图在虚轴上的剖面图，并将其与幅频曲线相比较，分析频域与复频域的对应关系。 (1)实验程序： a=-0:0.1:5; b=-20:0.1:20; [a,b]=meshgrid(a,b); s=a+i*b; %确定绘图区域 s=1./(s+3); Fs=abs(s); %计算拉普拉斯变换 subplot(211); mesh(a,b,Fs); %绘制曲面图 surf(a,b,Fs); view(-60,20) %调整观察视角 axis([-0,5,-20,20,0,0.5]); title('f (t) = e-3tu(t)拉普拉斯变换（S 域像函数）'); colormap(hsv); w=-20:0.1:20; %确定频率范围 Fw=1./(i.*w+3); %计算傅里叶变换 subplot(212); plot(w,abs(Fw)) %绘制信号振幅频 实验结果： 谱曲线 title('f (t) = e-3tu(t)傅里叶变换（振幅频谱曲线）') xlabel('频率w') 3、已知信号的拉普拉斯变换如下所示，试用MATLAB绘制曲面图，观察拉普拉斯变换零极点分布对曲面图的影响。 实验程序： function [p,q]=sjdt(A,B) %绘制连续系统零极点图程序 %A:系统函数分母多项式系数向量 %B:系统函数分子多项式系数向量 %p:函数返回的系统函数极点位置行向量 %q:函数返回的系统函数零点位置行向量 p=roots(A); %求系统极点 q=roots(B);; %求系统零点 实验结果： 拉普拉斯变换曲面图76543210420-2-4-5-4-3-2-10123p=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,'虚轴'); 连续系统零极点图43210-1-2-3-4-4-3-2-101234实轴虚轴