Page 64 - Algorithm Collections for Digital Signal Processing Applications using MATLAB

52 Chapter 1

for i=1:1:8
temp=P_cur(:,position(i));
temp=temp/sum(temp);
[value,index]=min(temp);
tempant(position(i))=notation(index);
P_cur(index,:)=[];
notation(index)=[];
end
ANTS=[ANTS;tempant];
end
c1=computecost(ANTS(1,:),matA,matB);
c2=computecost(ANTS(2,:),matA,matB);
c3=computecost(ANTS(3,:),matA,matB);
c4=computecost(ANTS(4,:),matA,matB);
end
_____________________________________________________________
computecost.m
function [s]=computecost(A,P,Q)
s=0;
for i=1:1:8
s=s+P(max(i,A(i)),min(i,A(i))).*Q(max(i,A(i)),min(i,A(i)));
end