三参数、四参数曲线拟合..培训讲学.ppt

正弦正弦(zhngxin)曲线拟合的三曲线拟合的三参数法与四参数法参数法与四参数法第一页，共32页。正弦正弦(zhngxin)曲线拟合的意义曲线拟合的意义 由正弦波形的采样序列获得(hud)其拟合正弦曲线函数，是一种基本信号处理方法，在许多场合下获得(hud)了应用，如评价数据采集系统的有效位数、采集速率、交流增益、通道间延迟、触发特性等，在调制信号的数字化解调和失真度测量中，也有应用。第二页，共32页。曲线拟合的一般曲线拟合的一般(ybn)过程过程正弦信号(xnho)采样A/D变换信号(xnho)处理拟合正弦曲线第三页，共32页。数学上，幅度、频率、相位和直流偏移4个参数可以唯一确定一条正弦曲线。曲线拟合的目的就是通过分析输入(shr)的正弦信号，得到正弦波形的四个参数值，从而得到拟合曲线。第四页，共32页。在已知输入正弦波形的前提(qint)下，怎样确定它的4个参数呢？第五页，共32页。正弦正弦(zhngxin)曲线拟合的总体思路曲线拟合的总体思路 主要是通过改变(gibin)拟合正弦函数的幅度、频率、相位和直流偏移，使拟合函数和采样序列各点的残差平方和最小，从而获得正弦波形序列最小二乘拟合结果。第六页，共32页。正弦曲线拟合的总体正弦曲线拟合的总体(zngt)思路思路假设采样点数是L，采样数据是D(I)，I:0,1,L-1拟合函数(hnsh)是S(t)=Asin(2ft+p)+C则残差的平方和为 为采样时间间隔拟合的目的就是找到让E最小的四个参数A、f、p、C第七页，共32页。三参数三参数(cnsh)法简介法简介 三参数正弦曲线拟合，特指信号频率已知时获取幅度、相位和直流偏移的波形拟合方法，它是一种闭合算法，无须迭代(di di)即能获得结果，没有收敛问题，具有良好的实用性。第八页，共32页。三参数三参数(cnsh)法的算法法的算法 在标准IEEE std1057-2007 IEEE Standard for Digitizing Waveform Recorders 的 Annex A 中给出了一种三参数正弦(zhngxin)拟合的算法。第九页，共32页。三参数拟合三参数拟合(n h)算法示例算法示例设理想正弦信号为三参数(cnsh)正弦波曲线拟合过程，即为输入信号的数字角频率已知，选取或寻找A，B，D，使下式所述残差平方和最小：则，参数(cnsh)A，B，D即为A0，B0，D0的最小二乘拟合值。为寻找出A，B，D，构造矩阵第十页，共32页。三参数三参数(cnsh)拟合算法示例拟合算法示例 残差平方和用矩阵表示为：当式E最小时可得x0 的最小二乘解为：拟合函数(hnsh)的幅度和相位表达形式为：其中：第十一页，共32页。三参数拟合三参数拟合(n h)算法示例算法示例 拟合残差为：拟合残差有效值为：其中：由于这是一种闭合算法,因而(yn r)收敛是肯定的。第十二页，共32页。四参数四参数(cnsh)法法 当正弦信号的四个参数都不知道时，一般采用四参数法进行拟合。四参数法也是最常用的一种正弦波拟合方法。与三参数正弦曲线拟合不同，四参数正弦曲线拟合是一个非线性迭代过程，没有解析公式可以直接应用(yngyng)获得结果，需要计算初始值进行迭代。第十三页，共32页。初始值的重要性初始值的重要性 初始值的精确度对于(duy)迭代结果有着很重要的影响。较大的初始误差将导致迭代发散，或收敛到局部最优值而非总体最优值上。第十四页，共32页。获取初始值的基本获取初始值的基本(jbn)方法方法频率f：(1)fft/dft (2)通过分析信号(xnho)过零点的时间间隔估计频率幅值A：峰峰值除以2直流偏移C：(1)计算信号(xnho)一个周期的平均值 (2)信号(xnho)最大值与最小值之和除以2相位p：第十五页，共32页。四参数四参数(cnsh)拟合的算法拟合的算法 四参数拟合有很多种算法。IEEE学会在标准IEEE std1057-2007 IEEE Standard for Digitizing Waveform Recorders 的 Annex A 中给出了一种方法，包括(boku)两种基本算法：一种通过矩阵运算，另一种通过迭代过程，二者均需要良好的初始条件估计。第十六页，共32页。四参数拟合的经典算法四参数拟合的经典算法(sun f)简介简介牛顿法：该方法是基于一阶泰勒展开与误差修正技术相结合的产物，搜索终止(zhngzh)的判据可以是参数增量，或残差平方和。顺序搜索法：顺序对每一个参数在初始值上使用增量搜索法寻找其最优点。第十七页，共32页。牛顿牛顿(ni dn)法简介法简介 牛顿法是对方程四个参数求偏微分，得到E对给定系数的增量的泰勒级数展开式。用增量对初始值进行校正，以此方法进行多次迭代，直到相关系数不再增大，或者设定一个迭代的次数，就可以(ky)得出四个值的最终结果。第十八页，共32页。四参数拟合的算法四参数拟合的算法(sun f)简介简介 顺序搜索法有一种算法是将四参数拟合过程拆分成两步走，可以避免四参数非线性迭代带来的收敛问题。该算法使用一种非线性迭代方法获得信号(xnho)频率估计值，然后在已知频率情况下，使用三参数最小二乘拟合算法获得最终结果。本质上是一种三参数方法。第十九页，共32页。四参数四参数(cnsh)顺序搜索算法示例顺序搜索算法示例()令 i=1,确定估计信号频率的大致区间.对于常见的等间隔采样,转步骤();对于非等间隔采样,直接转步骤().()利用 D F T 或 F F T 计 算信 号频率,设 为d,令迭代区间频率下限 ，迭代区间频率上限 (其中,c 为时钟频率,N为 D F T 或 F F T 的长度),转步骤().()观察采样序列过零点时刻,设第 m 个过“零点”(零点指采样序列的均值位置)时刻在区间tkm,tkm+1中,而第L(LM)个过“零点”时刻在区间 tkl,tkl+1 中,令 ，,其中m,l 为整数,转步骤().()令 ,从区间0l,0h中等(zhngdng)间距的取 2 M+1 个点(比如 M=5),利用三参数法分别计算出这些点对应的 A1j,B1j,C1j 和残差平方和 E1j(j=1,2,3,2 M+1).第二十页，共32页。()比较()中 2 M+1 个残差平方和,并找出最小残差平方和对应频率(记为 1)、正弦幅度(记为 A 1)、余弦幅度(记为 B1)以及直流偏移(记为C1).这就是正弦信号四参数的第 1 次估计值,其中,频率估计的最大偏差=0/M.()令 i=i+1,()从区间 il,ih 中等间距地取 2 M+1 个点,分别计算出这些(zhxi)点对应的 A ij,B ij,Cij 和误差平方和 E ij.第二十一页，共32页。()比较(bjio)()中 2 M+1 个误差平方和,并找出最小误差平方和对应的四个参数(分别记为 i,Ai,Bi 和 Ci),这就是正弦信号四参数第 i 次的估计值,其中,频率估计的最大偏差=0/Mi.()重复()(),直到找到满足精度要求的信号频率,将其记为,同时将与它对应的其他三个参数记为 A,B,C,那么,A,B,C 这四个参数就是正弦信号四参数的估计值.其中,频率估计的最大误差为max=0/Mi步骤()()给出了四参数估计法的一般步骤和频率估计的最大误差 max,可知,可以通过增大估计次数 i 来提高估计精度第二十二页，共32页。四参数四参数(cnsh)拟合的算法简介拟合的算法简介 还有学者使用遗传算法实现总体最优估计，以此实现四参数正弦参数的最小二乘估计，由于遗传算法原理本身可保证实现全局(qunj)最优逼近，可避免收敛到局部最优点上，从而具有良好的收敛性。第二十三页，共32页。拟合拟合(n h)误差误差 有很多因素会影响到拟合参数的精确度。序列长度、采样序列中含有的波形周期个数、采样量化误差、非线性误差等条件，都限制和影响了正弦(zhngxin)参数的估计。第二十四页，共32页。影响不确定影响不确定(qudng)度的因素度的因素(1)波形(b xn)采集速率。(2)波形(b xn)测量通道间延迟时间差。(3)采样序列的噪声及非谐波失真。(4)采样序列的抖动。(5)4参数正弦波拟合软件造成的测量不确定度，主要由于软件收敛判据、舍入误差、累积误差等造成；(6)另外，采集序列长度的变化、采集序列中所含信号的周期个数的变化，也将给测量带来影响，它们将体现在上述各项不确定度的分量中，不单独列出。第二十五页，共32页。拟合拟合(n h)参数的最小误差界参数的最小误差界 John P.Deyst1995年给出了正弦波四参数最小二乘拟合算法获得参数的误差界，使用蒙特卡罗搜索仿真法等对于各种可以想象的条件变化进行了细致研究，并分别以经验公式、误差界曲线等形式，给出了4个拟合参数随谐波次数和幅度、噪声、抖动、序列长度、序列所含信号周期个数等条件参量变化的规律。基本结论是：拟合获得的4个参数的误差界随着谐波阶次、序列长度、序列所含信号周期个数增大而变窄，随着谐波幅度、噪声、波形抖动的降低(jingd)而变窄。每个参数的误差界应该在一个确定区间内变化，最小误差界即是其Cramer-Rao界。第二十六页，共32页。拟合参数拟合参数(cnsh)的最小误差界的最小误差界Cramer-Rao界：在四参数(cnsh)正弦波拟合中，误差界的指数表达式给出了拟合参数(cnsh)误差随着信号周期个数和谐波阶次变化而变化的一个公式。其中：N为采样的周期个数，n为采样点数，h为谐波(xi b)次数（整数），A为幅值，p为相位，C为直流量。第二十七页，共32页。拟合参数拟合参数(cnsh)的不确定度的不确定度 利用Cramer-Rao界，可以估计出谐波失真造成的4个测量参数的不确定度。以相位为例，假设 h次谐波的幅值 造成的参数误差分别(fnbi)在各自的误差界内均匀分布，则 给 带来的测量不确定度为：由于三角函数基之间的正交性，不同谐波互不相关(xinggun)，所有谐波带来的总的不确定度为：第二十八页，共32页。拟合拟合(n h)参数的不确定度参数的不确定度同理也可以得到幅值和直流分量(fn ling)的不确定度。幅值：直流量(liling)：第二十九页，共32页。拟合拟合(n h)参数的不确定度参数的不确定度对于频率(pnl)不确定度的估计，可以做一下变换：同理可以(ky)计算得到频率的不确定度：第三十页，共32页。降低误差降低误差(wch)的手段的手段 使用滤波等手段(shudun)对正弦序列进行预处理后再进行拟合将是降低拟合误差界的一个有效手段(shudun)。有学者专门针对正弦波采样序列提出了一种滤波器，可用于正弦拟合的预处理，主要用来消除谐波因素影响，以便降低拟合误差。其特点是理论上可以滤除全部偶次谐波和任意指定的奇次谐波，且对4个待拟合的正弦波模型参数没有影响。第三十一页，共32页。谢谢 谢！谢！第三十二页，共32页。