信号与线性系统实验三：连续LTI时间系统的复频域分析

---

一、 实验目的

  1.通过该实验，掌握连续时间信号傅里叶变换和傅里叶逆变换的实现方法，以及傅立叶变换的时移特性、傅立叶变换的频移特性的实现方法。
  2.通过该实验，掌握连续时间信号拉普拉斯变换和拉普拉斯逆变换的实现方法及编程思想，掌握系统函数的零极点图的绘制方法。
  3.了解函数fourier() 和函数ifourier()的调用格式及作用。
  4.了解函数laplace()和函数ilaplace()的调用格式及作用。
  5.掌握傅里叶变换的数值计算方法，以及绘制信号频谱图的方法。
  6.了解连续系统零极点图的绘制及利用零极点图判断系统稳定性的方法。

二、实验内容与原理

第一部分：实验内容
（1）验证性实验
  A、连续系统的频域分析
  B、连续系统的复频域分析
（2）程序设计实验。
第二部分：实验原理
（1）函数调用格式
  A、函数fourier()和函数ifourier()的调用格式
  B、使用函数laplace()可以实现信号f(t)的拉普拉斯变换
  C、使用函数ilaplace()可以实现信号F(s)的拉普拉斯逆变换
（2）连续时间信号傅里叶变换的数值计算方法
（3）由系统函数的零极点图判断系统的稳定性

三、实验器材（设备、元器件、软件工具、平台）

  计算机，Matlab软件

四、实验步骤

第一部分：
验证性实验一：连续LTI系统频域分析
  1.傅里叶变换
  2.傅里叶逆变换
  3.傅里叶变换的数值计算
  4.傅里叶变换的时移特性
  5.傅立叶变换的频移特性
验证性实验二：连续LTI系统S域分析
  1.拉普拉斯变换
  2.拉普拉斯逆变换
  2.绘制连续系统零极点图
第二部分：程序设计实验
  题目在第五部分——实验数据及结果分析中，这里不在赘述。

五、实验数据及结果分析

1.已知某连续LTI系统的系统函数为：$H(s)=\frac{s^2+3s+2}{8s^4+2s^3+3s^2+s+5}$，试用Matlab编程确定系统的零极点，并绘出零极点分布图，同时判断系统的稳定性。
解：

>> a=[8,2,3,1,5];
b=[1,3,2];
lxljdt(a,b)
结果如下(p:函数返回的系统函数极点位置行向量;q:函数返回的系统函数零点位置行向量):
p =       
-0.6155 + 0.6674i
  -0.6155 - 0.6674i
   0.4905 + 0.7196i
   0.4905 - 0.7196i
q =    
    -2
    -1

图形分析：图中的‘X’代表极点。从图中可以看出，在s平面的右半平面存在极点，该系统不稳定。

2.通过编程绘出下列系统函数的零极点图。
$$H(s)=\frac{2s+1}{s^3+2s^2+2s+1}$$
解：

>> a=[1,2,2,1];
b=[2,1];
lxljdt(a,b);

图形分析：图中的‘X’代表极点。从图中可以看出，系统函数H(s)的极点均在s平面的左半开平面，该系统稳定。

3.分别确定下列信号的拉普拉斯变换，并绘制在s平面的三维曲面图。
（1）$f(t)=4e^{-5t}\epsilon (t)$
（2）$f(t)=3cos(5t)\epsilon (t)$
（3）$f(t)=8\delta (t-5)$
（4）$f(t)=3e^{-6t}cos(2t)\epsilon (t)$

(1)解：

>> syms t positive；             
f=4*exp(-5*t);
F=laplace(f)
F =
4/(s + 5)
 >> x1=-5.2:0.03:-4.8;
y1=-0.2:0.03:0.2;
[x,y]=meshgrid(x1,y1);
s=x+i*y;
Fs=abs(4./(s + 5));
mesh(x,y,Fs)
surf(x,y,Fs)
title('拉氏变换曲面图')

结果分析：$F(s)=4/(S+5)$,它的极点为-5。由$s=x+iy$,可知画图三维曲面时，x应该在-5左右取值，y应该在0左右取值。（上课抽查时，有的同学设定x和y时取到了极点，因为极点处是不存在值的，所以MATLAB中极点处会出现一个分叉，因此我x和y取值时有意跳过了极点，这样就不会出现分叉。）

