ht.asv

这个一个遗传算法的实例

clear
v1=0;a1=0;v6=0;a6=0;QBJ=0;QEJ=0;
q1=-10;q3=60;q4=20;q6=55;
h1=0.8139;h2=2.0387;h3=2.1821;h4=1.8392;h5=0.8910;%100次的优化结果
t1=0;t2=t1+h1;t3=t2+h2;t4=t3+h3;t5=t4+h4;t6=t5+h5;
%t1=0,t2=t1+h1,t3=t2+h2,t4=t3+h3,t5=t4+h4,t6=t5+h5
a11=3*h1+2*h2+h1^2/h2;a12=h2;
a21=h2-h1^2/h2;a22=2*(h2+h3);a23=h3;
a32=h3;a33=2*(h3+h4);a34=h4-h5^2/h4;
a43=h4;a44=3*h5+2*h4+h5^2/h4;
A=[a11 a12 0 0;
   a21 a22 a23 0;
   0 a23 a33 a34;
   0 0 a43 a44];
b1=6*(q3/h2+q1/h1)-6*(1/h1+1/h2)*(q1+h1*v1+(h1^2/3)*a1)-h1*a1;
b2=(6/h2)*(q1+h1*v1+h1^2*a1/3)+6*q4/h3-6*(1/h2+1/h3)*q3;
b3=(6/h4)*(q6-v6*h5+h5^2*a6/3)-6*(1/h4+1/h3)*q4+(6/h3)*q3;
b4=-6*(1/h5+1/h4)*(q6-v6*h5+h5^2*a6/3)+6*q6/h5+6*q4/h4-a6;
B=[b1;b2;b3;b4];
Q=A\B;
q2=q1+h1*v1+h1^2*a1/3+h1^2*Q(1)/6;
q5=q6-h5*v6+h5^2*a6/3+h5^2*Q(4)/6;
T1=t1:0.01:t2;T2=t2:0.01:t3;T3=t3:0.01:t4;T4=t4:0.01:t5;T5=t5:0.01:t6;

Q1=QBJ/(6*h1)*(t2-T1).^3+Q(1)/(6*h1)*(T1-t1).^3+(q2/h1-h1*Q(1)/6)*(T1-t1)+(q1/h1-h1*QBJ/6)*(t2-T1);
Q2=Q(1)/(6*h2)*(t3-T2).^3+Q(2)/(6*h2)*(T2-t2).^3+(q3/h2-h2*Q(2)/6)*(T2-t2)+(q2/h2-h2*Q(1)/6)*(t3-T2);
Q3=Q(2)/(6*h3)*(t4-T3).^3+Q(3)/(6*h3)*(T3-t3).^3+(q4/h3-h3*Q(3)/6)*(T3-t3)+(q3/h3-h3*Q(2)/6)*(t4-T3);
Q4=Q(3)/(6*h4)*(t5-T4).^3+Q(4)/(6*h4)*(T4-t4).^3+(q5/h4-h4*Q(4)/6)*(T4-t4)+(q4/h4-h4*Q(3)/6)*(t5-T4);
Q5=Q(4)/(6*h5)*(t6-T5).^3+QEJ/(6*h5)*(T5-t5).^3+(q6/h5-h5*QEJ/6)*(T5-t5)+(q5/h5-h5*Q(4)/6)*(t6-T5);

V1=-QBJ/(2*h1)*(t2-T1).^2+Q(1)/(2*h1)*(T1-t1).^2+(q2/h1-h1*Q(1)/6)-(q1/h1-h1*QBJ/6);
V2=-Q(1)/(2*h2)*(t3-T2).^2+Q(2)/(2*h2)*(T2-t2).^2+(q3/h2-h2*Q(2)/6)-(q2/h2-h2*Q(1)/6);
V3=-Q(2)/(2*h3)*(t4-T3).^2+Q(3)/(2*h3)*(T3-t3).^2+(q4/h3-h3*Q(3)/6)-(q3/h3-h3*Q(2)/6);
V4=-Q(3)/(2*h4)*(t5-T4).^2+Q(4)/(2*h4)*(T4-t4).^2+(q5/h4-h4*Q(4)/6)-(q4/h4-h4*Q(3)/6);
V5=-Q(4)/(2*h5)*(t6-T5).^2+QEJ/(2*h5)*(T5-t5).^2+(q6/h5-h5*QEJ/6)-(q5/h5-h5*Q(4)/6);
m=-Q(4)/(2*h5)*(t6-7.44).^2+QEJ/(2*h5)*(7.44-t5).^2+(q6/h5-h5*QEJ/6)-(q5/h5-h5*Q(4)/6)

A1=QBJ/h1*(t2-T1)+Q(1)/h1*(T1-t1);
A2=Q(1)/h2*(t3-T2)+Q(2)/h2*(T2-t2);
A3=Q(2)/h3*(t4-T3)+Q(3)/h3*(T3-t3);
A4=Q(3)/h4*(t5-T4)+Q(4)/h4*(T4-t4);
A5=Q(4)/h5*(t6-T5)+QEJ/h5*(T5-t5);

J1=(T1*(Q(1)-QBJ)/h1)./(T1+eps);
J2=T2*(Q(2)-Q(1))/h2./T2;
J3=T3*(Q(3)-Q(2))/h3./T3;
J4=T4*(Q(4)-Q(3))/h4./T4;
J5=T5*(QEJ-Q(4))/h4./T5;

plot(T1,Q1,'b:',T2,Q2,'b:',T3,Q3,'b:',T4,Q4,'b:',T5,Q5,'b:')
hold on
plot(T1,V1,'-b',T2,V2,'-b',T3,V3,'-b',T4,V4,'-b',T5,V5,'-b')
hold on
plot(T1,A1,'-r',T2,A2,'-r',T3,A3,'-r',T4,A4,'-r',T5,A5,'-r')
plot(T1,J1,'r:',T2,J2,'r:',T3,J3,'r:',T4,J4,'r:',T5,J5,'r:')
hold off
%grid