syserror.m

来自「matlab 电子测量 包括误差分析的很多函数」· M 代码 · 共 31 行

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