离散型函数拐点算法及应用.docx

1、离散型函数拐点算法及应用离散型函数拐点算法及应用 导读：就爱阅读网友为您分享以下“离散型函数拐点算法及应用”资讯，希望对您有所帮助，感谢您对的支持!离散型函数拐点算法及应用 李涛田晓君(黄石理工学院机电工程系, 湖北黄石,435003) 摘要在数学定义及物理意义的基础上，研究了利用计算机来计算数学拐点的一种快速有效的 方法，提出了对于离散数据以非拟合曲线的方式计算拐点的算法。该方法已成功地应用于漆包线生产过程质量监控。 关键词拐点算法离散数据非拟合曲线 文章编号：中图分类号：TP30文献标识码：A The Algorithm and application of discrete inflec

2、tion point Li Tao, Tian Xiaojun (Dept. of Mechanic and Electronic Engineering, HuangShi Institute of Technology, HuangShi, 435003) Abstract: Based on the mathematical definition and physical meaning, a fast effective method was studied to calculate mathematical inflection point using computer, and a

3、 algorithms was proposed tocalculate the mathematical inflection point of discrete data using non-fitting curve way. The methodhas been successfully applied in the quality control of enamel-insulated wire production. Key words: inflection point; algorithm; discrete data; non-fitting curve 1 引言 由于拐点所

4、反映的是相应曲线上升(或下降) 趋势发生变化的点，它在企业的生产管理等 方面比如销售分析、企业经营状况分析、产品质量的监控、生产成本的分析等各个方面都有 重要的作用。 在数学上，拐点的计算一般是通过计算函数二阶导数的零点及其不存在的点的方法来实 现的。这种方法在实际应用中往往遇到计算量非常大的问题。这是因为实际应用过程中的函 数及其二阶导数的形式往往非常复杂，这样就会造成在求解零点时计算量非常庞大，或者二 阶导数不存在的点的判别十分困难。此外，在实际应用中所得到数据的常常是以离散的点的 形式给出的，由此所构成的离散型函数在利用上述方法求解拐点时，还必须事先根据离散数 据点求解相应拟合曲线，这样

