专利名称：一种消除铸铁热分析曲线干扰的滤波方法
技术领域：
本发明属于钢铁铸造、信号处理技术领域，主要涉及在铸铁过程中，为了实现液铁质量在线监测，提高铸铁中碳、硅含量的分析精度，用于铁液温度曲线分析的平滑滤波，为一种消除铸铁热分析曲线干扰的滤波方法。
背景技术：
铁液成分的在线检测，要求迅速、准确提供铁液成分检测结果，以便在线调整铁液的化学成分，一般化学分析法难以适应在线检测的速度要求，热分析法是通过测定伴随物质的温度变化和化学反应的进行而发生的热力学性质、物理性质的变化，进而确定物质结构的热变化及化学反应常数的一种方法。其基本原理是在冷却曲线上的任何变化直接与凝固体的热量变化有关，即与相变有关。因此，冷却曲线实际上是特定合金凝固过程的一次记录，凝固过程的所有信息都包含在冷却曲线上，这样，只要求得合金在加热或冷却过程中的温度变化曲线，就可间接求出该合金相应的特性参数，如化学成分、组织及性能等。早期的热分析法是根据合金的冷却曲线定量的获得合金的固相、液相、转变相信息的一种方法，一般称为直接热分析法，但是冷却曲线常常不像人们所期望的那样清楚地显示出它的特征值，因而如何正确的分析热分析曲线成了直接热分析法测试中的关键。碳和硅两种元素的含量是铸铁过程中非常重要的两个元素，通过热分析的方法测量铸铁中的碳、硅含量通常的方法就是根据铁水冷却曲线，获得液相线温度和固相线温度两个温度点，经过经验公式计算获得。参考文献王利华，石德全，李大勇.热分析技术在铸造生成质量检测上的应用(一).金属铸锻焊技术，2009，Vol38(21) :72-76中描述了热分析法测量碳、硅含量的方法，但是对于如何高精度的提取液相线温度和固相线温度并未研究。为了使碳、硅分析的精度达到0. 1%以上，必须对冷却过程中的温度测量精度达到0. 1%以上，经过热电偶温度传感器将温度转换成电压信号后，在满量程50mv的条件下，1° C对应0.04mV左右，精度必须达到0. 1%以内。在如此高的精度要求下，获得的热分析曲线由于环境、电磁干扰特别是大用电量的设备功率变化等影响会导致产生纹波和毛刺，这些纹波和毛刺的存在会对液相线温度和固相线温度的判断带来影响。因此需要对热分析曲线进行一定的滤波处理。同时这种处理还要满足在线快速分析的要求。现场实验表明，如果采用传统的滤波方法在对曲线去除毛刺的同时，曲线在毛刺附近会形成一定起伏，造成局部台阶，并且会使得温度曲线整体变化趋于平滑，使得原有温度曲线的一阶导数和二阶导数都趋于缓和，从而掩盖一些热分析曲线原有的温度突变信息和曲线变化形状。而热分析曲线分析中最关键技术就是通过温度曲线变化捕捉液相线温度和固相线温度。因此，传统的滤波方法很难满足铸铁热分析曲线的滤波处理要求，本发明采用的方法可以在去除曲线纹波和毛刺的同时，保留曲线剧烈变化的温度点，既提高了温度变化曲线的平滑程度，也为液相点温度和固相线温度的高精度提取提供了条件。 发明内容
本发明要解决的问题是为了实现高精度的铁水成分在线分析，去除铸铁热分析曲线中存在的纹波和毛刺等干扰，需要对热分析处理进行处理，既要准确保留曲线原来的变化规律，又要去除干扰的影响，由于传统的滤波方法在去除毛刺的同时，对曲线整体进行了平滑，处理后的曲线会使得温度变化变得缓和，不利于对于液相温度和固相温度这两个关键温度点的提取。本发明的技术方案为一种消除铸铁热分析曲线干扰的滤波方法，包括以下步骤I)以小于0. 25秒的采样周期，对铸铁温度进行采集；2)取一段连续的温度采集点，个数为N，N取值为5 10的整数，将此N个温度采集点所在区域称之为滤波窗口；
3)根据滤波窗口参数N和已采集的温度数据，构成原始滤波序列IT1J2,
TnI ;4)在这个序列的前后各取M点，M < N/2，共2M个点构成子序列，gp IT1，…，T , Tn_m+1，…，Tn}，用最小二乘法对这个子序列IT1,…，Tm，Tn_m+1，…，Tn}进行曲线拟合，获得二次拟合曲线T=AkWa1I^atl，其中％、&1、&2为拟合曲线得到的曲线参数，根据拟合曲线计算获得N个采样点在拟合曲线上的的拟合数据点，构成拟合数据序列{T i，T' 2，…，T',，…，T' J5)根据原始滤波序列IT1, T2，…，Tk，一,T1J和拟合数据序列{T'2，…,T' k, * *, T' N}获得误差序列(A1, A2，…，Ak,…，AJ :
权利要求
1.一种消除铸铁热分析曲线干扰的滤波方法，其特征是包括以下步骤 1)以小于O.25秒的采样周期，对铸铁温度进行采集； 2)取一段连续的温度采集点，个数为N，N取值为5 10的整数，将此N个温度采集点所在区域称之为滤波窗ロ； 3)根据滤波窗ロ參数N和已采集的温度数据，构成原始滤波序列IT1,T2,…Tim, TJ ； 4)在这个序列的前后各取M点，M< N/2，共2M个点构成子序列，即(T1, - ,TM, T腿，…，TJ，用最小二乗法对这个子序列IT1, - ,TM, Tn_m+1，- ,TJ进行曲线拟合，获得二次拟合曲线Τ=#2+#+^其中为拟合曲线得到的曲线參数，根据拟合曲线计算获得N个采样点在拟合曲线上的的拟合数据点，构成拟合数据序列IT' dT' 2，…，T' k，···，!1' J 5)根据原始滤波序列IT1,T2, -,Tk, -,TJ和拟合数据序列IT' !,Ti 2，…，T' k,…，T' J 获得误差序列(A1, Δ 2) ···，Δ k, ···，Δν} δ: = |y - T2 ^ =|W| ’户 a- = Tn-Tn 6)根据误差序列W1,A,,···, Ak,..., Λ N}，计算得到该序列的平均误差 _I入 Δニ丄アΔ,； Ntt 7)设定单点误差阈值Λτ，两个平均误差阈值平均误差上限^和平均误差下限 其中毛く τΗ ,根据误差序列及平均误差与这三个阈值之间的关系，将曲线分成三种情况 (a)满足I< I并且Vi, I彡i彡N，Ai < Λ τ，则认为曲线平滑或存在小毛刺，不存在温度突变； (b)满足1<又并且ョもI彡i彡N，Ai> Λτ，则认为曲线中存在尖峰干扰，不存在温度突变； (c)除(a)(b)外的其他情况，认为曲线中有可能存在温度突变； 8)在步骤7)的(a)(b)情况下， 81)移动滤波窗ロ，将误差序列中的误差最大值对应的数据点移动至滤波窗ロ的中间，滤波窗ロ中仍为N个温度采集点，得到新的滤波序列ITR1, TR2,…TIV1, TRJ，对新的滤波序列进行步骤4)的曲线拟合，用拟合后得到的N个曲线拟合数据作为温度估计数据:f.; = TRi,i=l, “·Ν ； 82)如果误差序列中的误差最大值对应的数据点在滤波窗口中间，则用拟合数据序列IT'2，…，T' ,，···，Τ' Ν}作为温度估计数据，i=(i=l，…N; 9)在步骤7)的(c)情况下，保留原来N个温度測量点原始測量值，不予处理，以保留原有的曲线形状和可能的温度突变信息，即=Ι；“=1，…N;步骤8)及步骤9)得到的温度序列4 i = 1，···Ν即为温度序列1\パ=1，···Ν滤波后的结果； 对步骤I)采集的所有温度数据依次进行步骤2)-9)的滤波处理，得到完整的铸铁热分析曲线的滤波曲线。
2.根据权利要求I所述的ー种消除铸铁热分析曲线干扰的滤波方法，其特征是步骤7)中设定单点误差阈值Λτ=10，平均误差上限^ =4，平均误差下限毛=2。
3.根据权利要求I所述的ー种消除铸铁热分析曲线干扰的滤波方法，其特征是步骤7)中设定单点误差阈值Λτ=3，平均误差上限4 =3,平均误差下限毛=2。
全文摘要
一种消除铸铁热分析曲线干扰的滤波方法，设定滤波窗口，对窗口内的曲线分情况进行滤波，最终完成对热分析曲线的滤波。本发明滤波方法热分析曲线上的每个温度点首先判断属于毛刺、尖峰还是本身的温度突变。对于毛刺和尖峰干扰进行滤除，而本身的温度突变予以保留。因此，经过该滤波处理后，在去除毛刺的同时，保留了温度的突变信息，平滑了热分析曲线，提高了碳、硅含量分析精度。
文档编号G01N25/20GK102706921SQ20121019447
公开日2012年10月3日 申请日期2012年6月13日 优先权日2012年6月13日
发明者姚敏, 杨绍雨, 赵敏 申请人:南京第四分析仪器有限公司