本发明涉及工业无损检测领域，特别是涉及一种各向异性金属薄板感应式磁声图像的仿真方法及系统。
背景技术：
感应式磁声(magnetoacoustictomographywithmagneticinduction，mat-mi)无损检测技术是一种绘制被测金属试件表面电导率分布图的无损检测新方法。其原理是将被测金属试件置于静磁场中，通过线圈引入交变磁场，试件中产生与交变磁场同频率的感应涡流；在感应涡流与静磁场的共同作用下，材料质点由于受到洛伦兹力的作用而引起周期性振动，并以超声波的形式向外传播，即磁声信号；利用电磁超声换能器在试件周围接收磁声信号并送入计算机内，可重建出试件表面的电导率空间分布，据此判断被测试件的缺陷情况并对其进行定位。
目前，针对mat-mi无损检测技术的研究大多是建立在被测材料电导率分布各向同性的前提下。但是在实际情况中，金属材料通常是各向异性介质。电导率各向异性对被测试件中感应涡流密度的分布、声源分布以及声压分布均有影响。
技术实现要素：
本发明的目的是提供一种各向异性金属薄板感应式磁声图像的仿真方法及系统，通过改变各个方向上的电导率，精确地仿真电导率分布各向异性金属薄板表面的力分布和表面波位移。
为实现上述目的，本发明提供了如下方案：
一种各向异性金属薄板感应式磁声图像的仿真方法，所述仿真方法包括：
建立电导率各向异性金属薄板仿真结构模型，所述仿真结构模型中包括仿真被测试件和仿真折线线圈；
根据所述仿真结构模型，设置所述仿真被测试件不同方向的电导率；
将所述仿真结构模型置于静磁场中进行仿真实验，得到感应涡流；
根据所述感应涡流，确定所述被测试件受到的静态洛伦兹力和交变洛伦兹力；
根据所述静态洛伦兹力和所述交变洛伦兹力得到表面波；
根据各不同方向的电导率，计算所述表面波沿对应方向上的波速；
根据表面波沿不同方向上的波速，得到对应方向上的表面波位移；
根据各所述表面波位移，得到表面波位移分布图。
可选的，所述根据所述仿真结构模型，设置所述被测试件不同方向的电导率，具体包括：
基于所述被测试件建立三维笛卡尔直角坐标系，其中，以所述被测试件的表面中心为原点，与所述被测试件水平面平行的为x轴，与x轴垂直的为y轴，与所述xoy面垂直的为z轴；
设置被测试件在x轴、y轴和z轴三个方向上的电导率。
可选的，所述将所述仿真结构模型置于静磁场中，计算所述被测试件受到的静态洛伦兹力和交变洛伦兹力，具体包括：
根据
其中，fs(h,t)为被测试件受到的静态洛伦兹力,σ为电导率，μ0为空间磁导率，j是虚部单位，bs为静磁场强度，z0为折线线圈与被测试件的垂直距离，k为径向积分变量，i为激励电流的幅值，ω为激励电流的角频率，j0(k)为自变量为k的第一类零阶bessel函数，j0(kr)为自变量为kr的第一类零阶bessel函数。
根据
其中，fd(h,t)为被测试件受到的交变洛伦兹力，μ0为空间磁导率，i为激励电流的幅值，k为径向积分变量，z0为折线线圈与被测试件的垂直距离，j是虚部单位，j0(k)为自变量为k的第一类零阶bessel函数，j0(kr)为自变量为kr的第一类零阶bessel函数。
可选的，所述根据所述电导率计算所述表面波沿不同方向的波速，具体包括：
表面波的传播仅位于yoz平面，与x轴无关；
根据
其中，vry为所述表面波沿y轴传播的波速，σy为y轴方向的电导率，μm为所述被测试件的lame常数，ρm为所述被测试件的密度；
根据
其中，vrz为所述表面波沿z轴传播的波速，σz为z轴方向的电导率，μm为所述被测试件的lame常数，ρm为所述被测试件的密度。
可选的，所述根据所述波速得到所述表面波位移，具体包括：
根据
其中，kry为所述表面波沿y轴方向传播的面内波数，ωry为所述表面波沿y轴正方向的振动角频率，vry为所述表面波沿y轴传播的波速；
根据
其中，u(y,t)为被测试件中点h处的表面波在时刻t沿y轴的面内位移,ωry为表面波沿y轴正方向的振动角频率,α和β是表面波沿z轴正方向的衰减系数，kt为横波波数，kry为所述表面波沿y轴正方向的波数，μm为所述被测试件的lame常数；
根据
其中，krz为所述表面波沿z轴方向传播的面外波数，ωrz为所述表面波沿z轴正方向的振动角频率，vrz为所述表面波沿z轴传播的波速；
根据
其中，u(z,t)为被测试件中点h处的表面波在时刻t沿z轴的面外位移,ωrz为表面波沿z轴正方向的振动角频率，α和β为表面波沿z轴正方向的衰减系数，kt为横波波数，krz是所述表面波沿y轴正方向的波数，μm为所述被测试件的lame常数。
可选的，所述表面波位移分布图包括表面波面内位移分布图和表面波面外位移分布图；
所述根据所述表面波位移得到表面波位移分布图，具体包括：
根据公式
根据各所述灰度值gy(m,n)得到表面波面内位移分布图gy；
其中，gy(m,n)为所述表面波面内位移灰度图gy中像素(m,n)处的灰度值；min{u(y,t)|y＝0,1,...,w；t＝0,1,...,t}为所述表面波面内位移的最小值，max{u(y,t)|y＝0,1,...,w；t＝0,1,...,t}为所述表面波面内位移的最大值，u(y,t)为被测试件中点h处的表面波在时刻t沿y轴的面内位移；
根据公式
根据各所述灰度值gz(m,n)得到表面波面外位移分布图gz；
其中，gz(m,n)为所述表面波面外位移灰度图gz中像素(m,n)处的灰度值；min{u(z,t)|z＝0,1,...,h；t＝0,1,...,t}为所述表面波面外位移的最小值，max{u(z,t)|z＝0,1,...,h；t＝0,1,...,t}为所述表面波面外位移的最大值，u(z,t)为被测试件中点h处的表面波在时刻t沿z轴的面外位移。
为实现上述目的，本发明还提供了如下方案：
一种各向异性金属薄板感应式磁声图像的仿真系统，所述仿真系统包括：
仿真结构模型建立模块，用于建立电导率各向异性金属薄板仿真结构模型，所述仿真结构模型中包括仿真被测试件和仿真折线线圈；
电导率设置模块，用于根据所述仿真结构模型，设置所述被测试件不同方向的电导率；
感应涡流获取模块，用于将所述仿真结构模型置于静磁场中进行仿真实验，得到感应涡流；
洛伦兹力确定模块，用于根据所述感应涡流，确定所述被测试件受到的静态洛伦兹力和交变洛伦兹力；
表面波获取模块，用于根据所述静态洛伦兹力和所述交变洛伦兹力得到表面波；
波速计算模块，用于根据各不同方向的电导率，计算所述表面波沿对应方向上的波速；
表面波位移获取模块，用于根据表面波沿不同方向上的波速，得到对应方向上的表面波位移；
分布图获取模块，用于根据各所述表面波位移，得到表面波位移分布图。
可选的，所述电导率设置模块，具体包括：
直角坐标系单元，用于基于所述被测试件建立三维笛卡尔直角坐标系，其中，以所述被测试件的表面中心为原点，与所述被测试件水平面平行的为x轴，与x轴垂直的为y轴，与所述xoy面垂直的为z轴；
电导率设置单元，用于设置被测试件在x轴、y轴和z轴三个方向上的电导率。
可选的，所述表面波位移获取模块，具体包括：
面内波数获取单元，用于根据
面内位移获取单元，用于根据
其中，u(y,t)为被测试件中点h处的表面波在时刻t沿y轴的面内位移，ωry为表面波沿y轴正方向的振动角频率，α和β是表面波沿z轴正方向的衰减系数，kt为横波波数，kry为所述表面波沿y轴正方向的波数，μm为所述被测试件的lame常数；
面外波数获取单元，用于根据
面外位移获取单元，用于根据
其中，u(z,t)为被测试件中点h处的表面波在时刻t沿z轴的面外位移，ωrz为表面波沿z轴正方向的振动角频率，α和β为表面波沿z轴正方向的衰减系数，kt为横波波数，krz是所述表面波沿y轴正方向的波数，μm为所述被测试件的lame常数。
可选的，所述表面波位移分布图获取模块，具体包括：
第一灰度值获取单元，用于根据公式
面内位移分布图获取单元，用于根据各所述灰度值gy(m,n)得到表面波面内位移分布图gy；
其中，gy(m,n)为所述表面波面内位移灰度图gy中像素(m,n)处的灰度值；min{u(y,t)|y＝0,1,...,w；t＝0,1,...,t}为所述表面波面内位移的最小值，max{u(y,t)|y＝0,1,...,w；t＝0,1,...,t}为所述表面波面内位移的最大值，u(y,t)为被测试件中点h处的表面波在时刻t沿y轴的面内位移；
第二灰度值获取单元，用于根据公式
面外位移分布图获取单元，用于根据各所述灰度值gz(m,n)得到表面波面外位移分布图gz；
其中，gz(m,n)为所述表面波面外位移灰度图gz中像素(m,n)处的灰度值；min{u(z,t)|z＝0,1,...,h；t＝0,1,...,t}为所述表面波面外位移的最小值，max{u(z,t)|z＝0,1,...,h；t＝0,1,...,t}为所述表面波面外位移的最大值，u(z,t)为被测试件中点h处的表面波在时刻t沿z轴的面外位移。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明公开一种各向异性金属薄板感应式磁声图像的仿真方法，建立电导率各向异性金属薄板仿真结构模型；根据模型设置仿真被测试件不同方向的电导率；将模型置于静磁场中进行仿真实验，得到感应涡流；根据感应涡流，确定被测试件受到的静态洛伦兹力和交变洛伦兹力；根据静态洛伦兹力和交变洛伦兹力得到表面波；根据各不同方向的电导率，计算表面波沿对应方向上的波速；根据表面波沿不同方向上的波速，得到对应方向上的表面波位移；根据各表面波位移，得到表面波位移分布图，本发明的方法能够通过改变各个方向上的电导率，精确地仿真电导率分布各向异性金属薄板表面的力分布和表面波位移。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例各向异性金属薄板感应式磁声图像的仿真方法流程图；
图2为本发明实施例仿真结构模型图；
图3为本发明实施例直角坐标系图；
图4为本发明实施例被测试件与直角坐标系的立体图；
图5为本发明实施例各向异性金属薄板感应式磁声图像的仿真系统结构图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种各向异性金属薄板感应式磁声图像的仿真方法及系统，通过改变各个方向上的电导率，精确地仿真电导率分布各向异性金属薄板表面的力分布和表面波位移。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明实施例各向异性金属薄板感应式磁声图像的仿真方法流程图。如图1所示，一种各向异性金属薄板感应式磁声图像的仿真方法，所述仿真方法包括：
步骤101：建立电导率各向异性金属薄板仿真结构模型，所述仿真结构模型中包括仿真被测试件和仿真折线线圈；
步骤102：根据所述仿真结构模型，设置所述被测试件不同方向的电导率；
步骤103：将所述仿真结构模型置于静磁场中进行仿真实验，得到感应涡流；
步骤104：根据所述感应涡流，确定所述被测试件受到的静态洛伦兹力和交变洛伦兹力；
步骤105：根据所述静态洛伦兹力和所述交变洛伦兹力得到表面波；
步骤106：根据各不同方向的电导率，计算所述表面波沿对应方向上的波速；
步骤107：根据表面波沿不同方向上的波速，得到对应方向上的表面波位移；
步骤108：根据各所述表面波位移，得到表面波位移分布图。
其中，步骤102具体包括：
基于所述被测试件建立三维笛卡尔直角坐标系，其中，以所述被测试件的表面中心为原点，与所述被测试件水平面平行的为x轴，与x轴垂直的为y轴，与所述xoy面垂直的为z轴；
设置被测试件在x轴、y轴和z轴三个方向上的电导率。
步骤103具体包括：
根据
其中，fs(h,t)为被测试件受到的静态洛伦兹力，σ为电导率，μ0为空间磁导率，j是虚部单位，bs为静磁场强度，z0为折线线圈与被测试件的垂直距离，k为径向积分变量，i为激励电流的幅值，ω为激励电流的角频率，j0(k)为自变量为k的第一类零阶bessel函数，j0(kr)为自变量为kr的第一类零阶bessel函数。
根据
其中，fd(h,t)被测试件受到的交变洛伦兹力，μ0为空间磁导率，i为激励电流的幅值，k为径向积分变量，z0为折线线圈与被测试件的垂直距离，j是虚部单位，j0(k)为自变量为k的第一类零阶bessel函数，j0(kr)为自变量为kr的第一类零阶bessel函数。
步骤106具体包括：
表面波的传播仅位于yoz平面，与x轴无关；
根据
其中，vry为所述表面波沿y轴传播的波速，σy为y轴方向的电导率，μm为所述被测试件的lame常数，ρm为所述被测试件的密度；
根据
其中，vrz为所述表面波沿z轴传播的波速，σz为z轴方向的电导率，μm为所述被测试件的lame常数，ρm为所述被测试件的密度。
步骤107具体包括：
根据
其中，kry为所述表面波沿y轴方向传播的面内波数，ωry为所述表面波沿y轴正方向的振动角频率，vry为所述表面波沿y轴传播的波速；
根据
其中，u(y,t)为被测试件中点h处的表面波在时刻t沿y轴的面内位移，ωry为表面波沿y轴正方向的振动角频率，α和β是表面波沿z轴正方向的衰减系数，kt为横波波数，kry为所述表面波沿y轴正方向的波数，μm为所述被测试件的lame常数；
根据
其中，krz为所述表面波沿z轴方向传播的面外波数，ωrz为所述表面波沿z轴正方向的振动角频率，vrz为所述表面波沿z轴传播的波速；
根据
其中，u(z,t)为被测试件中点h处的表面波在时刻t沿z轴的面外位移，ωrz为表面波沿z轴正方向的振动角频率，α和β为表面波沿z轴正方向的衰减系数，kt为横波波数，krz是所述表面波沿y轴正方向的波数，μm为所述被测试件的lame常数。
步骤108中，所述表面波位移分布图包括表面波面内位移分布图和表面波面外位移分布图；
根据公式
根据各所述灰度值gy(m,n)得到表面波面内位移分布图gy；
其中，gy(m,n)为所述表面波面内位移灰度图gy中像素(m,n)处的灰度值；min{u(y,t)|y＝0,1,...,w；t＝0,1,...,t}为所述表面波面内位移的最小值，max{u(y,t)|y＝0,1,...,w；t＝0,1,...,t}为所述表面波面内位移的最大值，u(y,t)为被测试件中点h处的表面波在时刻t沿y轴的面内位移；
根据公式
根据各所述灰度值gz(m,n)得到表面波面外位移分布图gz；
其中，gz(m,n)为所述表面波面外位移灰度图gz中像素(m,n)处的灰度值；min{u(z,t)|z＝0,1,...,h；t＝0,1,...,t}为所述表面波面外位移的最小值，max{u(z,t)|z＝0,1,...,h；t＝0,1,...,t}为所述表面波面外位移的最大值，u(z,t)为被测试件中点h处的表面波在时刻t沿z轴的面外位移。
本发明提供的仿真方法，通过建立电导率各向异性金属薄板仿真结构模型；根据模型设置仿真被测试件不同方向的电导率；将模型置于静磁场中进行仿真实验，得到感应涡流；根据感应涡流，确定被测试件受到的静态洛伦兹力和交变洛伦兹力；根据静态洛伦兹力和交变洛伦兹力得到表面波；根据各不同方向的电导率，计算表面波沿对应方向上的波速；根据表面波沿不同方向上的波速，得到对应方向上的表面波位移；根据各表面波位移，得到表面波位移分布图，能够通过改变各个方向上的电导率，精确地仿真电导率分布各向异性金属薄板表面的力分布和表面波位移。
图2为本发明实施例仿真结构模型图。如图2所示，被测试件1是长方体金属薄板，其厚度大于4倍的超声波长，将其置于由永磁体2(如钕铁硼磁铁)提供的静磁场中，同时将通入正弦交流电流的折线线圈3置于待测试件上方，试件相对于线圈足够大，以避免产生边缘效应，由电磁超声换能器4在被测试件周围接收磁声信号。
图3为本发明实施例三维笛卡尔直角坐标系图。如图3所示，对被测试件建立三维笛卡尔坐直角标系oxyz，坐标原点o位于被测试件表面中心。被测试件的电导率张量为：
其中d是旋转矩阵：
式中，γ是各向异性走向角，ψ是各向异性倾角，χ是各向异性偏角。
式中，σx、σy、σz分别是x、y、z三个正交方向上的电导率。
对被测试件施加沿z轴负方向的静磁场5，强度为bs。折线线圈3共有n折且相邻导线之间的间距为l，折线线圈3与被测试件1之间的垂直距离是z0。将折线线圈3简化为相互平行且相互不连接的导线，其中通入与z轴垂直的正弦交流电流作为激励电流：
jc(t)＝iejωt
其中，h＝(x,y,z)是三维直角坐标系oxyz中一点的位置矢量；t是时间；jc(t)是时间t的激励电流强度；i是激励电流的幅值；ω是激励电流的角频率；j是虚部单位。激励电流在被测试件周围产生交变磁场，其在时间t、位置h处的强度为bd(h,t)。
线圈下方的作用区域包括ra和rb两部分。ra中的点坐标满足0＜z＜z0，ra中的磁矢势主要由线圈中的激励电流和被测试件的表面涡流引起，其计算公式为
其中，aa(h,t)是ra内位置h处在时刻t的磁矢势；
rb中的点坐标满足z≤0，rb中的磁矢势由被测试件的表面涡流产生，其计算公式为
其中，ab(h,t)是rb内位置h处在时刻t的磁矢势。
单折线圈在待测试件表面h处在时刻t产生的感应涡流为：
则n折线圈在待测试件表面h处在时刻t产生的感应涡流为：
je(h,t)＝nj1e(h,t)
感应涡流je分别与静磁场和交变磁场相互作用产生静态洛伦兹力和交变洛伦兹力，即为声源。具体计算方法是：je(h,t)与bs相互作用在时刻t位置h处产生的静态洛伦兹力为：
je(h,t)与bd(h,t)在时刻t位置h处的轴向分量相互作用产生的交变洛伦兹力为：
被测试件1中的移动电荷在静态洛伦兹力和交变洛伦兹力的共同作用下产生周期性的振动，且以表面波的形式向外传播。表面波位移仅位于yoz平面中，与x无关。在时刻t表面波沿y和z两个方向传播的面内位移和面外位移分别是：
其中，u(y,t)是被测试件1中点h处的表面波在时刻t的面内位移；u(z,t)是被测试件中点h处的表面波在时刻t的面外位移；
vl是纵波波速：
vry和vrz分别是表面波沿y轴正方向和z轴正方向的波速：
ρm是被测试件的密度；λm和μm是被测试件的lame常数。
将被测试件1表面各点处的表面波面内位移和面外位移曲线，分别转化成灰阶为256的灰度图像gy和gz。具体方法如下：
被测试件1表面各点处的表面波面内位移为u(y,t)，其中y＝0,1,...,w，t＝0,1,...,t，图4为本发明实施例被测试件1与直角坐标系的立体图。如图4所示，w是被测试件沿y方向的宽度，t是测量时间界限。则表面波面内位移灰度图gy中像素(m,n)处的灰度值为
其中，gy(m,n)是表面波面内位移灰度图gy中像素(m,n)处的灰度值；min{u(y,t)|y＝0,1,...,w；t＝0,1,...,t}和max{u(y,t)|y＝0,1,...,w；t＝0,1,...,t}分别是表面波面内位移的最小值和最大值。
被测试件表面各点处的表面波面外位移为u(z,t)，其中z＝0,1,...,h，t＝0,1,...,t，图4为本发明实施例被测试件与直角坐标系的立体图。如图4所示，h是被测试件沿z方向的厚度，t是测量时间界限。则表面波面外位移灰度图gz中像素(m,n)处的灰度值为
其中，gz(m,n)是表面波面外位移灰度图gz中像素(m,n)处的灰度值；min{u(z,t)|z＝0,1,...,h；t＝0,1,...,t}和max{u(z,t)|z＝0,1,...,h；t＝0,1,...,t}分别是表面波面外位移的最小值和最大值。
图5为本发明实施例一种各向异性金属薄板感应式磁声图像的仿真系统结构图。如图5所示，一种各向异性金属薄板感应式磁声图像的仿真系统，所述仿真系统包括：
仿真结构模型建立模块201，用于建立电导率各向异性金属薄板仿真结构模型，所述仿真结构模型中包括仿真被测试件和仿真折线线圈；
电导率设置模块202，用于根据所述仿真结构模型，设置所述被测试件不同方向的电导率；
感应涡流获取模块203，用于将所述仿真结构模型置于静磁场中进行仿真实验，得到感应涡流；
洛伦兹力确定模块204，用于根据所述感应涡流，确定所述被测试件受到的静态洛伦兹力和交变洛伦兹力；
表面波获取模块205，用于根据所述静态洛伦兹力和所述交变洛伦兹力得到表面波；
波速计算模块206，用于根据各不同方向的电导率，计算所述表面波沿对应方向上的波速；
表面波位移获取模块207，用于根据表面波沿不同方向上的波速，得到对应方向上的表面波位移；
分布图获取模块208，用于根据各所述表面波位移，得到表面波位移分布图。
其中，所述电导率设置模块202，具体包括：
直角坐标系单元，用于基于所述被测试件建立三维笛卡尔直角坐标系，其中，以所述被测试件的表面中心为原点，与所述被测试件水平面平行的为x轴，与x轴垂直的为y轴，与所述xoy面垂直的为z轴；
电导率设置单元，用于设置被测试件在x轴、y轴和z轴三个方向上的电导率。
所述表面波位移获取模块207，具体包括：
面内波数获取单元，用于根据
面内位移获取单元，用于根据
其中，u(y,t)为被测试件中点h处的表面波在时刻t沿y轴的面内位移，ωry为表面波沿y轴正方向的振动角频率，α和β是表面波沿z轴正方向的衰减系数，kt为横波波数，kry为所述表面波沿y轴正方向的波数，μm为所述被测试件的lame常数；
面外波数获取单元，用于根据
面外位移获取单元，用于根据
其中，u(z,t)为被测试件中点h处的表面波在时刻t沿z轴的面外位移，ωrz为表面波沿z轴正方向的振动角频率，α和β为表面波沿z轴正方向的衰减系数，kt为横波波数，krz是所述表面波沿y轴正方向的波数，μm为所述被测试件的lame常数。
所述表面波位移分布图获取模块208，具体包括：
第一灰度值获取单元，用于根据公式
面内位移分布图获取单元，用于根据各所述灰度值gy(m,n)得到表面波面内位移分布图gy；
其中，gy(m,n)为所述表面波面内位移灰度图gy中像素(m,n)处的灰度值；min{u(y,t)|y＝0,1,...,w；t＝0,1,...,t}为所述表面波面内位移的最小值，max{u(y,t)|y＝0,1,...,w；t＝0,1,...,t}为所述表面波面内位移的最大值，u(y,t)为被测试件中点h处的表面波在时刻t沿y轴的面内位移；
第二灰度值获取单元，用于根据公式
面外位移分布图获取单元，用于根据各所述灰度值gz(m,n)得到表面波面外位移分布图gz；
其中，gz(m,n)为所述表面波面外位移灰度图gz中像素(m,n)处的灰度值；min{u(z,t)|z＝0,1,...,h；t＝0,1,...,t}为所述表面波面外位移的最小值，max{u(z,t)|z＝0,1,...,h；t＝0,1,...,t}为所述表面波面外位移的最大值，u(z,t)为被测试件中点h处的表面波在时刻t沿z轴的面外位移。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。