grubbsmain.m

用格拉布斯坏值提出理论于剔除统计样本样本中坏值

%格拉布斯剔除准则
format long
fid=fopen('E:\lunwendata\singlex39km.txt','r');
[X,count]=fscanf(fid,'%g');
count;
lambda=input('请输入N=50时候的lambda值:\n');
fclose all;
fidhao=fopen('E:\毕业论文程序temp文件\singlex39kmhao.txt','w');
fidhuai=fopen('E:\毕业论文程序temp文件\singlex39kmhuai.txt','w');
bianhao=1:50:count;
for i=1:length(bianhao)-1
    [hao,huai]=grubsn(lambda,X(bianhao(i):bianhao(i)+49,1));
    for j=1:length(hao)
        fprintf(fidhao,'%4.4f\n',hao(j));
    end
    for j=1:length(huai)
        fprintf(fidhuai,'%4.4f\n',huai(j));
    end
end
count;
rem(count,50);
fprintf('还有N=%g个数值没有进行剔除:\n',rem(count,50));
lambdan=input('请输入时候的lambda值:\n');
    [hao,huai]=grubsn(lambdan,X(count-rem(count,50)+1:count,1));
    for j=1:length(hao)
        fprintf(fidhao,'%4.4f\n',hao(j));
    end
    for j=1:length(huai)
        fprintf(fidhuai,'%4.4f\n',huai(j));
    end
fclose all;