美赛数模论文 - 图文

Team # 35565 Page 45 of 47

237_0.shtml. 2014

[13] Li Guo, The difficult locating wreckage: Satellite Monitoring of Chinese Academy of Sciences. The 21st century business herald. 2014

Appendix

? Matlab Source Code for Submarine Topography

clc,clear

I=imread('map.jpg'); imshow(I); [x,y]=size(I); [x,y]=ginput(36);

xlswrite('book.xls',[x,y])

clc, clear

a=xlsread('book.xls'); x=a(:,1)';

y=a(:,2)'; z=-a(:,3)';

xmm=minmax(x); ymm=minmax(y); xi=xmm(1):xmm(2); yi=ymm(1):ymm(2);

zi1=griddata(x,y,z,xi,yi','cubic'); zi2=griddata(x,y,z,xi,yi','nearest'); zi=zi1;

zi(isnan(zi1))=zi2(isnan(zi1)); figure(1), plot(x,y,'*') figure(2), mesh(xi,yi,zi)

? Track of the Object In The Water

clc,clear m=242000; g=9.8;

Team # 35565 Page 46 of 47

p=1030;

k=-138.32815; v=200*0.28;

g_h=(m*g-p*g*v)/m; a=0.04; h=3400;

v0=200*0.28;

t_end=sqrt(2*h/g_h); t=0:0.1:t_end; h=-1/2.*g_h.*t.*t; x=v0.*t-1/2.*a.*t.*t;

plot(x,h); hold on;

x_end=v0*t_end-1/2*a*t_end*t_end h_end=-3400;

v0_end=v0-a*t_end

a_down=1-0.308/2420*(v0_end/2)^2 t2_end=v0_end/a_down; t2=0:0.1:t2_end;

l=v0_end.*t2-1/2.*a_down.*t2.*t2+x_end; [k len]=size(t2);

h2=zeros(1,len)-3400; plot(l,h2); figure(2) plot(x,h);

? Algorithm Flow of NonconvexDivide

input:CIj,j?(1,?,i)L?(w1,wk)S(IL)?SiwhileArea(ILr)?AreaRequired(S(CILr))andLe?wmdo

if Area(IL)?AreaRequired(S(CIL))then if Le?Si then

rrr

Team # 35565 Page 47 of 47

k1?1;

Ls?wk1;elseL1?(Ls,wk?1);T?(t1,t2,t3)r

if Area(IL)?AreaRequired(S(CIL)) then

rT??triangle(t1,t2,t3); DetachAndAssign(I?T?Predpoly(CI,(t1,t)));

DetachAndAssign(I?T);Area(CIL1?PredPoly(CI,(t1,t2)))?AreaReqired(S(CI))thenT?triangle(t1,t2,t3)DetachAndAssign(I?T);DetachAndAssign(I?T?PredPoly(CI,(t1,t)));elseT?(triangle(t1,IS,t3)); Else if AreaRequired(IS)?AreaRequired(S(CI))?Area(CI?T);

rDetachAndAssign(PredPoly(CI,(t1,t2)));DetachAndAssign(P?T);DetachAndAssign(P?T);elseL?(t,Si);DetachAndAssign(IL1);DetachAndAssign(IL1)lr