grubsn.m

来自「用格拉布斯坏值提出理论于剔除统计样本样本中坏值」· M 代码 · 共 22 行

function [hao,huai]=grubsn(lambda,yuandata)%参数为lambda和样本
[hang,lie]=size(yuandata);
D=var(yuandata,1);
u=mean(yuandata,1);
delta=sqrt(D);
zhunze=lambda*delta;%判断准则

nhao=find(abs(yuandata(:,1)-u)<=zhunze);%好值序号
len=length(nhao);
hao=zeros(len,1);
for i=1:len
   hao(i,1)=yuandata(nhao(i),1);
end
mhuai=find(abs(yuandata(:,1)-u)>zhunze);%坏值序号
huai=zeros(hang-len,1);
for i=1:hang-len
  huai(i,1)=yuandata(mhuai(i),1);
end