资源简介

详尽的PCA算法对TE过程的故障诊断分析,内含有部分故障数据,可供调试使用,对做故障诊断的非常有帮助。

资源截图

代码片段和文件信息

%%TE过程的传统主元分析在Matlab中的仿真程序
%建立模型:
%载入模型数据
Xtrain = load(‘d0_te.dat‘); 
Xtrain = double(Xtrain);

%载入测试数据
Xtest = load(‘d01.dat‘);
Xtest = double(Xtest);

%标准化处理
X_mean = mean(Xtrain);                             
X_std = std(Xtrain);                               
[X_rowX_col] = size(Xtrain);                   
% for i = 1:X_col
%     Xtrain(:i) = (Xtrain(:i) - X_mean(i)./X_std(i));
%     Xtest(:i) = (Xtest(:i) - X_mean(i)./X_std(i));
% end     
Xtrain = (Xtrain - repmat(X_meanX_row1))./repmat(X_stdX_row1);

%求协方差矩阵,并对协方差矩阵进行特征分解
sigmaXtrain = cov(Xtrain);
[Tlamda] = eig(sigmaXtrain);      
% disp(‘特征根(由小到大)‘);
% disp(lamda);
% disp(‘特征向量:‘);
% disp(T);                                            

%取对角元素,即lamda值,并上下反转使其从大到小排列,主元个数初值为1,若累计贡献率小于85%则增加主元个

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件      416052  2006-10-27 21:22  d00.dat
     文件      799680  2006-10-27 21:36  d00_te.dat
     文件      399840  2006-10-27 21:34  d01.dat
     文件      799680  2006-10-27 21:41  d01_te.dat
     文件      225140  2014-02-25 10:40  d02.dat
     文件      799680  2006-10-27 21:47  d02_te.dat
     文件      399840  2006-10-27 21:49  d03.dat
     文件      799680  2006-10-27 21:52  d03_te.dat
     文件      399840  2006-10-27 21:54  d04.dat
     文件      799680  2006-10-27 22:12  d04_te.dat
     文件      399840  2006-10-27 22:14  d05.dat
     文件      799680  2006-10-27 22:22  d05_te.dat
     文件      399840  2006-10-27 22:24  d06.dat
     文件      799680  2006-10-27 22:27  d06_te.dat
     文件      399840  2006-10-27 22:27  d07.dat
     文件      799680  2006-10-28 08:19  d07_te.dat
     文件      399840  2006-10-28 08:19  d08.dat
     文件      799680  2006-10-28 08:19  d08_te.dat
     文件      399840  2006-10-28 08:19  d09.dat
     文件      399840  2006-10-28 08:20  d10.dat
     文件        3578  2014-09-10 16:58  PCA_TE_myself.m

评论

共有 条评论