重建储粮温度场的方法与流程

本发明涉及温度场重建技术领域，具体涉及到自适应鸟群方法用于粮食温度场的重建。

背景技术：

中国人口众多，对粮食有很大的需求，因此，粮食的存储是头等大事。霉变、虫害、水分异常等都可能引起储粮温度的变化，储粮区域的温度直接反应该区域的粮食存储状况，因此，储粮区域的温度监测非常重要。目前国内外仓储粮食的温度检测多采用接触测温法，接触式测温方法简单直观，但测量时安全性较差，精度低，而且维护成本也高。非接触式测温是被测物体不需要与测量物体进行接触，通过一定的测量方式得到被测物体的温度。

王交锋等人利用声学法测温研究航空发动机的燃烧室出口的温度场，将16个声波收发器均匀放置在燃烧室的外壁，对数据进行采集，利用最小二乘法重建温度场，并对重建后的数据利用双三次插值法进行插值，获得更加细致的图像，取得了较好的温度场重建效果；沈国清等人提出将代数迭代重建法应用于温度场的重建，在实验条件下，利用重建法得到每个像素中心的温度值，利用双三次插值法进行插值，得到整个测量区域的温度场分布，得到的重建结果较为理想；田丰等人针对温度场重建法进行了深入研究，提出将高斯函数和正则化法相结合进行温度场重建，仿真实验表明，该法具有较强的抗噪能力；颜华等人提出正则化参数自适应选取的温度场重建法，该方法是一种像素数远多于声波路径数、且正则化参数可自适应调整，兼顾温度场细节重建和噪声抑制的声学ct温度场重建法，实验条件下，与常用的确定正则化参数的方法相比，该文中所使用的方法对应着更小的温度场重建误差，用于仓储粮食温度分布监测等对重建质量有较高要求的应用场合。

上述的温度场重建中所使用的最小二乘法、代数迭代重建法等传统方法重建速度慢。

技术实现要素：

本发明所要解决的主要技术问题在于克服上述传统重建方法的不足，提供一种重建精度高、收敛速度快的重建储粮温度场的方法。

解决上述技术问题所采用的技术方案由下述步骤组成：

(1)设置声波传感器

在被测储粮区域的四周设置声波传感器，形成声波路径，采用快速互相关结合抽样率变换的时延法得到有效路径的声波传播时间。

(2)用自适应鸟群方法重建储粮温度场

1)建立适应度函数

其中，是系数矩阵，代表m×n的维权因子，amn表示第m条声波路径穿过第n个像素的长度，m、n为有限的正整数，x＝(v1v2…vn)^t在温度场求解中代表n维声速分布向量，vn代表声波穿过第n个像素的平均速度，b＝(t1t2…tm)^t在温度场求解中代表m维的各有效传播路径上的声波传播时间向量，tm代表第m条声波的声波传播时间。

2)根据迭代次数k和飞行频率fq的关系确定鸟处于觅食、警戒状态或者处于飞行状态,k∈[1,...,m]，最大迭代次数m取值1000，fq为有限的正整数；

如果满足：k对fq取余不为0，鸟处于觅食或警戒状态，执行步骤3)，否则处于飞行状态，执行步骤4)。

3)在(0,1)内产生一个随机数，觅食概率r∈[0,1]，当产生的随机数小于r时，鸟处于觅食状态，否则处于警戒状态。

鸟处于觅食状态时，其行为满足公式(2)：

p＝1+exp(-(1-pfitid/pfiti))(3)

q＝1+exp(-(1-sumfit/pfiti))(4)

其中i∈[1,...,n]，鸟群数目n的值取为有限的正整数，j∈[1,...,d]，d为有限的正整数，rand(0,1)表示在[0,1]范围内独立均匀分布的随机数，表示第i只鸟在第j维时第k次迭代的位置，pi,j表示鸟群中第i只鸟在第j维的最佳位置，p是认知加速因子，q是社会加速因子，gj表示鸟群在第j维时的最佳位置，pfiti表示第i只鸟的最佳适应度值,pfitid表示第i只鸟局部极值的最佳适应度值，sumfit表示整个鸟群的最佳适应度值之和。

鸟处于警戒状态时，其行为满足公式为：

其中，h∈[1,...,n]且h≠i，rand(-1,1)表示在[-1,1]范围内独立均匀分布的随机数，meanj表示整个鸟群在第j维的平均位置，ε是一个极小的常数，用来避免分母出现0的情况，a1为一只鸟向鸟群中心移动过程中由环境引发的间接作用，a2为一只鸟向中间位置移动的直接影响因素，间接影响因子a1的值取0.001，直接影响因子a2的值取1。

