遗传算法-matlab代码

clear; clc; % x1=-2.687069; % x2=5.361653; % y=code(x1,x2) % disp(y); % [a1,a2]=codeback(y); % vpa(a1) % vpa(a2) [a,b]=init() for i=1:1:10 v(i,:)=code(a(i),b(i)); end f=fun(a,b); for i=1:1:1000 [p,q]=possible(f); r=rand(1,10); num=generation(r,q); C=cross_over(v,0.30,num) x=Mutation(C,0.01); for j=1:1:10 [x1(j),x2(j)]=codeback(x(j,:)); end f=fun(x1,x2); value(i)=max(f); if i==1 max_value(i)=value(i); else if value(i)>max_value(i-1) max_value(i)=value(i); else max_value(i)= max_value(i-1); end end end t=1:1:1000; plot(t,value,'r-'); hold on; plot(t,max_value,'g-'); vpa(max_value(1000))