第7章系统频率响应仿真精选文档.ppt

《第7章系统频率响应仿真精选文档.ppt》由会员分享，可在线阅读，更多相关《第7章系统频率响应仿真精选文档.ppt（38页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第7章系统频率响应仿真本讲稿第一页，共三十八页function y=f(x)sym xif x10 y=2*xelse y=(1-x)/4endx=input(x=)if x10 y=2*xelse y=(1-x)/4end2023/4/112本讲稿第二页，共三十八页functionsysm1,sysm2=modelm4(sysp)m1=sysp(1);m2=sysp(2);k=sysp(3);c=sysp(4);A=0 0 1 0;0 0 0 1;-k/m1 k/m1-c/m1 c/m1;k/m2-k/m2 c/m2-c/m2;B=0 0 0 1/m2;C=0 1 0 0;M=1 0 0 0

2、;D=0;sys2=ss(A,B,C,D);sys1=ss(A,B,M,D);sysm2=tf(sys2);sysm1=tf(sys1);sysp=12 38 1000 0.1;X1_F,X2_F=modelm4(sysp)2023/4/113本讲稿第三页，共三十八页试用试用MATLAB语言表示如图所示系统。当分别以语言表示如图所示系统。当分别以x2和和 f 为系统输出、输入时的传为系统输出、输入时的传递函数模型和状态空间模型。递函数模型和状态空间模型。Transfer function:0.1786 s2+0.02551 s+0.3571-s4+0.2679 s3+2.005 s2+0.25

3、 sfunctionsysm=modelm3(sysp)k=sysp(1);c1=sysp(2);c2=sysp(3);m1=sysp(4);m2=sysp(5);A=0 1 0 0;-k/m1-c1/m1 0 c1/m1;0 0 0 1;0 c1/m2 0-(c1+c2)/m2;B=0 0 0 1/m2;C=0 0 1 0;D=0;sys1=ss(A,B,C,D);sysm=tf(sys1);sysp=7 0.5 0.2 3.5 5.6;X2_F=modelm3(sysp)2023/4/114本讲稿第四页，共三十八页2023/4/115本讲稿第五页，共三十八页1/(m2s2)c1s1/(m1

4、s2)kc2sFx2x1FkFc1Fc1Fc2clcsys0=tf(1,3.5 0 0)sys1=feedback(sys0,7,-1)sys05=tf(0.5 0,1)sys2=series(sys1,sys05)sys3=feedback(sys2,1,-1)sys4=tf(1,5.6 0 0)sys02=tf(0.2 0,1)sys5=feedback(sys4,sys02,-1)n,d=tfdata(sys1,v)sys7=tf(d,n)sys8=series(sys7,sys3)sys=feedback(sys5,sys8,-1)2023/4/116本讲稿第六页，共三十八页2023/

5、4/117本讲稿第七页，共三十八页第第7章章 系统频率响应及其仿真系统频率响应及其仿真频率特性的一般概念频率特性的一般概念频率特性的图解方法频率特性的图解方法频率特性的频率特性的MATLAB函数函数离散系统的频域仿真离散系统的频域仿真系统分析图形用户界面系统分析图形用户界面系统的频域分析系统的频域分析2023/4/118本讲稿第八页，共三十八页7.1 频率响应与频率特性频率响应与频率特性频率特性是控制工程最基本的概念之一频率特性是控制工程最基本的概念之一何谓频率响应？何谓频率响应？频率响应有什么特性？频率响应有什么特性？如何求系统的频率响应？如何求系统的频率响应？为什么可以用频率响应来研究系统

6、的性能？为什么可以用频率响应来研究系统的性能？何谓频率特性？何谓频率特性？演示程序：演示程序：pltxys.mdl2023/4/119本讲稿第九页，共三十八页7.1.1 频率特性的概念与定义频率特性的概念与定义频率响应的特点：频率响应的特点：稳定系统稳定系统Xi=AisintXo=Ao()sint+()幅值比幅值比它是频率它是频率的函数，的函数，称为幅频特性；称为幅频特性；相位差相位差它是频率它是频率的函数，的函数，称为相频特性。称为相频特性。结论结论幅值比、相位差随频率幅值比、相位差随频率变化变化频率响应的这种特性，用频率特性描述。频率响应的这种特性，用频率特性描述。2023/4/1110本

