设传感器系统静态特性为:
下图为该传感器静态特性曲线及用最小二乘法与端基线法拟合的直线,其中
y=1+2x+0.02x^2
通过最小二乘法拟合的直线线性度为14.8%,通过端基线法拟合的直线线性度为19.24%。图2为MATLAB脚本编辑界面、线性度运算结果及相关曲线显示。
maltab代码如下,仅供参考
clear;clc
%原曲线
x=linspace(0,10,100);
S=[1,2,0.02,0];
y=S(1)+x*S(2)+S(3)* (x.^3)+S(4)*(x.^5);
plot(x,y,'Color','black','LineWidth',1)
%最小二乘拟合
N=length(x);
k=(sum(y.*x)-N*mean(y)*mean(x))/(sum(x.^2)-N*mean(x)^2);
b=mean(y)-k*mean(x);
y1_line=k*x+b;
hold on
plot(x,y1_line,'Color','red','LineWidth',1)
%线性度
Lm1=max(abs(y-y1_line));
yFS1=max(max(y))-min(min(y));
L1=Lm1/yFS1;
disp(L1);
%端基线法拟合
y2_line=y(1)+((y(100)-y(1))*x-x(1))/(x(100)-x(1));
hold on
plot(x,y2_line,'Color','blue','LineWidth',1)
%线性度
Lm2=max(abs(y-y2_line));
yFS2=max(max(y))-min(min(y));
L2=Lm2/yFS2;
disp(L2);
%显示图片名称
legend({'原曲线','最小二乘法拟合直线','端基线法拟合直线'})