function y=syserror(x)
% 系統誤差判別
y1=mean(x);  % 均值
y2=std(x);   % 試驗方差
vk=x-y1;     % 殘差
s=0;
[m n]=size(x);
for k=1:n-1
    s=s+vk(k)*vk(k+1);
end
s=abs(s);
t=sqrt(n-1)*y2.^2;
if s<=t       % Abbe-Helmert 判據
    disp('測量中不存在變值系差'); % y=0;
else
    s1=0;s2=0;
    for k=1:n/2
        s1=s1+vk(k);
        s2=s2+vk(n+1-k);
    end
    ss=abs(s1-s2);
    if ss>=max(vk)  %馬爾科夫判據
        disp('存在累進性系差'); % y=10;
    else
        disp('存在周期性系差'); % y=11;
    end
end