sol_qa.m

这个一个遗传算法的实例

function QA=SOL_QA(HL,VL,AL,PL,NL,nL)
%solve the QA,that is a N*n metrixs
B(1,1)=3*HL(1)+2*HL(2)+HL(1)^2/HL(2);B(1,2)=HL(2);
B(2,1)=HL(2)-HL(1)^2/HL(2);B(2,2)=2*(HL(2)+HL(3));B(2,3)=HL(3);
B(3,2)=HL(3);B(3,3)=2*(HL(3)+HL(4));B(3,4)=HL(4);
B(4,3)=HL(4);B(4,4)=2*(HL(4)+HL(5));B(4,5)=HL(5);
B(5,4)=HL(5);B(5,5)=2*(HL(5)+HL(6));B(5,6)=HL(6)-HL(7)^2/HL(6);
B(6,5)=HL(6);B(6,6)=3*HL(7)+2*HL(6)+HL(7)^2/HL(6);
QAL=zeros(6,3);
for i=1:NL;
    C(1,i)=6*(PL(i,3)/HL(2)+PL(i,1)/HL(1))-6*(1/HL(1)+1/HL(2))*(PL(i,1)+HL(1)*VL(i,1)+(HL(1)^2/3)*AL(i,1))-HL(1)*AL(i,1);
    C(2,i)=(6/HL(2))*(PL(i,1)+HL(1)*VL(i,1)+HL(1)^2*AL(i,1)/3)+6*PL(i,4)/HL(3)-6*(1/HL(2)+1/HL(3))*PL(i,3);
    C(3,i)=6*((PL(i,5)-PL(i,4))/HL(4)-(PL(i,4)-PL(i,3))/HL(3));
    C(4,i)=6*((PL(i,6)-PL(i,5))/HL(5)-(PL(i,5)-PL(i,4))/HL(4));
    C(5,i)=-(6/HL(6))*(PL(i,8)-VL(i,8)*HL(7)+HL(7)^2*AL(i,8)/3)-6*(1/HL(6)+1/HL(5))*PL(i,6)+(6/HL(5))*PL(i,5);%有争论
    %C(3,i)=(6/HL(4))*(PL(i,6)-VL(i,6)*HL(5)+HL(5)^2*AL(i,6)/3)-6*(1/HL(4)+1/HL(3))*PL(i,4)+(6/HL(3))*PL(i,3);
    %C(4,i)=-6*(1/HL(5)+1/HL(4))*(PL(i,6)-VL(i,6)*HL(5)+HL(5)^2*AL(i,6)/3)+6*PL(i,6)/HL(5)+6*PL(i,4)/HL(4)-AL(i,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;