m_linearity.m

来自「一组仪器性能分析程序」· M 代码 · 共 21 行

function linearity=M_linearity(Xu,Yu,Xd,Yd)
%输入为拟合直线的正反行程数组；输出为线性度
%设定输入正反行程数组为顺序对应的！！Xu=[1.1  2.1 4.9]  Xd=[1 1.9  5]
%计算方法按照书p16-17
x=(Xu+Xd)/2                         
y=(Yu+Yd)/2
p1=polyfit(x,y,1)
DltYl=y-polyval(p1,x)
MaxDltYl=max(abs(DltYl))

Xmin=min([Xu,Xd])
Xmax=max([Xu,Xd])
Yfs=abs(polyval(p1,Xmax)-polyval(p1,Xmin))
linearity=MaxDltYl/Yfs


%另外函数fliplr可以翻转数组
%添加的，该删除
%plot([Xu Xd],[Yu Yd],'.')
%plot(x,y,'*r')
%M_plotPolyVal(p1,Xmin,Xmax,1,1,1)