发布时间 : 星期日 文章北京理工大学信号与系统实验报告4 LTI系统的频域分析更新完毕开始阅读967fcd45a32d7375a41780d4
实验4 LTI系统的频域分析
（综合型实验）
一、实验目的
1） 加深对LTI系统频率响应的基本概念的掌握和理解。 2） 学习和掌握LTI系统频率特性的分析方法。 二、实验原理与方法
1. 连续时间系统的频率响应
系统的频率响应定义为系统单位冲击响应h(t)的傅里叶变换，即
??H(?)????h(?)e?j??d? （1）
若LTI连续时间系统的单位冲激响应为h(t)，输入信号为x(t)，根据系统的时域分析可知系统的零状态响应为
y(t)?x(t)*h(t) （2）
对上式两端分别求傅里叶变换，由时域卷积定理可得
Y(?)?X(?)H(?) （3）
因此系统的频率响应还可以由系统的零状态响应和输入的傅里叶变换之比得到：
H(?)?Y(?)/X(?) （4）
H(?)反映了LTI连续时间系统对不同频率信号的响应特性，是系统内在的固有特性，
与外部激励无关。H(?)又可以表示成： H(?)?|H(?)|ej?(?) （5）
j?t其中|H(?)|成为系统的幅度响应，?(?)成为系统的相位响应。当虚指数信号eLTI系统时，系统的零状态响应y(t)仍然是同频率的虚指数信号，即
作用
y(t)?ej?tH(?) （6）
由此还可以推导出正弦信号作用在系统上的响应如下表所示：
输入信号 响应 sin(?0t??),???t??? cos(?0t??),???t??? 对于下述微分方程描述的LTI连续时间系统
|H(?0)|sin(?0t??(?0)??) |H(?0)|cos(?0t??(?0)??) ?aynn?0N(n)(t)??bmx(m)(t) （7）
m?0M
其频率响应H(j?)可表示为（8）式所示的j?的有理多项式。
Y(?)bM(j?)M?bM?1(j?)M?1?...?b1j??b0 （8） H(?)??NN?1X(?)aN(j?)?aN?1(j?)?...?a1j??a0MATLAB的信号处理工具箱提供了专门的函数freqs，用来分析连续时间系统的频率响应，该函数有下列几种调用格式：
[h,w]?freqs(b,a)计算默认频率范围内200个频率点上的频率响应的取样值，这200
个频率点记录在w中。
h?freqs(b,a,w) b、a分别为表示H(j?)的有理多项式中分子和分母多项式的系数向
量，w为频率取样点，返回值h就是频率响应在频率取样点上的数值向量。
[h,w]?freqs(b,a,n)计算默认频率范围内n个频率点上的频率响应的取样值，这n个
频率点记录在w中。freqs(b,a,...) 这种调用格式不返回频率响应的取样值，而是以对数坐标的方式绘出系统的幅频响应和相频响应。 2. 离散时间系统的频率响应
LTI离散时间系统的频率响应定义为单位抽样响应h(n)的离散时间傅里叶变换。
H(e)?j?n????h(n)e???j?n （9）
对于任意输入信号x(n)，输入与输出信号的离散时间傅里叶变换有如下关系
Y(ej?)?H(ej?)X(ej?) （10）
因此，系统的频率响应还可以表示为
H(ej?)?Y(ej?)/X(ej?) （11）
当系统输入信号为x(n)?ej?n时，系统的输出为
y(n)?ej?n*h(n)?k????e??j?(n?k)h(k)?ej?nH(ej?) （12）
由（12）式可知，虚指数信号通过LTI离散时间系统后信号的频率不变，信号的幅度由系统的频率响应的幅度值确定，所以H(e)表示了系统对不同频率信号的衰减量。 一般情况下离散系统的频率响应H(e)是复值函数，可用幅度和相位来表示。
j?j?H(ej?)?|H(ej?)|ej?(?) （13）
其中|H(e)|称为系统的幅度响应，?(?)称为系统的相位响应。
j?
若LTI离散系统可以由如下差分方程描述。
?ay(n?i)??bx(n?j) （14）
iji?0j?0NM则由（11）描述的离散时间系统的频率响应H(ej?)可以表示为ej?的有理多项式。
Y(ej?)b0?b1e?j??...?bMe?jM? （15） H(e)??X(ej?)a0?a1e?j??...?aNe?jN?j?MATLAB的信号处理工具箱提供了专门的函数freqz，用来分析连续时间系统的频率响应，该函数有下列几种调用格式：
[H,w]?freqz(b,a,n) b、a分别为有理多项式中分子和分母多项式的系数向量，返回
值H是频率响应在0~pi范围内n个频率等分点上的数值向量，w包含了这n个频率点。
[H,w]?freqz(b,a,n,'whole')计算0~2?n个频率点上的频率响应的取样值，这n个频
率点记录在w中。
H?freqz(b,a,w) w为频率取样点，计算这些频率点上的频率响应的取样值。 freqz(b,a,...)这种调用格式不返回频率响应的取样值，而是直接绘出系统的幅频响应和
相频响应。 三、实验内容
（1）已知一个RLC电路构造的二阶高通滤波器如下图所示，其中R?L， 2C L=0.4H，C=0.05F
1） 计算该电路系统的频率响应及高通截止频率；
Y(?)?X(?)1111??11j?CRj?L?Rj?L经整理
?1
(j?)2H(?)?Y(?)/X(?)?(j?)2?10j??50 ，使H(?)?0.707得高通截止频率??7.07
2）利用MATLAB绘制幅度响应和相位响应曲线，比较系统的频率特性与理论计算的结果是否一致。
b=[1 0 0]; a=[1 10 50];
[H,w]=freqs(b,a); subplot(211); plot(w,abs(H));
xlabel('\\omega(rad/s)'); ylabel('Magnitude');
set(gca,'ytick',[0 0.4 0.707 1]); title('|H(j\\omega)|');
|H(j?)|1grid on; subplot(212); plot(w,angle(H)); xlabel('\\omega(rad/s)'); ylabel('Phase'); title('\\phi(\\omega)'); grid on
Magnitude0.7070.400102030405060708090100?(rad/s)?(?)43Phase2100102030405060708090100?(rad/s)
set(gca,'xtick',[0:10]); axis([0 10 0 1]);
set(gca,'ytick',[0 0.4 0.707 1]); title('|H(j\\omega)|'); grid on;
为更准确观察，对以上程序略作修改
b=[1 0 0]; a=[1 10 50]; [H,w]=freqs(b,a); plot(w,abs(H));
xlabel('\\omega(rad/s)'); ylabel('Magnitude');
|H(j?)|1由左图可知，H(?)?0.707时，频率约
为7.07，实验结果与理论值近似相符。
0.707Magnitude0.40012345678910?(rad/s)
（2）已知一个RC电路如下图所示。