确定鸟的个体适应度函数值，如果当前位置优于前次保留位置，则保留当前鸟的个体位置为最佳位置，并对鸟的位置进行更新。

4)鸟处于飞行行为时，鸟群分成生产者和乞讨者两类，生产者积极寻觅食物，乞讨者随机跟随一个生产者寻找食物，飞行行为中生产者和乞讨者的位置更新公式分别为：

其中，randn(0,1)代表产生服从均值为0，标准差为1的高斯分布的一个随机数，step表示移动步长，step∈[1,5]，fl(fl∈[0,2])表示乞讨者跟随生产者寻找食物。

确定鸟的个体适应度函数值，如果当前位置优于前次保留位置，则保留当前鸟的个体位置为最佳位置，并对鸟的位置进行更新。

5)当达到最大迭代次数m时，获得鸟群的最佳位置，得到声波穿过温度场的声速分布向量x，根据粮食中声波传播速度cm与温度t的关系式：

其中z为气体介质的声音常数，λ为声速转换因子，获得测量区域温度场每个像素中心的温度值。

(3)确定被测储粮区域的温度场

用四格点样条插值法进行插值，得到被测储粮区域的温度场分布。

在本发明的自适应鸟群方法重建储粮温度场步骤(2)的子步骤3)的公式(3)、(4)中，其中i∈[1,...,n]，鸟群数目n的值取为40。

在本发明的自适应鸟群方法重建储粮温度场步骤(2)的子步骤4)的公式(8)中，step表示移动步长，step∈[1,5]，step取值为2。

本发明采用自适应鸟群方法重建黄豆温度场，能准确地重建出所测黄豆区域的温度场分布，在鸟群的飞行行为中加入移动步长，并根据鸟的局部及全局极值，动态地改变学习因子的大小，使其在前期具有较强的自学能力，提高了全局搜索能力，后期强化鸟的社会学习能力，加快了收敛速度。本发明具有重建精度高、收敛速度快等优点，可用于重建储粮区域的温度场。

附图说明

图1是实施例1重建的被测黄豆区域的温度场的灰度图。

图2是实施例1重建的被测黄豆区域的温度场的等高线图。

具体实施方式

下面结合附图和实施例对本发明进一步详细说明，但本发明不限于这些实施例。

实施例1

以测量的黄豆区域为1.2m×1.2m的正方形区域，热点位置在区域中心，且为一个半径为20cm，中心点坐标为(0,0)的圆，温度为298k，被测黄豆区域的温度场重建方法如下：

1、设置声波传感器

在正方形被测储粮区域的四周均匀设置8个声波传感器，形成24条声波路径，被测储粮区域划分成16个像素，采用快速互相关结合抽样率变换的时延法得到有效路径的声波传播时间，快速互相关结合抽样率变换的时延法为已知方法，在声学法测温系统及声传播时间测量方法研究，沈阳工业大学,2009中公开。

2、自适应鸟群方法重建储粮温度场

1)建立适应度函数

其中，是系数矩阵，在温度场求解中代表m×n的维权因子，amn表示第m条声波路径穿过第n个像素的长度，m的值取24，n的值取16，x＝(v1v2…vn)^t在温度场求解中代表n维声速分布向量，vn代表声波穿过第n个像素的平均速度，b＝(t1t2…tm)^t在温度场求解中代表m维的各有效传播路径上的声波传播时间向量，tm代表第m条声波的声波传播时间。

2)根据迭代次数k和飞行频率fq的关系确定鸟处于觅食、警戒状态或者处于飞行状态。k∈[1,...,m]，最大迭代次数m取值1000，fq取值10。

如果满足：k对fq取余不为0,鸟处于觅食或警戒状态，执行步骤3)，否则处于飞行状态，执行步骤4)；

3)在(0,1)内产生一个随机数，觅食概率r∈[0,1]，当产生的随机数小于r时，鸟处于觅食状态，否则处于警戒状态。

鸟处于觅食状态时，其行为满足公式：

p＝1+exp(-(1-pfitid/pfiti))(3)

q＝1+exp(-(1-sumfit/pfiti))(4)

其中i∈[1,...,n]，鸟群数目n的值取40，j∈[1,...,16]，rand(0,1)表示在[0,1]范围内独立均匀分布的随机数，表示第i只鸟在第j维时第k次迭代的位置，pi,j表示鸟群中第i只鸟在第j维的最佳位置，p是认知加速因子，q是社会加速因子，gj表示鸟群在第j维时的最佳位置，pfiti表示第i只鸟的最佳适应度值,pfitid表示第i只鸟局部极值的最佳适应度值，sumfit表示整个鸟群的最佳适应度值之和。