7、讲稿第十页，共三十八页频率特性的定义频率特性的定义 线性稳定系统在正弦信号作用下，当频率从零变化到无穷线性稳定系统在正弦信号作用下，当频率从零变化到无穷时，稳态输出与输入的幅值比、相位差随频率变化的特性，称时，稳态输出与输入的幅值比、相位差随频率变化的特性，称为频率特性。频率特性由幅频特性和相频特性两部分组成。为频率特性。频率特性由幅频特性和相频特性两部分组成。幅频特性幅频特性相频特性相频特性 系统稳态正弦输出信号与相应的正弦输人信号的幅值之比随输入频率的变系统稳态正弦输出信号与相应的正弦输人信号的幅值之比随输入频率的变比而变化的特性称为幅频特性，它描述了系统对输入信号幅值的放大、衰减特比而变

8、化的特性称为幅频特性，它描述了系统对输入信号幅值的放大、衰减特性。性。系统稳态正弦输出信号与相应的正弦输入信号的相位之差随输入频率的变系统稳态正弦输出信号与相应的正弦输入信号的相位之差随输入频率的变化而变化的特性称为相频特性，它描述了系统输出信号相位对输入信号相位的化而变化的特性称为相频特性，它描述了系统输出信号相位对输入信号相位的超前、迟后特性。超前、迟后特性。2023/4/1111本讲稿第十一页，共三十八页频率特性的求取方法频率特性的求取方法两两个个复复数数积积或或商商的的模模等等于于两两个个复复数数的的模模的的积积或或商商；两两个个复复数数积或商的相位等于两个复数的相位的和（加）积或商的

9、相位等于两个复数的相位的和（加）或差（减）。或差（减）。注意复数的运算：注意复数的运算：2023/4/1112本讲稿第十二页，共三十八页求系统的稳态输出求系统的稳态输出已知系统传递函数和输入，求该系统的稳态输出。已知系统传递函数和输入，求该系统的稳态输出。解：求稳态输出，首先必须求得幅频特性和相频特性。解：求稳态输出，首先必须求得幅频特性和相频特性。求频率特性求频率特性求稳态输出求稳态输出2023/4/1113本讲稿第十三页，共三十八页7.1.2 频率特性的图解方法频率特性的图解方法设系统的频率特性为设系统的频率特性为2023/4/1114本讲稿第十四页，共三十八页Bode图坐标系图坐标系 2

10、023/4/1115本讲稿第十五页，共三十八页7.1.3 稳定裕量稳定裕量相位裕度相位裕度 在幅值交界频率上，使系统达到临界稳定状态所需附加的相位滞后量。在幅值交界频率上，使系统达到临界稳定状态所需附加的相位滞后量。幅值裕度幅值裕度Kg 在相位交界频率上，使开环幅值达到在相位交界频率上，使开环幅值达到1所需放大的倍数，用所需放大的倍数，用Kg或或GM表示。表示。定义式定义式定义式定义式计算式计算式计算式计算式2023/4/1116本讲稿第十六页，共三十八页在在BODEBODE图上的图上的幅值、相位裕度幅值、相位裕度2023/4/1117本讲稿第十七页，共三十八页7.2 频率特性的频率特性的MA

11、TLAB函数函数频率响应的计算频率响应的计算频率特性的图示方法频率特性的图示方法本讲稿第十八页，共三十八页7.2.1 频率响应的计算频率响应的计算Y=polyval（P，X）Abs（Y）Angle（Y）幅频特性幅频特性相频特性相频特性例例71多项式计算指令多项式计算指令多项式系数向量jh=freqs（b,a,w）h,w=freqs（b,a）h,w=freqs（b,a,f）freqs（b,a,w）2023/4/1119本讲稿第十九页，共三十八页例例7-12023/4/1120本讲稿第二十页，共三十八页例例7-22023/4/1121本讲稿第二十一页，共三十八页例例7-2的图的图图图72图图732

