sol_qasy.asv

这个一个遗传算法的实例

function QASY=SOL_QASY(HL,VL,AL,PL,NL,nL)
%solve the QA,that is a N*n metrixs
B(1,1)=0.5*HL(1)+HL(1)^2/(6*HL(2))+HL(2)/3;B(1,2)=HL(2)/6;
B(2,1)=(HL(1)+HL(2))/2;B(2,2)=HL(2)/2+HL(3)/3;B(2,3)=HL(3)/6;
B(3,1)=(HL(1)+HL(2))/2;B(3,2)=(HL(2)+HL(3))/2;B(3,3)=HL(3)/2+HL(4)/3;B(3,4)=HL(4)/6;
B(4,1)=(HL(1)+HL(2))/2;B(4,2)=(HL(2)+HL(3))/2;B(4,3)=(HL(3)+HL(4))/2;B(4,4)=HL(4)/2+HL(5)/3;B(4,5)=HL(5)/6;
B(5,1)=(HL(1)+HL(2))*(HL(5)+HL(6))/2;B(5,2)=(HL(2)+HL(3))*(HL(5)+HL(6))/2;B(5,3)=(HL(3)+HL(4))*(HL(5)+HL(6))/2;B(5,4)=(HL(4)+HL(5))*(HL(5)+HL(6))/2;B(5,5)=HL(5)*HL(6)/2+(HL(5)+HL(6)*)HL(7)/2+HL(6)^2/3;B(5,6)=HL(6)*HL(7)/2+HL(6)^2/6+HL(7)^2/3;
B(6,1)=(HL(1)+HL(2))/2;B(6,2)=(HL(2)+HL(3))/2;B(6,3)=(HL(3)+HL(4))/2;B(6,4)=(HL(4)+HL(5))/2;B(6,5)=(HL(5)+HL(6))/2;B(6,6)=(HL(6)+HL(7))/2;
QAL=zeros(6,3);
for i=1:NL;
    C(1,i)=(PL(i,3)-PL(i,1)-(HL(1)+HL(2))*VL(i,1)-(HL(1)^2/3+HL(1)*HL(2)/2)*AL(i,1))/HL(2);
    C(2,i)=(PL(i,4)-PL(i,3))/HL(3)-HL(1)*AL(i,1)/2-VL(i,1);
    C(3,i)=(PL(i,5)-PL(i,4))/HL(4)-HL(1)*AL(i,1)/2-VL(i,1);
    C(4,i)=(PL(i,6)-PL(i,5))/HL(5)-HL(1)*AL(i,1)/2-VL(i,1);
    C(5,i)=(PL(i,8)-PL(i,6))-HL(7)*AL(i,8)/6;
    C(6,i)=
    -6*(1/HL(7)+1/HL(6))*(PL(i,8)-VL(i,8)*HL(7)+HL(7)^2*AL(i,8)/3)+6*PL(i,8)/HL(7)+6*PL(i,6)/HL(6)-AL(i,8)*HL(7);
    QAL(:,i)=B\C(:,i);
    
    
    
end;
QA=AL;
for j=1:NL;
    for i=2:nL-1;
        QA(j,i)=QAL(i-1,j);
    end;
end;