最小二乘法在二传感器信息融合中的应用

1、最小二乘法在传感器信息融合中的应用摘 要 本文用多维回归方程建立被测目标参量与传感器输出量之间的对应关系。并进行多维标定/校准试验，然后，按最小二乘法原理由试验标定/校准数据计算出均方误差最小条件的回归方程中的系数。用已知系数的多维回归方程计算出相应的输入被测目标参数。关键词 最小二乘法 信息融合 传感器1引言通常传感器都存在交叉灵敏度，表现在传感器的输出值不仅决定于一个参量，当其他参量变化时输出值也要发生变化。传感器信息融合技术就是通过对多个参数的监测并采用一定的信息处理方法达到提高每一个参量测量精度的目的。在只要求测量一个目标参量的场合，为达到提高被测目标参量的测量精度的目的，其他参量都是

2、干扰量，其影响应被消除，既然检测了多个参量，每一个参量测量精度都获得提高。压阻式压力传感器存在对静压、温度的交叉灵敏度，尤其是对温度的敏感成为它最大的缺点。人们为了消除温度对它的影响付出了长期的努力和高昂的代价。近来此类传感器采用了信息融合处理技术使得温度附加误差小于0.25%FS/55，测量精度达到（0.1-0.075）%FS。本文用多维回归方程建立被测目标参量与传感器输出量之间的对应关系。并进行多维标定/校准试验，然后，按最小二乘法原理由试验标定/校准数据计算出均方误差最小条件的回归方程中的系数。这样，测量时当测得了传感器的输出值，就可用已知系数的多维回归方程计算出相应的输入被测目标参数。

3、2二传感器信息融合智能传感器已知压力传感器输出电压U，且存在温度灵敏度。因此只对压力传感器进行一维标定实验，并由此获得输入（压力P）与输出（电压U）特性曲线来求取被测压力值会有较大误差。因为被测压力P不是输出值U的一元函数。现在由另一温度传感器输出电压Ut代表温度信息t，则压力参量P可以用U及Ut二元函数来表示才较完备，即 同理，可将压力传感器输出电压U描述为压力参量P和温度传感器输出Ut的二元函数，即 由二维坐标(Uk,Utk)决定的Pk在一平面上，可利用曲面拟合方程，即二维回归方程来描述，同样，由二维坐标（Pk，Utk）决定的Uk也在一个平面上，也可由二维回归方程来描述。如果回归方程中的各

4、个系数已知，那么用于检测压力P和输出U的二元输入-输出特性，即曲面拟合方程就确定了。为此，首先要进行二维标定实验，然后根据标定的输入、输出有最小二乘法原理确定常系数。2.1实验标定数据硅压阻式压力传感器采用恒流源供电，电路原理如图1所示。电桥B、D两端输出电压U为输入参量压力P的输出信号；A、C两端输出电压Ut为输入参量温度t的输出信号。实验标定数据如表1所示。压力标定点取6个，即n=6；温度标定点的数量也取6个，即m=6。 表1 压力传感器二维实验标定数据P/104Pa01.02.03.04.05.021.5U/mV-13.8410.6928.8847.0565.1983.36Ut/mV27

5、.6426.9526.4325.9225.4524.9428U/mV-13.499.3226.3443.1259.9976.82Ut/mV34.4133.9333.4732.9332.4731.9134U/mV-10.807.5424.8442.0559.2576.38Ut/mV37.7636.9236.4435.9735.3935.0944U/mV-9.726.5623.8741.2158.5875.87Ut/mV54.8853.9752.8752.4151.9351.5550U/mV-8.624.8621.8438.7056.3273.75Ut/mV65.7764.7963.8462.91

6、61.9961.0670U/mV-7.723.7221.2538.655.5673.28Ut/mV86.1284.9483.7882.6581.5580.45UtUAR4R3R2R1BDC图1 硅压阻式压力传感器的两个输出信号由表所列实验标定数据可见，传感器的压力输出信号随工作温度的升高而减小，在工作温度21.5至70范围内，在满量程压力值（5.0104Pa）输入时，输出值随温度的最大变化值为10.08mV，压力传感器的温度灵敏度系数，压力传感器存在的温度附加误差在满量程时的相对误差为，压力传感器的零位温度系数，温度传感器压力灵敏度系数。2.2数据融合处理为便于与未处理前各性能指标加以比较，本

7、文根据输入压力P和温度输出Ut来拟合压力传感器输出电压U。 根据温度输出Ut及传感器输出电压U拟合输入温度t。 为确定式所表征的二次曲面拟合方程式的常系数，依最小二乘法原理，求得系数值满足均方误差最小条件。由二次曲面拟合方程计算得到的值与标定值之间存在误差2k为： k=1，2，mn 总计有mn个标定点，其均方误差R1应最小 有式可见，其均方误差R1是常系数的函数，根据多元函数求极值条件，令下列各偏导数为零，即；则得下列方程：式中：L= mn。同理，有可得到如下方程：解方程组得常系数： 于是曲面拟合方程就确定为 2.3数据融合处理效果二传感器信息融合智能传感器系统框图如图2所示。系统输入信号压力

8、为P，PtUUtU或PUt或t二功能传感器数 据融 合图2二传感器信息融合智能传感器系统框图温度为t，它们的相应电压输出为U、Ut。数据融合后后的输出可以是电压P、温度t，或者是它们相应的电压U、Ut。选用何种输出形式依所选用的拟合方程式而定。采用上式的数据融合输出值为U和t，列于表2。表2 按式数据融合的输出值U与tP/104Pa01.02.03.04.05.021.5修正温度U/mV-11.537.6423.8943.0360.5878.06t/23.1923.7423.7623.8624.1224.3628修正温度U/mV-11.407.5123.8242.8960.5577.81t/2

9、7.5627.6227.7427.6127.9627.6734修正温度U/mV-11.367.3223.7642.5959.5877.18t/30.2730.5831.2130.5731.0331.0444修正温度U/mV-11.277.2523.3642.4458.9276.95t/41.8341.9142.0442.4941.5842.6150修正温度U/mV-11.206.9622.8742.2858.7376.23t/53.3152.7353.1353.8353.2352.8770修正温度U/mV-11.126.6422.8441.6857.7875.87t/69.2169.0168.3569.3869.1669.16由表2融合输出U可见，其零位温度漂移最大为0.14Mv，可得零位温度系数0=.10-4/，温度灵敏度系数s=5.7810-4/，与其融合前相比，性能提高一个数量级，传感器输出值受温度的影响大大减小。同时由表修正温度值可以看出，同样在70时，压力变化5104Pa，温度改变仅为0.05，故其压力灵敏度系数p=7.010-3/5.0104Pa，温度值受压力的影响也减小了一个数量级。3结论用最小二乘法对二功能传感器数