大学交通运输工程学院《机械工程控制基础实验报告》上机题及答案解析

0第3题............................................11s2kk-forj=1:6gtext('k=0.1','linewidth',1.5,'fontsize',10)%用鼠标放置文本gtext('k=0.2','linewidth',1.5gtext('k=0.5','linewidth',1.5gtext('k=0.8','linewidth',1.5gtext('k=1.0','linewidth',1.5gtext('k=2.4','linewidth',1.523forj=1:6gtext('k=0.1','linewidth',1.5gtext('k=0.2','linewidth',1.5gtext('k=0.5','linewidth',1.5gtext('k=0.8','linewidth',1.5gtext('k=1.0','linewidth',1.5gtext('k=2.4','linewidth',1.5331sys=zpk([],[-1+3*i-finalvalue=polyval(n,0)/polyvalwhile(y(l)>0.98*finalvalue)&(y(l)<1.02*finalval45EQ\*jc3\*hps34\o\al(\s\up13(0),5)2s2650.2s2s(s50.2s2s(s2ess1=polyval(n1,0)/ess2=polyval(n2,0)/7杆铰接处无摩擦。假定系统中所有变量变化范围很小，则可以将该系统视ssssX(s)s10(s1s1s10(s10(s5ttitle('thestepresponseofthesysxlabel('time(second8xkc54m707my2x(t)2.452t92222sX(s)x2=0.707*ones(size(ttitle('ResponseofHangingxlabel('Time(second)');ylabel('inputx(t)andoutputy(t)EQ\*jc3\*hps34\o\al(\s\up17(0),5)EQ\*jc3\*hps34\o\al(\s\up17(K),s)s2s222n2nw=[0,logspace(-2,2,20forj=1:5sys=tf([wn*wn],[1,2*e(j)*K,并判断闭环系统的稳定性。K,并判断闭环系统的稳定性。EQ\*jc3\*hps34\o\al(\s\up17(16),s)Wc=interp1(mag1,W,0,'splineW_3db=interp1(mag1,W,-3,'splineW_90=interp1(pha1,W,-90,'spline第5题已知单位反馈系统开环传递函数5(s2)5(s2)title('NyquistDiagraKK221xlabel('Frequency/(rad/title('feedbackbodediatitle('feedbackbodediafinalvalue=polyval(num,0)/polyval(dewhile(y(l1)>0.98*finalvalue)&(y(l1)<1.02*finalvalwhile(y(l2)>0.95*finalvalue)&(y(l2)<1.05*finalval1+ccoo_ylabel('x_o(t)')title('StepResponseofGK(s)=k/(s+1)(s+2)(s+5)')1H(s)为单位反馈，对系统采用比例微分控制，比例系title('StepResponseofGK(s)=k/(s+1)(s+2)(s+5)')1title('StepResponseofGK(s)=k/(s+1)(s+2)(s+5)')K,试用MATLAB,试用MATLAB设计PID校正装置，使得4s－1。若要校正更大的相位角，则值要更大，故选β=10。于是相位滞后校正装置另一转角频deng=conv([1,0],conv([0title(['校正后系统的相位裕度=',num2str(PmEQ\*jc3\*hps34\o\al(\s\up18(K),5s)解：首先设计超前校正环节，并绘制校正前后的开环伯德[mag,phase,w]=bode(sys_oT=1/(wcg*sqrt(alpha));Tz=title(['校正后系统的相位裕度=',num2str(Pmstep(sys_old,'r',sys解：设计滞后校正环节，并绘制校正前后的开环伯德den=conv([1,0],conv([1wcg=spline(phase(1,:),w',dPm-18step(Gy_close,'--b',Gx_ctitle(['校正后系统的相位裕度=',num2str(Pmfori=find((mag1<=-19.8)&(mag1>=-w_1=w1(ii)fori=find((mag1<=0.2)&(mag1>=jj=input('enterindexfordesiredmians=0.67810.68000.68180.68370.68560.6875在上述数据组中选择数值308，得到转角频率之一ians=在上述数据组中选择数值808，得到转角频率之一sysc=tf(numc1,denc1)*tf(nu,设计要求：相位裕度Wc=13.5;%rad/seconddisp('Modeloforigdisp('Result:modeloflag-leadisp('Result:modelofnewopendisp('Thesystem