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;