x0=x(1)-(T0-12)*(x(21)-x(1)); y0=x0*p(1)+p(2);
l=(x0^2+y0^2)^(1/2); end
%主程序：本程序采用遗传算法接力进化，
%将上次进化结束后得到的最终种群作为下次输入的初始种群 clc;
close all; clear all; %进化的代数 T=100;
optionsOrigin=gaoptimset('Generations',T/2); x=zeros(1,3);
[x,fval]=ga(@ch14_1f,3);
%,optionsOrigin||reason,output,finnal_pop %进行第二次接力进化
% options1=gaoptimset('Generations',T/2,'InitialPopulation',finnal_pop,... % 'PlotFcns',@gaplotbestf);
[x,fval]=ga(@ch14_1f,3,[],[],[],[],[0;-90;0],[180;90;365]); Bestx=x
BestFval=fval
function y=distance(address,order) nmb=size(address,1); y=0;
for i=1:nmb-1
y=y+sqrt((address(order(i+1),1)-address(order(i),1))^2+(address(order(i+1),2)-address(order(i),2))^2); end
y=y+sqrt((address(order(i+1),1)-address(order(1),1))^2+(address(order(i+1),2)-address(order(1),2))^2);
function y=exhgpath(order) while 1
b=size(order,1); r=unidrnd(b,1,2); if r(1)-r(2)~=0 break end end
b=order(r(2));
order(r(2))=order(r(1));
全国大学生数学建模竞赛一等奖论文
order(r(1))=b; y=order;
function [theta1,rough]=ch(jingdu1,weidu1,day) % load Qb.mat % x=Qb(:,2); % y=Qb(:,3);
[x,y]=yizhi(jingdu1,weidu1,day); g=length(weidu1); t1=12.7:1/20:13.7; n=length(t1);
beta0=chiwei(2015,day); beta=beta0; ?ta=18; for i=1:g
b(i)=xianxin(jingdu1,x,y);
h(i)=abs(tand(90-(beta+weidu1(i)))*b(i)); end
for i=1:g
t=dangdishi(jingdu1,120,t1); % b=xianxin(jingdu1(i),x,y); % h=tand(90-(beta+weidu1(i)))*b; % h=fun2(jingdu1(i),x,y); alpha=weidu1(i); for j=1:n
theta(j)=(12-t(j))*15;
x(j)=-cosd(beta)*sind(theta(j))/(cosd(alpha)*cosd(theta(j))*cosd(beta)+sind(alpha)*sind(beta))*h(i);
y(j)=(sind(alpha)*cosd(beta).*cosd(theta(j))-cosd(alpha)*sind(beta))/(cosd(alpha)*cosd(theta(j))*cosd(beta)+sind(alpha)*sind(beta))*h(i); end
[theta1(i,:),rough(i,:)]=ji(x,y); end
%子函数：适应度函数同时也是目标函数,函数存储名称为ch14_2f.m function f=ch14_1f(x) jingdu=x(1); weidu=x(2); day=x(3);
f=ff(jingdu,weidu,day); end
全国大学生数学建模竞赛一等奖论文
问题三
clc,clear; N=41;
for i=1:41
b=huidu(i);
data(:,:,i)=rgb2gray(b); end
for i=1:41
a(:,:,i)=data(200:1000,600:1800,i); end
% a(:,:,2)=a(:,:,2)-a(:,:,1); %
for i=1:41
for j=1:801
for k=1:1201
c(j,k,i)=a(j,k,i)-a(j,k,1); end end end
imshow(c(:,:,2))
load c.mat for i=1:41
d(:,:,i)=c(600:700,1000:1100,i); end
for i=2:41 if i<=10
[j,k]=find(d(:,:,i)>=10); [A,B]=max(k); x(i-1)=A+1000; y(i-1)=k(B)+600; else
[j,k]=find(d(:,:,i)>=15); [A,B]=max(k); x(i-1)=A+1000; y(i-1)=k(B)+600; end end
x=x-277; y=y-678;
k1=2/((878-320)^2+9^2)^(1/2); x=k1*x; y=k1*y;
全国大学生数学建模竞赛一等奖论文
p4=k1*(960^2+540^2)^(1/2);
clc,clear; e11=[-62,0,1]; e12=[43,0,1]; e21=[43,0,1]; e22=[43,21,1]; e31=[43,0,1]; e32=[59,4,1]; a=cross(e11,e12); b=cross(e21,e22); c=cross(e31,e32);
m1=[a(1)*a(1),a(1)*a(2),a(1)*a(3); a(2)*a(1),a(2)*a(2),a(2)*a(3); a(3)*a(1),a(3)*a(2),a(3)*a(3)]; m2=[b(1)*b(1),b(1)*b(2),b(1)*b(3); b(2)*b(1),b(2)*b(2),b(2)*b(3); b(3)*b(1),b(3)*b(2),b(3)*b(3)]; m3=[c(1)*c(1),c(1)*c(2),c(1)*c(3); c(2)*c(1),c(2)*c(2),c(2)*c(3); c(3)*c(1),c(3)*c(2),c(3)*c(3)]; m=m1+m2+m3; [x,y]=eig(m);
eigenvalue=diag(y); y_lamda = x(:, 3);
问题二的实测照片：
全国大学生数学建模竞赛一等奖论文
问题三的实测xy坐标： X Y X Y X Y 0.269204 0.189235 0.32719 0.186698 0.393107 0.183429 0.277091 0.188911 0.336062 0.186281 0.403334 0.18289 0.2851 0.188575 0.345101 0.185848 0.413794 0.18233 0.293237 0.188226 0.354317 0.185398 0.424499 0.181749 0.301509 0.30992 0.187865 0.18749 0.363716 0.373308 0.184933 0.18445 0.435461 0.446692 0.181146 0.18052 0.318479 0.187101 0.383101 0.183949 0.458207 0.179869
问题四的实际坐标： X Y X Y X Y X Y X Y X Y -3.43325 -40.5397 -3.36515 -43.4685 -3.43325 -40.5397 -3.36515 -43.4685 -3.36157 -43.6226 -3.36515 -43.4685 -3.42608 -40.848 -3.36515 -43.4685 -3.43325 -40.5397 -3.36157 -43.6226 -3.36157 -43.6226 -3.36515 -43.4685 -3.42608 -40.848 -3.42966 -40.6939 -3.43325 -40.5397 -3.36515 -43.4685 -3.36157 -43.6226 -3.36515 -43.4685 -3.36157 -43.6226 -3.35799 -43.7768 -3.42608 -40.848 -3.42966 -40.6939 -3.42966 -40.6939 -3.43325 -40.5397 -3.36515 -43.4685 -3.43325 -40.5397 -3.42608 -40.848 -3.36157 -43.6226 -3.36157 -43.6226 -3.43683 -40.3856 -3.42608 -40.848 -3.43325 -40.5397 -3.36515 -43.4685 -3.42966 -40.6939 -3.35799 -43.7768 -3.35799 -43.7768 -3.43325 -40.5397 -3.36157 -43.6226 -3.36157 -43.6226 -3.36157 -43.6226
全国大学生数学建模竞赛一等奖论文