m_bessel.m

一组仪器性能分析程序

function S=M_bessel(Yuij,Ydij)
%输入p21页的公式(2.3.32)中的Yuij以及相应的Ydij
%其中，第一维 i=1 to N 个测点； 第二维 j=1 to M 次实验  (N行M列)
%输出p21页的公式(2.3.34)中的S
  if any(size(Yuij)~=size(Ydij))    %size(Yuij)~=size(Ydij)返回一向量，各个位代表相应的维是否长度相同
      return                        %ANY    True if any element of a vector is a nonzero number
  end                               %只要有一维不等长度，就返回true
  %上面的必须要any(),否则判断含有0值(虽然有非0值）的向量返回false
  S=sqrt(Si(Yuij)+Si(Ydij))
end

function Si=Si(Yij)
  n=size(Yij,1)
  m=size(Yij,2)
  YiP=mean(Yij,2)
  dlt=zeros(n,m)
  %两种方式：按行相减和按列相减。
  %若按行相减程序如下（按照测点个数）
  for i=1 :1: n
      dlt(i,:)=Yij(i,:)-YiP(i,:)          %相当于p21页Yuij-YuiP的矩阵
  end
  %若按列相减程序如下（按照实验次数）
  %for i=1 :1: m
  %    dlt(:,i)=Yij(:,i)-YiP
  %end  
  dlt2=dlt.^2
  Sui2=sum(dlt2,2)/(m-1)              %相当于p21页公式（2.3.32）
  Si=sum(Sui2)/(2*n)                  %相当于p21页公式（2.3.34）中根号下的前半部分
end