北理工信号与系统实验(4)

实验4LTI系统的频域分析一、实验目的1.加深对LTI系统频率响应基本概念的掌握和理解2.学习和掌握LTI系统频率特性的分析方法二、实验原理1.连续时间系统的频率响应系统的频率响应定义为系统单位冲激响应h(t)的傅里叶变换，即()()jHhed若LTI连续时间系统的单位冲激响应为h(t)，输入信号为x(t)，根据系统的时域分析可知系统的零状态响应为(t)x(t)h(t)y，对等式两边分别求傅里叶变换，根据时域卷积定理可以得到()X()H()Y。因此，系统的频率响应还可以由系统的零状态响应和输入的傅里叶变换之比得到H()Y()/X()。H()反映了LTI连续时间系统对不同频率信号的响应特性，是系统内在固有的特性，与外部激励无关。H()又可以表示为()H()=|H()|ej，其中|H()|称为系统的幅度响应，()称为系统的相位响应。当虚指数信号ejt作用于LTI系统时，系统的零状态响应y(t)仍为同频率的虚指数信号，即y(t)=eH()jt。对于由下述微分方程描述的LTI连续时间系统(n)(m)00(t)(t)NMnmnmaybx，其频率响应H(j)可以表示为有理多项式11101110(j)(j)...j()()()(j)(j)...jMMMMNNNNbbbbYHXaaaaMATLAB的信号处理工具箱提供了专门的函数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(t)的离散时间傅里叶变换。(e)(n)ejjnHh对于任意输入信号x(n)，输入与输出信号的离散时间傅里叶变换有如下关系Y(e)H(e)X(e)jjj因此，系统的频率响应还可以表示为H(e)=Y(e)/X(e)jjj当系统输入信号为(n)ejnx时，系统的输出(nk)y(n)e(n)(k)e(e)jnjjnjkhehH虚指数信号通过LTI离散时间系统后信号的频率不变，信号的幅度由系统频率响应的幅度值确定，所以H(e)j表示了系统对不同频率信号的衰减量。一般情况下离散系统的频率响应H(j)是复值函数，可用幅度和相位表示。()H(e)=|H(e)|ejjj，其中|H(e)|j称为系统的幅度响应，()称为系统的相位响应。若LTI离散系统可以由如下差分方程描述00(n)(n)NMijijayibxj，其频率响应H(e)j可以表示为有理多项式0101e...e(e)(e)(e)e...ejjMjjMjjjNNbbbYHXaaa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电路构造的二阶高通滤波器，其中FCHLCLR05.0,4.0,2（1）计算该电路系统的频率响应及高通截止频率解：由电路图50)(10)()(1111111)()()(22jjjCjLjRLjRtxtyH高通截止频率为7.07（2）利用MATLAB绘制幅度响应和相位响应曲线，比较系统的频率特性与理论计算的结果是否一致。解：b=[100];a=[11050];[H,w]=freqs(b,a);subplot(211);plot(w,abs(H));set(gca,'xtick',[0;10]);set(gca,'ytick',[00.40.7071]);xlabel('\omega(rad/s)');ylabel('Magnitude');title('|H(j\omega)|');gridon;subplot(212);plot(w,angle(H));set(gca,'xtick',[0;10]);xlabel('\omega(rad/s)');ylabel('Phase');title('|phi(\omega)|');gridon;01000.40.7071(rad/s)Magnitude|H(j)|01001234(rad/s)Phase|phi()|由图像，系统的频率特性与理论计算的结果一致。2.已知一个RC电路（1）对不同的RC值，用MATLAB画出系统的幅度响应曲线|H()|，观察实验结果，分析RC电路具有什么样的频率特性（高通、低通、带通或带阻）？系统的频率特性随着RC值的改变，有何变化规律？解：电路的频率响应为1111)()()(RCjRCjCjtxtyHsymsrck;r=input('r=');c=input('c=');k=r*c;b=[1];a=[k1];[H,w]=freqs(b,a);subplot(211);plot(w,abs(H));set(gca,'xtick',[0,10]);set(gca,'ytick',[00.40.7071]);xlabel('\omega(rad/s)');ylabel('Magnitude');title('|H(j\omega)|');gridon;subplot(212);plot(w,angle(H));set(gca,'xtick',[0,10]);xlabel('\omega(rad/s)');ylabel('Phase');title('|phi(\omega)|');gridon;r=1c=0.101000.40.7071(rad/s)Magnitude|H(j)|010-1.5-1-0.50(rad/s)Phase|phi()|r=10c=0.101000.40.7071X:1.012Y:0.703(rad/s)Magnitude|H(j)|010-1.5-1-0.50(rad/s)Phase|phi()|r=100c=0.1000.40.7071X:0.1012Y:0.703(rad/s)Magnitude|H(j)|0-1.5-1-0.50(rad/s)Phase|phi()|（2）系统输入信号(t)cos(100t)cos(3000t),t0~0.2sx，该信号包含了一个低频分量和一个高频分量。试确定适当的RC值，滤除信号中的高频分量，并绘出滤波前后的时域信号波形及系统的频率响应曲线。解：b=[1];a=[0.00151];h=tf(b,a);[H,w]=freqz(b,a);t=0:0.001:0.2;x=cos(100*t)+cos(3000*t);plot(t,x);lsim(h,x,t);plot(w,abs(H));xlabel('\omega(rad/s)');ylabel('Magnitude');title('|H(j\omega)|');gridon;RC=0.001500.020.040.060.080.10.120.140.160.180.2-2-101200.511.522.533.50.9980.99911.0011.0021.003(rad/s)Magnitude|H(j)|LinearSimulationResultsTime(seconds)AmplitudeRC=0.00100.020.040.060.080.10.120.140.160.180.2-2-101200.511.522.533.50.9980.99911.0011.0021.003(rad/s)Magnitude|H(j)|LinearSimulationResultsTime(seconds)Amplitude3.已知离散系统的系统框图（1）写出M=8时系统的差分方程和系统函数解：差分方程：y(n)=x(n)+x(n-1)+x(n-2)+x(n-3)+x(n-4)+x(n-5)+x(n-6)+x(n-7)系统函数：8765432111)(zzzzzzzzzH（2）利用MATLAB计算系统的单位抽样响应b=[111111111];a=[1];impz(b,a,-10:10);-10-8-6-4-2024681000.10.20.30.40.50.60.70.80.91n(samples)AmplitudeImpulseResponse（3）试利用MATLAB绘出其系统零极点分布图、幅频和相频特性曲线，并分析该系统具有怎样的频率特性。解：系统的极点分布图程序及图像如下：b=[111111111];a=[1];zplane(b,a,-10:10);-1-0.500.51-1-0.8-0.6-0.4-0.200.20.40.60.818RealPartImaginaryPart系统的幅频和相频特性曲线如下：b=[111111111];a=[1];[H,w]=freqz(b,a);subplot(211);plot(w/pi,abs(H));xlabel('\omega(\pi)');ylabel('Magnitude');title('|H(e^j^\Omega)|');gridon;subplot(212);plot(w/pi,angle(H)/pi);xlabel('\omega(\pi)');ylabel('Phase(\pi)');title('\theta(\Omega)');gridon;00.10.20.30.40.50.60.70.80.910510()Magnitude|H(ej)|00.10.20.30.40.50.60.70.80.91-1-0.500.5()Phase()()3.已知一离散时间LTI系统的频率响应H(e)j，输入信号为(n)cos(0.3n)0.5cos(0.8n)x。试分析正弦信号0sin(t)通过频率响应为H(e)j的离散时间系统的响应，并根据分析结果计算系统对于x(n)的响应y(n)，用MATLAB绘出系统输入与输出波形。观察实验结果，分析系统具有什么样的频率特性（高通、低通、带通或带阻）?从输入输出信号上怎么反映出系统的频率特性？解：n=0:1:50;x=cos(0.3*pi*n)+0.5*cos(0.8*pi*n);y=sinc(0.5*pi*n);z=conv(x,y);subplot(211);stem(x,'filled');axis([050-22]);subplot(212);stem(x,'filled');axis([050-22]);05101520253035404550-2-101205101520253035404550-2-1012四、实验总结通过本次实验，掌握和理解了LTI系统频率响应的