gajisuankekaoduqiaoshi.m

来自「粒子群算法用于桥式系统的可靠性优化。算法简单可靠方便。」· M 代码 · 共 26 行

function [R,eval]=GAjisuankekaoduqiaoshi(R,options)
% kekaodu=R(1)*R(4)+R(2)*R(5)+R(2)*R(3)*R(4)+R(1)*R(3)*R(5)+2*R(1)*R(2)*R(3)*R(4)*R(5)...
%     -R(2)*R(3)*R(4)*R(5)-R(1)*R(3)*R(4)*R(5)-R(1)*R(2)*R(4)*R(5)-R(1)*R(2)*R(3)*R(5)-R(1)*R(2)*R(3)*R(4);
% num=length(R);cost=0;
% for i=1:5,
%     cost=cost+exp(0.0003/(1-R(i)));
% end
% if kekaodu<0.99
%     eval=-10;
% else
%     eval=100-cost;
% end
kekaodu=R(1)*R(4)+R(2)*R(5)+R(2)*R(3)*R(4)+R(1)*R(3)*R(5)+2*R(1)*R(2)*R(3)*R(4)*R(5)...
    -R(2)*R(3)*R(4)*R(5)-R(1)*R(3)*R(4)*R(5)-R(1)*R(2)*R(4)*R(5)-R(1)*R(2)*R(3)*R(5)-R(1)*R(2)*R(3)*R(4);
while(kekaodu<0.99)
    for i=1:5,
        R(i)=rand;
    end
    kekaodu=R(1)*R(4)+R(2)*R(5)+R(2)*R(3)*R(4)+R(1)*R(3)*R(5)+2*R(1)*R(2)*R(3)*R(4)*R(5)...
    -R(2)*R(3)*R(4)*R(5)-R(1)*R(3)*R(4)*R(5)-R(1)*R(2)*R(4)*R(5)-R(1)*R(2)*R(3)*R(5)-R(1)*R(2)*R(3)*R(4);
end
cost=0;
for i=1:5,
    cost=cost+exp(0.0003/(1-R(i)));
end
eval=100-cost;