鸟处于警戒状态时，其行为满足公式为：

其中，h∈[1,...,40]且h≠i，rand(-1,1)表示在[-1,1]范围内独立均匀分布的随机数，meanj表示整个鸟群在第j维的平均位置，ε是一个极小的常数，用来避免分母出现0的情况，a1为一只鸟向鸟群中心移动过程中由环境引发的间接作用，a2为一只鸟向中间位置移动的直接影响因素，间接影响因子a1的值取0.001，直接影响因子a2的值取1。

确定鸟的个体适应度函数值，如果当前位置优于前次保留位置，则保留当前鸟的个体位置为最佳位置，并对鸟的位置进行更新。

4)鸟处于飞行行为时，鸟群分成生产者和乞讨者两类，生产者积极寻觅食物，乞讨者随机跟随一个生产者寻找食物，飞行行为中生产者和乞讨者的位置更新公式分别为：

其中，randn(0,1)代表产生服从均值为0，标准差为1的高斯分布的一个随机数，step表示移动步长，step的值为2,fl(fl∈[0,2])表示乞讨者跟随生产者寻找食物。

确定鸟的个体适应度函数值，如果当前位置优于前次保留位置，则保留当前鸟的个体位置为最佳位置，并对鸟的位置进行更新。

5)当达到最大迭代次数m时，获得鸟群的最佳位置，得到声波穿过温度场的声速分布向量x，根据粮食中声波传播速度cm与温度t的关系式：

粮食中声波传播速度cm与温度t的关系式为已知关系式，在声学法仓储粮食温度检测关键技术的研究，沈阳工业大学,2012中公开，其中z为气体介质的声音常数，λ为声速转换因子，获得测量区域温度场每个像素中心的温度值。

3、确定被测储粮区域的温度场

在matlab软件中调用四格点样条插值法进行插值，得到被测储粮区域的温度场分布，重建效果以伪彩图和等高线图的形式展现出来，见图1、图2。由图1、2可见本发明能有效地重建出所测黄豆区域的温度场，准确地检测出所测区域的热点区域。

实施例2

以测量的黄豆区域为1.2m×1.2m的正方形区域，热点位置在区域中心，且为一个半径为20cm，中心点坐标为(0,0)的圆，温度为298k，被测黄豆区域的温度场重建方法如下：

1、设置声波传感器

设置声波传感器与实施例1相同。

步骤1与实施例1相同。

2、自适应鸟群方法重建储粮温度场

在自适应鸟群方法重建储粮温度场2的子步骤4)中，鸟处于飞行行为时，鸟群分成生产者和乞讨者两类，生产者积极寻觅食物，乞讨者随机跟随一个生产者寻找食物，飞行行为中生产者和乞讨者的位置更新公式分别为：

其中，randn(0,1)代表产生服从均值为0，标准差为1的高斯分布的一个随机数，step表示移动步长，step的值为1,fl(fl∈[0,2])表示乞讨者跟随生产者寻找食物。

确定鸟的个体适应度函数值，如果当前位置优于前次保留位置，则保留当前鸟的个体位置为最佳位置，并对鸟的位置进行更新。

该步骤的其它步骤与实施例1相同。

其它步骤与实施例1相同。

实施例3

以测量的黄豆区域为1.2m×1.2m的正方形区域，热点位置在区域中心，且为一个半径为20cm，中心点坐标为(0,0)的圆，温度为298k，被测黄豆区域的温度场重建方法如下：

1、设置声波传感器

设置声波传感器与实施例1相同。

步骤1与实施例1相同。

2、自适应鸟群方法重建储粮温度场

在自适应鸟群方法重建储粮温度场2的子步骤4)中，鸟处于飞行行为时，鸟群分成生产者和乞讨者两类，生产者积极寻觅食物，乞讨者随机跟随一个生产者寻找食物，飞行行为中生产者和乞讨者的位置更新公式分别为：

其中，randn(0,1)代表产生服从均值为0，标准差为1的高斯分布的一个随机数，step表示移动步长，step的值为5,fl(fl∈[0,2])表示乞讨者跟随生产者寻找食物。

确定鸟的个体适应度函数值，如果当前位置优于前次保留位置，则保留当前鸟的个体位置为最佳位置，并对鸟的位置进行更新。

该步骤的其它步骤与实施例1相同。

其它步骤与实施例1相同。