(综合型实验)
一、实验目的
1) 加深对LTI系统频率响应的基本概念的掌握和理解。 2) 学习和掌握LTI系统频率特性的分析方法。 二、实验原理与方法
1. 连续时间系统的频率响应
系统的频率响应定义为系统单位冲击响应h(t)的傅里叶变换,即
H()h()ejd (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)
jt其中|H()|成为系统的幅度响应,()成为系统的相位响应。当虚指数信号eLTI系统时,系统的零状态响应y(t)仍然是同频率的虚指数信号,即
作用
y(t)ejtH() (6)
由此还可以推导出正弦信号作用在系统上的响应如下表所示:
输入信号 响应 sin(0t),t cos(0t),t 对于下述微分方程描述的LTI连续时间系统
|H(0)|sin(0t(0)) |H(0)|cos(0t(0)) aynn0N(n)(t)bmx(m)(t) (7)
m0M其频率响应H(j)可表示为(8)式所示的j的有理多项式。
Y()bM(j)MbM1(j)M1...b1jb0 (8) H()NN1X()aN(j)aN1(j)...a1ja0MATLAB的信号处理工具箱提供了专门的函数freqs,用来分析连续时间系统的频率响应,该函数有下列几种调用格式:
[h,w]freqs(b,a)计算默认频率范围内200个频率点上的频率响应的取样值,这200
个频率点记录在w中。
a分别为表示H(j)的有理多项式中分子和分母多项式的系数向hfreqs(b,a,w) b、
量,w为频率取样点,返回值h就是频率响应在频率取样点上的数值向量。
[h,w]freqs(b,a,n)计算默认频率范围内n个频率点上的频率响应的取样值,这n个
频率点记录在w中。freqs(b,a,...) 这种调用格式不返回频率响应的取样值,而是以对数坐标的方式绘出系统的幅频响应和相频响应。
2. 离散时间系统的频率响应
LTI离散时间系统的频率响应定义为单位抽样响应h(n)的离散时间傅里叶变换。
H(e)jnh(n)ejn (9)
对于任意输入信号x(n),输入与输出信号的离散时间傅里叶变换有如下关系
Y(ej)H(ej)X(ej) (10)
因此,系统的频率响应还可以表示为
H(ej)Y(ej)/X(ej) (11)
当系统输入信号为x(n)ejn时,系统的输出为
y(n)ejn*h(n)kej(nk)h(k)ejnH(ej) (12)
由(12)式可知,虚指数信号通过LTI离散时间系统后信号的频率不变,信号的幅度由系统的频率响应的幅度值确定,所以H(e一般情况下离散系统的频率响应H(ejj)表示了系统对不同频率信号的衰减量。
)是复值函数,可用幅度和相位来表示。
H(ej)|H(ej)|ej() (13)
其中|H(ej)|称为系统的幅度响应,()称为系统的相位响应。
若LTI离散系统可以由如下差分方程描述。
ay(ni)bx(nj) (14)
iji0j0NM则由(11)描述的离散时间系统的频率响应H(ejj)可以表示为ej的有理多项式。
Y(ej)b0b1ej...bMejM (15) H(e)X(ej)a0a1ej...aNejNMATLAB的信号处理工具箱提供了专门的函数freqz,用来分析连续时间系统的频率响应,该函数有下列几种调用格式:
[H,w]freqz(b,a,n) b、a分别为有理多项式中分子和分母多项式的系数向量,返回
值H是频率响应在0~pi范围内n个频率等分点上的数值向量,w包含了这n个频率点。
[H,w]freqz(b,a,n,'whole')计算0~2n个频率点上的频率响应的取样值,这n个频
率点记录在w中。
Hfreqz(b,a,w) w为频率取样点,计算这些频率点上的频率响应的取样值。
而是直接绘出系统的幅频响应和freqz(b,a,...)这种调用格式不返回频率响应的取样值,相频响应。 三、实验内容
(1)已知一个RLC电路构造的二阶高通滤波器如下图所示,其中RL, 2C L=0.4H,C=0.05F
1) 计算该电路系统的频率响应及高通截止频率;
Y()X()111111jCRjLRjL经整理
1
(j)2H()Y()/X()(j)210j50 ,使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电路如下图所示。
1)对不同的RC值,用MATLAB画出系统的幅度响应曲线|H()|,观察实验结果,分析图中所示的RC电路具有什么样的频率特性(高通、低通、带通或带阻)?系统的频率特性随着RC值得改变,有何变化规律? 系统微分方程为RCx'(t)x(t)y(t)
RC=input('RC='); b=[1]; a=[RC,1];
[H,w]=freqs(b,a); plot(w,abs(H));
set(gca,'ytick',[0.1 0.3 0.5 0.707 1]);
xlabel('\\omega(rad/s)'); ylabel('Magnitude'); title('|H(j\\omega)|'); grid on;
>> DFTfourth_2_2 RC=0.25
|H(j)|1>> DFTfourth_2_2 RC=1
|H(j)|10.7070.707Magnitude0102030405060708090100Magnitude0.50.50.30.30.10.1012345678910(rad/s)
>> DFTfourth_2_2 RC=16
(rad/s)
>> DFTfourth_2_2 RC=4
|H(j)|1|H(j)|10.7070.707Magnitude0.50.300.10.20.30.40.50.60.70.80.91Magnitude0.50.30.100.10.20.30.40.50.60.70.80.91(rad/s)
(rad/s)
综上分析:RC电路具有低通特性。随着RC值的增大,系统的频率响应衰减加快,低通截止频率逐渐降低。
2) 系统输入信号x(t)cos(100t)cos(3000t),t0~0.2s,该信号包含了一个低频分
量和一个高频分量。试确定适当的RC值,滤除信号中的高频分量。并绘出滤波前后的时域信号波形及系统的频率响应曲线。
x(t)cos(100t)cos(3000t)
X()[(100)(100)(3000)(3000)]
要滤掉高频分量,则要求低通截止频率10003000。
当RC取0.005时,低通截止频率为200rad/s,频率响应曲线如下图:
RC=input('RC='); w=linspace(-4000,4000); t=0:0.001:0.2;
x=cos(100*t)+cos(3000*t); subplot(221); plot(t,x); title('x(t)'); xlabel('t');
X=x*exp(-j*t'*w)*0.001; subplot(223); plot(w,X); xlabel('w'); title('X(w)');
subplot(224); b=[1]; a=[RC 1]; H=freqs(b,a,w); Y=H.*X; plot(w,Y); xlabel('w'); title('Y(w)'); subplot(222); sys=tf(b,a); lsim(sys,x,t) xlabel('t'); title('y(t)')
取RC=0.005
x(t)22y(t)10-1-200.050.1tX(w)0.150.2Amplitude10-1-200.050.1t (seconds)0.150.2Y(w)0.010-0.01-0.020.060.040.020-0.02-4000-20000w20004000-0.03-4000-20000w20004000
(3)已知离散系统的系统框图如图所示。
x(n) 第1个 z-1 第2个 z-1 第M个 z-1 y(n) 1)写出M=8时系统的差分方程和系统函数; 差分方程:y(n)x(n)x(n1)...x(n8)
系统函数:H(z)1zzzzzzzz 2)利用MATLAB计算系统的单位抽样响应。 >> b=ones(1,9); >> a=1;
>> impz(b,a);
Impulse Response10.90.80.7
12345678 Amplitude0.60.50.40.30.20.1001234n (samples)5678
3)试利用MATLAB绘出其系统的零极点分布图、幅频和相频特性曲线,并分析该系统具有怎样的频率特性。 接上小题 >> zplane(b,a)
10.80.60.4Imaginary Part0.20-0.2-0.4-0.6-0.8-1-1-0.50Real Part0.518
>> [H,w]=freqz(b,a); >> subplot(211); >> plot(w,abs(H)); >> grid on;
>> title('|H(e^{j\\Omega})|') 幅频、相频特性曲线如下:
|H(ej)|864200.511.522.53>> axis tight >> subplot(212); >> plot(w,angle(H)); >> grid on
>> title('\heta(\\Omega)')
()20-2-400.511.522.533.5
所以具有低通频率特性
(4)已知一离散时间LTI系统的频率响应H(ej)如下图所示,输入信号
x(n)cos(0.3n)0.5cos(0.8n)。根据式(12)分析正弦信号sin(0t)通过频率响应
为H(ej)的离散时间系统的响应,并根据分析结果计算系统对于x(n)的响应y(n),用MATLAB
绘出系统输入与输出波形。
x(n)cos(0.3n)0.5cos(0.8n)ej0.3nej0.3nej0.8nej0.8n24由式(12):
y(n)ejn*h(n)kej(nk)h(k)ejnH(ej)
y(n)ej0.3nH(ej0.3)ej0.3nH(ej0.3)ej0.8nH(ej0.8)ej0.8nH(ej0.8)24 j0.3nj0.3nee222cos(0.3n)
n=-30:30;
x=cos(0.3*pi*n)+0.5*cos(0.8*pi*n);
y=2*cos(0.3*pi*n); subplot(211); stem(n,x,'filled'); xlabel('n'); title('x(n)'); subplot(212); stem(n,y,'filled'); xlabel('n'); title('y(n)')
x(n)210-1-2-30-20-100ny(n)102030210-1-2-30-20-100n102030
观察实验结果,分析图中所示的系统具有什么样的的频率特性(高通、低通、带通或带阻)?从输入输出信号上怎么反映出系统的频率特性? 答:
由实验结果,图中所示的系统具有低通的频率特性。
将输入、输出信号二者对比可知,输入信号中的高频分量在输出信号中被滤掉了,即系统具有低通性。
四、实验收获体会
这次实验加深了我对LTI系统频率响应的基本概念的掌握和理解,学会了利用MATLAB对LTI系统进行频域分析的方法。
因篇幅问题不能全部显示,请点此查看更多更全内容