资源描述:
《控制系统matlab常用函数.doc》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、1.tfnum=[013];den=[1221];tf(num,den)2.zpkk=10;z=[-5];p=[-0.5-2-3];sys=zpk(z,p,k)3.ord2已知二阶系统的自然频率和阻尼比,建立传递函数[num,den]=ord2(1,0.5);G=tf(num,den)4.residue多项式模型转化为部分分式展开式模型,或反过来,根据输入参数不同而异5.zp2tf6.tf2zp7.series串联连接8.parallel并联连接9.feedback反馈连接numg=[251];deng=[123];numh=[510];denh=[110];G
2、=tf(numg,deng)H=tf(numh,denh)Gm=feedback(G,H)10.ss2tf将状态空间模型转化成多项式模型11.tf2ss传递函数到状态空间转换12.eig求系统的特征根num=10;den=[1210];roots(den)sys=tf(num,den);eig(sys)13.damp求系统阻尼和自然频率num=10;den=[114];H=tf(num,den);damp(H)14.stepnum=100;den=[150];sys=tf(num,den);sysc=feedback(sys,1);step(sysc)由step
3、函数生成的图形,可以右击直接得到各种动态性能指标,如图1所示图1num=100;i=0;forsigma=0:0.25:1.25den=[12*sigma*10100];damp(den)sys=tf(num,den);i=i+1;step(sys,2)holdon;endgridholdoff;title('阻尼比不同时的阶跃响应曲线')lab1='ζ=0';text(0.3,1.9,lab1),lab2='ζ=0.25';text(0.3,1.5,lab2),lab3='ζ=0.5';text(0.3,1.2,lab3),lab4='ζ=0.75';text
4、(0.3,1.05,lab4),lab5='ζ=1';text(0.35,0.9,lab5),lab6='ζ=1.25';text(0.35,0.8,lab6)输出结果如图2所示图2可见当一定时,系统随着阻尼比ζ的增大,闭环极点的实部在s左半平面的位置更加远离原点。sigma=0.25;i=0;forwn=10:20:50num=wn^2;den=[12*sigma*wnwn^2];sys=tf(num,den);i=i+1;step(sys,2)holdon;grid;endholdoff;title('wn变化时的阶跃响应曲线')lab1='wn=10';t
5、ext(0.35,1.4,lab1),lab2='wn=30';text(0.12,1.3,lab2),lab3='wn=50';text(0.05,1.2,lab3)输出结果如图3所示图3可见,当ζ一定时,随着增大,系统响应加速,振荡频率增大,系统调整时间缩短,但是超调量没变化。1.lsim用于求出系统任意输入激励的响应num=10;den=[1210];G=tf(num,den);t=0:0.1:10;u=5+2*t+8*t.^2;lsim(G,u,t)gridon;2.impulse用于求系统的单位脉冲响应3.nyquistnum=0.5;den=[121
6、0.5];H=tf(num,den)nyquist(H)4.bode5.margin计算系统的增益裕度和相位裕度s=tf('s');k=1;G=k/(s*(0.5*s+1)*(0.1*s+1))margin(G)6.pzmap绘制系统的零极点图num=[21];den=[132];sys=tf(num,den)pzmap(sys)1.rolcus绘制系统的根轨迹图k=1;z=[];p=[0-1-2];sys=zpk(z,p,k)rlocus(sys),grid2.dcgain求连续系统的稳态增益。3.residue部分分式展开。4.linspace构造线性分布的
7、向量。5.roots求多项式根。6.limit求极限,可用于求解系统的稳态位置、速度与加速度误差系数。例:已知某稳定单位负反馈系统的闭环传递函数为。试求系统的稳定位置、速度与加速度误差系数Kp、Kv和Ka。clearP=[556];roots(P)symssphibGbKpKvKa;phib=5/(5*s^2+5*s+6);[Gb]=solve('5/(5*s^2+5*s+6)=Gb/(1+Gb)',Gb)Kp=limit(Gb,s,0,'right')Kv=limit(s*Gb,s,0,'right')Ka=limit(s^2*Gb,s,0,'right')7
8、.conv两个多项式的乘