遗传算法matlab源代码

%% 子函数 % %题 目：Caldata % %输 入： % pop 种群 % bestSeita 最佳seita角 % bestTime 最佳时间 %输 出： % data 各种参量 %历 史： % 2015-8-22 编写程序 %% function data = Caldata(pop,bestSeita,bestTime) %% 数据 g = 1.633;%49.3247*10^11;%月球重力加速度 m = 2400*ones(1,1);%嫦娥3号质量2400Kg V_e = 2940;%冲比2940m/s R = 1737013;%平均月球半径1737013m，着陆点高度-2640m h_0 = 15000;%变轨起点高度15000m h_end = 3000;%变轨终点高度3000m V_r = zeros(1,1);%向心速度0m/s V_st = 1692.46*ones(1,1);%切向速度1692.46m/s r = (R+h_0)*ones(1,1);%距月心距离 aerfa = zeros(1,1);%极角 T = 7500; %% data = {[]};%初始化cell for k=1:size(pop,2)%分别计算每个个体的适应度 %读入seita多项式系数和T for t=1:size(bestSeita,2) seita(t) = -8*10^-9*t*(t-bestTime)*(t-pop(1,k))+bestSeita(t); %求夹角seita m(t+1) = m(t)-T/V_e; a(t) = T/m(t); %求加速度 V_r(t+1) = V_r(t)-g/r(t)^2+V_st(t)^2/r(t)-a(t)*sin(seita(t)); V_st(t+1) = V_st(t)-V_r(t)*V_st(t)/r(t)-a(t)*cos(seita(t)); r(t+1) = r(t)-V_r(t); aerfa(t+1) = aerfa(t)+V_st(t)/r(t); end data{k,1} = {seita a V_r V_st aerfa seita r [t] m [T*t/V_e]}; end