《试验二用MATLAB建立传递函数模型》由会员分享，可在线阅读，更多相关《试验二用MATLAB建立传递函数模型（74页珍藏版）》请在金锄头文库上搜索。
1、精选优质文档-倾情为你奉上自动控制原理实验指导书实验指导书北京科技大学自动化学院控制科学与工程系 2013年4月 专心-专注-专业目 录实验一 典型系统的时域响应和稳定性分析1实验二 用MATLAB建立传递函数模型5实验三 利用MATLAB进行时域分析13实验四 线性定常控制系统的稳定分析25实验五 利用MATLAB绘制系统根轨迹29实验六 线性系统的频域分析37实验七 基于MATLAB控制系统频域法串联校正设计52附录1 MATLAB简介59附录2 SIMULINK简介68实验一 典型系统的时域响应和稳定性分析一、实验目的1研究二阶系统的特征参量 (、n) 对过渡过程的影响。2研究二阶对象的
2、三种阻尼比下的响应曲线及系统的稳定性。3熟悉Routh判据，用Routh判据对三阶系统进行稳定性分析。二、实验设备PC机一台，TD-ACC+教学实验系统一套。三、实验原理及内容 1典型的二阶系统稳定性分析(1) 结构框图：如图1-1所示。 图1-1(2) 对应的模拟电路图：如图1-2所示。图1-2(3) 理论分析系统开环传递函数为：G(s)=?开环增益： K=?先算出临界阻尼、欠阻尼、过阻尼时电阻R的理论值，再将理论值应用于模拟电路中，观察二阶系统的动态性能及稳定性，应与理论分析基本吻合。在此实验中由图1-2，可以确地1-1中的参数。， ， 系统闭环传递函数为：其中自然振荡角频率：；阻尼比：。
3、2典型的三阶系统稳定性分析(1) 结构框图：如图1-3所示。 图1-3(2) 模拟电路图：如图1-4所示。图1-4(3) 理论分析系统的开环传函为: 系统的特征方程为:。 (4) 实验内容实验前由Routh判断得Routh行列式为：S3 S2 S1 S0 为了保证系统稳定，第一列各值应为正数，因此可以确定 系统稳定 K值的范围 系统临界稳定K 系统不稳定K值的范围四、实验步骤1）将信号源单元的“ST”端插针与“S”端插针用“短路块”短接。由于每个运放单元均设置了锁零场效应管，所以运放具有锁零功能。将开关分别设在“方波”档和“500ms12s”档，调节调幅和调频电位器，使得“OUT”端输出的方波
4、幅值为1V，周期为10s左右。2.）典型二阶系统瞬态性能指标的测试(1) 按模拟电路图1-2接线，将1中的方波信号接至输入端，取R = 10K。(2) 用示波器观察系统响应曲线C(t)，测量并记录超调MP、峰值时间tp和调节时间tS。(3) 分别按R = 20K；40K；100K；改变系统开环增益，观察响应曲线C(t)，测量并记录性能指标MP、tp和tS，及系统的稳定性。并将测量值和计算值进行比较 (实验前必须按公式计算出)。将实验结果填入表1-1中。3）典型三阶系统的性能(1) 按图1-4接线，将1中的方波信号接至输入端，取R = 30K。(2) 观察系统的响应曲线，并记录波形。(3) 减小
5、开环增益 (R = 41.7K；100K)，观察响应曲线，并将实验结果填入表1-3中。表1-3中已填入了一组参考测量值，供参照。五、实验现象分析1）典型二阶系统瞬态性能指标实验测试值表1-1参数项目R(K)KnC(tp)C()MP (%)tP (s)tS (s)响应情况理论值测量值理论值测量值理论值测量值0 1过阻尼200 其中， ，2）典型三阶系统在不同开环增益下的响应情况表1-2R(K)开环增益K稳定性表13R(K)开环增益K稳定性30167不稳定发散41 .712临界稳定等幅振荡1005稳定衰减收敛注意：在做实验前一定要进行对象整定，否则将会导致理论值和实际测量值相差较大。实验二 用MA
6、TLAB建立传递函数模型一、实验目的（1）熟悉MATLAB实验环境，掌握MATLAB命令窗口的基本操作；（2）掌握MATLAB建立控制系统数学模型的命令及模型相互转换的方法；（3）掌握使用MATLAB命令化简模型基本连接的方法；（4）学会使用Simulink模型结构图化简复杂控制系统模型的方法。二、实验原理及内容控制系统常用的数学模型有四种：传递函数模型（tf对象）、零极点增益模型（zpk对象）、结构框图模型和状态空间模型（ss对象）。经典控制理论中数学模型一般使用前三种模型，状态空间模型属于现代控制理论范畴。1、有理函数模型线性连续系统的传递函数模型可一般地表示为： (1)将系统的分子和分母
7、多项式的系数按降幂的方式以向量的形式输入给两个变量和，就可以轻易地将传递函数模型输入到MATLAB环境中。命令格式为：; （2）;（3）用函数tf ( ) 来建立控制系统的传递函数模型，该函数的调用格式为：Gtf(num，den);（4）例1一个简单的传递函数模型:可以由下面的命令输入到MATLAB工作空间中去。num=1，2;den=1，2，3，4，5;G=tf(num，den)运行结果：Transfer function: s + 2-s4 + 2 s3 + 3 s2 + 4 s + 5这时对象G可以用来描述给定的传递函数模型，作为其它函数调用的变量。例2一个稍微复杂一些的传递函数模型：
8、num=6*1，5; den=conv(conv(1，3，1，1，3，1)，1，6); tf(num,den)运行结果Transfer function: 6 s + 30-s5 + 12 s4 + 47 s3 + 72 s2 + 37 s + 6其中conv()函数（标准的MATLAB函数）用来计算两个向量的卷积，多项式乘法也可以用这个函数来计算。该函数允许任意地多层嵌套，从而表示复杂的计算。【自我实践1】建立控制系统的传递函数模型：2、零极点模型线性系统的传递函数还可以写成极点的形式： 将系统增益K、零点-zi和极点-pj以向量的形式输入给三个变量、Z和P，命令格式为：(6) (7) (8
9、)用函数命令zpk ( ) 来建立系统的零极点增益模型，其函数调用格式为：G=zpk(Z,P,KGain) (9)例3某系统的零极点模型为： KGain=6;z=-1.9294;-0.0353+0.9287j;-0.0353-0.9287j;p=-0.9567+1.2272j;-0.9567-1.2272j;0.0433+0.6412j;0.0433-0.6412j;G=zpk(z,p,KGain)运行结果:Zero/pole/gain: 6 (s+1.929) (s2 + 0.0706s + 0.8637)-(s2 - 0.0866s + 0.413) (s2 + 1.913s + 2.42
10、1)注意：对于单变量系统，其零极点均是用列向量来表示的，故z、p向量中各项均用分号（;）隔开。【自我实践2】建立控制系统的零极点模型：3、控制系统模型间的相互转换零极点模型转换为传递函数模型：num , den = zp2tf ( z, p, k ) 传递函数模型转化为零极点模型：z , p , k = tf2zp ( num, den ) 例4 给定系统传递函数为：对应的零极点模型可由下面的命令得出num=6.8, 61.2, 95.2;den=1, 7.5, 22, 19.5, 0;G=tf(num,den);G1=zpk(G)显示结果：Zero/pole/gain: 6.8 (s+7)
11、(s+2)-s (s+1.5) (s2 + 6s + 13)可见，在系统的零极点模型中若出现复数值，则在显示时将以二阶因子的形式表示相应的共轭复数对。例5给定零极点模型：可以用下面的MATLAB命令立即得出其等效的传递函数模型。输入程序时要注意大小写。Z=-2,-7;P=0,-3-2j,-3+2j,-1.5;K=6.8;G=zpk(Z,P,K);G1=tf(G)结果显示：Transfer function: 6.8 s2 + 61.2 s + 95.2-s4 + 7.5 s3 + 22 s2 + 19.5 s【自我实践3】已知系统传递函数，求其等效的零极点模型。【自我实践4】建立控制系统的多项
12、式模型：。4、反馈系统结构图模型设反馈系统结构图如图2- 1所示。图2-1反馈系统结构图两个环节反馈连接后，其等效传递函数可用feedback ( ) 函数求得，其调用格式为：sys = feedback（G1, G2, sign）其中sign是反馈极性，sign缺省时，默认为负反馈，sign-1；正反馈时，sign1；单位反馈时，G21，且不能省略。series( )函数：实现两个模型的串联；多于两个必须嵌套使用。parallel( )函数：实现两个模型的并联；多于两个必须嵌套使用。例6若反馈系统如图2-1中的两个传递函数分别为： ，则反馈系统的传递函数： G1=tf(1,1,2,1); G2=tf(1,1,1); G=feedback(G1,G2)运行结果:Transfer function: s + 1-s3 + 3 s2 + 3 s + 2若采用正反馈连接结构输入命
