[转帖](来自罗卜驿站)模拟太阳系运动

huright

1. function solar_system
   
2. % 模拟太阳系运动
   
3. t=linspace(0,2*pi,100);
   
4. fill(cos(t),sin(t),'r');
   
5. hold on;
   
6. plot(4*cos(t),sin(t)*4,'k');
   
7. set(gca,'position',[0 0.11 0.775 0.815])
   
8. a=0.1;b=0;
   
9. xe=4*cos(a)+cos(t)*0.6;
   
10. ye=4*sin(a)+sin(t)*0.6;
    
11. He=fill(xe,ye,'b');
    
12. xm=4*cos(a)+cos(b);
    
13. ym=4*sin(a)+sin(b);
    
14. set(gcf,'doublebuffer','on');
    
15. Hm=plot(xm,ym,'c.','markersize',24);
    
16. aa=gca;
    
17. axis([-6,6,-6,6]);
    
18. axis square;
    
19. k=1;da=0.1;db=0.5;
    
20. xlabel('Please press "space" key and stop this program!',...
    
21.   'fontsize',12,'color','r');
    
22. title('simulate solar system')
    
23. axes('position',[0.75,0.11,0.25,0.8]);
    
24. fill(0.2+cos(t)*0.18,0.75+sin(t)*0.08,'r');
    
25. ylim([0,1]);xlim([0,0.9]);
    
26. text(0.5,0.75,'Sun');hold on;
    
27. fill(0.2+cos(t)*0.11,0.5+sin(t)*0.05,'b');
    
28. text(0.5,0.5,'Earth');
    
29. plot(0.2,0.3,'c.','markersize',24);
    
30. text(0.5,0.3,'Moon');
    
31. axis off
    
32. axes(aa);
    
33. while k;
    
34.    s=get(gcf,'currentkey');
    
35.    if strcmp(s,'space');
    
36.        clc;k=0;
    
37.    end
    
38.    a=a+da;
    
39.    b=b+db;
    
40.    xe=4*cos(a)+cos(t)*0.6;
    
41.    ye=4*sin(a)+sin(t)*0.6;
    
42.    xm=4*cos(a)+cos(b);
    
43.    ym=4*sin(a)+sin(b);
    
44.    set(He,'xdata',xe,'ydata',ye);
    
45.    set(Hm,'xdata',xm,'ydata',ym);
    
46.    pause(0.1);
    
47.    if a<80;
    
48.        plot(xm,ym);
    
49.    end
    
50. end
    
51. figure(gcf);

复制代码

[ 本帖最后由 suffer 于 2006-10-9 19:42 编辑 ]

hunter_009

呵呵，这也不能叫太阳系吧,感觉简单了一点。<BR>楼主能否做一个更好一点的。

bb20040291

楼主 我一运行就是错误啊 怎么回事啊？

suffer

原帖由 bb20040291 于 2006-7-25 01:04 发表
楼主 我一运行就是错误啊 怎么回事啊？

请说明错误提示

suffer

原帖由 bb20040291 于 2006-7-25 01:04 发表
楼主 我一运行就是错误啊 怎么回事啊？

有可能是版本问题造成的，我这里在6.5下运行没有什么错误

puyen

我在6.5下运行没有什么错误

frinim

可以进一步进行修改吗？好像有一点粗糙啊。

suffer

原帖由 frinim 于 2006-10-28 21:16 发表
可以进一步进行修改吗？好像有一点粗糙啊。

你想修改什么？

geoer

那就改掉吧，好好跟他学点。

[ 本帖最后由 geoer 于 2007-1-11 08:30 编辑 ]

心灯

原帖由 geoer 于 2007-1-10 22:10 发表
萝卜 :
...

没必要把萝卜的个人真实资料贴这儿来吧？...

刚才到萝卜驿站上看了一下，发现他公布的自己的信息包括上面的内容了，估计geoer也是从上面粘贴过来的....

[ 本帖最后由 心灯 于 2007-1-10 22:34 编辑 ]