matlab模擬(word文檔良心出品)

1、MATLAB模擬太陽系實(shí)驗(yàn)報告小組成員：鄒智康 U201313766 王思懋 U201313475 陳桃斌 U201313531楊飛U201313530陳鐵U201313522班級：信卓 1301指導(dǎo)老師：文灝一、 實(shí)驗(yàn)?zāi)康?、熟悉 MATLAB 的基本功能，了解 MATLAB 的編譯環(huán)境。2、掌握 MATLAB 的基本繪圖函數(shù)，熟悉繪圖函數(shù)的用法、簡單圖形標(biāo)注、簡單顏色鑒定。3、通過組內(nèi)分工、協(xié)商合作，增強(qiáng)組員自主學(xué)習(xí)的能力，加強(qiáng)組員間團(tuán)結(jié)合作的能力。二、 預(yù)期目標(biāo)1、繪制太陽系三維動態(tài)圖形。2、模擬八大行星和彗星繞太陽運(yùn)行動態(tài)軌道。3、模擬月球繞地球運(yùn)動.小星星帶及土星環(huán)。4、模擬出大致的

2、行星大小.軌道大小及運(yùn)行周期的比例三、 實(shí)驗(yàn)前的準(zhǔn)備工作1、對 MATLAB的初步了解： MATLAB（矩陣實(shí)驗(yàn)室）是 MATrix LABoratory 的縮寫，是一款由美國 TheMathWorks 公司出品的商業(yè)數(shù)學(xué)軟件。MATLAB是一種用于算法開發(fā)、數(shù)據(jù)可視化、數(shù)據(jù)分析以及數(shù)值計算的高級技術(shù)計算語言和交互式環(huán)境。除了矩陣運(yùn)算、繪制函數(shù)/ 數(shù)據(jù)圖像等常用功能外，MATLAB還可以用來創(chuàng)建用戶界面及與調(diào)用其它語言（包括C， C+和 FORTRAN）編寫的程序。2、組內(nèi)進(jìn)行多次討論，并明確組員分工：陳鐵、王思懋兩人負(fù)責(zé)前期資料 ( 包括太陽系的各種所需數(shù)據(jù)、 MATLAB的編譯環(huán)境等 )

3、的收集，陳桃斌負(fù)責(zé) ppt 的制作，楊飛負(fù)責(zé)主體代碼，鄒智康負(fù)責(zé)實(shí)驗(yàn)報告的書寫。由楊飛將主題代碼寫好后，再把其余代碼分塊分給其余組員各自書寫，以保證每個組員的參與度。四、 實(shí)驗(yàn)成果展示1、全部代碼2 、3 、4 、5 、6 、7 、8 、9 、10 、11 、12 、%3?ê ?éè ?figure('name', 'D? ?Yê ?' ); %éè ?±êì a?× ?pausetime=.01;%éè ?Yíê±

4、; ?set(gca,'xlim',-50 50,'ylim',-50 30,'zlim',-5050);set(gcf,'doublebuffer', 'on' )%?3y?xlabel('x?á'),ylabel('y? á'),zlabel('z?á');axisequal;gridon ;view(3 5 2);holdona=8.5 12.5 20 30 50 60 80 100 90 1.5;%3¤ ?á1

5、3 、14 、15 、16 、17 、18 、b=8 12 18 26 45 55 70 90 30 1.5;%?ì ?áomga=4 1.25 1 0.5 0.1 0.05 0.1 0.09 0.1 5;%?ù ?èr=0.35 0.8 0.8 0.5 3 2.5 1.5 1.5 0.5 0.35;%?òì ?° ?°° ?a? ó |DDD? ?ú ?a?D? ?úê ?a? òc=sqrt(a.2-b.2);h=pi/18;h1=pi/10;f=pi

6、/9;g=pi/8;g1=pi/6;aby=h h 0;h1 h 0;h h 0;h h 0;h h 0;h h 0;h g 0;h h h;g0 g;0 g1 0;%?1ìà ?ò ?éè ?19 、20 、21 、22 、23 、24 、25 、%D?DDD? ?ó ?íá D?·éè ?runu=35:0.5:40;theta=(0:0.05*pi:2*pi)'runa=2.8:0.4:5.6;xx=cos(theta)*runu+20;yy=0.9*sin(theta)

