auncertainty.m

来自「matlab进行不确度计算」· M 代码 · 共 22 行

function  [auncertainA,auncertainB,auncertain]=Auncertainty(X,Syserr,alpha)
%auncertain=Auncertainty(X,alpha)
%本函数用来求A类不确定度
%X为样本数据
%alpha为置信度

[m,n]=size(X);%输入样本数据X的大小
if (m==1||n==1)%X必须为一维行向量
 %  average=mean(X);%样本均值
%      ss=0;%样本与均值的差的平方和
%     for i=1:n
%         ss=ss+(X(i)-average)^2;
%     end
%     sigmahat=sqrt(ss/(n-1));%样本方差
    sigmahat = std(X);%样本方差
    auncertainA=sigmahat/sqrt(n)*tinv(1-alpha/2,n-1);%A类不确定度
    auncertainB = Syserr;       %B类不确定度
    auncertain = sqrt(auncertainA.^2+auncertainB.^2);   %综合不确定度
    
else
    disp('输入样本数据必须为一维行向量，请重试');     %输入不正确
end