自动控制原理实验报告

自动控制原理实验报告姓名学号时间地点实验楼B院系专业实验一系统的数学模实验二控制系统的时域分析实验三控制系统的频域分析实验一系统的数学模、实验目的和任务1、 学会使用MATLA啲命令；2、 掌握MATLAB有关传递函数求取及其零、极点计算的函数。3、 掌握用MATLAB求取系统的数学模型二、 实验仪器、设备及材料1、 计算机2、 MATLAB^件三、 实验原理1、 MATLAB^件的使用2、 使用MATLAB软件在计算机上求取系统的传递函数四、 实验内容1、特征多项式的建立与特征根的求取在命令窗口依次运行下面命令，并记录各命令运行后结果>>p=[1,3,0,4];P=13 0 4>>r=roots(p)-3.3553+O.OOOOi0.1777+1.0773i0.1777-1.0773i>>p=poly(r)P=1.0000 3.0000 -0.0000 4.00002、求单位反馈系统的传递函数：在命令窗口依次运行下面命令，并记录各命令运行后结果>>numg=[1];deng=[500,0,0];>>numc=[1,1];denc=[1,2];>>[num1,den1]=series(numg,deng,numc,denc);>>[num,den]=cloop(num1,den1,-1)num=den=500>>printsys(num,den)num/den=10001500sA3+1000sA2+s+13、传递函数零、极点的求取在命令窗口依次运行下面命令，并记录各命令运行后结果>>num1=[6,0,1];den1=[1,3,3,1];>>z=roots(num1);>>p=roots(den1);>>n仁[1,1];n2=[1,2];d仁[1,2*i];d2=[1,-2*i];d3=[1,3];>>num2=conv(n1,n2)num2=1 3 2>>den2=conv(d1,conv(d2,d3))den2=1 3 4 12>>printsys(num2,den2)sA2+3s+2sA3+3sA2+4s+12>>num=conv(num1,den2);den=conv(den1,num2);>>printsys(num,den)6sA5+18sA4+25sA3+75sA2+4s+12sA5+6sA4+14sA3+16sA2+9s+2>>pzmap(num,den),title(极点一零点图')2概点一專点图(spuaosm)嗖<ajeu^e--2' >—-3 卫_5-2 ^1.5 -1RealAxis(seconds1)―1 ft-05 04、求反馈联接系统的传递函数:在命令窗口依次运行下面命令，并记录各命令运行后结果>>numg=[1];deng=[500,0,0];>>numh=[1,1];denh=[1,2];>>[num,den]=feedback(numg,deng,numh,denh)num=0012den=500>>printsys(num,den)500>>printsys(num,den)s+210001500sA3+1000sA2+s+15、自行利用MATLAB命令求取以下系统传递函数，并记录下结果。X(s) ►+>>g2=tf([1,2],[1,3]);>>gg仁feedback(g1,g2,+1);>>g3=tf(10,[1,1]);>>gg2=series(gg1,g3);>>g4=tf([5,0],[1,6,8]);>>sys=feedback(gg2,g4)sys=20sA3+180sA2+520s+480sA6+11sA5+43sA4+67sA3+118sA2+252s-32Continuous-timetransferfunction.思考题思考题1、 实验内容中的第5项系统的反馈是什么反馈?有正反馈和负反馈。2、 闭环系统的零、极点个数相等吗？各为多少？不等。>>pzmap(sys)21.5Pole-ZeroMapRealAxis(seconds)5

o.5515>>[p,z]=pzmap(sys)-4.4999+1.9707i-4.4999-1.9707i-2.9863+O.OOOOi0.4332+1.8760i0.4332-1.8760i0.1198+O.OOOOi-4.OOOO-3.OOOO-2.OOOO实验二控制系统的时域分析、实验目的和任务1、 掌握用MATLAB对系统进行时间响应分析；2、 掌握一阶惯性系统以及二阶系统的时间响应特征以及系统性能与系统参数之间的关系。二、 实验仪器、设备及材料1、 计算机2、 MATLAB软件三、 实验原理1、利用计算机对控制系统进行时域分析。2、使用MATLAB^件在计算机上对一、二阶系统进行时域分析。使用MATLAB求一阶惯性系统的单位阶跃响应曲线。系统传递函数如下:G(s)=G(s)=100.5s1在命令窗口依次运行下面命令，并记录各命令运行后结果>>t=[0:.5:5];>>y=1-exp(-2*t);>>plot(t,y,'r');>>axis([0501.1]);>>set(gca,'ytick',0:.1:1.1);

>>title(‘y-t)xp(-2t)');>>xlabel('t');>>ylabel('y(t)');>>gridy(t)=1-exp(-2t)y(t)=1-exp(-2t)若系统传递函数为:G（s）二弋自行编制在命令窗口运行命令,报s+1自行编制在命令窗口运行命令,报求其单位阶跃响应，并与上面的结果进行比较，把结论写入实验告>>y=1-exp(-1*t);>>plot(t,y,'r');>>axis([0501.1]);>>set(gca,'ytick',0:.1:1.1);>>title('y(t)=1-exp(-2t)');>>xlabel('t');>>ylabel('y(t)')>>title('y(t)=1-exp(-1t)');>>gridy(t)=1-exp(-1t)比较结果，可得结论： T值的大小反映系统的惯性。 T值小，惯性就小，响应速度快;响应速度慢。2•使用MATLAB求二阶系统的单位阶跃响应曲线。系统传递函数如下：T值大,2-'nG(s)2 2(其中：n=0.4s+2加nS+国n'从0变化到2)在命令窗口如下运行命令，，并记录各命令运行后结果。>>symssforzeta=[0:0.2:0.8,1:0.5:2]wn=0.4;wn=sym(num2str(wn));zet=sym(num2str(zeta));ifzeta==0figure(1)ezplot(ilaplace(wnA2/s/(sA2+wnA2)),[080]);gridontitle('\xi=0' )elseifzeta==1figure(2)ezplot(ilaplace(wnA2/s/(s+wn)A2),[080]);holdon;elsefigure(2)ezplot(ilaplace(wnA2/s/(sA2+2*zet*wn*s+w门人2)),[080]);holdon;endendgridon;title( 'xi:0.2,0.4,0.6,0.8,1.0,1.5,2.0' )axis([O8001.8])gtext( '0.4')gtext( '1.0')gtext( '2.0')^=0比较不同值系统的响应曲线，有何结论?当时，振荡越来越小都趋向于1O四、思考题一阶系统的时间常数T越大，系统阶跃响应越慢。2、二阶系统阻尼比对系统阶跃响应的影响。二阶系统阻尼比越大，系统阶跃响应越慢。实验三控制系统的频域分析一、 实验目的和任务1、 禾U用MATLAB画一阶和二阶控制系统的伯德图2、 利用MATLAB计算所给系统的相角裕量和幅值裕量二、 实验仪器、设备及材料0、计算机1、MATLAB软件三、 实验原理0、利用计算机对控制系统进行频域分析。1、使用MATLAB^件在计算机对系统进行频域分析。四、 实验步骤1、 使用MATLAB画一阶惯性系统的伯德图传递函数如下：1G(s)：4s+1C(s)>>num=1;C(s)>>den=[41];>>g=tf(num,den);>>bode(g,'r')□Sod«Di^granifffip)lallpntiumEs305£-EQioJ t&'1 thFrequencyiradZsec]54fswE富ELImior2、 使用MATLAB画二阶系统的伯德图。传递函数如下：2亠©nG(s)2 2s+2%ns+con其中：-0.8 =0.1,0.5,1,1.5,2>>w=[0,logspace(-2,2,200)];>>wn=0.8;>>forzeta=[0.10.511.52];G=tf([wn*wn],[wnT2*zeta*wnwn*wn]);bode(Gw);holdon;EndBodeDiagrammp)epwEcn罗5010210*110°101102Frequency(rad/s)so5493_l>1N(Bap)as佰匸亠-1B03、 使用MATLAB画如下系统的伯德图并计算出系统的相角裕量和幅值裕量。500(0.0167s1G(s)500(0.0167s1G(s)s(0.05s1)(0.0025s1)(0.001s1)>>num=500*[0.01671];>>den仁conv([10],[0.051]);>>den2=conv([0.00251],[0.0011]);>>den=conv(den1,den2);>>Go=tf(num,den);>>w=logspace(0,4,50);>>bode(Go,w)

1D0BodeDiagram500500090-1-Ifmp)epmEMN(map)as佰匸亠-27010°W1Frequency(rad/s)3&60JIJImargin(Go)mp)epnlEm啰BodeDiagramGm=17J1D0BodeDiagram500500090-1-Ifmp)epmEMN(map)as佰匸亠-27010°W1Frequency(rad/s)3&60JIJImargin(Go)mp)epnlEm啰BodeDiagramGm=17JdB(at587rad/s),Pm=45.5deg(at162r胡⑸100mep)asm匸亠Frequency(rad/s)>>[Gm,Pm,Wcg,Wcp]=margin(Go)Gm=7.1968Pm=