从数据作手吧,修改一下数据,把明显不服从正态分布的数据作为异常剔除掉
对于样本X,(X为一个n*1的样本矩阵)执行下面的代码,剔除坏值:mark=0;while mark==1 n=length(x); x=x(find(x<3*std(x))); %超过3倍sigma的剔除 if n==length(x) mark=1; endend