wucha.m

峰值检测

%莱特准则判断粗大误差(适合次数较多时)
 a=[1.1 1.2 1.3 0.9 0.96 160];b=mean(a);
sum=0;
for i=1:length(a)
v(i)=abs(a(i)-b);
sum=sum+v(i).^2;
ruang=sqrt(sum/(length(a)-1));
if v(i)>3*ruang
a(i)=[];
end
end
%肖维捏准则
format short
n=10;
a=[1.1 1.2 1.3 0.9 0.96 1.60 2.3 1.9 1.2 1.3];b=mean(a);
sum=0;
fine=(2*n-1)/(4*n);
z=solve('1/2*erf(z/sqrt(2))=0.4750');
ruang=std(a);
V=ruang*z;
for i=1:length(a)
v(i)=abs(a(i)-b);
if v(i)>V
    a(i)=[];
end
end