m_repeatability.m

一组仪器性能分析程序

function r=M_repeatability(Yuij,Ydij,X)
%输入p21页的公式(2.3.32)中的Yuij以及相应的Ydij
%其中，第一维 i=1 to N 个测点； 第二维 j=1 to M 次实验  (N行M列)

%这里暂时采用书上的方法，即使用一维的X（如书P23表中X=[2 4 6 8 10]）
%如果以后更改为输入Xuij和Xdij则直接用原始数据拟和直线，可能更合理。
%得到输入数据对，并拟和直线，yfs采用最小二乘直线求得
y=[Yuij Ydij]
for i=1:size(y,2)           %之前用length(y)返回最长维的长度，而我需要第二维m，实验的次数
    xij(:,i)=X
end
p=polyfit(xij,y,1)
mx=polyval(p,max(X))        %或max(max(xij))
mi=polyval(p,min(X))
yfs=mx-mi           
s=M_bessel(Yuij,Ydij)       %调用Bessel公式求s（书第21页公式2.3.34）
r=(3*s)/yfs                 %计算重复性r（书第21页公式2.3.36）