a.m

求解离散问题的pso aco 算法

tic

for j=1:345%%二折交叉验证开始
        dm=10000;%%给距离比较系数赋一个大值，便于比较
        for b=346:690%%对前345个样本计算与后345个样本的距离
            d=sqrt(pso(i,1)*(feature(b,1)-feature(j,1))*(feature(b,1)-feature(j,1))+pso(i,2)*(feature(b,2)-feature(j,2))*(feature(b,2)-feature(j,2))+pso(i,3)*(feature(b,3)-feature(j,3))*(feature(b,3)-feature(j,3))+pso(i,4)*(feature(b,4)-feature(j,4))*(feature(b,4)-feature(j,4))+pso(i,5)*(feature(b,5)-feature(j,5))*(feature(b,5)-feature(j,5))+pso(i,6)*(feature(b,6)-feature(j,6))*(feature(b,6)-feature(j,6))+pso(i,7)*(feature(b,7)-feature(j,7))*(feature(b,7)-feature(j,7))+pso(i,8)*(feature(b,8)-feature(j,8))*(feature(b,8)-feature(j,8))+pso(i,9)*(feature(b,9)-feature(j,9))*(feature(b,9)-feature(j,9))+pso(i,10)*(feature(b,10)-feature(j,10))*(feature(b,10)-feature(j,10))+pso(i,11)*(feature(b,11)-feature(j,11))*(feature(b,11)-feature(j,11))+pso(i,12)*(feature(b,12)-feature(j,12))*(feature(b,12)-feature(j,12))+pso(i,13)*(feature(b,13)-feature(j,13))*(feature(b,13)-feature(j,13))+pso(i,14)*(feature(b,14)-feature(j,14))*(feature(b,14)-feature(j,14))+pso(i,15)*(feature(b,15)-feature(j,15))*(feature(b,15)-feature(j,15)));
            if d<dm%%最短距离算法
                x(j,2)=x(b,1);
                dm=d;
            end
        end
        pso_rmse(i,1)=pso_rmse(i,1)+(x(j,1)-x(j,2))*(x(j,1)-x(j,2));%%离差在循环中累加，共345次；
        if x(j,1)==1&x(j,2)==1%%判定分类正误
            pso_rmse(i,3)=pso_rmse(i,3)+1;
        end
        if x(j,1)==1&x(j,2)==2%%判定分类正误
            pso_rmse(i,4)=pso_rmse(i,4)+1;
        end
        if x(j,1)==2&x(j,2)==1%%判定分类正误
            pso_rmse(i,5)=pso_rmse(i,5)+1;
        end
        if x(j,1)==2&x(j,2)==2%%判定分类正误
            pso_rmse(i,6)=pso_rmse(i,6)+1;
        end
    end
    for j=346:690
        dm=10000;%%给距离比较系数赋一个大值，便于比较
        for b=1:345%%对后345个样本计算与前345个样本的距离
            d=sqrt(pso(i,1)*(feature(b,1)-feature(j,1))*(feature(b,1)-feature(j,1))+pso(i,2)*(feature(b,2)-feature(j,2))*(feature(b,2)-feature(j,2))+pso(i,3)*(feature(b,3)-feature(j,3))*(feature(b,3)-feature(j,3))+pso(i,4)*(feature(b,4)-feature(j,4))*(feature(b,4)-feature(j,4))+pso(i,5)*(feature(b,5)-feature(j,5))*(feature(b,5)-feature(j,5))+pso(i,6)*(feature(b,6)-feature(j,6))*(feature(b,6)-feature(j,6))+pso(i,7)*(feature(b,7)-feature(j,7))*(feature(b,7)-feature(j,7))+pso(i,8)*(feature(b,8)-feature(j,8))*(feature(b,8)-feature(j,8))+pso(i,9)*(feature(b,9)-feature(j,9))*(feature(b,9)-feature(j,9))+pso(i,10)*(feature(b,10)-feature(j,10))*(feature(b,10)-feature(j,10))+pso(i,11)*(feature(b,11)-feature(j,11))*(feature(b,11)-feature(j,11))+pso(i,12)*(feature(b,12)-feature(j,12))*(feature(b,12)-feature(j,12))+pso(i,13)*(feature(b,13)-feature(j,13))*(feature(b,13)-feature(j,13))+pso(i,14)*(feature(b,14)-feature(j,14))*(feature(b,14)-feature(j,14))+pso(i,15)*(feature(b,15)-feature(j,15))*(feature(b,15)-feature(j,15)));
            if d<dm%%最短距离算法
                x(j,2)=x(b,1);
                dm=d;
            end
        end
        pso_rmse(i,1)=pso_rmse(i,1)+(x(j,1)-x(j,2))*(x(j,1)-x(j,2));%%离差在循环中累加，共345次，与上面
        %%的345次累加完成对690个样本的离差和累加
        if x(j,1)==1&x(j,2)==1%%判定分类正误
            pso_rmse(i,3)=pso_rmse(i,3)+1;
        end
        if x(j,1)==1&x(j,2)==2%%判定分类正误
            pso_rmse(i,4)=pso_rmse(i,4)+1;
        end
        if x(j,1)==2&x(j,2)==1%%判定分类正误
            pso_rmse(i,5)=pso_rmse(i,5)+1;
        end
        if x(j,1)==2&x(j,2)==2%%判定分类正误
            pso_rmse(i,6)=pso_rmse(i,6)+1;
        end
    end

toc