粗大误差处理.doc

莱以特准则
load a.txt
while(1)
i=1:length(a);
n=length(a);
v(i)=a(i)-mean(a);
bzc=sqrt(sum(v(i).^2)/(length(a)-1));
d=3*bzc;
[maxv,I]=max(abs(v(i)));
if maxv>d
fprintf('cdw is %f\n',a(I));
a(I)=[];
else break;
end
end
cdw is 29.520000
cdw is 28.400000
罗曼诺夫斯基准则
load a.txt
n=input('please input n:\n');
xzd=input('please input xzd:\n');
switch xzd
case xzd==0.05
x=1;
otherwise
x=2;
end
b=a(n);
a(n)=[];
while(1)
c=mean(a);
i=1:length(a);
n=length(a);
v(i)=a(i)-mean(a);
bzc=sqrt(sum(v(i).^2)/(length(a)-1));
k=[4.97 3.56 3.04 2.78 2.62 2.51 2.43 2.37 2.33 2.29 2.26 2.24 2.22 2.20 2.18 2.17 2.16 2.15 2.14 2.13 2.12 2.11 2.10 2.10 2.09 2.09 2.08;11.46 6.53 5.04 4.36 3.96 3.71 3.54 3.41 3.31 3.23 3.17 3.12 3.08 3.04 3.01 3.00 2.95 2.93 2.91 2.90 2.88 2.86 2.85 2.84 2.83 2.82 2.81];
g=k(x,n-2);
f=g*bzc;
e=abs(b-c);
if e>f
fprintf('cdw is %f\n',b);
else fprintf('wcdw\n');
end
break;
end
please input n:
4
please input xzd:
0.05
cdw is 29.520000

格罗布准则
load a.txt
xzd=input('please input xzd:\n');
switch xzd
case xzd==0.05
x=1;
otherwise
x=2;
end
while(1)
i=1:length(a);
n=length(a);
v(i)=ai)-mean(a);
bzc=sqrt(sum(v(i).^2)/(length(a)-1));
a=sort(a);
c=min(a);
d=max(a);
g1=(mean(a)-c)/bzc;
gn=(d-mean(a))/bzc;
T=[1.15 1.46 1.67 1.82 1.94 2.03 2.11 2.18 2.23 2.28 2.33 2.37 2.41 2.44 2.48 2.50 2.53 2.56 2.58 2.60 2.62 2.64 2.66 2.74 2.81 2.87 2.96 3.17;1.16 1.49 1.75 1.94 2.10 2.22 2.32 2.41 2.48 2.55 2.61 2.66 2.70 2.75 2.78 2.82 2.85 2.88 2.91 2.94 2.96 2.99 3.01 3.10 3.18 3.24 3.34 3.59];
g=T(x,n-2);
if g1>g
fprintf('cdw is %f\n',c);
a(1)=[];
end
if gn>g
fprintf('cdw is %f\n',d);
a(n)=[];
else break;
end
end
please input xzd:
0.05
cdw is 29.520000
cdw is 28.400000

狄克松准则
load a.txt
xzd=input('please input xzd:\n');
switch xzd
case xzd==0.05
x=1;
otherwise
x=2;
end
while(1)
a=sort(a);
i=1:length(a);
n=length(a);
T=[0.988 0.889 0.780 0.698 0.637 0.683 0.635 0.597 0.679 0.642 0.651 0.641 0.616 0.595 0.577 0.561 0.547 0.535 0.524 0.514 0.505 0.497 0.489;0.341 0.765 0.642 0.560 0.507 0.554 0.512 0.477 0.576 0.546 0.521 0.546 0.525 0.507 0.490 0.475 0.462 0.450 0.440 0.430 0.421 0.413 0.406];
if n<=7
r101=(a(n)-a(n-1))/(a(n)-a(1));
r102=(a(1)-a(2))/(a(1)-a(n));
g=T(x,n-2);
if r101>g
fprintf('cdw is %f\n',a(n));
a(n)=[];
else fprintf('wcdw\n');
end
if r102>g
fprintf('cdw is %f\n',a(1));
a(1)=[];
else break;
end
end
if (n<=8)&&(n<=10)
r111=(a(n)-a(n-1))/(a(n)-a(2));
r112=(a(1)-a(2))/(a(1)-a(n-1));
g=T(x,n-2);
if r111>g
fprintf('cdw is %f\n',a(n));
a(n)=[];
else fprintf('wcdw\n');
end
if r112>g
fprintf('cdw is %f\n',a(1));
a(1)=[];
else break;
end
end
if (n>=11)&&(n<=13)
r211=(a(n)-a(n-2))/(a(n)-a(2));
r212=(a(1)-a(3))/(a(1)-a(n-1));
g=T(x,n-2);
if r211>g
fprintf('cdw is %f\n',a(n));
a(n)=[];
else fprintf('wcdw\n');
end
if r212>g
fprintf('cdw is %f\n',a(1));
a(1)=[];
else break;
end
end
if n>=14
r221=(a(n)-a(n-2))/(a(n)-a(3));
r222=(a(1)-a(3))/(a(1)-a(n-2));
g=T(x,n-2);
if r221>g
fprintf('cdw is %f\n',a(n));
a(n)=[];
else fprintf('wcdw\n');
end
if r222>g
fprintf('cdw is %f\n',a(1));
a(1)=[];
else break;
end
end
end
please input xzd:
0.05
cdw is 29.520000
cdw is 28.400000

THANKS !!!

致力为企业和个人提供合同协议，策划案计划书，学习课件等等
打造全网一站式需求

欢迎您的下载，资料仅供参考

-可编辑修改-