基于MATLAB的格拉布斯准则

%a=[12345420212345432123484321];
%a=[12345678910987654321234567891098765432123456789109876520321234567891098765432123456789409876543212345678910]; 
a=[12345678910987654321234567891098765432123456789109876520321234567891098765432123456789409876543212345678910]; 
alpha=input（‘（请输入alpha的值0.01或0.05）‘）;
subplot（121）;stem（a‘filled‘）;
xlabel（‘序列个数n‘）;
ylabel（‘幅值‘）;
title（‘原序列‘）;
%axis（[0 25 0 30]）;
axis（[0 100 0 50]）;
flag_1=1;
flag_2=0;
while flag_1==1
n=length（a）;
avg=mean（a）;
st=std（a1）; 
reg=abs（（a-avg）./st）;
reg_max=max（reg）;
T=[0 0 0 1.153 1.463 1.672 1.822 1.938 2.032 2.11 2.176 2.234 2.285 2.331 2.371 2.409 2.443 2.475 2.501 2.532 2.557 2.58 2.603 2.624 2.644 2.663 2.681 2.698 2.714 2.73 2.745 2.759 2.773 2.786 2.799 2.811 2.823 2.835 2.846 2.857 2.866 2.877 2.887 2.896 2.905 2.914 2.923 2.931 2.940 2.948 2.956 2.964 2.971 2.978 2.986 2.992 3.000 3.006 3.013 3.019 3.025 3.032 3.037 3.044 3.049 3.055 3.061 3.066 3.0