ICA故障检测与诊断

训练集采用TE过程正常状态参数，测试集采用TE过程故障10状态参数，采用ICA方法对TE过程进行故障检测与诊断

%对已知的密度函数求取控制限
function ConInt=ComCon（densityxmeshalpha）	%density 密度矩阵 
							%xmesh为密度对应的x值的行向量（或矩阵）
							%alpha 置信度
ConInt=zeros（12）;			%输出ConInt为控制限
n=length（density）;			%计算每一列的控制限

x1=min（xmesh）;
x2=max（xmesh）;
dx=（x2-x1）/（n-1）;

%黄金分割搜索从开始到某一点的面积为（1-alpha）/2以求取下分位数;
while（1）
     x3=x1+（sqrt（5）-1）/2*（x2-x1）;
     t=x1:dx:x3;
%    s=trapz（tdensity（1:length（t）））;
     s=trapz（density（1:length（t）））*dx;
     if abs（x2-x1）<1.0e-3 
%       if s-alpha<1.0e-4
        ConInt（1）=x3; 
     break;
%    elseif s     elseif s<（1-alpha）/2
            x1=x3;
     else
          x2=x3;
     end
end
%     求上分为数，从开始到某一点的面积为1-（1-alpha）/2;
x1=min（xmesh）;
x2=max（xmesh）;
while（1）
      x3=x1+（sqrt（5）-1）/2*（x2-x1）;
      t=x1:dx:x3;
%     s=trapz（tdensity（1:length（t）））;
      s=trapz（density（1:length（t）））*dx;
      if abs（x2-x1）<1.0e-3 
         ConInt（2）=x3; 
         break;
      elseif s<（1+alpha）/2
             x1=x3;
      else
           x2=x3;
      end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2014-05-27 17:11  ICA\
     目录           0  2014-05-04 20:38  ICA\ICA\
     文件        1112  2014-05-27 15:53  ICA\ICA\ComCon.m
     文件        1343  2014-06-04 22:49  ICA\ICA\ICA_monitor.m
     文件        1462  2014-06-04 22:49  ICA\ICA\ICA_normal.m
     文件         199  2014-04-23 20:35  ICA\ICA\control_limit.m
     文件      399840  2002-04-23 00:20  ICA\ICA\d10.dat
     文件      799680  2006-10-28 08:20  ICA\ICA\d10_te.dat
     文件        1651  2014-06-04 22:36  ICA\ICA\main_n.m
     文件         293  2014-04-28 20:30  ICA\ICA\sort_IC.m
     文件         191  2014-05-27 16:07  ICA\ICA\variable_c.m
     目录           0  2014-05-27 17:11  ICA\ICA_modified\
     文件        1112  2014-05-27 15:53  ICA\ICA_modified\ComCon.m
     文件        1166  2014-05-27 21:27  ICA\ICA_modified\ICA_monitor.m
     文件        1592  2014-05-27 21:20  ICA\ICA_modified\ICA_normal.m
     文件         199  2014-04-23 20:35  ICA\ICA_modified\control_limit.m
     文件      399840  2002-04-23 00:20  ICA\ICA_modified\d10.dat
     文件      799680  2006-10-28 08:20  ICA\ICA_modified\d10_te.dat
     文件        1496  2014-05-27 21:27  ICA\ICA_modified\main_n.m
     文件         293  2014-04-28 20:30  ICA\ICA_modified\sort_IC_VRE.m
     文件         313  2014-05-27 19:28  ICA\ICA_modified\sort_IC_cum.m
     文件         293  2014-04-28 20:30  ICA\ICA_modified\sort_IC_negentropy.m
     文件         187  2014-05-27 21:08  ICA\ICA_modified\variable_c.m