elcf.m

matlab 电子测量 包括误差分析的很多函数

function y=elcf(index,n,p)
% index 1莱特准则 2肖维纳准则 
% 3格布拉斯准则 p=1 0.95;p=2 0.99
% eliminate coefficient
if index==1  % 莱特准则
    y=3;
else if index==2 % 肖维纳准则
ch=[1.65 1.73 1.79 1.86 1.92 1.96 2.00 2.04 2.07 2.10 2.13 2.16 2.18 2.20 2.22 2.24 2.26 2.28 2.30 2.32 2.33 2.34 2.35 2.37 2.38 2.39 2.45 2.50 2.58 2.64 2.74 2.81 3.02];
if n<=30
    y=ch(n-4);
else  switch n
        case 35
            y=2.45;
        case 40
            y=2.50;
        case 50
            y=2.58;
        case 60
            y=2.64;
        case 80
            y=2.74;
        case 100
            y=2.81;
        case 200
            y=3.02;
    end
end

    else if index==3 %格布拉斯准则
gb1=[1.15 1.46 1.67 1.82 1.94 2.03 2.11 2.18 2.23 2.29 2.33 2.37 2.41 2.44 2.47 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];
gb2=[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.71 2.75 2.79 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];
gb=[gb1;gb2];
if n<=25
    y=gb(p,n-2);
else if n==30
        y=gb(p,24);
        else if n==35
                y=gb(p,25);
            else if n==40
                    y=gb(p,26);
                    else if n==50 
                            y=gb(p,27);
                        else if n==100
                                y=gb(p,28);
                        end
                end
            end
    end
end
    end
    end
    end
end