连续系统的复频域分析.pdfVIP

连续系统的复频域分析 一、拉普拉斯变换及其曲面图 用MATLAB 绘制拉普拉斯变换的曲面图 拉普拉斯变换是分析连续时间信号的有效手段，对于当 t 时信号幅度不衰减或增长 的时间信号，其傅里叶变换不存在，但我们可以用拉普拉斯变换来分析它们。 连续时间信号 f (t ) 的拉普拉斯变换定义为： 其中 s σ jw ，若以 s 为横坐标（实轴）， jw 为纵坐标（虚轴），复变量 s 就 构成了一个复平面，称为 s 平面。 显然， F(s) 是复变量 s 的复函数，为了便于理解和分析 F(s) 随s 的变化规律，我们可以 将 F(s) 写成 其中 | F(s)| 为复信号 F(s) 的模，而 j (s) 为 F(s) 的相角。 从三维几何空间的角度来看， 和 φ 对应着复平面上的两个曲面，如果我们能绘出 F(s) (s) 它们的三维曲面图，我们就可以直观地分析连续信号的拉普拉斯变换 F(s) 随复变量 s的变化。 上述过程，我们可以利用 MATLAB 的三维绘图功能来实现。现在考虑如何用 MATLAB 来绘制 s 平面的有限区域上连续时间信号 f (t ) 的拉普拉斯变换 F(s) 的曲面图，我们以单位 阶跃信号 e (t ) 为例来说明实现过程。 我们知道，对单位阶跃信号 f (t ) e (t ) ，其拉普拉斯变换为 。 首先，我们用两个向量来确定绘制曲面图的 s 平面的横、纵坐标的范围。例如，我们可定 义绘制曲面图的横坐标范围向量 x1和纵坐标范围向量 y1分别为： x1=::; y1=::; 然后再调用前面介绍过的 meshgrid() 函数来产生矩阵 S ，并用该矩阵来表示绘制曲面图的复 平面区域，对应的 MATLAB 命令如下： [x,y]=meshgrid(x1,y1); s=x+i*y; 上述命令产生的矩阵 s 包含了复平面 σ 、 jw 范围内以间隔 取样的所有样点。 最后我们再计算出信号拉普拉斯变换在复平面的这些样点上的值，即可用函数 mesh绘出其 曲面图，对应命令为： fs=abs(1./s); %计算拉氏变换在复平面上的样点值 mesh(x,y,fs); %绘制的氏变换曲面图 surf(x,y,fs); title( 单位阶跃信号拉氏变换曲面图 ); colormap(hsv); axis([,,,,0,60]); rotate3d; 执行上述命令后，绘制的单位阶跃信号拉普拉斯曲变换面图如图所示 例： 已知连续时间信号 f (t ) sin( t )e (t ) ，求出该信号的拉普拉斯变换，并用 MATLAB 绘制拉普拉斯变换的曲面图。 解：该信号的拉普拉斯变换为 我们可以用上面介绍的方法来绘出单边正弦信号的拉普拉斯变换的曲面图，实现这一 过程的程序如下： %绘制单边正弦信号拉普拉斯变换曲面图程序 clf; a=::; b=::; [a,b]=meshgrid(a,b); d=ones(size(a)); c=a+i*b; %确定绘制曲面图的复平面区域 c=c.*c; c=c+d; c=1./c; c=abs(c); %计算拉普拉斯变换的样值 mesh(a,b,c); %绘制曲面图 surf(a,b,c); axis([,,-2,2,0,15]); title( 单边正弦信号拉氏变换曲面图 ); colormap(hsv); 上述程序绘制的曲面图如图所示。 二、 由拉普拉斯曲面图观察频域与复频域的关系 我们