sov_ht.asv

这个一个遗传算法的实例

function HT=SOV_HT(PHTL,QAHTL,VL,HL,TL)
%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;
for i=1:5;
    T(1,:)=TL(i):0.01:TL(i+1);
end;
%T1=t1:0.01:t2;T2=t2:0.01:TL(3);T(3,:)=TL(3):0.01:TL(4);T(4,:)=TL(4):0.01:TL(5);T(5,:)=TL(5):0.01:TL(6);

Q1=QAHTL(1)/(6*HL(1))*(TL(2)-T(1,:)).^3+QAHTL(2)/(6*HL(1))*(T(1,:)-TL(1)).^3+(PHTL(2)/HL(1)-HL(1)*QAHTL(2)/6)*(T(1,:)-TL(1))+(PHTL(1)/HL(1)-HL(1)*QAHTL(1)/6)*(TL(2)-T(1,:));
Q2=QAHTL(2)/(6*HL(2))*(TL(3)-T(2,:)).^3+QAHTL(3)/(6*HL(2))*(T(2,:)-TL(2)).^3+(PHTL(3)/HL(2)-HL(2)*QAHTL(3)/6)*(T(2,:)-TL(2))+(PHTL(2)/HL(2)-HL(2)*QAHTL(2)/6)*(TL(3)-T(2,:));
Q3=QAHTL(3)/(6*HL(3))*(TL(4)-T(3,:)).^3+QAHTL(4)/(6*HL(3))*(T(3,:)-TL(3)).^3+(PHTL(4)/HL(3)-HL(3)*QAHTL(4)/6)*(T(3,:)-TL(3))+(PHTL(3)/HL(3)-HL(3)*QAHTL(3)/6)*(TL(4)-T(3,:));
Q4=QAHTL(4)/(6*HL(4))*(TL(5)-T(4,:)).^3+QAHTL(5)/(6*HL(4))*(T(4,:)-TL(4)).^3+(PHTL(5)/HL(4)-HL(4)*QAHTL(5)/6)*(T(4,:)-TL(4))+(PHTL(4)/HL(4)-HL(4)*QAHTL(4)/6)*(TL(5)-T(4,:));
Q5=QAHTL(5)/(6*HL(5))*(TL(6)-T(5,:)).^3+QAHTL(6)/(6*HL(5))*(T(5,:)-TL(5)).^3+(PHTL(6)/HL(5)-HL(5)*QAHTL(6)/6)*(T(5,:)-TL(5))+(PHTL(5)/HL(5)-HL(5)*QAHTL(5)/6)*(TL(6)-T(5,:));

V1=-QAHTL(1)/(2*HL(1))*(TL(2)-T(1,:)).^2+QAHTL(2)/(2*HL(1))*(T(1,:)-TL(1)).^2+(PHTL(2)/HL(1)-HL(1)*QAHTL(2)/6)-(PHTL(1)/HL(1)-HL(1)*QAHTL(1)/6);
V2=-QAHTL(2)/(2*HL(2))*(TL(3)-T(2,:)).^2+QAHTL(3)/(2*HL(2))*(T(2,:)-TL(2)).^2+(PHTL(3)/HL(2)-HL(2)*QAHTL(3)/6)-(PHTL(2)/HL(2)-HL(2)*QAHTL(2)/6);
V3=-QAHTL(3)/(2*HL(3))*(TL(4)-T(3,:)).^2+QAHTL(4)/(2*HL(3))*(T(3,:)-TL(3)).^2+(PHTL(4)/HL(3)-HL(3)*QAHTL(4)/6)-(PHTL(3)/HL(3)-HL(3)*QAHTL(3)/6);
V4=-QAHTL(4)/(2*HL(4))*(TL(5)-T(4,:)).^2+QAHTL(5)/(2*HL(4))*(T(4,:)-TL(4)).^2+(PHTL(5)/HL(4)-HL(4)*QAHTL(5)/6)-(PHTL(4)/HL(4)-HL(4)*QAHTL(4)/6);
V5=-QAHTL(5)/(2*HL(5))*(TL(6)-T(5,:)).^2+QAHTL(6)/(2*HL(5))*(T(5,:)-TL(5)).^2+(PHTL(6)/HL(5)-HL(5)*QAHTL(6)/6)-(PHTL(5)/HL(5)-HL(5)*QAHTL(5)/6);
%m=-QAHTL(5)/(2*HL(5))*(TL(6)-7.44).^2+QAHTL(6)/(2*HL(5))*(7.44-TL(5)).^2+(PHTL(6)/HL(5)-HL(5)*QAHTL(6)/6)-(PHTL(5)/HL(5)-HL(5)*QAHTL(5)/6)

A1=QAHTL(1)/HL(1)*(TL(2)-T(1,:))+QAHTL(2)/HL(1)*(T(1,:)-TL(1));
A2=QAHTL(2)/HL(2)*(TL(3)-T(2,:))+QAHTL(3)/HL(2)*(T(2,:)-TL(2));
A3=QAHTL(3)/HL(3)*(TL(4)-T(3,:))+QAHTL(4)/HL(3)*(T(3,:)-TL(3));
A4=QAHTL(4)/HL(4)*(TL(5)-T(4,:))+QAHTL(5)/HL(4)*(T(4,:)-TL(4));
A5=QAHTL(5)/HL(5)*(TL(6)-T(5,:))+QAHTL(6)/HL(5)*(T(5,:)-TL(5));

J1=(T(1,:)*(QAHTL(2)-QAHTL(1))/HL(1))./(T(1,:)+eps);
J2=T(2,:)*(QAHTL(3)-QAHTL(2))/HL(2)./T(2,:);
J3=T(3,:)*(QAHTL(4)-QAHTL(3))/HL(3)./T(3,:);
J4=T(4,:)*(QAHTL(5)-QAHTL(4))/HL(4)./T(4,:);
J5=T(5,:)*(QAHTL(6)-QAHTL(5))/HL(4)./T(5,:);

plot(T(1,:),Q1,'b:',T(2,:),Q2,'b:',T(3,:),Q3,'b:',T(4,:),Q4,'b:',T(5,:),Q5,'b:')
hold on
plot(T(1,:),V1,'-b',T(2,:),V2,'-b',T(3,:),V3,'-b',T(4,:),V4,'-b',T(5,:),V5,'-b')
hold on
plot(T(1,:),A1,'-r',T(2,:),A2,'-r',T(3,:),A3,'-r',T(4,:),A4,'-r',T(5,:),A5,'-r')
plot(T(1,:),J1,'r:',T(2,:),J2,'r:',T(3,:),J3,'r:',T(4,:),J4,'r:',T(5,:),J5,'r:')
hold off
%grid