matlab校正系统程序

《matlab校正系统程序》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、clear;Gl=tf(20,[1ZO.4,0])%系统的开环传递函数G=feedback(Gl,lf-l)殆原系统的闭环传递函数figure(1)step(G)%原系统的单位阶跃响应figure(2)rlocus(G1)$原系统的根轨迹figure(3)margin(G1)$原系统的bode图[GmO,PmO,WgO,WcO]=margin(G1)deta=5%补偿角Phim=45-Pm0+deta匕校正装置最大超前角a=sin(Phim*pi/180)alpha=(1+a)/(1-a)Wm=W

2、cO*alpha/s(1/4)T=1/(Wm*sqrt(alpha));alpha_T=alpha*TGc=tf([alpha_Tzl]z[T,l])%校正装置Ga=Gl*Gc%校正后的传递函数Gll=feedback(Glf1);$原系统闭环传递函数Gaa=feedback(Ga,1);$校正后系统闭环传递函数figure(4)rlocus(Ga)%校正后系统的根轨迹legend('Ga')figure(5)margin(Ga)衛校正后系统的频域特性[Gm,Pm,Wcm,Wpm]=margin(

3、Ga)figure(6)step(Gil,'-.k'zGaa)%原系统与校正后系统的时域特性的比较legend('Gl1,'Ga')figure(7)bode(Glz'-.k'zGc,1:Ga)%原系统与校正后系统的频域特性的比较legend('G1'Gc)G1=Continuous-timetransferfunction•sA2+0.4s+20Continuous-timetransferfunction•GmO=InfPmO=5.1214WgO=InfWcO=4.4631deta=5P

4、him=44.87860.7056alpha=5.7936Wm=6.9243alpha_T=0.3476Gc=0.3476s+10.06s+1function.0•4sfunction•Continuous-timetransferGa=6.952s+200.06sA3+1.024s入2+Continuous-timetransferGm=InfE6・9HUIdMJUIHUIOMSZ(spuooes)eluH5Z0Aspuooes)8sfiLoinleAleuizomelsxsL6L(spuogs

5、)oE-?-£asoWSSASTT7h~cooBsuodsaad£s0945089-0976RootLocus068030.05谕-0.05(「spuooas)曲xvAJeuoew-l0994-0994r0976094503503089.45114-0.4-035015005OJ10680503-0.3-025-02-0.15-0.1RealAxis(seconds1)00500.05(8P)epnuu6esBodeDiagramGm=InfdB(atInfrad⑸，Pm=5.12deg(at4.

6、46rad/s)351•(Bap)aseqd10110°Frequency(racVs)101校正后T""oSoRea-Axis(seconds.s_epResponseSystemGaseeing5®(seconds)090806Time(seconds)>SystemGa58Ume(seconds^0」62SystemGaF-EvaErcx'1100500510090(8P)0-90BodeDiagramGm=InfdB(atInfrad/s),Pm=48.2deg(at6.94rad/s)1

7、0'1101100500-50100-90(8P)9pn_c6es35(Bop)heqd103Frequency(rad/s)