资源描述：
《实验4 LTI系统的频域分析.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、一，实验目的针对LTI系统频率响应，加深了对于基本概念的掌握和理解，学习并掌握了关于LTI系统频率特性的分析方法。二，实验原理1.连续时间系统的频率响应调用函数freqs:[h,w]:freqs(b,a)计算默认频率范围内200个频率点上的频率响应的取样值，这200个频率点记录在w中。h=freqs(b,a,w)b,&分别为表示H(jw)的有理多项式中分子和分母多项式的系数向量，w为频率取样点，返冋值h就是频率响应在频率取样点上的数值向量。Fh,w]:freqs(b,a,n)计算默认频率范围内n个频率点上的频率
2、响应的取样值，这n个频率点记录在w中。freqs(b,a,…)这种调用格式不返冋频率响应的取样值，而是以对数坐标的方式绘出系统的频率响应和相频响应。2.离散时间系统的频率响应调用函数freqz:[H,w]:freqz(b,a,"whole')计算0〜2nn个频率点上的频率响应的取样值，这n个频率点记录在w中。H=freqz(b,a,n)b,a分别为有理多项式中分子和分母多项式的系数向量,返回值H就是频率响应在0到pi范围内n个频率等分点上的数值向量，w包含了这n个频率响应。[H,w]:freqz(b,a,w)w
3、为频率取样点，计算这些频率点上的频率响应的取样值。freqz(b,a,…)这种调用格式不返回频率响应的取样值，而是直接绘出系统的频率响应和相频响应。三,实验内容(1)已知一个RLC电路构造的二阶高通滤波器如图，其中厶R=—,L=0.4乩C=0.05F①计算该电路系统的频率响应及高通截止频率。H((0)=(閻(“y)‘+10Je+50②利用MATLAB绘制幅度响应和相位响应曲线，比较系统的频率特性与理论计算的结果是否一致。MATLAB程序如下：b=[0.0400]a=[0.040.42][H,w]=freqs(b
4、za)subplot(211)plot(w,abs(H))set(gc“1xtick')set(gca,'ytick'z[00.40.7071])xlabel(1omega(rad/s)T)ylabel(1Magnitude1)title('
5、H(jomega)I')gridonsubplot(212)plot(wzangle(H))set(gcaz1xtick1)xlabel(1omega(rad/s)1)ylabel('Phase')title('phi(omega)')gridon①prnc6e
6、乏
7、H(P)
8、0.7070.4010203040506070800403(rad/s)讽3)囲2Ql01020304050607080901003(rad/s)(2)已知一个RC系统电路如图。%1对不同的RC值，用MATLAB画岀系统的幅度响应曲线
9、H(w)
10、,观察实验结果，分析如图所示电路具有什么样的频率特性？系统的频率特性随着RC值的改变，有何变化规律？MATLAB程序如下：A=input('A=,)b=[1]a=[A1][H,w]=freqs(b,a)plot(w,abs(H))set(gca,ytick
11、JO0.40.7071])xlabel(,omega(rad/s),)ylabel('Magnitude')title('
12、H(jomega)「)gridon程序执行如下：A=0.01
13、HS01002003004005006007008009001000cu(rad/s)7▲o①7Co①pruc6osA=1Magnitudeos§020.020.030.040・050・06067008OS0・1E(rad、s)A巴00-IsMagnitudeo由程序执行结果可以看出，RC电路具有带通特性，随着RC值的减小，
14、带通频率增加。%1系统输入信号x(t)=cos(100t)+cos(3000t),t二0~0.2s,该信号包含了一个低频分最和一个高频分立，试确定适当的RC值，滤岀信号中的高频分量，并绘出滤波前后的时域信号波形及系统的频率响应曲线。MATLAB程序如下：A=input('A=');b=[1]a=[A1][H,w]=freqs(b,a)plot(w,abs(H))set(gca,'ytick',[O0.40.7071])xlabel('omega(rad/s)')ylabel('Magnitude')title
15、('
16、H(jomega)
17、,)gridonMATLAB程序执行如下：A=0.0004
18、H(M
19、0.7070.4010002000300040005000600070008000900010000•(rad/s)①prqc6w2t=0：0.0001：0.2x=cos(100*t)+cos(3000*t)plot(t,x)xlabel('t')title('x(t)')可得：输入：x(t
