基于粒子群算法的投影寻踪法matlab代码

function [ Q ] = touying（a）
%touying为投影目标函数
load X.txt;
n=size（X2）;
z=a*X;  %指标数为列数
for k=1:100  %100为种群规模
    for i=1:n
        Z（ki）=z（ki）-mean（z（k:））;
    end
end
for k=1:100  %100为种群规模
    S（k）=norm（Z（k:）2）/sqrtm（n-1）;
    D（k）=0;
    for i=1:n