?? pcathrd.m
字號:
%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);
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -