MATLAB 格拉布斯準則 剔除異常數(shù)據(jù) Grubbs criterion

%%Grubbs criterion %%輸入數(shù)據(jù)

x=input('輸入數(shù)據(jù)');

alpha=input('輸入alpha值(0.01或0.05)'); flag_1=1; flag_2=0; %% 計算

while flag_1==1 x=sort(x); n=length(x); avg=mean(x); st=std(x,1);

reg=abs((x-avg)./st); reg_max=max(reg); %% 確定g值

T=[1.153 1.463 1.672 1.822 2.285 2.331 2.371 2.409 2.58 2.603 2.624 2.644 2.745;...

1.155 1.492 1.749 1.944 2.55 2.607 2.659 2.705 2.912 2.939 2.963 2.987 3.103]; switch alpha case 0.05

g=T(1,n); case 0.01

g=T(2,n); otherwise

disp('輸入了錯誤的alpha值'); end

%% 比較確定異常值 ifreg_max>g flag_2=1; abn=x(reg>g);

disp('被剔除的數(shù)據(jù)為'); disp(abn);

x=x(reg<=g); else flag_1=0; end; end;

if flag_2==0

disp('沒有異常數(shù)據(jù)'); end;

1.938 2.032 2.443 2.475 2.663 2.681 2.097 2.22 2.747 2.785 3.009 3.029 2.11 2.176 2.501 2.532 2.698 2.714 2.323 2.41 2.821 2.954 3.049 3.068 2.2342.5572.732.4852.8843.085

總結

以上是生活随笔為你收集整理的grubbs准则 matlab_MATLAB-格拉布斯准则(MATLAB-Grubbs criterion)-M的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。