机电系统动力学仿真matlab课后答案 刘白雁

(2(2)(2(2)机电系统动力学仿真实验报告习题一机工1004班201003130150宋康习题一1.3Abewu_2004两个是合法的。1.4⑴XCommandWindow>>(12+2^(.7-4)):辺呂=2>>(2)?回HowtoAdd[2What'sS[S'艮丫IdersS[S'艮丫Idersholdersortbj;Shaw.-rsi;nt-5.mactionsthumlShownow,.Jfunction/2es._會(3)clear;x=-8:0.5:8;y=x';X=ones(size(y))*x;Y=y*ones(size(x));Directoiy-HAM...XCommandWindow»A=[1,.^3:4,5-'68,9]A=1234567S9■»R=sqrt(X.A2+Y.A2)+eps;Z=sin(R)./R;mesh(X,Y,Z);colormap(hot)xlabel('x'),ylabel('y'),zlabel('z')1.7(1)JWhart's^NewXCommandWindow：>>^l=；2'xsinC0.3*piV-Cl+sQrt：.C5))；yi二0.5000rSIthuml4window寧》y1=2*5111（0.3*pi）/（1+sqrt⑸）■71=O.EOOO>>y2=2^cos.C0.1+sqtt⑸）y2=0^30331.11习题二2.1».-k=1!0..5<2K=1.00001.-2000»y=2':0J：Zi1y-二EmptymatrixK=1.00001.-2000»y=2':0J：Zi1y-二Emptymatrix:l~by-0»1.40001.60001.80002.00002.2»3e=lLH5pA£e[£l7D.l^reE=GolunnEIthicojch1.60fl.12BS□,2£flE0.3Q47□.E]£90.$4L]n?caaD.897SL.02SQ].L54LL2S23].4]0G1.53BTL.甌M3.T9E2L9254ColuriiuJTthrnu^h322.051.72,17932.300L2bA363£.2.6325&B2]02rSH953,07?5■3.EUSr3.333P3.^6223.53Ma.fl跖乳3751C£>1lUi»£31tkrou^.h4Q■I.L临du2315d.^LQSD4.6162乩7WILB7275uODQ9乩12615.257-i5.3Q565.5j395■创20氐77D35u3M5乩0K7ColunnE4.9tkTDU^JisnE.L5E0L2832v讣u=O!：2hpi/43：2^pi5E=0fl-12620.25d50fl-12620.25d5D.saa?0.6329CLfiJH.0.泗AD.B97fldeluruuL『3£2.Q5Lr2-]?£O2.303]2.43632.G^3g2-692S2-D2JD3.舸93SLumu33thlCQJitia1D334.23LEL359S4.4SQO4.6362L7J444.B72?5.DQ09C^l-unais->19tLzoufh5D6ul55D6u2B32Lthe日口ghi咅LOSEB].L54]3.2323L.430EL.E3B7LflfifQ1-^9621.92313-OF7E3.2DE73.S3393.<10223.5904a.?L8fi3-646&3.9?6i区璋谄5-129J5.25745.38血5.G13®E.6420G.?7035-B9S5fi.oaer2.4Comm-andWindow»[1,2;:0=[5.6;7,8]:C=A*Bc二5...?XCommandWindow»A=[X.2.^4;1,.^5,0]:■-B-ll,h酸礬谶:MB,A|B,扎卫二二B,A>Bk:iI_LS=0011SIiioakI笺袞丄I习题三3.1clf9clc,cieart=(0:18/100:18)';xi=0.2:0.2:0.8;beita=sqrt(1-xi.A2);sita=atan(beita./xi);y=1-exp(-t*xi).*sin(t*beita+ones(101,1)*sita)./(ones(101,1)*beita);figure(1)piot(t,y(:,1),'b:',t,y(:,2),'r-.',t,y(:,3),'g-',t,y(:,4),'ko')legend('\xi=0.2','\xi=0.4','\xi=0.6','\xi=0.8')3.2clf,clearx=-50:1:50;y=-40:1:60;z=1./(sqrt((1-x)•八2+y.A2)+sqrt((1+x)・八2+『・八2));figure(1)plot3(x,y,z),boxon[X,Y]=meshgrid(x,y);Z=1・/(sqrt((1-X).A2+Y.A2)+sqrt((1+X)・八2+丫・八2));figure(2)mesh(X,Y,Z)figure(3)surf(Z)10050501005050100505010050505050■000.50.2-50■000.50.2-50-500.5-0.5-0.5-0.5-0.4-0.3-0.2-0.1-0.4-0.3-0.2-0.1-1500.150003.3clc;clearall;closeall;t=0:pi/50:2*pi;x=sin(t);y=cos(t);z=t;plot3(x,y,z,'g-')8-8-8-8-习题四4.1XCommandWindow总i=0；i=0；A.virilei<=63/k二k+2"i；i=i+l;endksk=jrrill<-S44Fe+0191>>■—xICommandWindow>；^：k=0:fork=k+2；'''I；endxxCommandWindowxxCommandWindowX1.044fe-FO19X1.044fe-FO19»i=(o：-：i■■jg'a)?；^=2.K=sum聶逹)Kl^srnCOKI=1.044re-bO19»习题五5.1>>a二[13529110931700];>>b=[l2892542541468-41700]:>>c=tfb)Transferfuric丈丄口门；s4+35事3+291s"2+1093s+1700s^B+289s^4+254s^3+2541s"2+4684s+1F0Q>>z=-3;>>p=[-1—5—15]；>>k=15;>>sys=zpkPjk)Zero/pole/gain:IE(s+3)(s+1)(s+5)(s+15)>>>>k二lOCi;>>==[0,-2,-2];»P=[-l1]；sysl=zpk(zjpjk);>>a=[l32];>>b=[l252];>>sys2=tf(a3b);>>sys3=series(sysljsys2)Zero/pole/gain:100eCe+2)"3(e+O(s+O(e+0.46680(e-1)+1.533e+4.284)7、I5.4»A=[321;046;0-33];»B=[l23厂;»C=[l25];»D=0;>>Eye=SE(A3E3C,D);>>SYS=tf(sys)Transfer■function:205^2-123s+234sa3-10sa2+51s-90>>sys2=zpk(sysOZero/pole/gain:20-6.15s+11.7)(s-3)伝"2-7s+30)»5.5>>nun=[l1550500]>>deri=[l20]:>>3ys=tf(derijnun)Transferfundion:s"2十2e+15s^2+50s+500>>3=[1364];>>b=[40];>>sysl=t£(b,a)Transferfunciion:4as^3+3s"2+6s+45.6a=KlKJk4si0100k2-?-n.nn.0001k-40□.089290-0.125TOC\o"1-5"\h\zulk10k20k30x40.1786c=xlk2k3k4ylUU1Ud=ulyl0Cc»ntinuLou出一七imumod_e1.Tr：ai_Lsferfui-LCtion:LL1786s_'2+U.02551s十U.3571s■'■4十匚i・2679s■'■3十2■匚山厅s"2十0・2厅s5.7»ml=12;»m2=38;k=1000;c=0.1;>>sysl=tf([cjk]j[mlck]}Transferfunction:0.1s+100012『2+Cl.1s+1000»£ys2=tf([12□.11000],[12*1112,0.1*^1+112)^1000*J0,0]^Transferfun匚tion:12『2十0.1s十1DDD456护4+5s"3+50000>>sys3=series(sysl^sys2)TransferfunctzltirL：1.2e：+1.2e004十200e+leOOfi5472日花十105.6srt5十l.D罚罚血sn4十10000日“3十5e007事2»l习题六6.26.3三］三］Figure1□回Ma-/Zero/pole/gam:20s伝+115)s(s+0.1)(s+2)(s+10)>>6.4FileEdit\/iewInsgrtToolsDesktopWindowHelp□Sfe題&貯曷«UkJh□OJpn七-dLUVStepResponse-2040Time(.sec)s6QQ4.2O&o.ap_Ty-CLuuv6042OJpn七-dLUVStepResponse-2040Time(.sec)s6QQ4.2O&o.ap_Ty-CLuuv6042oDEJFigure1〔=I回FileEditViewInsertToolsDesktopWindowHelp□ua£典Qc®4c□3B□StspResponse1-Bjiiiiin1.6一习题七7.1(1)Figure1FileEditViewInsertToolsDesktopWindowHelp□HSfe題目的倒退□Bb□Gm=050Bo-dsdiagram(at1rad/sec),Pm=0deg(st1rsd/sae)Nyquist匚iagramo5»fl〔暑)IbwELIdGm=050Bo-dsdiagram(at1rad/sec),Pm=0deg(st1rsd/sae)Nyquist匚iagramo5»fl〔暑)IbwELId-x10°Jn—1Q15642-2RealRealAsisFrequency(rad/sec)临界稳定(2)Figure1Figure1*2—Figure1Figure1*2—FileEditVievwInsertToolsDesktopWindowHelpHa®^5□Es□OpenFileGm=-19.2dS(a100BDoeDiagram0.31Srad/sac>,Pm=-34.3deg(atO.Gfifirad/sec)2000.-NyquistCiagrammp)jjpn七UBE2ooooooo

o5T—00

o02101Q£-22(BOJP)2W时工dscVP匚」3T—w_xv亡闵匚_m旳Lu-nu-1000-15Q0-270iaFrequency(radJsec)不稳定(3)-2000」1-150--1GO-500RealAxisFigureFigure1日回|-£WFileEditVievi.i1InsertToolsDesktopWindovu'Help□HSfe汶泳曹@霰□SB□Cm1QQ500-50-1QQ-15Q-2GG-18Q-225-270-3-15-360BodeDiagramCm1QQ500-50-1QQ-15Q-2GG-18Q-225-270-3-15-360BodeDiagram=In1,Pm=-461deg(at2.SErad/sec)NyquistDiagram15-1QQ-5010E匚030o3一XVLm匸一口mE_RealAx.is不稳定(4)

FileEditViewInsertToolsDesktopWindowHeIpDAkQC@□Sh□BodeDiagramInf,Pm=-33.6deg(atO.SS2rad/sec)3O-QQajpnl!LlmE左ooooo(GwpoQcfl閒L-d10-am-6-1CCC20-Q010-Q0-20-QQNyquistDiagram-4-202RealAxis.oFrequency(rad/secj不稳定7.2K=36OOFigure1IG凹b—OOFigure1IG凹b—FileEditViewInsertToolsDesktopWindowHelpDSa|feC®®□目画0BodeDiagram50・50-1000(署)畫工QQQFigure1|=■丨回QQFigure1|=■丨回FileEditViewInsertToolsDesktopWindowHelp*DQak汶quB□LinearSimulationRsbuhbT—S.5Mpnll-duJV3a2520Zern/pale/gain:LinearSimulationRsbuhbT—S.5Mpnll-duJV3a252020s(s+O.5〕5(s+0.1)(s+2)(s+10)FileEditViexyInsertToolsDesktopWindowHelp□QAteQC®00Yd5-tepResponseBpn-t±_clsLinearSimulationResults5-tepResponseBpn-t±_clsLinearSimulationResults7.3nyure丄nyure丄22f-(dEd_f<-evwInsertTQ0-5DeskfopXAZ-ndowHe-PAmplitude吕bjuseTIRewpon器TimeCD7.4FileEditViewInsertTaalsDesktopWindowHelp□HSk⑥晏DinVn〔邕oBodeDiagram5aa〔詈jOJSBLIQ.101匚m3T-j-10〔邕oBodeDiagram5aa〔詈jOJSBLIQ.101匚m3T-j-10-io-1ioaFrequency(rad/sec)Filfi上ditViftwr[nsfirT丄ccIf;UAskrnpWinrinwrHrIp□Ur电氏零⑥T4TT□冒「-□MyquisiDiegram100□oo55SIXy/JBU-mELU--21.4StepResponse■1.2453Q205FileEditViewInsertToolsDesktopWindQVii'Help□话日彎fe愆色鋼®晏ns■口.s勇a.a'.gpn七-duuvo.2oa.oa453riyure峠FileEditViewInsertToolsDesktopWindowIHelp□QO氏愆色侨⑥运□0■口BodeDiagramSooo42mp*pnc^w1■11■111I111111111111111I111■111111111111111I111■111111111111111I111I11111111111LJ-——L-_l-l-l111——1——1--1-1-11L1111_IIIII11111111111111111111111111I--1-—X_l111111111111r-r-1-t-|-tT-r*-_^r