2009年全国大学生数学建模大赛C题优秀论文 下载本文

end

附件2:卫星或飞船地面轨迹图Matlab程序

load('topo.mat','topo','topomap1'); whos topo topomap1

contour(0:359,-89:90,topo,[0 0],'b') axis equal box on

set(gca,'XLim',[0 360],'YLim',[-90 90], ... 'XTick',[0 60 120 180 240 300 360], ... 'Ytick',[-90 -60 -30 0 30 60 90]); hold on

x=linspace(0,500*pi,1000); y1=58*sin(0.05*x+0.25*pi); y2=58*sin(0.05*x+0.5*pi); y3=58*sin(0.05*x+0.75*pi); y4=58*sin(0.05*x+pi);

y5=58*sin(0.05*x+1.25*pi); y6=58*sin(0.05*x+1.5*pi); y7=58*sin(0.05*x+1.75*pi); y8=58*sin(0.05*x+2*pi);

plot(x,y1,'k',x,y2,'g',x,y3,'r',x,y4,'y',x,y5,'b',x,y6,'g',x,y7,'m',x,y8,'r') hold on z1=0*x+58; z2=-z1;

plot(x,z1,'k',x,z2,'k')

附件3:圆内接六边形全覆盖图Matlab程序 load iris_dataset

net = newsom(irisInputs,[16 3]); plotsomtop(net); hold on

x=-1:0.05:16.5;

y1=0.85+1.15*sin(x);

y2=0.85+1.15*sin(x+0.25*pi); y3=0.85+1.15*sin(x+0.5*pi); y4=0.85+1.15*sin(x+0.75*pi); y5=0.85+1.15*sin(x+pi);

y6=0.85+1.15*sin(x+1.25*pi); y7=0.85+1.15*sin(x+1.5*pi); y8=0.85+1.15*sin(x+1.75*pi); z1=0*x+2; z2=0*x-0.3;

plot(x,y1,'k',x,y2,'g',x,y3,'r',x,y4,'y',x,y5,'b',x,y6,'g',x,y7,'m',x,y8,'r',x,z1,'k',x,z2,'k')

15

附件4:function y=cekongzhan(R,H,b) b=b./180.*pi;

s1=2.*pi.*(R+H).*(1-sin(b)); s=4.*pi.*(R+H).^2; s2=s-s1; e=0.827;

r=(R+H).*sin(87./180.*pi-asin(R.*(sin(93./180.*pi))./(R+H))); s3=pi.*r.^2; s4=s3.*e;

y=s2./s4; y=ceil(y); end

附件表1 神舟七号测控站地理位置

1 2 3 4 5 6 7 8 9 10 11 12 13 14

卫星观测站 北京航天指挥控制中心 酒泉卫星发射中心 西安卫星控制中心 山西兴县卫星测控站 陕西渭南卫星测控站 厦门卫星测控站 喀什卫星测控站 远望一号测控船 远望三号测控船 卡拉奇卫星测控站 和田卫星测控站 南非卫星测控站 远望二号测控船 愿望四号测控船 经度 120.25°E 100.50°E 108.91°E 111.22°E 109.50°E 118.01°E 75.99°E 131.83°E 13.72°W 66.99°E 79.94°E 18.42°E 75.84°W 151.00°E 纬度 48.07°N 41.83°N 34.23°N 38.47°N 34.50°N 24.07°N 39.45°N 37.23°N 09.50°N 24.82°N 37.12°N 33.92°N 19.84°N 33.00°S 16