一种噪声背景下准确估计超声流量计传播时间方法
一种噪声背景下准确估计超声流量计传播时间方法
【专利说明】一种噪声背景下准确估计超声流量计传播时间方法 所属技术领域
[0001] 本发明属于流量测量技术领域，涉及超声流量计。
【背景技术】
[0002] 超声流量计作为新兴的流量仪表近十几年发展迅速，在大型水利水电项目以及天 然气贸易结算领域已有广泛应用，其中时差法超声流量计作为计量仪表应用最多。精确确 定超声换能器接收信号的起始位置，得到精确的超声波传播时间是时差法超声流量计能够 准确测量流量的关键。目前，传统测量方法中普遍采用的是双门限电平法（简称双门限 法），首先通过非过零比较，确定使用哪个周期的信号作为计时起始位置，再通过过零比较 确定超声波传播时间。但这种方法的缺点是当有噪声存在时，过零比较得到的时间会有较 大波动，严重影响了流量计的精度，而且当噪声大到一定程度时，非过零比较无法准确确定 信号位置，从而导致错波发生，此时由于误差太大，流量计已不能使用。
[0003] 为了在有噪声时仍能获得准确的超声波传播时间，研宄者提出基于超声脉冲接收 信号模型对超声波传播时间进行估计的方法。目前超声脉冲接收信号的经验模型主要有两 种，高斯模型和混合指数模型。高斯模型用于描述比较对称的宽带高频信号。周方等学者 将此模型与高斯一牛顿法、模拟退火法、蚁群算法等相结合，证明了蚁群算法是一种有效的 回波参数估计方法，具有较高精度[1];针对混合指数模型，姚振静等人将此模型与卡尔曼 滤波及其优化算法进行结合求取超声波传播时间，用来提高超声测距精度[2]。虽然已有学 者基于超声脉冲接收信号模型提出了一些超声波传播时间估计算法，但研宄仅针对算法本 身，应用也仅限于超声测距领域。对于超声流量测量来说，时间测量精度要求远高于超声测 距，已有方法不能直接用来解决噪声背景下超声流量计的准确测时问题。
[0004] [1]周方，张小凤，张光斌，超声回波参数的蚁群算法估计，陕西师范大学学报（自 然科学版），2012,40 (2)，35 ~40
[0005] [2]Zhen-JingYao,Qing-HaoMeng,MingZeng,Improvementintheaccuracyof estimatingthetime-〇f-flightinanultrasonicrangingsystemusingmultiple square-rootunscentedKalmanfilters,ReviewofScientificInstruments，2010,81， 104901-1 ~104901-7
【发明内容】
[0006] 本发明为解决噪声背景下超声流量计中超声波传播时间准确测量问题，基于超声 脉冲接收信号混合指数模型，提供一种超声波传播时间准确估计方法，达到降低测量误差， 提高超声流量计测量精度的目的。
[0007] 本发明的技术方案如下：
[0008] -种噪声背景下准确估计超声流量计传播时间方法，包括下列步骤：：
[0009] 步骤一：进行超声流量计实验，根据试验结果，建立超声脉冲接收信号离散模型如 下：
[0010] Am (x) =A(kts)sin[2 3ifc (kts-t) + 0 ] (1)
[0011] 其中，
[0012] 式中，Am(X)为超声脉冲接收离散模型信号，k为采样点数序号（k= 1,2，…N)，N 为采样点数，ts为采样时间间隔，f。为超声换能器的中心频率，0为初始相位角，A。为接收 信号幅值，T和m是超声换能器的两个特性参数，t为超声波传播时间，U(kts-T)为单位 阶跃信号，当采样频率、采样点数、超声换能器及其发射频率确定后，模型中的k、N、ts、f。和 9即为定值，随着被测流量改变，模型信号Am(X)中的和t会相应变化，其中，参数 向量x=[心mTt]，对x估计得越准确，模型信号越逼近实验获得的真实接收信号，利用 下面的遗传一蚁群算法对x中这四个参数进行估计，最终t的最优估计值即为超声波传播 时间；
[0013] 步骤二：根据实验采样获得的超声接收信号与模型信号Am(X)确定目标函数 f(x):基于最小二乘思想确定目标函数，将4中各采样点电压与Am(X)中采样点电压之差的 平方和作为目标函数f(x)，当目标函数f(x)取得最小值时，x为最优估计值；
[0014] 步骤三：基于目标函数f(x)构造适应度函数eval(x)，目标函数值越小，适应度函 数值就越大：
[0016] 其中，a是使得a_f(Xi)为正的任意实数，n为种群数目，n值设定越大，收敛精度 越高，收敛速度越慢，角标i代表第i个种群个体，在参数边界内随机产生；
[0017] 步骤四：设定x中四个参数边界，根据适应度函数值确定每个个体被复制的次数： 首先确定每个个体的预复制次数，它等于相应个体的适应度函数值与种群数目乘积的整数 部分，设R为总的预复制次数与种群数目之差，为了保持种群数目不变，从最优个体到第R 个优秀个体，每个个体的复制数目加1，即适应度函数值较大的个体，被复制的概率越大；
[0018] 步骤五：交叉、变异迭代获得优化解：从最优解到最差解，将种群中每两个相邻个 体分为一组，选定一个固定的判断概率〇. 9,并随机生成一个在0到1之间的实数P。，如果 P。小于判断概率〇. 9则进行交叉，否则不交叉，为获得最优个体，每次交叉、变异之后都用上 一代的最优个体代替下一代的最差个体；步骤二至步骤五为遗传算法，通过迭代获得优化 解作为后续蚁群算法的信息素初始条件；
[0019] 步骤六：利用上述的遗传算法得到的优化解初始化蚁群信息素；
[0020]步骤七：蚂蚁移动步长a初始化，根据实验获得的超声接收信号确定参数向量X 的变化范围；
[0021] 步骤八：判断蚂蚁是否转移及得到新的优化解：在生成新的参数解时，分别对参 数向量X中的每个参数单独进行，首先对Ac!进行，令
，如果/?) 小于/(41)，则
，如果/(〇小于
;其中，；^为第t代第s只蚂蚁的参数向量 值，xf1为第t_l代迭代之后按照信息素更新第s只蚂蚁最好的参数向量值，(1\在[_a，a] 之间随机产生；同理，对T、m和t三个参数进行优化，得到新的优化解；
[0022] 步骤九：信息素挥发和更新：按照下列规则设定信息素挥发系数：当迭代次数较 小时，信息素挥发系数设定较大，有利于全局寻优；当迭代次数较大时，信息素挥发系数设 定较小，有利于提高局部收敛精度、加快收敛速度；重新计算信息素，蚂蚁重新分配，返回步 骤七进行迭代计算，迭代结束判据设为：在算法迭代t代后进行判断，如果第t代的最优目 标函数值与第t-j(j= 1，2, 3, 4, 5)代的最优目标函数值之差的绝对值均小于1(T7,则迭代 结束，此时参数向量x为最优解，其中的t即为超声波传播时间最优估计值。
[0023] 本发明将遗传算法和蚁群算法相结合，利用遗传算法的全局寻优能力使得在噪声 背景下仍能找到全局最优解，且遗传算法收敛速度比蚁群算法快。同时结合蚁群算法的高 度局部寻优能力，提高参数估计的精确度。
【具体实施方式】
[0024] 以下将以实验获得的气体超声接收信号为目标对象，结合技术方案中的步骤方 法，给出各步的参数设置和实施方法：
[0025] 步骤一 ：DN100单声道气体超声流量计的超声换能器中心频率为120KHZ，实验采 用泰克TBS1102示波器采集超声换能器接收到的原始信号，时基选择lOus/div，整个示波 器窗口采样时间为l〇us/div*10格=lOOus，在2500个数据点存储深度下，相邻采样点的时 间间隔为lOOus/2. 5k= 40ns，则超声脉冲接收信号混合指数模型[2]中（式（1) (2))的相 关参数设置为：采样点序号k= 1，2,…2500,采样时间间隔ts= 40ns，换能器中心频率 = 120KHz，初始相位角0=0rad。其余参数AQ、T、m和t为待估计参数，定义为参数向量x= [A0mTt] 0
[0026] 步骤二：根据实验采样获得的超声接收信号与模型信号Am(X)确定目标函数 f(x) 〇
[0028] 其中，是示波器采样获得的真实超声接收信号，Am(X)是基于超声脉冲接收信号 混合指数模型重构的信号，x= [A。mTt]为待估计参数向量，采样点数N= 2500,角标 k代表在第k个采样点处的信号。
[0029] 步骤三：基于目标函数f(x)构造适应度函数eval(x)。式⑶中，设定a= 100， 种群数目n= 20,种群个体由Matlab中rand函数随机生成。
[0030] 步骤四：通过分析示波器采样超声接收信号特点，将x中四个参数边界分别设为 A0G[0? 025,0. 045]V，mG[3.4,3. 8]，TG[12. 3, 12.7]，tG[384, 439]us。根据适应度 函数值确定每个个体被复制的次数乂：
[0031] Mj=abs(eval(xj) *n)
[0033] Mi是第i个种群个体的预复制数目，R是种群数目与总的预复制数目之差。为了 保持种群数目不变，从最优个体到第R个优秀个体，每个个体的复制数目加1。
[0034] 步骤五：交叉、变异迭代获得优化解：从最优解到最差解，将种群中每两个相邻个 体分为一组，选定一个固定的判断概率〇. 9,采用Matlab中rand函数随机生成一个 在0到 1之间的实数P。，如果P。小于判断概率〇. 9则进行交叉，否则不交叉。
[0035] 假设存在两个相邻个体，编码值为&和x2,t为迭代次数，以概率P。做如下交叉：
[0038] 为获得最优个体，每次交叉、变异之后都用上一代的最优个体代替下一代的最差 个体。
[0039] 遗传算法迭代结束判据与步骤十中蚁群算法迭代结束判据相同。该部分迭代结束 获得的优化解将作为蚁群算法的信息素初始条件。
[0040] 步骤六：利用遗传算法得到的优化解初始化蚁群信息素，信息素公式为：
[0041] Dj=fi-fmin
[0044] 其中，Di为第i个种群个体结果与最好个体结果的欧几里得距离，X是标准偏差 取0. 005，ni为蚁群中第i个种群个体参数向量值的信息素值。
[0045] 蚁群算法通过信息素值更新种群的方法与前述遗传算法中通过适应度函数值实 现种群复制的方法相同。
[0046] 步骤七：蚂蚁移动步长a初始化：根据示波器采样信号变化特点，确定心、T和m 的最大变化范围分别为〇. 02,0. 4,0. 4,实验所用的超声换能器中心频率为120KHz，即周期 为8. 33ys，故设置蚂蚁移动步长a= [0. 02 0. 4 0. 4 8. 33]。
[0047] 步骤八：判断蚂蚁是否转移及得到新的优化解：在生成新的参数解时，分别对参 数向量x中的每个参数单独进行。首先对怂进行，令
，如果/(<) 小于 _
_ ，如果/(x〈）小于
=其中，（为第t代第s只蚂蚁的参数向 量值，xf1为第t-1代迭代之后按照信息素更新第S只蚂蚁最好的参数向量值。（1\在 [_a，a]之间随机产生。
[0048] 同理，对T、m和t三个参数进行优化，得到新的优化解。
[0049] 步骤九：信息素挥发、更新：信息素挥发系数直接用于控制蚂蚁移动步长，进而决 定了算法的全局寻优能力和收敛速度。信息素挥发系数设定如下：
[0050] 如果迭代次数t〈 = 3
[0051]a(1)=a(l)*〇. 92，a(2)=a(2)*0. 95，a(3)=a(3)*0. 95，a(4)= a⑷ *〇? 8
[0052] 否则
[0053] a(1)=a(l)*〇. 7,a(2) =a(2)*0. 8,a(3) =a(3)*0. 8,a(4) =a(4)*0. 7
[0054] 根据信息素挥发系数设定方法重新计算信息素，蚂蚁重新分配，同步骤六。
[0055] 步骤十：为了避免算法过早停滞、收敛于非最优解，将迭代结束判据设为：在算法 迭代t代后进行判断，如果第t代的最优目标函数值与第t_j(j= 1，2, 3, 4, 5)代的最优目 标函数值之差的绝对值均小于1〇_7,则迭代结束。此时参数向量x为最优解，其中的t即 为超声波传播时间最优估计值。
[0056] 为验证本发明在噪声背景下估计超声波传播时间的准确性，首先开展了仿真实 验。将本发明提供的遗传一蚁群算法与传统双门限电平法和同样基于超声脉冲接收信号模 型的遗传算法、蚁群算法分别进行比较，不同信噪比下对超声波传播时间估计结果对比如 表1所示。
[0057] 表1不同信噪比下各方法对比结果
[0059]
[0060] 结果表明：无噪声时，双门限法的平均相对误差均小于模型法，说明测时精度优于 模型法；而在有噪声时，模型法优于双门限法。在15dB时双门限法错波次数已达5个，说明 此时该方法已不适用于超声波传播时间的测量。从基于模型法的三种算法对比来看，本发 明提供的遗传一蚁群算法无论是在收敛精度、收敛速度、错波次数和抗噪声能力方面，均优 于遗传算法和蚁群算法，验证了本发明方法在噪声背景下准确估计超声波传播时间的可行 性。
[0061] 为进一步验证本发明，对DN100,声路角为45°的单声道气体超声流量计开展实 流实验，比较采用遗传一蚁群算法与传统双门限法进行超声流量测量的准确性，结果如表2 所示。
[0062] 表2实流实验结果
[0063]
[0065]由于气体流动的不稳定性使得超声换能器原始接收信号噪声较大，易造成双门限 法错波，从而导致平均相对误差和重复性均较差。而本发明提供的遗传一蚁群算法对噪声 的适应性较强，平均相对误差和重复性均优于双门限法，特别是在低流速测量时效果明显。 实验结果再次表明本发明方法用于解决噪声背景下精确测时问题的优越性。
【主权项】
1. 一种噪声背景下准确估计超声流量计传播时间方法，包括下列步骤：： 步骤一：进行超声流量计实验，根据试验结果，建立超声脉冲接收信号离散模型如下： Am (x) = A (kts) sin [2 π fc (kts- τ ) + θ ] (I) 其中，式中，A111(X)为超声脉冲接收离散模型信号，k为采样点数序号（k = 1，2，"·Ν)，N为 采样点数，ts为采样时间间隔，f。为超声换能器的中心频率，Θ为初始相位角，Atl为接收信 号幅值，T和m是超声换能器的两个特性参数，τ为超声波传播时间，u(kts-〇为单位阶 跃信号，当采样频率、采样点数、超声换能器及其发射频率确定后，模型中的k、N、ts、f^P0 即为定值，随着被测流量改变，模型信号Am(X)中的ApKm和τ会相应变化，其中，参数向 量X = [Atl m T τ]，对X估计得越准确，模型信号越逼近实验获得的真实接收信号，利用下 面的遗传一蚁群算法对X中这四个参数进行估计，最终τ的最优估计值即为超声波传播时 间； 步骤二：根据实验采样获得的超声接收信号Ae与模型信号Am(X)确定目标函数f(x): 基于最小二乘思想确定目标函数，将Ae中各采样点电压与A m (X)中采样点电压之差的平方 和作为目标函数f (X)，当目标函数f (X)取得最小值时，X为最优估计值； 步骤三：基于目标函数f (X)构造适应度函数eval (X)，目标函数值越小，适应度函数值 就越大：其中，a是使得a-f (Xi)为正的任意实数，η为种群数目，η值设定越大，收敛精度越高， 收敛速度越慢，角标i代表第i个种群个体，在参数边界内随机产生； 步骤四：设定X中四个参数边界，根据适应度函数值确定每个个体被复制的次数：首 先确定每个个体的预复制次数，它等于相应个体的适应度函数值与种群数目乘积的整数部 分，设R为总的预复制次数与种群数目之差，为了保持种群数目不变，从最优个体到第R个 优秀个体，每个个体的复制数目加1，即适应度函数值较大的个体，被复制的概率越大； 步骤五：交叉、变异迭代获得优化解：从最优解到最差解，将种群中每两个相邻个体分 为一组，选定一个固定的判断概率0. 9,并随机生成一个在O到1之间的实数P。，如果P。小 于判断概率0. 9则进行交叉，否则不交叉，为获得最优个体，每次交叉、变异之后都用上一 代的最优个体代替下一代的最差个体；步骤二至步骤五为遗传算法，通过迭代获得优化解 作为后续蚁群算法的信息素初始条件； 步骤六：利用上述的遗传算法得到的优化解初始化蚁群信息素； 步骤七：蚂蚁移动步长α初始化，根据实验获得的超声接收信号确定参数向量X的变 化范围； 步骤八：判断蚂蚁是否转移及得到新的优化解：在生成新的参数解时，分别对参数 向量X中的每个参数单独进行，首先对Atl进行，令< =.(1+[叫,，〇,〇,〇]，如果/〇〇 小于/(彳V则·<=·4-^[^,ο,ο,ο]，否则令λ-χ-[料,〇,〇,〇]，如果/(〇小于 /(#)，则-[#，〇,〇,〇]，否则χχ;其中，< 为第t代第S只蚂蚁的参数向 量值，为第t-ι代迭代之后按照信息素更新第S只蚂蚁最好的参数向量值，（1\在 [-α，α ]之间随机产生；同理，对T、m和τ三个参数进行优化，得到新的优化解； 步骤九：信息素挥发和更新：按照下列规则设定信息素挥发系数：当迭代次数较小时， 信息素挥发系数设定较大，有利于全局寻优；当迭代次数较大时，信息素挥发系数设定较 小，有利于提高局部收敛精度、加快收敛速度；重新计算信息素，蚂蚁重新分配，返回步骤七 进行迭代计算，迭代结束判据设为：在算法迭代t代后进行判断，如果第t代的最优目标函 数值与第t-j (j = 1，2, 3, 4, 5)代的最优目标函数值之差的绝对值均小于ΚΓ7,则迭代结束， 此时参数向量X为最优解，其中的τ即为超声波传播时间最优估计值。
【专利摘要】本发明涉及一种噪声背景下准确估计超声流量计传播时间方法，包括：进行超声流量计实验，根据试验结果，建立超声脉冲接收信号离散模型，确定包括超声波传播时间的参数向量x；根据实验采样获得的超声接收信号与模型信号确定目标函数：确定目标函数；构造适应度函数；进行遗传算法，通过迭代获得优化解作为后续蚁群算法的信息素初始条件；利用上述的遗传算法得到的优化解初始化蚁群信息素；根据实验获得的超声接收信号确定参数向量x的变化范围；判断蚂蚁是否转移及得到新的优化解；信息素挥发和更新，进行迭代计算，得到超声波传播时间最优估计值。本发明能够提高参数估计的精确度。
【IPC分类】G01F25/00, G06N3/12
【公开号】CN104897248
【申请号】CN201510324809
【发明人】郑丹丹, 侯惠让, 张涛
【申请人】天津大学
【公开日】2015年9月9日
【申请日】2015年6月12日
【专利说明】一种噪声背景下准确估计超声流量计传播时间方法 所属技术领域
[0001] 本发明属于流量测量技术领域，涉及超声流量计。
【背景技术】
[0002] 超声流量计作为新兴的流量仪表近十几年发展迅速，在大型水利水电项目以及天 然气贸易结算领域已有广泛应用，其中时差法超声流量计作为计量仪表应用最多。精确确 定超声换能器接收信号的起始位置，得到精确的超声波传播时间是时差法超声流量计能够 准确测量流量的关键。目前，传统测量方法中普遍采用的是双门限电平法（简称双门限 法），首先通过非过零比较，确定使用哪个周期的信号作为计时起始位置，再通过过零比较 确定超声波传播时间。但这种方法的缺点是当有噪声存在时，过零比较得到的时间会有较 大波动，严重影响了流量计的精度，而且当噪声大到一定程度时，非过零比较无法准确确定 信号位置，从而导致错波发生，此时由于误差太大，流量计已不能使用。
[0003] 为了在有噪声时仍能获得准确的超声波传播时间，研宄者提出基于超声脉冲接收 信号模型对超声波传播时间进行估计的方法。目前超声脉冲接收信号的经验模型主要有两 种，高斯模型和混合指数模型。高斯模型用于描述比较对称的宽带高频信号。周方等学者 将此模型与高斯一牛顿法、模拟退火法、蚁群算法等相结合，证明了蚁群算法是一种有效的 回波参数估计方法，具有较高精度[1];针对混合指数模型，姚振静等人将此模型与卡尔曼 滤波及其优化算法进行结合求取超声波传播时间，用来提高超声测距精度[2]。虽然已有学 者基于超声脉冲接收信号模型提出了一些超声波传播时间估计算法，但研宄仅针对算法本 身，应用也仅限于超声测距领域。对于超声流量测量来说，时间测量精度要求远高于超声测 距，已有方法不能直接用来解决噪声背景下超声流量计的准确测时问题。
[0004] [1]周方，张小凤，张光斌，超声回波参数的蚁群算法估计，陕西师范大学学报（自 然科学版），2012,40 (2)，35 ~40
[0005] [2]Zhen-JingYao,Qing-HaoMeng,MingZeng,Improvementintheaccuracyof estimatingthetime-〇f-flightinanultrasonicrangingsystemusingmultiple square-rootunscentedKalmanfilters,ReviewofScientificInstruments，2010,81， 104901-1 ~104901-7
【发明内容】
[0006] 本发明为解决噪声背景下超声流量计中超声波传播时间准确测量问题，基于超声 脉冲接收信号混合指数模型，提供一种超声波传播时间准确估计方法，达到降低测量误差， 提高超声流量计测量精度的目的。
[0007] 本发明的技术方案如下：
[0008] -种噪声背景下准确估计超声流量计传播时间方法，包括下列步骤：：
[0009] 步骤一：进行超声流量计实验，根据试验结果，建立超声脉冲接收信号离散模型如 下：
[0010] Am (x) =A(kts)sin[2 3ifc (kts-t) + 0 ] (1)
[0011] 其中，
[0012] 式中，Am(X)为超声脉冲接收离散模型信号，k为采样点数序号（k= 1,2，…N)，N 为采样点数，ts为采样时间间隔，f。为超声换能器的中心频率，0为初始相位角，A。为接收 信号幅值，T和m是超声换能器的两个特性参数，t为超声波传播时间，U(kts-T)为单位 阶跃信号，当采样频率、采样点数、超声换能器及其发射频率确定后，模型中的k、N、ts、f。和 9即为定值，随着被测流量改变，模型信号Am(X)中的和t会相应变化，其中，参数 向量x=[心mTt]，对x估计得越准确，模型信号越逼近实验获得的真实接收信号，利用 下面的遗传一蚁群算法对x中这四个参数进行估计，最终t的最优估计值即为超声波传播 时间；
[0013] 步骤二：根据实验采样获得的超声接收信号与模型信号Am(X)确定目标函数 f(x):基于最小二乘思想确定目标函数，将4中各采样点电压与Am(X)中采样点电压之差的 平方和作为目标函数f(x)，当目标函数f(x)取得最小值时，x为最优估计值；
[0014] 步骤三：基于目标函数f(x)构造适应度函数eval(x)，目标函数值越小，适应度函 数值就越大：
[0016] 其中，a是使得a_f(Xi)为正的任意实数，n为种群数目，n值设定越大，收敛精度 越高，收敛速度越慢，角标i代表第i个种群个体，在参数边界内随机产生；
[0017] 步骤四：设定x中四个参数边界，根据适应度函数值确定每个个体被复制的次数： 首先确定每个个体的预复制次数，它等于相应个体的适应度函数值与种群数目乘积的整数 部分，设R为总的预复制次数与种群数目之差，为了保持种群数目不变，从最优个体到第R 个优秀个体，每个个体的复制数目加1，即适应度函数值较大的个体，被复制的概率越大；
[0018] 步骤五：交叉、变异迭代获得优化解：从最优解到最差解，将种群中每两个相邻个 体分为一组，选定一个固定的判断概率〇. 9,并随机生成一个在0到1之间的实数P。，如果 P。小于判断概率〇. 9则进行交叉，否则不交叉，为获得最优个体，每次交叉、变异之后都用上 一代的最优个体代替下一代的最差个体；步骤二至步骤五为遗传算法，通过迭代获得优化 解作为后续蚁群算法的信息素初始条件；
[0019] 步骤六：利用上述的遗传算法得到的优化解初始化蚁群信息素；
[0020]步骤七：蚂蚁移动步长a初始化，根据实验获得的超声接收信号确定参数向量X 的变化范围；
[0021] 步骤八：判断蚂蚁是否转移及得到新的优化解：在生成新的参数解时，分别对参 数向量X中的每个参数单独进行，首先对Ac!进行，令
，如果/?) 小于/(41)，则
，如果/(〇小于
;其中，；^为第t代第s只蚂蚁的参数向量 值，xf1为第t_l代迭代之后按照信息素更新第s只蚂蚁最好的参数向量值，(1\在[_a，a] 之间随机产生；同理，对T、m和t三个参数进行优化，得到新的优化解；
[0022] 步骤九：信息素挥发和更新：按照下列规则设定信息素挥发系数：当迭代次数较 小时，信息素挥发系数设定较大，有利于全局寻优；当迭代次数较大时，信息素挥发系数设 定较小，有利于提高局部收敛精度、加快收敛速度；重新计算信息素，蚂蚁重新分配，返回步 骤七进行迭代计算，迭代结束判据设为：在算法迭代t代后进行判断，如果第t代的最优目 标函数值与第t-j(j= 1，2, 3, 4, 5)代的最优目标函数值之差的绝对值均小于1(T7,则迭代 结束，此时参数向量x为最优解，其中的t即为超声波传播时间最优估计值。
[0023] 本发明将遗传算法和蚁群算法相结合，利用遗传算法的全局寻优能力使得在噪声 背景下仍能找到全局最优解，且遗传算法收敛速度比蚁群算法快。同时结合蚁群算法的高 度局部寻优能力，提高参数估计的精确度。
【具体实施方式】
[0024] 以下将以实验获得的气体超声接收信号为目标对象，结合技术方案中的步骤方 法，给出各步的参数设置和实施方法：
[0025] 步骤一 ：DN100单声道气体超声流量计的超声换能器中心频率为120KHZ，实验采 用泰克TBS1102示波器采集超声换能器接收到的原始信号，时基选择lOus/div，整个示波 器窗口采样时间为l〇us/div*10格=lOOus，在2500个数据点存储深度下，相邻采样点的时 间间隔为lOOus/2. 5k= 40ns，则超声脉冲接收信号混合指数模型[2]中（式（1) (2))的相 关参数设置为：采样点序号k= 1，2,…2500,采样时间间隔ts= 40ns，换能器中心频率 = 120KHz，初始相位角0=0rad。其余参数AQ、T、m和t为待估计参数，定义为参数向量x= [A0mTt] 0
[0026] 步骤二：根据实验采样获得的超声接收信号与模型信号Am(X)确定目标函数 f(x) 〇
[0028] 其中，是示波器采样获得的真实超声接收信号，Am(X)是基于超声脉冲接收信号 混合指数模型重构的信号，x= [A。mTt]为待估计参数向量，采样点数N= 2500,角标 k代表在第k个采样点处的信号。
[0029] 步骤三：基于目标函数f(x)构造适应度函数eval(x)。式⑶中，设定a= 100， 种群数目n= 20,种群个体由Matlab中rand函数随机生成。
[0030] 步骤四：通过分析示波器采样超声接收信号特点，将x中四个参数边界分别设为 A0G[0? 025,0. 045]V，mG[3.4,3. 8]，TG[12. 3, 12.7]，tG[384, 439]us。根据适应度 函数值确定每个个体被复制的次数乂：
[0031] Mj=abs(eval(xj) *n)
[0033] Mi是第i个种群个体的预复制数目，R是种群数目与总的预复制数目之差。为了 保持种群数目不变，从最优个体到第R个优秀个体，每个个体的复制数目加1。
[0034] 步骤五：交叉、变异迭代获得优化解：从最优解到最差解，将种群中每两个相邻个 体分为一组，选定一个固定的判断概率〇. 9,采用Matlab中rand函数随机生成一个 在0到 1之间的实数P。，如果P。小于判断概率〇. 9则进行交叉，否则不交叉。
[0035] 假设存在两个相邻个体，编码值为&和x2,t为迭代次数，以概率P。做如下交叉：
[0038] 为获得最优个体，每次交叉、变异之后都用上一代的最优个体代替下一代的最差 个体。
[0039] 遗传算法迭代结束判据与步骤十中蚁群算法迭代结束判据相同。该部分迭代结束 获得的优化解将作为蚁群算法的信息素初始条件。
[0040] 步骤六：利用遗传算法得到的优化解初始化蚁群信息素，信息素公式为：
[0041] Dj=fi-fmin
[0044] 其中，Di为第i个种群个体结果与最好个体结果的欧几里得距离，X是标准偏差 取0. 005，ni为蚁群中第i个种群个体参数向量值的信息素值。
[0045] 蚁群算法通过信息素值更新种群的方法与前述遗传算法中通过适应度函数值实 现种群复制的方法相同。
[0046] 步骤七：蚂蚁移动步长a初始化：根据示波器采样信号变化特点，确定心、T和m 的最大变化范围分别为〇. 02,0. 4,0. 4,实验所用的超声换能器中心频率为120KHz，即周期 为8. 33ys，故设置蚂蚁移动步长a= [0. 02 0. 4 0. 4 8. 33]。
[0047] 步骤八：判断蚂蚁是否转移及得到新的优化解：在生成新的参数解时，分别对参 数向量x中的每个参数单独进行。首先对怂进行，令
，如果/(<) 小于 _
_ ，如果/(x〈）小于
=其中，（为第t代第s只蚂蚁的参数向 量值，xf1为第t-1代迭代之后按照信息素更新第S只蚂蚁最好的参数向量值。（1\在 [_a，a]之间随机产生。
[0048] 同理，对T、m和t三个参数进行优化，得到新的优化解。
[0049] 步骤九：信息素挥发、更新：信息素挥发系数直接用于控制蚂蚁移动步长，进而决 定了算法的全局寻优能力和收敛速度。信息素挥发系数设定如下：
[0050] 如果迭代次数t〈 = 3
[0051]a(1)=a(l)*〇. 92，a(2)=a(2)*0. 95，a(3)=a(3)*0. 95，a(4)= a⑷ *〇? 8
[0052] 否则
[0053] a(1)=a(l)*〇. 7,a(2) =a(2)*0. 8,a(3) =a(3)*0. 8,a(4) =a(4)*0. 7
[0054] 根据信息素挥发系数设定方法重新计算信息素，蚂蚁重新分配，同步骤六。
[0055] 步骤十：为了避免算法过早停滞、收敛于非最优解，将迭代结束判据设为：在算法 迭代t代后进行判断，如果第t代的最优目标函数值与第t_j(j= 1，2, 3, 4, 5)代的最优目 标函数值之差的绝对值均小于1〇_7,则迭代结束。此时参数向量x为最优解，其中的t即 为超声波传播时间最优估计值。
[0056] 为验证本发明在噪声背景下估计超声波传播时间的准确性，首先开展了仿真实 验。将本发明提供的遗传一蚁群算法与传统双门限电平法和同样基于超声脉冲接收信号模 型的遗传算法、蚁群算法分别进行比较，不同信噪比下对超声波传播时间估计结果对比如 表1所示。
[0057] 表1不同信噪比下各方法对比结果
[0059]
[0060] 结果表明：无噪声时，双门限法的平均相对误差均小于模型法，说明测时精度优于 模型法；而在有噪声时，模型法优于双门限法。在15dB时双门限法错波次数已达5个，说明 此时该方法已不适用于超声波传播时间的测量。从基于模型法的三种算法对比来看，本发 明提供的遗传一蚁群算法无论是在收敛精度、收敛速度、错波次数和抗噪声能力方面，均优 于遗传算法和蚁群算法，验证了本发明方法在噪声背景下准确估计超声波传播时间的可行 性。
[0061] 为进一步验证本发明，对DN100,声路角为45°的单声道气体超声流量计开展实 流实验，比较采用遗传一蚁群算法与传统双门限法进行超声流量测量的准确性，结果如表2 所示。
[0062] 表2实流实验结果
[0063]
[0065]由于气体流动的不稳定性使得超声换能器原始接收信号噪声较大，易造成双门限 法错波，从而导致平均相对误差和重复性均较差。而本发明提供的遗传一蚁群算法对噪声 的适应性较强，平均相对误差和重复性均优于双门限法，特别是在低流速测量时效果明显。 实验结果再次表明本发明方法用于解决噪声背景下精确测时问题的优越性。
【主权项】
1. 一种噪声背景下准确估计超声流量计传播时间方法，包括下列步骤：： 步骤一：进行超声流量计实验，根据试验结果，建立超声脉冲接收信号离散模型如下： Am (x) = A (kts) sin [2 π fc (kts- τ ) + θ ] (I) 其中，式中，A111(X)为超声脉冲接收离散模型信号，k为采样点数序号（k = 1，2，"·Ν)，N为 采样点数，ts为采样时间间隔，f。为超声换能器的中心频率，Θ为初始相位角，Atl为接收信 号幅值，T和m是超声换能器的两个特性参数，τ为超声波传播时间，u(kts-〇为单位阶 跃信号，当采样频率、采样点数、超声换能器及其发射频率确定后，模型中的k、N、ts、f^P0 即为定值，随着被测流量改变，模型信号Am(X)中的ApKm和τ会相应变化，其中，参数向 量X = [Atl m T τ]，对X估计得越准确，模型信号越逼近实验获得的真实接收信号，利用下 面的遗传一蚁群算法对X中这四个参数进行估计，最终τ的最优估计值即为超声波传播时 间； 步骤二：根据实验采样获得的超声接收信号Ae与模型信号Am(X)确定目标函数f(x): 基于最小二乘思想确定目标函数，将Ae中各采样点电压与A m (X)中采样点电压之差的平方 和作为目标函数f (X)，当目标函数f (X)取得最小值时，X为最优估计值； 步骤三：基于目标函数f (X)构造适应度函数eval (X)，目标函数值越小，适应度函数值 就越大：其中，a是使得a-f (Xi)为正的任意实数，η为种群数目，η值设定越大，收敛精度越高， 收敛速度越慢，角标i代表第i个种群个体，在参数边界内随机产生； 步骤四：设定X中四个参数边界，根据适应度函数值确定每个个体被复制的次数：首 先确定每个个体的预复制次数，它等于相应个体的适应度函数值与种群数目乘积的整数部 分，设R为总的预复制次数与种群数目之差，为了保持种群数目不变，从最优个体到第R个 优秀个体，每个个体的复制数目加1，即适应度函数值较大的个体，被复制的概率越大； 步骤五：交叉、变异迭代获得优化解：从最优解到最差解，将种群中每两个相邻个体分 为一组，选定一个固定的判断概率0. 9,并随机生成一个在O到1之间的实数P。，如果P。小 于判断概率0. 9则进行交叉，否则不交叉，为获得最优个体，每次交叉、变异之后都用上一 代的最优个体代替下一代的最差个体；步骤二至步骤五为遗传算法，通过迭代获得优化解 作为后续蚁群算法的信息素初始条件； 步骤六：利用上述的遗传算法得到的优化解初始化蚁群信息素； 步骤七：蚂蚁移动步长α初始化，根据实验获得的超声接收信号确定参数向量X的变 化范围； 步骤八：判断蚂蚁是否转移及得到新的优化解：在生成新的参数解时，分别对参数 向量X中的每个参数单独进行，首先对Atl进行，令< =.(1+[叫,，〇,〇,〇]，如果/〇〇 小于/(彳V则·<=·4-^[^,ο,ο,ο]，否则令λ-χ-[料,〇,〇,〇]，如果/(〇小于 /(#)，则-[#，〇,〇,〇]，否则χχ;其中，< 为第t代第S只蚂蚁的参数向 量值，为第t-ι代迭代之后按照信息素更新第S只蚂蚁最好的参数向量值，（1\在 [-α，α ]之间随机产生；同理，对T、m和τ三个参数进行优化，得到新的优化解； 步骤九：信息素挥发和更新：按照下列规则设定信息素挥发系数：当迭代次数较小时， 信息素挥发系数设定较大，有利于全局寻优；当迭代次数较大时，信息素挥发系数设定较 小，有利于提高局部收敛精度、加快收敛速度；重新计算信息素，蚂蚁重新分配，返回步骤七 进行迭代计算，迭代结束判据设为：在算法迭代t代后进行判断，如果第t代的最优目标函 数值与第t-j (j = 1，2, 3, 4, 5)代的最优目标函数值之差的绝对值均小于ΚΓ7,则迭代结束， 此时参数向量X为最优解，其中的τ即为超声波传播时间最优估计值。
【专利摘要】本发明涉及一种噪声背景下准确估计超声流量计传播时间方法，包括：进行超声流量计实验，根据试验结果，建立超声脉冲接收信号离散模型，确定包括超声波传播时间的参数向量x；根据实验采样获得的超声接收信号与模型信号确定目标函数：确定目标函数；构造适应度函数；进行遗传算法，通过迭代获得优化解作为后续蚁群算法的信息素初始条件；利用上述的遗传算法得到的优化解初始化蚁群信息素；根据实验获得的超声接收信号确定参数向量x的变化范围；判断蚂蚁是否转移及得到新的优化解；信息素挥发和更新，进行迭代计算，得到超声波传播时间最优估计值。本发明能够提高参数估计的精确度。
【IPC分类】G01F25/00, G06N3/12
【公开号】CN104897248
【申请号】CN201510324809
【发明人】郑丹丹, 侯惠让, 张涛
【申请人】天津大学
【公开日】2015年9月9日
【申请日】2015年6月12日
最新回复(0)