1、莱以特准则load a.txtwhile(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 maxvdfprintf(cdw is %fn,a(I);a(I)=;else break;endendcdw is 29.520000cdw is 28.400000罗曼诺夫斯基准则load a.txtn=input(please input n:n);xzd=input(please input xzd:n); switch x
2、zdcase xzd=0.05x=1;otherwise x=2;endb=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
3、.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 effprintf(cdw is %fn,b);else fprintf(wcdwn);endbreak;endplease input n:4please input xzd:0.05cdw is 29.520000格罗布准则load a.txtxzd=input(please input x
4、zd:n); switch xzdcase xzd=0.05x=1;otherwise x=2;endwhile(1)i=1:length(a);n=length(a);v(i)=a(i)-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
5、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 g1gfprintf(cdw is %fn,c);a(1)=;endif gngfprintf(cdw is %fn,d);a(n)=;else break;endendplease input xzd
6、:0.05cdw is 29.520000cdw is 28.400000狄克松准则load a.txtxzd=input(please input xzd:n); switch xzdcase xzd=0.05x=1;otherwise x=2;endwhile(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
7、.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 ngfprintf(cdw is %fn,a(n);a(n)=;else fprintf(wcdwn);endif r102gfprintf(cdw is %fn,a(1);a(1)=;else break;endendif (ngfprintf(cdw is %fn,a(n);a(n)=;el
8、se fprintf(wcdwn);endif r112gfprintf(cdw is %fn,a(1);a(1)=;else break;endendif (n=11)a(n)=;else fprintf(wcdwn);endif r212gfprintf(cdw is %fn,a(1);a(1)=;else break;end endif n=14r221=(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 r221gfprintf(cdw is %fn,a(n);a(n)=;else fprintf(wcdwn);endif r222gfprintf(cdw is %fn,a(1);a(1)=;else break;endendendplease input xzd:0.05cdw is 29.520000cdw is 28.400000
