实验五连续系统频域分析.doc

实验五 连续系统频域分析 1 实验目的 (1) 掌握系统频率响应特性的概念及其物理意义。 (2) 掌握利用MATLAB计算系统频率响应特性及曲线绘制方法，理解具有不同频率响应特性的滤波器对信号的滤波作用。 (3) 掌握利用MATLAB进行系统频响特性分析的方法。 2 实验原理及方法 2.1 连续时间LTI系统的频率响应 所谓频率特性，也称为频率响应特性，简称频率响应，是指系统在正弦信号激励下的稳态响应随频率变化的情况，包括响应的幅度随频率的变化情况和响应的相位随频率的变化情况两个方面。 图5-1 LTI系统的时域及频域分析图 图5-1中、分别为系统的时域激励信号和响应信号，是系统的单位冲激响应，它们三者之间的关系为：，由傅里叶变换的时域卷积定理可得到： 3.1 或者： 5-1 为系统的频域数学模型，它实际上就是系统的单位冲激响应h(t)的傅里叶变换。即： 5-2 由于实际上是系统单位冲激响应的傅里叶变换，如果是收敛的，或者说是绝对可积的话，那么一定存在，而且通常是复数，因此，也可以表示成复数的不同表达形式。在研究系统的频率响应时，更多的是把它表示成极坐标形式： 5-3 上式中，称为幅度频率响应，反映信号经过系统后，信号各频率分量的幅度发生变化的情况；称为相位频率响应，反映信号经过系统后，信号各频率分量在相位上发生变化的情况。和都是频率(的函数。 对于一个系统，其频率响应为，其幅度响应和相位响应分别为和，如果作用于系统的信号为，则其响应信号为： 5-4 若输入信号为正弦信号，即x(t)=sin((0t)，则系统响应为： 5-5 可见系统对某一频率分量的影响表现为两个方面，一是信号幅度要被加权，二是信号相位要被移相。 由于和都是频率(的函数，所以系统对不同的频率分量造成的幅度和相位上的影响是不同的。 2.2 用MATLAB计算系统频率响应 MATLAB提供了专门对连续系统频率响应进行分析的函数freqs()。该函数可以求出系统频率响应的数值解，并可绘出系统的幅频及相频响应曲线。freqs()函数有如下三种调用格式： (1) [H,w] = freqs(b,a)：b，a分别为连续时间LTI系统的微分方程右边的和左边的系数向量w中。 (2) H = freqs(b,a,w) 在指定的频率范围内计算系统的频率响应特性。在使用这种形式的函数时，要在前面先指定频率变量w的范围。 例如在语句H = freqs(b,a,w)之前加上语句： w = 0:2*pi/256:2*pi; (3) [H,w] = freqs(b,a,n) 调用该格式将计算默认频率范围内n个频率点上系统频率响应的样值，并赋给返回变量H，n个频率点存放在w中。 下面举例说明调用freqs()函数求系统频率响应的方法。 例5-1：理想低通滤波器在物理上是不可实现的，但传输特性近似于理想特性的电路却能找到。图5-2是常见的用RLC元件构成的二阶低通滤波器(一般说来，阶数越高，实际滤波器的特性越接近于理想特性)。该电路的频率响应为： 设，绘出该电路频率响应。 图5-2 RLC二阶低通滤波器电路 解：令截止频率。 将R、L、C值代入的表达式，得： 其中： 当时，。 实现该系统响应的MATLAB程序如下： % Program5_1 b = [1];a = [0.08 0.4 1]; [H,w] = freqs(b,a,100); Hm = abs(H);phi = angle(H); subplot(211),plot(w,Hm),grid on title(幅度响应),xlabel( w ),ylabel(幅度) subplot(212),plot(w,phi),grid on title(相位响应),xlabel( w ),ylabel(相位(弧度) ) 运行结果如图5-3所示： 图5-3 RLC二阶低通滤波器的幅度及相位响应 从图5-3可见，当从0增大时，该低通滤波器的幅度从1降到0，约为3.5，而从0o降到-180o，与理论分析的结果一致。 例的极点位于左半平面，零点位于右半平面，且零点与极点对于轴与互为镜像对称网络。它可保证不影响传输信号的幅频特性，只改变信号的相频特性。一全通网络的系统函数为 试画出和。% Program5_2 b=[-2 10];a=[2 10]; w=0:0.01:100; H=freqs(b,a,w);