利用最小二乘法减少铂热电阻测温非线性误差

利用最小二乘法减少铂热电阻测温的非线性误差1引言各种智能仪表中，广泛使用钳热电阻作为温度传感器。典型的用法是前端采用不平衡电桥测量钳热电阻随温度变化的毫伏信号输出,再经过放大和A/D转换，送到单片机中进行运算。这种用法中，钳热电阻的非线性和不平衡电桥的非线性，给最后的温度测量带来一定的误差。解决钳热电阻不平衡电桥测温法中非线性误差的方案有许多种，我们采用的是最小二乘法，只需要一个简单的多项式即可以获得满意的测量精度。2钳热电阻的非线性R=R0(1+at+bt2) (1)R=R0(1+at+bt2) (1)t0、显然，随着温度的升高，钳热电阻的非线性越来越严重。以钳电阻Pt100为例，式中的常数为:RO=100Qa=3.90802X10-3°C-1b=-5.802X10-7C-2如果采用折线法拟合，则在100〜200C区间段做简单的直线拟合，即设:R1=A（t-100）Ro+R100-^00 100= 100xR 100x100A=R200一R100=75.84一38.5=3.734x10--^00 100= 100xR 100x100则可以得到在此区间的误差，见表1。表1Pt100折线法测温100〜200C区间误差温度T/温度T/C电阻R/Q100138.50120146.06140153.58160161.04180168.46200175.84电阻R/Q误差e/%138.500.00145.960.07153.440.08160.900.09168.370.05175.840.00采用折线法固然可以减少误差，但是如果对于0〜800C测温方案，则需要8个折线公式，而且需要准确地选取折线点。3电桥的非线性当钳热电阻阻值变化幅度较大时，不平衡电桥即存在严重的非线性。由图1可知Uo和鸟的关系是：U=U(—^t—— —)0R+R]R3+R2若令Ri=R2,R3=R0,设Rt=R0+R&，则有：

URUo_-(R+R)(R+R+R)R

0 10 1 NtUR式中R=R0(at+bt2)UI (R0+R「2RR1N其误差为:UI (R0+R「2RR1N其误差为:U=U"U0= RNt0 U0 R0+R]图1不平衡电桥线性处理的误差见表2。表2不平衡电桥线性化处理误差温度T/°C电阻RAt/Q误差e/%5019.40.3710038.50.7415057.311.1020075.841.45300112.022.15400147.042.83500180.903.48

600213.594.15800275.515.30由此可见，在实际应用中，系统的误差主要是不平衡电桥的非线性，同时，当测温范围＞100°C时，这种线性化处理的误差将随着量程的增加而增加，最终达到不能允许的程度。4解决方案铂热电阻不平衡电桥法测温方案中非线性误差的常用解决方案有插值法、折线法和迭代法。但是插值法和折线法对内存空间有一定要求，同时要求选取合适的插值点，而迭代法则存在运算复杂,运算周期长的缺点。在设计智能热量计的过程中，我们采用最小二乘法，很好地解决了铂热电阻不平衡电桥法测温方案中的非线性误差。经过推导，可知送到A/D转换电路的电压Ut为：KUR(RKUR(R0+R)(R0+q+R&)RAt(5)U= 气"*1 R(at+bt2)t (R0+R「(R0+R]+R0(at+bt2))0而实际应用中是已知电压值U求温度值t。只要有足够的数据点，就可以用最小二乘法进行曲线拟合。式(5)可以简化为：(6)U= K R(6)t5.2x(5.2+R『At式中RyR0(at+bt2)KrK0UR!对于使用Pt100，量程为0〜800C，对应输入电压为0〜5V的测温方案，R=0.27551kQ，AtU=5.000V，则可求得Kt=516.726289。再将各温度对应下的R^值代入公式(6)，即可得到一组数据点，见表3。表3温度、电压对应表电压U/V温度T/C电压U/V温度T/C0.00003.0404500.369503.3415000.7301003.6345501.0831503.9216001.4282004.200650

1.7662504.4737002.0963004.7407502.4183505.0008002.733400在实际计算过程中，为防止运算溢出，拟合多项式的形式为：t=a0+a1(U-U)+a2(U-U)2+a3(U-U)3式中— 147 一一U=tU=2.64544i=1用最小二乘法进行拟合得：a0=385.949a/=159.635a2=6.06577a=0.336545即对于使用Pt100，量程为0〜800°C，对应输入电压为0〜5V的测温方案，其多项式拟合公式为:t=385.949+159.635(U-2.64544)+6.0677(U-2.64544)2+0.336545(U-2.64544)3其误差见表4。表4Pt100测温拟合公式0〜800C误差实际温度/C对应电压/V计算温度/C误差绝对值/C00.000-0.140.1433.740.25033.730.0167.970.50068.060.09102.760.750102.870.11138.081.000138.200.12174.031.250174.090.06210.511.500210.550.04247.641.750247.630.01285.402.000285.350.05323.642.250323.750.11362.942.500362.860.08402.792.750402.710.08443.413.0000443.330.08484.793.250484.750.04527.003.500527.010.01570.303.750570.130.17614.064.000614.150.09659.004.250659.100.10

704.904.500705.010.11751.904.750751.910.01800.005.000799.840.16此公式与数据点误差绝对值的最大值为0.17°C。同理，对于使用Pt100,量程为0〜200°C,对应输入电压为0〜5V的测温方案，其多项式拟合公式为：z=99.097+39.9947