PCA主成分分析代码，函数形式

森之张卫东

1. function [contribution,T2,SPE,t2cl,s_cl] =  PCA_model(Xtrain,Xtest)
   
2. X_mean = mean(Xtrain);               
   
3. X_std = std(Xtrain);                  
   
4. [X_row ,X_col]= size(Xtrain);
   
5. for i = 1:X_col
   
6.     Xtrain(:,i) = (Xtrain(:,i)-X_mean(i))./X_std(i);
   
7.     Xtest(:,i) = (Xtest(:,i)-X_mean(i))./X_std(i);
   
8. end
   
9. [U,S,V]=svd(Xtrain./sqrt(size(Xtrain,1)-1),0);        
   
10. D= S^2;
    
11. lamda=diag(D);
    
12. num_pc=1;
    
13. while sum(lamda(1:num_pc))/sum(lamda)<0.9
    
14.     num_pc=num_pc+1;
    
15. end
    
16. D=diag(lamda);
    
17. P=V(:,1:num_pc);
    
18. &#91;a,b&#93;=size(Xtest);
    
19. &#91;r,y&#93;=size(P*P');
    
20. I=eye(r,y);
    
21. e=Xtest*(I-P*P');
    
22. for i=1:a  
    
23. T2(i)=Xtest(i,:)*P*inv(D(1:num_pc,1:num_pc))*P'*Xtest(i,:)';
    
24. end
    
25. for l=1:a
    
26.     SPE(l)=e(l,:)*e(l,:)';
    
27. end
    
28. 
    
29. for j=1:b
    
30.     contribution(j)=(norm(e(:,j)))^2;
    
31. end
    
32. 
    
33. t2cl=num_pc*(X_row-1)*(X_row+1)*icdf('f',0.99,num_pc,X_row-num_pc)/(X_row*(X_row-num_pc));
    
34. for i=1:3
    
35.      theta(i)=trace((D(num_pc+1:X_col,num_pc+1:X_col))^i);        
    
36. end
    
37. % 另一种SPE控制线算法
    
38. % h=(theta(1)^2)/theta(2);
    
39. % g=theta(2)/theta(1);
    
40. % conf=0.95;                              
    
41. % df=round(h);                           
    
42. % delta2a1=g*pinv(df,2);
    
43. 
    
44. h0=1-2*theta(1)*theta(3)/(3*theta(2)^2);
    
45. ca=icdf('norm',0.99,0,1);
    
46. s_cl=theta(1)*(ca*sqrt(2*theta(2)*h0^2)/theta(1)+1+theta(2)*h0*(h0-1)/theta(1)^2)^(1/h0);

复制代码

森之张卫东

上传一个测试程序。