在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日