grubsn.asv

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

function [hao,huai]=grubsn(lambda,yuandata)%参数为lambda和样本
[hang,lie]=size(yuandata);
D=var(X,1);
delta=sqrt(D);
zhunze=lambda*delta
nhao=find(yuandata(:,1)<=zhunze);%好值序号
for i=1:
hao(i,1)=yuandata(hao(i),1);

mhuai=find( yuandata(:,1)>zhunze);%坏值序号
huai(i,1)=yuandata(mhuai(i),1);

hao=nhao;