投影寻踪法matlab.zip

投影寻踪( projection pursuit，PP) 方法属于直接由样本数据驱动的探索性数据分析方法，是美国科学家 Kruscal 于 20 世纪 70 年代提出的，在高维性、非线性、非正态数据分析处理方面有独到之处，运用matlab编程，其计算步骤如下：( 1) 指标体系无量纲化； （2） 构造投影指标函数； （3）构造投影目标函数；( 4) 确定最佳投影方向；( 5) 确定投影值

function Qa=Project_Pursuit（XaAlpha）
%% 输入参数列表
% X 本指标矩阵，n*p的矩阵，每一行为一个样本
% Xij表示第i行第j列指标，X是否已经均一化均可
% a 投影向量，1*p的矩阵，元素取值范围-1~1，要求其元素平方和等于1
% Alpha 窗口半径系数，典型取值0.1
X=[0.02 0.05 0.75 0.18;
   0.03 0.79 0.11 0.07;
   0.04 0.49 0.33 0.14;
   0.04 0.12 0.69 0.15;
   0.04 0.32 0.48 0.16;
   0.05 0.44 0.34 0.17;
   0.03 0.72 0.15 0.10;
   0.14 0.42 0.25 0.19;
   0.04 0.59 0.25 0.12;
   0.05 0.15 0.60 0.20;
   0.12 0.50 0.22 0.16;
   0.03 0.83 0.09 0.05;
   0.05 0.38 0.41 0.16;
   0.03 0.12 0.51 0.34;
   0.02 0.86 0.08 0.04;
   0.08 0.64 0.18 0.10;
   0.04 0.32 0.48 0.16;
   0.03 0.77 0.13 0.07;
   0.09 0.56 0.21 0.14]; 
a=[a1a2a3a4a5a6a7a8a9a10a11a12a13a14a15a16a17a18a19];
Alpha=0.1;
%% 输出参数列表
% Qa 投影指标函数
%% 第0步：对a的预处理
b=sqrt（sum（a.^2））;
a=a/b;
%% 第一步：指标体系无量纲化   % 缺少负指标的归一化处理
[np]=size（X）;
x=zeros（np）;
Xjmax=max（X）;
Xjmin=min（X）;
for i=1:n
    x（i:）=（X（i:）-Xjmin）./（Xjmax-Xjmin）;
end
%% 第二步：构造投影指标函数
Z=zeros（n1）;
for i=1:n
    Z（i）=sum（a.*x（i:））;
end
%%第三步：计算投影指标函数
%% 计算类间类内矩阵散度
meanZ=mean（Z）;
Sa=0;
for k=1:n
    sa=（Z（i）-meanZ（i））.^2;
    Sa=Sa+sa;
    Sa=sqrt（Sa/n）;
end
R=Alpha*Sa;   %窗口半径
Da=0;
for k=1:n
    rik=abs（Z（i）-Z（k））;
    if R>rik
        Da=Da+rik;
    end
end
Qa=Sa*Da;

    

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-05-16 11:49  投影寻踪法（matlab）\
     文件       43520  2019-05-12 19:08  投影寻踪法（matlab）\matlab和投影寻踪.doc
     目录           0  2019-05-16 11:49  投影寻踪法（matlab）\PP-matlab\
     文件         753  2010-04-20 15:22  投影寻踪法（matlab）\PP-matlab\diaoyong.m
     文件         118  2010-03-30 21:19  投影寻踪法（matlab）\PP-matlab\Feasibility.m
     文件        3877  2010-03-30 21:50  投影寻踪法（matlab）\PP-matlab\RAGA.m
     文件         441  2010-03-30 21:19  投影寻踪法（matlab）\PP-matlab\Target.m
     文件        1452  2019-05-12 19:16  投影寻踪法（matlab）\Project_Pursuit.m
     文件        1155  2019-05-12 19:27  投影寻踪法（matlab）\Untitled2.m