(2)

 >> syms t positive;             
f=3*cos(5*t);
F=laplace(f)
F =
(3*s)/(s^2 + 25)
>> x1=-5:0.3:5;
y1=-10:0.3:10;
[x,y]=meshgrid(x1,y1);
s=x+i*y;
Fs=abs((3*s)./(s.^2 + 25));
mesh(x,y,Fs)
surf(x,y,Fs)
title('拉氏变换曲面图')

结果分析：$F(s)=3s/(s^2+25)$,它的极点为-5i和5i。可知画图三维曲面时，x和y的取值应该囊括所以的极点。由$s=x+iy$,所以x取值在0左右，y取值时最小值小于-5，最大值大于5。

(3)

>> syms t positive;             
f=8*dirac(t-5);
F=laplace(f) 
F = 
8*exp(-5*s)
>> x1=-0.2:0.03:0.2;
y1=-0.2:0.03:0.2;
[x,y]=meshgrid(x1,y1);
s=x+i*y;
Fs=abs(8*exp(-5*s));
mesh(x,y,Fs)
surf(x,y,Fs)
title('拉氏变换曲面图')

结果分析：$F(s)=8e^{5S}$,它不存在极点。画图三维曲面时，x和y的取值适当即可。

(4)

>> syms t positive;             
f=3*exp(-6*t)*cos(2*t);
F=laplace(f)
F =
(3*(s + 6))/((s + 6)^2 + 4)
>> x1=-10:0.3:0;
y1=-5:0.3:5;
[x,y]=meshgrid(x1,y1);
s=x+i*y;
Fs=abs((3*(s + 6))./((s + 6).^2 + 4));
mesh(x,y,Fs)
surf(x,y,Fs)
title('拉氏变换曲面图')

结果分析：$F(s)=\frac{3(s+6)}{(s+6{)}^2+4}$,它的极点为-6+2i和-6-2i。可知画图三维曲面时，x和y的取值应该囊括所以的极点。由s=x+iy,所以x取值在-6左右，y取值时最小值小于-2，最大值大于2。

4.已知LTI系统的激励函数为$f(t)=2e^{-3t}\epsilon (t)$，单位冲激响应$h(t)=t{e}^{-t/4}\epsilon (t)$，试采用变换域分析方法确定该系统的零状态响应。
解：

>> syms t f h;  
t=sym('t','positive');
f=2*exp(-3*t);
h=t*exp(-t/4);
F=laplace(f);
H=laplace(h);
ys=F*H;
yt=ilaplace(ys)
yt =
(32*exp(-3*t))/121 - (32*exp(-t/4))/121 + (8*t*exp(-t/4))/11

结果分析：这是利用了卷积定理以及系统的零状态响应等于该激励函数与单位冲激响应的卷积。先从S域求出Y(s)，然后求拉普拉斯逆变换。

5.分别确定下列信号的拉普拉斯逆变换。
$（1）F(s)=\frac{1}{(s+2)(s+4)}$
$（2）F(s)=\frac{s^2+s+2}{s^3+4s^2+s-6}$
解：（1）

>> syms s;
F=1/(s+2)/(s+4);
f=ilaplace(F)
f =
exp(-2*t)/2 - exp(-4*t)/2

结果分析：F(s)的极点为-2和-4。

（2）

>> syms s;
F=(s*s+s+2)/(1*s*s*s+4*s*s+1*s-6);
f=ilaplace(F)
f =
2*exp(-3*t) - (4*exp(-2*t))/3 + exp(t)/3

结果分析：我尝试修改分母，使其得出的结果可以看懂。我的想法就是先确定一个极点为1，然后修改了F(s)的分母不同阶数前的系数，进过大量实验，变成了$s^3+4s^2+s-6$。此时，经由拉普拉斯逆变换求得的式子明显可以看懂。故F(s)的极点为-3，-2和1。

6．试确定下列信号的傅里叶变换的数学表达式。
$（1）f(t)=\epsilon (t+1)-\epsilon (t-1)$
$（2）f(t)=e^{-3t}\epsilon (t)$
$（3）f(t)=e^{-t}\epsilon (t)$
$（4）f(t)=\delta ″(t)$
(1)

>> syms t;
f=heaviside(t+1)-heaviside(t-1); 
F=fourier(f)
F =
- (- sin(w) + cos(w)*1i)/w + (sin(w) + cos(w)*1i)/w

