[工学]现代机械控制工程实验报告

《[工学]现代机械控制工程实验报告》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、华东理工大学现代机械控制工程仿真实验实验报告学院________________班级________________姓名________________学号________________2010年12月现代机械控制工程仿真实验实验一1.比例环节：G(s)=kk=12num=[12];num=[12];den=[1];den=[1];t=[0:0.1:10];t=[0:0.1:10];y=impluse(num,den,t);y=step(num,den,t);plot(t,y);plot(t,y);gridongridon积分环节：G(s)=1/sn

2、um=[1];num=[1];den=[1,0];den=[1,0];t=[0:0.1:10];t=[0:0.1:10];y=impluse(num,den,t);y=step(num,den,t);plot(t,y);plot(t,y);gridongridon微分环节：G(s)=snum=[1,0];den=[1]num=[1,0];num=[1,0];den=[1];den=[1];t=[0:0.1:10];t=[0:0.1:10];y=impluse(num,den,t);y=step(num,den,t);plot(t,y);plot(t,y

3、);gridongridon不存在惯性环节：G(s)=1/(Ts+1)T取2num=[1];num=[1];den=[2,1];den=[2,1];t=[0:0.1:10];t=[0:0.1:10];y=impluse(num,den,t);y=step(num,den,t);plot(t,y);plot(t,y);gridongridon二阶振荡系统：G(s)=50/((0.05s^2+(1+50T)s+50)T取1num=[50];num=[50];den=[0.05,51,50];den=[0.05,51,50];t=[0:0.1:10];t=[

4、0:0.1:10];y=impluse(num,den,t);y=step(num,den,t);plot(t,y);plot(t,y);gridongridonSimulink练习：实验二例1系统模型如下所示，判断系统的稳定性，以及系统是否为最小相位系统。num=[3,16,41,28];den=[1,14,110,528,1494,2117,112];G=tf(num,den);sys=zpk(G);pzmap(sys)clc;clear;closeall;num=[20];den=[1836400];[numc,denc]=cloop(num,d

5、en);t=0:0.1:10;y=step(numc,denc,t);[y1,x,t1]=step(numc,denc);plot(t,y,'r:',t1,y1);title('thestepresponce');xlabel('time-sec');disp(‘稳态值dc为');dc=dcgain(numc,denc)clc;clear;closeall;num=[20];den=[1836400];[numc,denc]=cloop(num,den);t=0:0.1:10;y=impulse(numc,denc,t);[y1,x,t1]=impul

6、se(numc,denc);plot(t,y,'r:',t1,y1);title('theimpulseresponce');xlabel('time-sec');disp(‘稳态值dc为');dc=dcgain(numc,denc)clc;clear;closeall;num=[25];den=[1625];t=0:0.1:10;y=step(num,den,t);[y1,x,t1]=step(num,den);plot(t,y,'r:',t1,y1);title('thestepresponce');xlabel('time-sec');disp(

7、‘稳态值dc为');dc=dcgain(num,den)clc;clear;closeall;num=[1025];den=[0.161.961025];[numc,denc]=cloop(num,den);t=0:0.1:10;y=step(numc,denc,t);[y1,x,t1]=step(numc,denc);plot(t,y,'r:',t1,y1);title('thestepresponce');xlabel('time-sec');disp(‘稳态值dc为');dc=dcgain(numc,denc)clearclccloseall%输

8、入期望的超调量及峰值时间pos=40('inputpleaseinputexpectpos(%