locate_main.m

genetic遗传算法

function [sol,val] = locate_main(sol,options)
global bounds
% t=[0 1738 1688 0 1762 1586 1754	1660 1666 1472 1670 1478 1256 1330 1386 1564]/1000;
t=[0	1164	1126	0	1308	1068	0	1044	1104	0	0	1168	1134	1048	1182	1168]/1000;
x_cor=[39487714.54 39487190.16 39486897.50 39486467.13 39485936.18 39486344.51... 
    39485491.82 39487555.29 39487938.09 39486000.43 39485251.45 39485786.05...
    39486147.64 39486620.20 39485806.34 39485849.05];
y_cor=[3939197.67 3939147.28 3938985.83 3939299.85 3939319.40 3938493.23 3939196.54...
    3938345.63 3937955.98 3937954.13 3938591.75 3937837.04 3937302.06... 
    3937329.38 3937405.75 3938350.12];
z_cor=[-904.30 -847.60 -847.70 -664.20 -590.20 -850.50 -505.00 -850.00 -850.00... 
    -850.00 -505.00 -758.40 -846.50 -850.00 -754.90 -760.60];
v=3600;
LB=bounds(:,1)';UB=bounds(:,2)';
X0=sol(1:(end-1))';
options1=optimset('Diagnostics','off','LargeScale','off','Display','off');
[source fval]=fmincon(@(x)locate(x,x_cor,y_cor,z_cor,t,v),X0,[],[],[],[],LB,UB,[],options1);
sol(1:(end-1))=source';
val=-fval;