实验六、系统的频域分析

《实验六、系统的频域分析》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、重庆三峡学院《信号与系统分析》实验实验六、系统的频域分析1实验目的1）学会利用MATLAB对连续系统进行频域分析；2）学会利用MATLAB分析离散系统函数的零极点分布与时域特性的关系；3）学会利用MATLAB进行离散系统的频率特性分析。2实验原理及实例分析（实验原理见教材的第五章、第六章及第七章）2.1连续LTI系统的频率特性例1：已知连续LTI系统的微分方程为求该系统的频率响应，并用MATLAB绘出其幅频特性和相频特性图。解：MATLAB程序如下：clccloseallclearallb=[137];a=[11085];w=-3

2、*pi:0.01:3*pi;H=freqs(b,a,w);subplot(211);plot(w,abs(H),'Linewidth',2);grid;xlabel('omega(rad/s)');title('

3、H(jomega)

4、');subplot(212);plot(w,angle(H),'Linewidth',2);grid;xlabel('omega(rad/s)');title('phi(omega)');程序产生的图形如图1所示。6重庆三峡学院《信号与系统分析》实验图1例1程序产生的波形图2.2连续LTI

5、系统的频域分析例2：设系统的频率响应为，若输入信号为，用MATLAB命令求其稳态响应。解：MATLAB程序如下：clccloseallclearallt=0:0.01:20;w1=1;w2=10;H1=1/(-w1^2+j*3*w1+2);H2=1/(-w2^2+j*3*w2+2);f=5*cos(t)+2*cos(10*t);yss=abs(H1)*cos(w1*t+angle(H1))+abs(H2)*cos(w2*t+angle(H2));subplot(211);plot(t,f,'Linewidth',2);grid;x

6、label('t(sec)');title('f(t)');subplot(212);plot(t,yss,'Linewidth',2);grid;xlabel('t(sec)');title('y_s_s(t)');程序产生的图形如图2所示。6重庆三峡学院《信号与系统分析》实验图2例2程序产生的波形图2.3离散系统函数的零极点分布与其时域特性的关系例3：用MATLAB命令绘出下列系统函数的零极点分布图以及对应的时域单位取样响应的波形图，并分析系统函数的极点对时域波形的影响。1）2）3）4）5）6）解：MATLAB程序如下：clc

7、closeallclearallb1=[10];a1=[10.8];subplot(221);zplane(b1,a1);subplot(222);impz(b1,a1,30);grid;b2=[10];a2=[1-1.20.72];subplot(223);zplane(b2,a2);subplot(224);impz(b2,a2,30);grid;figureb3=[10];a3=[1-1];subplot(221);zplane(b3,a3);subplot(222);impz(b3,a3,30);grid;6重庆三峡学院《

8、信号与系统分析》实验b4=[10];a4=[1-1.61];subplot(223);zplane(b4,a4);subplot(224);impz(b4,a4,30);grid;figureb5=[10];a5=[1-1.2];subplot(221);zplane(b5,a5);subplot(222);impz(b5,a5,30);grid;b6=[10];a6=[1-21.36];subplot(223);zplane(b6,a6);subplot(224);impz(b6,a6,30);grid;程序产生的图形如图3所示

9、。6重庆三峡学院《信号与系统分析》实验图3例3程序产生的图形2.4离散LTI系统的频率特性分析例4：用MATLAB命令绘制离散系统的频率响应曲线。解：MATLAB程序如下：clccloseallclearallb=[1-0.960.9028];a=[1-1.560.8109];[H,w]=freqz(b,a,600,'whole');subplot(211);plot(w-pi,abs(H),'Linewidth',2);grid;xlabel('omega(rad/s)');title('

10、H(jomega)

11、');subp

12、lot(212);plot(w-pi,angle(H),'Linewidth',2);grid;xlabel('omega(rad/s)');title('phi(omega)');程序产生的图形如图4所示。6重庆三峡学院《信号与系统分析》实