matlab上机练习5.docVIP

上机练习5 姓名 学号 设线性定常系统的传递函数为：，试用MATLAB建立系统的传递函数模型。num=[2,3,16,0,8]; den=[1,4,0,2,17,0]; sys=tf(num,den) Transfer function: 2 s^4 + 3 s^3 + 16 s^2 + 8 -------------------------- s^5 + 4 s^4 + 2 s^2 + 17 s 已知一系统的传递函数为：，试用MATLAB建立系统的零极点模型，并求其等效的分子分母多项式模型。 K=12; Z=[-5]; P=[-11,-2,-9]; SYS=zpk(Z,P,K) Zero/pole/gain: 12 (s+5) ------------------ (s+11) (s+9) (s+2) 已知一单位反馈系统的开环传递函数为：，试绘制该闭环系统的单位脉冲响应曲线，单位阶跃响应曲线，单位斜坡响应曲线和单位加速度响应曲线。 num=25; den=[1,3,25]; g=tf(num,den); sys=feedback(g,1) t=0:0.01:10; subplot(2,2,1);impulse (sys,t) subplot(2,2,2);step (sys,t) subplot(2,2,3);sys1=sys*tf(1,[1,0]);step(sys1,t) subplot(2,2,4);sys2=sys*tf(1,[1,0,0]);step(sys2,t) Transfer function: 25 -------------- s^2 + 3 s + 50 对典型二阶系统闭环传递函数，输入单位阶跃信号，且ωn = 1，ζ分别等于-0.1；-0.2；0；0.1；0.5；1；1.5。及ζ= 0.5，ωn 分别等于1；10；30；50的响应曲线。分析ζ、ωn 对响应曲线的影响。 t=0:0.01:10; subplot(7,1,1);wn=1;zeta=-0.1;num=wn^2;den=[1,2*zeta*wn,wn^2];closys=tf(num,den);step(closys,t) subplot(7,1,2);wn=1;zeta=-0.2;num=wn^2;den=[1,2*zeta*wn,wn^2];closys=tf(num,den);step(closys,t) subplot(7,1,3);wn=1;zeta=0;num=wn^2;den=[1,2*zeta*wn,wn^2];closys=tf(num,den);step(closys,t) subplot(7,1,4);wn=1;zeta=0.1;num=wn^2;den=[1,2*zeta*wn,wn^2];closys=tf(num,den);step(closys,t) subplot(7,1,5);wn=1;zeta=0.5;num=wn^2;den=[1,2*zeta*wn,wn^2];closys=tf(num,den);step(closys,t) subplot(7,1,6);wn=1;zeta=1;num=wn^2;den=[1,2*zeta*wn,wn^2];closys=tf(num,den);step(closys,t) subplot(7,1,7);wn=1;zeta=1.5;num=wn^2;den=[1,2*zeta*wn,wn^2];closys=tf(num,den);step(closys,t) t=0:0.01:10; subplot(4,1,1);wn=1;zeta=0.5;num=wn^2;den=[1,2*zeta*wn,wn^2];closys=tf(num,den);step(closys,t) subplot(4,1,2);wn=10;zeta=0.5;num=wn^2;den=[1,2*zeta*wn,wn^2];closys=tf(num,den);step(closys,t) subplot(4,1,3);wn=30;zeta=0.5;num=wn^2;den=[1,2*zeta*wn,wn^2];closys=tf(num,den);step(closys,t) subplot(4,1,4);wn=50;zeta=0.5;num=wn^2;den=[1,2*zeta*wn,wn^2];closys=tf(num,den);step(closys,t)