7、*runu;zz=-0.17*xx-0.17*yy; plot3(xx,yy,zz,':k');%D?DD ?éè ?hx=cos(theta)*runa;hy=sin(theta)*runa;hz=-0.1*hx-0.2*hy;%color='y','m','b','m','r','c','b','b'26 、27 、28 、29 、30 、31 、32 、33 、34 、35 、36 、37 、38 、39 、40 、41 、4

8、2 、43 、44 、%ì ?éè ?X,Y,Z=sphere(40);surf(5*X,5*Y,5*Z);colormap(autumn)%éè ?ì ?light ('position',1 0 2,'style', 'infinite')lightingphongmaterialshinyt=0:0.01*pi:50*pi; t'num=length(a);%1ìà· ?3ìforn=1:numx(:,n)=a(n)*cos(omga

9、(n)*t)+c(n);y(:,n)=b(n)*sin(omga(n)*t);z(:,n)=0*t;xuanz(:,:)=1 0 0;0 cos(aby(n,1) -sin(aby(n,1);0sin(aby(n,1) cos(aby(n,1)*cos(aby(n,2) 0sin(aby(n,2);0 1 0;-sin(aby(n,2) 0cos(aby(n,2)*cos(aby(n,3) -sin(aby(n,3)0;sin(aby(n,3) cos(aby(n,3) 0;0 0 1;xyz(:,:)=x(:,n) y(:,n) z(:,n)*xuanz(:,:);x(:,n)=xyz(:,

10、1);y(:,n)=xyz(:,2);z(:,n)=xyz(:,3);45 、46 、47 、48 、49 、50 、51 、52 、53 、ifn<=9p(n)=surf(r(n)*X+x(1,n),r(n)*Y+y(1,n),r(n)*Z+z(1,n);shadinginterpplot3(x(:,n),y(:,n),z(:,n),'-k');%?-3?ùó D1ì ? ?elsep(n)=surf(r(n)*X+x(1,n)+x(1,3),r(n)*Y+y(1,n)+y(1,3),r(n)*Z+z(1,n)+z(1,3);shading

11、interppmuw(n-9)=plot3(x(:,n)+x(1,3),y(:,n)+y(1,3),z(:,n)+z(1,3),'-k' );%?ò 1ìàéè ?endend54 、%íá D?·55 、forn=1:length(runa)56 、ph(n)=plot3(hx(:,n)+x(1,6),hy(:,n)+y(1,6),hz(:,n)+z(1,6), '-c');57 、end58 、59 、60 、61 、62 、63 、%DDD?é?set(p(1),&#

12、39;facecolor','y');set(p(2),'facecolor','m');set(p(3),'facecolor','b');set(p(4),'facecolor','m');set(p(5),'facecolor','r');set(p(6),'facecolor','c');set(p(7),'facecolor','b');set(p(8),'fac

13、ecolor','b');set(p(9),'facecolor','r');set(p(10),'facecolor', 'k');64 、 %?-3?í ?D?65 、 form=1:500066 、 forn=1:num67 、 ifn<=968 、 set(p(n),'xdata',r(n)*X+x(m,n),'ydata',r(n)*Y+y(m,n),'zdata' ,r(n)*Z+z(m,n);%?ùó D ?&

14、#234;± ?69 、else70 、 set(p(n),'xdata',r(n)*X+x(m,n)+x(m,3),'ydata',r(n)*Y+y(m,n)+y(m,3),'zdata',r(n)*Z+z(m,n)+z(m,3);71 、 set(pmuw(n-9),'xdata',x(:,n)+x(m,3),'ydata',y(:,n)+y(m,3),'zdata',z(:,n)+z(m,3);72 、end73 、end74 、75 、76 、77 、78 、forn=1:length(runa)set(ph(n),'xdata',hx(:,n)+x(m,6),'ydata',hy(:,n)+y(m,6), 'zdata',hz(:,n)+z(m,6);%1a?· ?ê± ?end79 、80 、81 、pause(pausetime);%?Yíò ?ádrawnowend2、分塊詮釋(1)這是初始的設(shè)置，包括一些必要數(shù)據(jù)的設(shè)置以及后面要用到的參數(shù)的設(shè)置。(2)這是小行星帶、土星帶以