结果分析：上面的式子不是最简式，可以写成F=2sin(w)/w=2Sa(w)

(2)

 >> syms t;
f=exp(-3*t)*heaviside(t);
F=fourier(f) 
F =
1/(3 + w*1i) 

结果分析：上面的式子和我们平时的写法不太一样，可以写成F=1/(3+jw)

(3)

>> syms t;
f=exp(-t)*heaviside(t);
F=fourier(f) 
F =
1/(1 + w*1i) 

结果分析：上面的式子和我们平时的写法不太一样，可以写成F=1/(1+jw)

(4)

 >> syms t;
f=diff(dirac(t),t,2);
F=fourier(f)
F =
-w^2

结果分析：这是冲激函数的二阶微分的傅里叶变换,F=-w^2.

7. 试画出信号$f(t)=e^{-3t}\epsilon (t)$的频谱图，并画出信号f(t-4)以及信号$f(t)e^{-j4t}$的频谱图。
   解：

>> syms t;
f=exp(-3*t)*heaviside(t);
f1=exp(-3*(t-4))*heaviside(t-4);
F1=fourier(f)
F2=fourier(f1)
F3=fourier(f*exp(-4*t*i))
subplot(3,2,1)
ezplot(abs(F1),[-30 30]),title('f(t)的幅频曲线')
subplot(3,2,2)
ezplot(angle(F1),[-30 30]),title('f(t)的相频曲线')
subplot(3,2,3)
ezplot(abs(F2),[-30 30]),title('f(t-4)的幅频曲线')
subplot(3,2,4)
ezplot(angle(F2),[-30 30]),title('f(t-4)的相频曲线')
subplot(3,2,5)
ezplot(abs(F3),[-30 30]),title('f(t)exp(-4t)的幅频曲线')
subplot(3,2,6)
ezplot(angle(F3),[-30 30]),title('f(t)exp(-4t)的相频曲线')
F1 =
1/(3 + w*1i)
F2 =
exp(-w*4i)/(3 + w*1i)
F3 =
1/(w*1i + 3 + 4i)  

结果分析：上面是我进行的符号计算（我认为符号计算相对简单），不是数值计算。从图形中可以分析的得出：1. f(t-4)在时域向右移动了4个单位（发生了时移），但是从频谱图中可以看出，它的幅频响应和f(t)一样，但是相频响应有明显差异（没有频移，一般频谱图对应的是幅频响应，故是否频移看的是幅频响应）
2. $f(t)e^{-j4t}$虽然没有时移，但是它的相频响应和f(t)一样，但是幅频响应有移动（发生了频移）。

六、实验结论（整个实验总结性的结论性陈述）

  本次实验对连续LTI时间系统进行了频域、复频域（S域）的验证性实验和程序设计实验，相关的总结如下：
  1.我认为本次的实验的难点就在于熟练掌握函数的调用格式，例如：laplace (f) 是符号函数f的拉普拉斯变换；ilaplace()可以实现信号F(s)的拉普拉斯逆变换；fourier(f) 是符号函数f的傅里叶变换；laplace (f) 是符号函数f的拉普拉斯变换
  2.知道零极点图的绘制方法，并学会如何由系统的零极点图来判断系统的稳定性。对于连续LTI系统，若其系统函数 的极点均在 平面的左半开平面，则系统稳定；若在虚轴上只存在一阶极点，其余极点全部位于 平面的左半开平面，则系统临界稳定；若在虚轴上存在二阶及二阶以上极点或在 平面的右半平面存在极点，则系统不稳定。

七、其他（主要是心得体会）

  本次实验相比于前两次实验而言，比较简单。但是我也遇到了一些问题：
  1.做程序设计实中的第三题时，我用MATLAB做出来的三维曲面图是一片黑色，没有网格及颜色的变化。我实验了很多次，发现了原来是x和y定义的步长太小，导致用MATLAB做出来的三维曲面图中网格很小很小，因此网格可极限成黑线。在视觉效果上，由黑线构成的三维图看不出网格及颜色的变化。
  2. 做程序设计实中的第七题时（符号计算做的），我开始的时候没有规定频率的范围，导致MATLAB做出来的频谱图中只是系统默认的一部分，因此只能看出频谱的局部变化。后来我划分了频率的不同范围，找到了一个我认为合适的范围，可以观察频谱图的整体变化。