C++作业设计.docVIP

基于单片机的传感器曲线拟合 **学院 学号：1010***** 陈** 摘要：传感器使用前要进行标定，标定时要进行曲线拟合，若用微型计算机处理很简单，但实际中用单片机标定时，因它无矩阵运算，只能进行一般的代数运算，处理就显得非常不方便。推导了二次多项式曲线拟合待定系数的代数计算公式，应用这些公式来处理数据非常方便。通过实例说明了其应用，实践证明，效果良好，可广泛用于基于单片机标定中的二次项式曲线拟合。 关键词：曲线拟合；最小二乘法；传感器；标定；单片机 引 言 在工程实践和科学实验中，经常利用曲线拟合来寻求实验数据与控制该实验物理规律的数学模型的相互关系。常见的曲线拟合有最小二乘法曲线拟合、指数函数曲线拟合、傅里叶级数曲线拟合和多项式曲线拟合等。多项式曲线拟合在实际中应用广泛，若用最小二乘法拟合，通常利用微机进行矩阵运算能快捷方便地予以处理，但在实践中，有时因有其它要求没有使用微机而是使用单片机，例如：在电容式微位移传感器中就是利用单片机进行测量标定，这时会因内存不够或无数组等原因而不能进行矩阵运算，而只能通过一般的代数运算予以处理，并且这时拟合的次数也不高。所以，就必须先求出拟合函数中的待定系数a0、a1、a2等的代数计算公式，然后，再进行数据处理。因此，本文推导出了用最小二乘法拟合应用较广的二次多项式曲线拟合的待定系数a0、a1、a2的计算公式，并举出实例予以说明其应用[1]。 1、公式推导 设有一组实测数据 i=1,2,3…..n,其拟合函数为 ， 实测值为y,则任意y 与拟合值的偏差为，则偏差的均方值为 为使偏差均方值为最小，则对上式分别对a0、a1、a2，求偏导并令为零，整理后得 此时可以转换成矩阵的形式 上式等价于 ， 令则 再令;即BA=C,则 （2） 矩阵B的伴随矩阵为 矩阵B 的逆矩阵为,所以 最后得到待定系数为 （3） 这样，求a0、a1、a2和的运算式均为代数计算，即可用单片机对实验数据进行拟合。 2 应用实例 一电容式微位移传感器测量系统是用单片机来处理数据。单片机的输入量x是传感器的电压信号；输出量y是位移量，送入液晶显示器显示。为建立x和y之间的数学关系，就必须对系统进行标定[2]。通过理论推导，可用该方法进行曲线拟合标定。 这里记录了部分数据，如表1所示。 表1 实验标定数据 利用（1）、（2）式得到处理数据如表2 表2 处理数据 再由式（3）即得到待定系数a0=2.988；a1=2.060；a2=0.972 其相应的拟合方程为 =2.988+2.060x+0.972 由此可得到的处理数据如表3 表3 处理数据 由表3可知，拟合最大相差误差为0.99%，拟合方程的神谕标准误差为 其中n为测量次数；q为待定系数的个数。 3 结论 利用本文推导的二次项式曲线拟合待定系数的代数计算公式求待定系数比较方便，然后，再进行曲线拟合，其拟合准确度也很高，尤其用单片机进行标定时优点更为显著。因此，它可以广泛用于单片机标定中的二次项式曲线拟合。 参考文献： ［1］王福明，贺正辉，索! 瑾 应用数值计算方法［4］ 北京：科学出版社1992.190—298 ［2］黄长征，李圣怡 超精密车床主轴回转误差运动的动态测试 航空精密制造技术 2002,38（4）：1—3