资源描述：
《《机械控制工程基础》实验报告2.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、《机械控制工程基础》实验报告二班级•学号201079250420姓名黄煌实验日期2012年10月30口任课教师李碧琼实验名称实验二基于MATLAB系统动态特验证型性分析%1.实验目的和要求：1•要求掌握在MATLAB中所给系统的单位脉冲响应，单位阶跃响应以及在任意输入下的响应2•要求掌握MATLAB控制系统工具箱中用于分析系统动态特性的函数，绘制二阶系统动态响应曲线%1.上机内容1•绘制出所给系统的单位脉冲响应曲线2•绘制出所给系统的单位阶跃响应曲线3•绘制出所给系统在任意输入下的响应曲线%1.实验平台Windows2000或XPMATLAB7.0以上版
2、本四.操作过程，源程序和计算结果:L=J源程序:t=[0:0.1:10];num=1；zetal=0.1；deni=[12*zetal1]；sysl=tf(num,deni)；zeta2=0.25；den2=[12*zeta21]；sys2=tf(num,den2)；zeta3=0.5；den3=[12*zeta31]；sys3=tf(num,den3)；zeta4=1•0；den4=[12*zeta41]；sys4=tf(num,den4)；impulse(sysl,sys2,sys3,sys4,t)；gridon；运行结果源程序:t=[0:0.1:1
3、0];num=1；zetal=0.1；deni=[12*zetal1]；sysl=tf(num,deni)；zeta2=0.25；den2=[12*zeta21]；sys2=tf(num,den2)；zeta3=0.5；den3=[12*zeta31]；sys3=tf(num,den3)；zeta4=1•0；den4=[12*zeta41]；sys4=tf(num,den4)；impulse(sysl,sys2,sys3,sys4,t)；gridon；运行结果题6.7源程序t=[0:0.1:12];num=[1]；zetal=0.I；deni=[l2*z
4、etall];sysl=tf(num,deni)；zeta2=0.2；den2=[l2*zeta2l];sys2=tf{num,den2)；zeta3=0.4；den3=[l2*zeta3l]；sys3=tf(num,den3)；zeta4=0.7；den4=[l2*zeta4l]；sys4=tf(num,den4)；zeta5=l•0；den5=[I2*zeta5I]；sys5=tf(num,den5)；zeta6=2.0；den6=[l2*zeta6l];sys6=tf(num,den6)；step(sysl,sys2,sys3,sys4,sys5,
5、sys6,t)；title(1zeta=0.10・20.40.7l・02・0°)gridon；运行结果丿FigureFileEditViewInsertToolsDesktopWindowHelp除致®礬®渥口匡
6、El"CM0.20.40.71.02.024681012Time(sec)6・o4.o1.8o.opnu-d$o.2oo8642••••1111题6.8源程序t=[0:0.01:5];numl=1；num2=3；den=[10.21]；sysl=tf(numl,den)；sys2=tf(num2,den)；u=sin(t)；lsim(sysl
7、,sys2,uzt)；gridon；运行结果-Inlxi•丿FigoreFileEditViewInsertToolsDesktopWindowHelp除®®卩®渥口匡1QTime(sec)题6.9源程序num=[0.1]；den=[10.240.1]；sys=tf(num,den)；[wn,zzp]=damp(sys)运行结果wn=0.31620.3162z=0.37950.3795P=-0.1200+0.2926i-0.1200-0.29261题3.3k=100；z=[-4];p=[0-0.5-50-50];[num,den]=zp2tf(zzp,k
8、)；bode(num,den)；title(1BodePlot')gridon运行结果（截图如下）0〔邕(DPWU豊-200BodePlot•90-135-180-225-27010°102103io"2io'1io110*Figure1FileEditViewInsertToolsDesktopWindowHelp-50-100-150Frequency(rad/sec)题3・4wn=6；kosi=[0.1:0.1:1.0]；w=logspace(-1,1,100)；figure(1)num=[wnA2]；forkos=kosiden=[1,2*kos
9、*wn,wnA2]；[mag,pha,wl]=bode(num,den,w)；s
