grubbs准则 matlab_MATLAB-格拉布斯准则(MATLAB-Grubbs criterion)_0

%剔除可疑数据

%DESCRIPTION OF VALUBLES

%x-原始实验数据

%alpha-危险度(可取 0.05、0.025 或 0.01)

%avg-平均值

%st-标准差

%tx-计算T值

%T-T值表

%t-所选取的T值

%a-剔除后数据

%b-被剔除的数据

%VERSION V1.0 2013\3\17

%% 获得实验数据及危险度

x=input('输入实验数据');

alpha=input('输入危险度');

%% 排序得到数据容量计算平均值、标准差、计算T值

x=sort(x);

n=length(x);

avg=mean(x);

st=std(x);

tx=abs((x-avg)./st);

%% 根据危险度选取比较T值

T=[1.15 1.46 1.67 1.82 1.94 2.03 2.11 2.18 2.23 2.29;...

1.15 1.48 1.71 1.89

2.02 2.13 2.21 2.29 2.36 2.41;...

1.15 1.49 1.75 1.94

2.10 2.22 2.32 2.41 2.48 2.55]; switch alpha

case 0.05

t=T(1,n-2);

case 0.025

t=T(2,n-2);

case 0.01

t=T(3,n-2);

otherwise

error('WRONG VALUE OF ALPHA');

end

%% 比较剔除可疑数据

a=x(txb=x(tx>=t);

disp('剔除后数据为');

disp(a);

disp('被剔除的数据为');

disp(b);