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