5、往往费时费力而且精度不高。因此，寻求一种快捷、精确的拐 点算法十分必要。 2 数学计算模型 2.1拐点判别式根据拐点特性，在拐点前后曲线的上升(或下降) 趋势不一样，如图 1所示。若P=(p0, f(p0)为该曲线的拐点，A 、B 点为P 前任意两点，C 、D 为P 后任意两点，则有判别式(kAB-kBP)(kPC-kCD)0 (1)其中kAB 、kBP 、kPC 、kCD 分别为直线AB 、直线BP 、直线PC 、直线CD 的斜率。 D C P A B x0=ax1 x2 xk xn-1 xn=b图1 曲线拐点图2 拐点计算 2.2拐点判别算法利用判别式(1)即可获得曲线拐点判别算法。对于区

6、间a, b上的函数 f(x)，为求其拐点，先将区间a, b分割为小区间，小区间长度最大者趋于0。设分点分别为x0, x1, xk,xn-1，对应纵座标为f(x0),f(x1),f(xk), f(xn-1)，如图2所示。如果拐点处于点xk-1与xk(k=3,4,5n-3)之间，根据判别式(1)有 f (x ).f (x ) f (x ).f (x ) f(x ) .f (x ) f (x ) .f (x ) k.1 k.2 k.2 k.3 k+2 k+1 k+1 k F(k) =( .)( .) 0 (2) x .xx .xx .xx .x k.1 k.2 k.2 k.3 k+2 k+1 k+1

7、 k 为简化算法，不失一般性，可将区间a,bn等分，则式(2)可变换为如下判别式F(k)=(f(xk-1)+ f(xk-3)-2 f(xk-2)(f(xk)+ f(xk+2)-2 f(xk+1)0 (3) 其中k=3,4,5,n-3。若判别式(2)或(3)在k=m处成立，则说明位于xm 与xm+1之间。显然，当0(或n ) 时，则xm p0，即m 充分大时，(xm,ym)P=(p0,f(p0)。 2.3算法终止判别算法开始时分别将各等分点代入(3)式，计算出F(k)，当F(k)小于0时算法即可终止。在实际应用过程中，n 取无穷大不现实也无必要。实际应用中 n 的取值由数据特性、区间a,b的宽度

8、，以及精度要求来决定。区间a,b的宽度取5cm ，计算机屏幕取1024768分辨率时，n 的取值对函数图形的影响如图3所示。 n=50 n=200 n=500 n=2000 图3n 的取值对函数图形的影响 由图3可知，在实际应用中，对精度的要求不是很高的情况下，一般 n 取3003000就能取得较好的效果。 3 实际应用计算模型 由上述理论分析知，对于连续函数而言，上述方法实际上就是通过离散取点的方式将原 曲线离散化后计算出拐点。对于由离散的点构成的函数，计算判别式(2)与(3)一样适用。在实际工程应用中，用于待求拐点的信号往往是一个或多个模拟信号，在送入计算机处 理之前先要进行一次A/D转换

9、。比较常用的一种做法是使用数据采集卡来实现被测信号的数字化。数据采集卡在工作的过程中，以一定的采样频率对被测信号取样转换。得到的数据 就是一组以时间xk 为横轴X ，以被测量大小f(xk)为纵轴Y 的离散点(xk,f(xk)。因而在实际应用中，是使用判别式(2)或(3)来求取拐点的。如果进行数据随机采样，则利用判别式(2)进行计算，如果进行数据等时间间隔周期采样，则利用判别式(3)进行计算。 4 数据处理 计算机端数据处理所要进行的工作主要包含数据预处理、绘制曲线、拐点计算等 3个方面。各相应步骤以算法在漆包线的绝缘状况的监控中的应用实例来说明。 4.1对采集到的数据进行预处理数据预处理主要包

10、含以下2个方面。(1)去掉重复点。由于温度场(恒温箱) 的温度变化比较缓慢，相应参量(介质损耗角正切)的变化也就比较缓慢，而数据采集卡的采样频率较高，这样就会采集到大量的相同数据，这 些相同数据在代入判别式(2)或(3)计算F(k)时会出现除零错，即使对判别式(2)或(3)进行改进，使它不出现除零错的情况，这些相同数据也会大大增加计算量，所以在将采集到的数据带入 公式前必须要将重复值去掉。 (2)去掉错误值。在数据采集过程中由于各种因素的影响(如车间环境恶劣，电气干扰较为严重) ，往往会造成数据采集误差，而出现错误点。这些错误点如果偏离真实值较远会严 重影响计算结果，因此必须去掉。实验结果表明

11、，通过去极值取平均的方式就可以较好对采集到的数据进行上述两个方面 的处理。 4.2 绘制曲线根据采集到的数据在计算机屏幕上绘制参量(介质损耗角正切)-温度曲线，以坐标横轴为 温度，坐标纵轴为参量(介质损耗角正切) 。由于恒温箱中温度是以与时间成线性的方式增加，这样在实际编程过程中可以用时间通过计算来取得温度，时间可以通过一个软件计数器取 得，得到时间后只需要进行一次简单的线性转换就可以得到温度了。根据采集到的数据点以 及计算出来的对应温度即可在坐标轴上描出相应的点。相邻点用直线连接起来即可得到近似 的介质损耗角正切-温度曲线。 4.3拐点计算将得到的各坐标点依次代入判别式(2)或(3)计算F(

12、k)，判断代入公式的点是否为拐点。拐点有可能存在多个，所以在计算拐点时一定要将所有的点均带入计算一遍。如果采集的数据不需要存储，绘制曲线和计算拐点可以同时进行，程序框图如图 4所示。 图4 程序框图 5 应用实例 在漆包线生产企业中，漆包线的质量监控非常重要，而漆包线的绝缘状况的监控又显得 尤为重要。测量介质损耗对判断电气设备的绝缘状况是一种传统的、十分有效的方法，绝缘 能力的下降直接反映为介质损耗增大。介质损耗即绝缘材料在电场作用下，由于介质电导和 介质极化的滞后效应，在其内部引起的能量损耗。通常在实际应用中采用测量介质损耗角正 切的方法来反映介质损耗。通过检测不同温度下的介质损耗量还可以反

13、应出漆包线在不同温度下的绝缘性能。同时，随着温度的升高，漆包线的绝缘漆还会出现由玻璃态到高弹态的转 变。转变点反应出了漆包线的耐高温特性，它在介质损耗角正切-温度曲线上即表现为一个 拐点。所以，漆包线的质量监控通常是通过绘制介质损耗角正切-温度曲线来描述其绝缘性 能，通过介质损耗角正切-温度曲线的拐点来描述其耐高温特性。检测的过程如图 5所示。 打印显示介质测量仪 数据采集卡 计算机 图5 检测过程 其中，试样被放置在一个恒温箱里，恒温箱的温度从某一起点温度开始，以与时间成线 性的方式增加，这样介质损耗测量仪测得的数据就是一组随时间变化而变化的介质损耗角正 切值。由于介损仪测得的数据是一组模拟

14、电压量，数据采集卡按照一定的采样频率对这组模 拟信号进行采样，最后送入计算机处理。 6 结论 根据拐点的特性，提出了一种新的简单快捷的采用非拟合曲线的方式计算离散型函数的近似拐点的方法。理论分析表明，通过非拟合曲线近似求拐点的方法可以大大的减轻计算 拐点的计算量，特别适合于工程实际中的应用。利用该方法同时还可以绘制近似曲线，这样 对于工程实际中以离散点形式构成的函数，比如销售分析、企业经营状况分析、产品质量的 监控、生产成本的分析等各个方面都有重要的作用。 本文作者创新点在于，采用非拟合曲线的方式计算离散型函数的近似拐点。该方法不 但可以极大地减少运算量，同时还可以绘制近似曲线。 参考文献 1

15、汤锋, 王章野, 刘钢等. 基于矢量化方法的套色图光滑变换算法J.中国图像图形学报. 2003, 8(z1):166-169.2 Shu Chien Huang, Yung Nien SunPolygonal approximation using genetic algorithmsC Proceedings ofCAIP99, Ljubljana, Slovenia, 1-3 September, 1999: 175-182.3 Vincent Beau, Mark Singer, Reduced resolution and scale space for dominant featur

16、e detectionin contours ,Pattern RecognitionJ, 2001, 34(2): 287-297.4朱亚云, 方承远. Windows环境下实现实时数据采集的可行性J.上海理工大学学报 .1999, 21(2):153-1585严蔚敏, 吴伟民. 数据结构M.清华大学出版社.1992.6董海卫, 江早, 王永军. 基于矢量化二值工程图得符号提取算法J.计算机辅助设计与图形学学报.2000,12(4):294-298.7张永强, 李学东, 刘富军. 应变实验数据采集与处理系统的硬件设计与实现J.微计算机信息.2005,1:65-668杨琨超, 毕荣. 电力电缆(电线) 及附件生产加工工艺与质量检测标准实用手册M.安徽文化音像出版社.2003.9电缆电线生产使用质量检测与技术标准应用手册M.北京科海电子出版.2003.作者简介：李涛(1980) ，男，湖北荆门人，湖北黄石理工学院机电工程学院教师，主要从事智能测控方面的研究工作，硕士。田晓君(1969) ，男，湖北浠水人，湖北黄石理工学院机电工程学院教师，主要从事智能测控方面的研究工作，硕士。XX搜索“就爱阅读”,专业资料,生活学习,尽在就爱阅读网,您的在线图书馆