2014年高教社杯A题《嫦娥三号软着陆轨道设计与控制策略》全国一等奖优秀论文
s=std(b(1:end),1); s=s*s; yb1(i,j)=s; end end yb2=[]
for i=1:989 for j=1:971 b=A(i:i+11,j:j+29); s=std(b(1:end),1); s=s*s; yb2(i,j)=s; end end
mb1=min(yb1) mb2=min(yb2) mmb1=min(mb1) mmb2=min(mb2) if mmb1<mmb2
[row,column]=find(yb1==mmb1) else
[row,column]=find(yb2==mmb2) end