12、023/4/1122本讲稿第二十二页，共三十八页7.2.2 频率特性图示法频率特性图示法Nyquist图的绘制图的绘制Bode图的绘制图的绘制计算幅值、相位裕度计算幅值、相位裕度例例7-3K=5，K=30nyquist(sys)nyquist(sys,w)nyquist(sys1,sys2,.,sysn)re,im,w=nyquist(sys)nyquist(sys1,sys2,.,sysn,w)2023/4/1123本讲稿第二十三页，共三十八页例例7-3的程序的程序2023/4/1124本讲稿第二十四页，共三十八页7.2.2 频率特性图示法频率特性图示法Bode图的绘制图的绘制计算幅值、相位

13、裕度计算幅值、相位裕度bode(sys)bode(sys,w)bode(sys1,sys2,.,sysn)bodemag(sys)mag,phase,w=bode(sys)margi(sys)Gm,Pm,Wcg,Wcp=margi(sys)Gm,Pm,Wcg,Wcp=margi(mag,phase,w)2023/4/1125本讲稿第二十五页，共三十八页例例7-4和和7-5的程序的程序2023/4/1126本讲稿第二十六页，共三十八页7.3 离散系统频域仿真离散系统频域仿真Z变换变换例例7-62023/4/1127本讲稿第二十七页，共三十八页例例7-67-6仿真结果仿真结果2023/4/1128

14、本讲稿第二十八页，共三十八页仿真结果仿真结果2023/4/1129本讲稿第二十九页，共三十八页7.4 系统分析图形用户界面系统分析图形用户界面在指令窗中建立系统模型在指令窗中建立系统模型在指令窗中输入：在指令窗中输入：ltiview点击菜单点击菜单File，选择，选择Import在在LTI浏览器中，选择浏览器中，选择sys，调入，调入图图7-14 LTI Viewer窗口窗口图图7-15 系统模型调入窗口系统模型调入窗口显示调入系统模型的仿真曲线显示调入系统模型的仿真曲线2023/4/1130本讲稿第三十页，共三十八页例例7-7 LTI Viewer使用演示使用演示在窗口中右击，弹出在窗口中右

15、击，弹出现场功能菜单。现场功能菜单。2023/4/1131本讲稿第三十一页，共三十八页现场菜单的主要功能现场菜单的主要功能Plot TypeCharacteristicsGridPropertiessystemsNormalizeFull View2023/4/1132本讲稿第三十二页，共三十八页多个图形窗口显示多个图形窗口显示在在LTI Viewer窗口中，打开窗口中，打开Edit菜单，选择菜单，选择Plot Configurations，弹出图形配置窗口，弹出图形配置窗口选择显示图形数和图形类型，点击选择显示图形数和图形类型，点击OK2023/4/1133本讲稿第三十三页，共三十八页多个图

16、形窗口显示多个图形窗口显示右击各图形右击各图形选择选择characteristics标注特征量标注特征量左双击各图形左双击各图形弹出弹出Property Editor编辑线型、字编辑线型、字体大小、色彩体大小、色彩等。等。2023/4/1134本讲稿第三十四页，共三十八页第第7 7章章 习题习题3 3 G1G2G3-试用试用LTI Wiewre分析系统的稳定性，并求出系统的稳态裕度和分析系统的稳定性，并求出系统的稳态裕度和单位阶跃响应峰值。单位阶跃响应峰值。2023/4/1135本讲稿第三十五页，共三十八页习题习题3 3 演示演示2023/4/1136本讲稿第三十六页，共三十八页2023/4/1137本讲稿第三十七页，共三十八页2023/4/1138本讲稿第三十八页，共三十八页