机械工程控制基础matlab报告

PAGEPAGE11《机械工程控制基础》MATLAB编程实例班级183234学号2114793姓名吴洲试验一：用MATLAB进行部分分式展开1、试验目的：（1）对MATLAB进行初步的了解；（2）掌握应用MATLAB对高阶函数进行部分分式的展开。2、试验学时：2学时3、试验方法：MATLAB有一个命令用于求B（S）/A（S）的部分分式展开式。设S的有理分式为F（S）=B（S）/A(S)=num/den=(b0Sn+b1Sn-1+…+bn)/(Sn+a1Sn-1+…+an)式中ai和bi的某些值可能是零。在MATLAB的行向量中，num和den分别表示F(S)分子和分母的系数，即num=[b0b1…bn]den=[1a1…an]命令[r,p,q]=residue(num,den)MATLAB将按下式给出F（S）部分分式展开式中的留数、极点和余项：有：r(1)、r(2)…r(n)是函数的留数；p(1)、p（2）、…p(n)是函数的极点；K（s）是函数的余项。4、试验内容（1）试求下列函数的部分分式的展开式，写出程序语句和部分分式的结果程序语句：>>num=[111395226];>>den=[110355024];>>[r,p,q]=residue(num,den)结果：r=1.00002.5000-3.00000.5000p=-4.0000-3.0000-2.0000-1.0000q=1（2）试求下列函数的部分分式的展开式，写出程序语句和部分分式的结果程序语句：>>num=[146];>>den=[1331];>>[r,p,q]=residue(num,den)结果r=1.00002.00003.0000p=-1.0000-1.0000-1.0000q=[]试验二：时域特性的计算机辅助分析1、试验目的：（1）对MATLAB进行进一步的了解；（2）掌握应用MATLAB求解高阶系统的各种时域响应的求法；（3）掌握系统零极点的求法，从而判断系统的稳定性。2、试验学时：2学时3、试验方法：（1）impluse(num,den)，得到传递函数G(S)=num(S)/den(s)表示的脉冲响应曲线；（2）impluse(num,den,t)，得到传递函数G(S)=num(S)/den(s)表示的脉冲响应曲线，并可指定时间t；（3）step(num,den)，得到以得到传递函数G(S)=num(S)/den(s)表示的脉冲响应曲线；（4）step(num,den,t)，得到以得到传递函数G(S)=num(S)/den(s)表示的脉冲响应曲线,并可指定时间t；(5)lsim(num,den,u,t)计算连续系统的任意输入响应,得到以得到传递函数G(S)=num(S)/den(s)表示的u输入的响应曲线,并可指定时间t；（6）roots(den)，得到多项式den的根；（7）pzmap(num,den)，得到传递函数G(S)=num(S)/den(s)表示的系统的零极点图，“×”表示极点，“O”表示零点。4、试验内容（1）已知某高阶系统的传递函数为求该系统的单位脉冲响应、单位阶跃响应和单位速度响应和单位加速度响应。获得单位脉冲响应程序语句及图形：>>num=[22050];>>den=[11584223309240100];>>impulse(num,den)获得单位阶跃响应程序语句及图形：>>num=[22050];>>den=[11584223309240100];>>step(num,den)获得单位速度响应程序语句及图形：>>num=[22050];>>den=[11584223309240100];>>t=[0:0.01:1];>>u=(t);>>lsim(num,den,u,t)获得单位加速度响应程序语句及图形：>>num=[22050];>>den=[11584223309240100];>>t=[0:0.01:1];>>u=(0.5*t.*t);>>lsim(num,den,u,t)（2）已知某高阶系统的传递函数为试求该系统的极点并判断系统的稳定性。程序及结果：>>den=[123456789];>>roots(den)ans=-1.2888+0.4477i-1.2888-0.4477i-0.7244+1.1370i-0.7244-1.1370i0.1364+1.3050i0.1364-1.3050i0.8767+0.8814i0.8767-0.8814i（3）已知某高阶系统的传递函数为试求该系统的零极点并判断系统的稳定性，绘制系统的零极点图。程序及结果：>>num=[3120006250000];>>den=[110080004400006240000];>>pzmap(num,den)试验三：频域特性的计算机辅助分析1、试验目的：（1）对MATLAB进行加深了解；（2）掌握应用MATLAB求解系统的Nyquist图；（3）掌握应用MATLAB求解系统的Bode图；2、试验学时：2学时3、试验方法：（1）nyquist(num,den,w)，绘制Nyquist图。num和den分别表示传递函数的分子和分母中中包含以S的降幂排列的多项式系数。命令Nyquist(num,den,w)可以按指定的频率点W，绘制系统的Nyquist图。（2）bode(num,den,w)，绘制Bode图。num和den分别表示传递函数的分子和分母中中包含以S的降幂排列的多项式系数。命令bode(num,den,w)可以按指定的频率点W，绘制系统的bode图。4、试验内容（1）考虑下列系