[转帖](来自罗卜驿站)模拟太阳系运动
function solar_system% 模拟太阳系运动
t=linspace(0,2*pi,100);
fill(cos(t),sin(t),'r');
hold on;
plot(4*cos(t),sin(t)*4,'k');
set(gca,'position',)
a=0.1;b=0;
xe=4*cos(a)+cos(t)*0.6;
ye=4*sin(a)+sin(t)*0.6;
He=fill(xe,ye,'b');
xm=4*cos(a)+cos(b);
ym=4*sin(a)+sin(b);
set(gcf,'doublebuffer','on');
Hm=plot(xm,ym,'c.','markersize',24);
aa=gca;
axis([-6,6,-6,6]);
axis square;
k=1;da=0.1;db=0.5;
xlabel('Please press "space" key and stop this program!',...
'fontsize',12,'color','r');
title('simulate solar system')
axes('position',);
fill(0.2+cos(t)*0.18,0.75+sin(t)*0.08,'r');
ylim();xlim();
text(0.5,0.75,'Sun');hold on;
fill(0.2+cos(t)*0.11,0.5+sin(t)*0.05,'b');
text(0.5,0.5,'Earth');
plot(0.2,0.3,'c.','markersize',24);
text(0.5,0.3,'Moon');
axis off
axes(aa);
while k;
s=get(gcf,'currentkey');
if strcmp(s,'space');
clc;k=0;
end
a=a+da;
b=b+db;
xe=4*cos(a)+cos(t)*0.6;
ye=4*sin(a)+sin(t)*0.6;
xm=4*cos(a)+cos(b);
ym=4*sin(a)+sin(b);
set(He,'xdata',xe,'ydata',ye);
set(Hm,'xdata',xm,'ydata',ym);
pause(0.1);
if a<80;
plot(xm,ym);
end
end
figure(gcf);
[ 本帖最后由 suffer 于 2006-10-9 19:42 编辑 ] 呵呵，这也不能叫太阳系吧,感觉简单了一点。<BR>楼主能否做一个更好一点的。 楼主 我一运行就是错误啊 怎么回事啊？ 原帖由 bb20040291 于 2006-7-25 01:04 发表
楼主 我一运行就是错误啊 怎么回事啊？
请说明错误提示 原帖由 bb20040291 于 2006-7-25 01:04 发表
楼主 我一运行就是错误啊 怎么回事啊？
有可能是版本问题造成的，我这里在6.5下运行没有什么错误 我在6.5下运行没有什么错误 可以进一步进行修改吗？好像有一点粗糙啊。 原帖由 frinim 于 2006-10-28 21:16 发表
可以进一步进行修改吗？好像有一点粗糙啊。
你想修改什么？ 那就改掉吧，好好跟他学点。
[ 本帖最后由 geoer 于 2007-1-11 08:30 编辑 ] 原帖由 geoer 于 2007-1-10 22:10 发表
萝卜 :
...
没必要把萝卜的个人真实资料贴这儿来吧？...
刚才到萝卜驿站上看了一下，发现他公布的自己的信息包括上面的内容了，估计geoer也是从上面粘贴过来的....
[ 本帖最后由 心灯 于 2007-1-10 22:34 编辑 ]
页:
[1]