%格拉布斯剔除準(zhǔn)則
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:\畢業(yè)論文程序temp文件\singlex39kmhao.txt','w');
fidhuai=fopen('E:\畢業(yè)論文程序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個數(shù)值沒有進(jìn)行剔除:\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;