《机电系统动力学仿真matlab课后答案》由会员分享，可在线阅读，更多相关《机电系统动力学仿真matlab课后答案（32页珍藏版）》请在金锄头文库上搜索。
1、机电系统动力学仿真matlab课 后答案-刘白雁机电系统动力学仿真实验报告机工 1004班 201003130150宋康习题一13Abe wu_2004两个是合法的。1.4(1)x Command Window;-(12+2iU-4nZ(32)H启=2(2)Wiats NewIdersholders4_|ort byShawr siactionsthumlhowShowfunction/蠶zes懿矗;I?回 Howto Add 0Direcloiy H證M ：0 圍 E- TXCommand WindowA A= 13;4, 5-6 ;7；_.8, 9A =1234567S9elear;x=-
2、8:0.5:8;y=x;X=ones(size(y)*x;Y=y*ones(size(x);-5R=sqrt(X.A2+Y.A2)+eps;Z=sin(R)./R;mesh(X,Y,Z);colormap(hot)xlabel(x),ylabel(y),zlabel(z)0.5-0.5 W10-51.7(1)J WhartsNew0XCommand Window： -l=；2xsinC0. 3*piVCl+sqrt；(5)；yi 二0.5000r Sthuml4 y 1=2*5111(0. 3*pi) f： (1+sqrt (5卜)O.EOOO y2=2*eos.(0.(1+sqtt )淀=0
3、.36331.110.5-0.5.-10-5-5习题二2.12.4 k x Comm-and Window肚厲玄乱4:0=5.6;7, 8: C=A*Bc=A.*B.C.-19 a 4350512汶1323i0athum323435a050a0ao5nCommand Window厂SIarts =ans A=EX 乙” 1爲 5, 班 0, 5,3;1, 5,0,5: MB, A|B, Aj A=B, AB:ELTLS 二:TL 二习题三3.1clf9clc,cieart=(0:18/100:18);xi=0.2:0.2:0.8;beita=sqrt(1-xi.A2);sita=atan(be
4、ita./xi);y=1-exp(-t*xi).*sin(t*beita+ones(101,1)*sita)./(ones(101,1)*beita);figure(1)plot(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)g=0.2g=0.4幣g=0.63.2clf,clearx=-50:1:50;y=-4O:1:60;z=1/(sqrt(1-x) 八2+y.A2)+sqrt(1+x)八2+八2); figure(1)plot3(x,y,z),box on X,Y=me
5、shgrid(x,y);Z=1/(sqrt(1-X) .A2+Y.A2)+sqrt(1+X)八2+丫八2); figuremesh(X,Y,Z) figure surf(Z)0.08 k0.06-0.04-0.020.5-r-0.2-0.1000.5-0.4 -P0.3 -0.2-0.1 -3.3clc; clear all; close all; t=0:pi/50:2*pi; x=sin(t); y=cos(t);z=t; plot3(x,y,z,g-)8642010.510.500.-0.5-0.5习题四4.1 冥 Command Windowx ICommandWindowSI jrr
6、il麴 i=0；i=0； labile i=63/k二k+2i； i=i+l;endkk =1；：k=0 :for i=0;方工k=k+2； I；endk =1.8447e+019Command Window 1= W： U63FJ.;K=2. 1；K=sum 徑，1)Kl = sm C)1.044fe-FO 191.044re-bO 19x习题五5.1 a二1 35 291 1093 1700; b=l 289 254 2541 468-4 1700: c=tf b)Transfer furic丈丄口门；s4 + 35 事3 + 291 s2 + 1093 s + 1700+ 289 s4
7、+ 254 s3 + 2541 s2 + 4684 s + 1F0Q z=-3; p= -1 5 15； k=15; sys=zpk Pj k)Zero/pole/gain:IE (s+3)(s+1)(s+5)(s+15) k二lOCi; =0,-2,-2; P=-l 1；sys l=zpk (zj pj k); a=l 3 2; b=l 2 5 2; sys2=tf (a3 b); sys3=series(syslj sys2)Zero/pole/gain:100 e Ce+2)3 (e+O(s+O (e+0. 46680 (e-1)+ 1.533e + 4. 284)7、 I5.4 A=3
8、 2 1 ;0 4 6;0 -3 3; B=l 2 3厂; C=l 2 5; D=0; Eye=SE (A3 E3 C, D); SYS=tf(sys)Transfer function:20 52 - 123 s + 234sa3 - 10 sa2 + 51 s - 90 sys2=zpk(sysOZero/pole/gain:20- 6. 15s + 11. 7)(s-3)伝2- 7s + 30)5.5 nun=l 15 50 500 deri=l 2 0: 3ys=tf (derij nun)Transfer fund ion:s2 十 2 e+ 15 s2 + 50 s + 500 3
9、= 1 3 6 4; b= 4 0; sysl=t (b, a)Transf er funci ion:4 as3 + 3 s2 + 6 s + 45.6a =KlKJk4si0100k2-?-n.nn.0001k-40 . 089290-0. 125ul k10k20k30x4 0. 1786c =xl k2 k3 k4yl U U 1 Ud =ulyl 0C cnt inuLou 出一七 imu mo d_e 1.Tr：ai_Lsf er f ui-LCt ion:LL 1786 s_2 + U. 0 2551 s 十 U.3 571s 4十匚i2679 s 3十2匚山厅s2十02厅s5.
10、7 ml=12; m2=38;k=1000;c=0.1; sysl=tf(cj kj ml c kTransfer function:0.1 s + 1000122 + Cl. 1 s + 1000 ys2=tf (12 . 1 1000, 12*1112, 0. 1*1+112) 1000*J 0, 0Transfer fun匚tion:122 十 0. 1 s 十 1DDD456 *4 + 5 s3 + 50000 $2 sys3=series(sysl sys2)Transfer functzltirL：1. 2 e： + 1. 2e004十 200 e + leOOfi5472日花 十
11、105. 6 srt5十l.D罚罚血sn4十10000日“3十5e007事2l习题六6.2Zero/pole/gam:20 s伝+11 5)s (s+0. 1)(s+2)(s+10)6.4Figure 1File Edit /iew Insert Tools Desktop Window Help S订題&貯曷直币韵苗口Time (sec)1.s.64 Q Q QQJpn 七-d Lu匚 oo8 c4 2 o DE3 Fgure 1File Edit View Insert Tools Desktop Window HelpauhD矗|典敢7| | 口目| 二I*吕一dLW习题七7.1（1）a20-2-60101010Nyquist 匚 iagram-1-00-270Real AsiaFrequency (rad/sec)临界稳定Bo-dls DiagramGm = 0 dB at 1 rad/see), Pm = 0 deg (At 1 rsd/sac)50.S Figure 1丨=| 回 |BelFile EditViewInsert ToolsDesktop WindowHelpQ? JS E
