格拉布斯准则判断异常数据代码.zip

clear clc G1=[1.15,1.46,1.67,1.82,1.94,2.03,2.11,2.18,2.23,2.29,2.33,2.37,2.41,2.44,2.47,2.50,2.53,2.56] G2=[1.15,1.49,1.75,1.94,2.10,2.22,2.32,2.41,2.48,2.55,2.61,2.66,2.70,2.74,2.78,2.82,2.85,2.85] grubbstab=input('Please input measuring datas:') g=input('please input confidence level ') if g==0.99 G=G2 else G=G1 end meangrubbstab =mean(grubbstab) stdgrubbstab=std(grubbstab) len=length(grubbstab) for i=1:len A(i)=abs(grubbstab(i)-meangrubbstab) end Max1=max(A) location1=find(A==Max1) if abs(grubbstab(location1)-meangrubbstab)< G(len-2)*stdgrubbstab disp('no wide error!') else grubbstab(location1)=[] A(location1)=[] meangrubbstab =mean(grubbstab) stdgrubbstab=std(grubbstab) Max2=max(A) location2=find(A==Max2) len=length(grubbstab) if abs(grubbstab(location2)-meangrubbstab)< G(len-2)*stdgrubbstab disp('no more wide error!') else grubbstab(location2)=[] end end