:y=initial(sys,x0,t)x0作用同上(3)综合应用例1:某系统的传递函数为G(s)=5s+100/s4+2s3+16s2+18s+18,输入为u=1+exp(-t)*cos(5t),当系统的初始条件为零时,求出其时域响应解:>>num=[5100];>>den=[12161818];>>G=tf(num,den);>>t=0:0.04:4;>>u=1+exp(-t).*cos(5*t);>>Y=lsim(G,u,t);>>plot(t,Y,'o',t,u,'-')例2:已知系统的闭环传递函数为Φ(s)=16/s2+8ζs+16,其中ζ=0.707.求二阶系统的单位脉冲响应、单位阶跃响应和单位斜坡响应。解:>>zeta=0.707;>>num=[16];den=[18*zeta16];>>sys=tf(num,den);>>p=root(den)p=-2.8280+2.8289i-2.8280-2.8289i>>t=0:0.01:3;>>figure(1)>>impulse(sys,t);grid>>impulse(sys,t);grid>>step(sys,t);grid>>u=t;>>lsim(sys,u,t,0);grid(3)线性系统的根轨迹1)绘制零、极点分布图———命令格式:[p,z]=pzmap(sys)2)绘制根轨迹图例1:G(s)=k/s(s2+2s+2)(s2+6s+13)的根轨迹绘制解:>>num=[1];>>den=[182738260];>>rlocus(num,den)>>rlocfind(num,den)%确定增益及其相应位置的闭环极点***说明***执行rlocfind命令后,根轨迹图上会出现“+”提示符,通过鼠标将提示符移到根轨迹相应的位置,然后按回车键,所对应闭环极点及其对应的参数K就会在命令行中显示。4.控制系统的频域分析
全文预览
