%PCA的T方控制限和SPE控制限求算函數
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%輸入X為標準化后的樣本矩陣（行向量樣本），latent為主元分析的，alpha為檢驗水平，k為主元個數
%輸出ts_ctrl為T方控制限，spe_ctrl為SPE控制限
function [ts_ctrl,spe_ctrl]=PCAThrd(X,latent,alpha,kp)
[n,m]=size(X);
%利用F分布求T方控制限
ts_ctrl=kp*(n-1)/(n-kp)*finv(alpha,kp,n);
% 利用正態分布求SPE控制限
theta=zeros(3,1);
for i=1:3
    for j=kp+1:m
        theta(i)=theta(i)+latent(j)^(i);
    end
end
h0=1-2*theta(1)*theta(3)/(3*theta(2)^2);
spe_ctrl=theta(1)*(norminv(alpha)*(2*theta(2)*h0^2)^0.5/theta(1)+1+theta(2)*h0*(h0-1)/theta(1)^2)^(1/h0);