基于粒子群优化单稳态自适应随机共振的自主式水下机器人信号处理方法技术

本发明专利技术提供的是一种基于粒子群优化单稳态自适应随机共振的自主式水下机器人信号处理方法。通过粒子群优化算法对单稳态随机共振系统结构参数进行优化，以实现AUV控制信号和状态信号的自适应随机共振，并提高其随机共振效果，通过随机共振现象最终达到抑制AUV控制信号和状态信号中所含外部随机干扰，增强故障信号特征的目的。本发明专利技术解决传统单参数固定步长自适应随机共振方法优于忽略了参数之间的交互作用而导致的结构参数选取缺乏合理的理论依据，以及难以实现最优随机共振效果的问题，经自适应随机共振系统处理后的AUV控制信号和状态信号，抑制了外部干扰且增强了故障信号特征，可用于AUV推进器故障诊断、容错控制等领域。

【技术实现步骤摘要】
基于粒子群优化单稳态自适应随机共振的自主式水下机器人信号处理方法
本专利技术涉及的是一种自主式水下机器人推进器故障诊断、容错控制方法。
技术介绍
随着陆地资源日渐减少，人类开发海洋的步伐越来越快。自主式水下机器人(AUV：AutonomousUnderwaterVehicle)是目前唯一能够在无人情况下在深海进行探测、开发的载体，一直受到国内外研究人员的高度重视。推进器是AUV最重要的执行部件且负荷最重，一旦其出现故障直接影响AUV的安全性。对AUV推进器状态进行故障诊断，是提高AUV安全性的重要技术手段。由于外部干扰的影响，直接采用传感器采集到的AUV控制信号和状态信号进行故障诊断效果不佳，存在容易误检和漏检的问题。目前的外部干扰抑制方法能够有效抑制AUV信号中所含外部干扰，但在抑制的同时也极易造成过抑制问题，即对于故障诊断极为重要的信号特征也被抑制，影响故障诊断的准确性。随机共振理论是在研究古气象冰川问题时提出的，它是指在一个非线性双稳系统中，当仅在周期信号的驱动下，不足以使系统输出在两个稳态间进行跃迁，但在噪声的协助下，系统输出会按信号频率在两稳态间进行切换，从而达到抑制噪声并强化周期信号的效果。为实现信号的随机共振，传统方法通常采用向原信号中添加噪声或固定系统一个结构参数并按固定步长调节另一个结构参数的自适应调整策略。但在研究AUV信号随机共振的过程中发现，AUV信号中外部干扰与真实故障信号耦合在一起难以剥离，添加噪声的方法难以实现其随机共振；固定一个参数调整另一个参数的调整策略能够实现AUV信号随机共振，但由于只对单一参数进行自适应调整，忽略了参数之间的交互作用，存在随机共振系统参数选取缺乏合理的理论依据以及难以实现最优随机共振效果的问题。
技术实现思路
本专利技术的目的在于提供一种能抑制外部干扰、增强故障信号特征、提高故障诊断的准确性的基于粒子群优化单稳态自适应随机共振的自主式水下机器人信号处理方法。本专利技术的目的是这样实现的：(1)首先，对自主式水下机器人采集到的数据进行滑动窗处理，当采集到数据长度为L的传感器和控制器信号后开始检测，当再次采集到新的数据后，舍弃原数组第一个数据并将新采集回来的数据放在原数组的末尾，始终保持数据长度为L＝250～500；(2)采用粒子群优化算法对单稳态随机共振系统结构参数进行优化，首先进行种群初始化，具体方法为：设置单稳系统结构参数a、b的搜索范围在[0,7]区间，Tmax＝100，随机初始化N＝50个粒子组成一个种群；(3)采用式所示的单稳态随机共振模型处理自主式水下机器人信号，式中x(t)为单稳随机共振系统输出信号，为中x(t)的导数，s(t)为步骤(1)所述的自主式水下机器人控制信号和状态信号，a、b为单稳随机共振系统结构参数；(4)采用互相关系数R作为适应度函数来评价粒子的个体适应度，R值越大则粒子的适应性越好；R由式求得，式中其中：s为自主式水下机器人纵向速度原始信号，x为单稳系统输出信号，L为自主式水下机器人信号数据长度，j代表取L个数据中的第j个数据；(5)采用式Vi(t+1)＝w×Vi(t)+c1×R1×(Pi(t)-Xi(t))+c2×R2×(G(t)-Xi(t))所示的方式更新所有粒子的分形速度和位置，其中：t为当前迭代次数，i代表种群中第i个粒子，Xi(t)、Vi(t)分别为第t次迭代时第i个粒子的位置和飞行速度，Vi(t+1)为第t+1次迭代时第i个粒子的飞行速度，Vi(t)为第t次迭代时第i个粒子的飞行速度，c1＝2和c2＝2为学习因子，R1和R2为介于[0,1]之间相互独立均匀分布的随机数，Pi(t)为经过t次迭代第i个粒子迄今为止搜索到的最优解，G(t)为经过t次迭代整个种群所有粒子迄今为止搜索到的最优解，w为惯性权重、由式w(t)＝wmax-(wmax-wmin)×t/Tmax求得，式中wmax＝25为惯性权重的上限，wmin＝1为惯性权重的下限，t为当前迭代次数，Tmax为最大迭代次数，得到更新后的粒子飞行速度Vi(t+1)后，带入式Xi(t+1)＝Xi(t)+Vi(t+1)，式中Xi(t+1)为第t+1次迭代时第i个粒子的位置，Xi(t)为第t次迭代时第i个粒子的位置。经上述计算后，得到更新后粒子的位置和飞行速度；(6)判断是否符合终止条件，当满足下列条件中的任何一个时，迭代终止：1.进化代数达到设定的最大阈值Tmax；2.两代之间的最优适应度值之差小于预设精度，满足上述两个条件任何一个则迭代优化过程终止，进入下一步骤；两个条件均不满足则回到步骤(3)继续进行迭代优化；(7)经粒子迭代更新优化后，得到最优结构参数a和b的值并输出给单稳随机共振系统；(8)采用最优结构参数a和b的单稳态随机共振：根据步骤(7)得到的最优结构参数a和b的值，固化步骤(3)所述单稳态随机共振系统的结构参数值，对自主式水下机器人控制信号和状态信号进行随机共振处理。本专利技术根据AUV信号的非周期特性，采用单稳态随机共振系统实现其随机共振，采用粒子群优化算法实现随机共振系统结构参数的自适应优化调整，构成一种新型的基于粒子群优化单稳态自适应随机共振的自主式水下机器人信号处理方法，有效解决传统方法难以实现AUV信号最佳随机共振效果问题，以抑制外部干扰，增强故障信号特征，进而提高故障诊断的准确性。本专利技术与现有技术相比的有益效果主要体现在：与现有传统双稳态自适应随机共振方法相比，本专利技术提出了一种新型的单稳态自适应随机共振实现方法，该方法有效解决了传统方法仅对单一参数进行调节，缺乏理论依据的问题；解决了传统方法由于忽略了参数之间的交互作用，难以实现最优随机共振效果的问题。经随机共振处理后的AUV信号在抑制外部干扰的同时增强了故障信号特征，为故障诊断提供更为准确的AUV状态信息。附图说明图1为本专利技术方法的流程图。图2为推进器推力损失程度为20％的AUV多普勒信号时域曲线。图3为本专利技术方法的最优结构参数收敛曲线。图4为推力损失程度为20％的AUV多普勒信号经本专利技术方法处理后时域曲线。图5为传统单参数调节自适应随机共振方法实现的随机共振效果。具体实施方式下面结合附图举例对本专利技术做更详细的描述。结合图1，本专利技术的基于粒子群优化单稳态自适应随机共振的自主式水下机器人信号处理方法，其具体实施步骤如下：1、首先，对原始数据采用滑动窗进行截取，原始数据如图2所示。当采集到数据长度为L＝300的多普勒数据后启动检测算法，当再次采集到新的数据后，舍弃原数组第一个数据并将新采集回来的数据放在原数组的末尾，始终保持数据长度为L。2.采用粒子群优化算法对单稳态随机共振系统结构参数进行优化，首先进行种群初始化，具体方法为：设置单稳系统结构参数a、b的搜索范围在[0,7]区间，Tmax＝100，随机初始化N＝50个粒子组成一个种群。3.采用如下式所示的单稳态随机共振模型处理AUV信号：式中x(t)为单稳随机共振系统输出信号，为其导数，s(t)为步骤(1)所述的自主式水下机器人控制信号和状态信号，a、b为单稳随机共振系统结构参数；4.采用互相关系数R作为适应度函数来评价粒子的个体适应度，R值越大则粒子的适应性越好。互相关系数R由下式求得：式中其中：R为所求互相关系数，s为AUV纵向速度原始信号本文档来自技高网...

【技术保护点】
一种基于粒子群优化单稳态自适应随机共振的自主式水下机器人信号处理方法，其特征是：(1)首先，对自主式水下机器人采集到的数据进行滑动窗处理，当采集到数据长度为L＝250～500的传感器和控制器信号后开始检测，当再次采集到新的数据后，舍弃原数组第一个数据并将新采集回来的数据放在原数组的末尾，始终保持数据长度为L；(2)采用粒子群优化算法对单稳态随机共振系统结构参数进行优化，首先进行种群初始化，具体方法为：设置单稳系统结构参数a、b的搜索范围在[0,7]区间，Tmax＝100，随机初始化N＝50个粒子组成一个种群；(3)采用式所示的单稳态随机共振模型处理自主式水下机器人信号，式中x(t)为单稳随机共振系统输出信号，为中x(t)的导数，s(t)为步骤(1)所述的自主式水下机器人控制信号和状态信号，a、b为单稳随机共振系统结构参数；(4)采用互相关系数R作为适应度函数来评价粒子的个体适应度，R值越大则粒子的适应性越好；R由式R=Σj=1L(s(j)-s‾)(x(j)-x‾)Σj=1L(s(j)-s‾)2Σj=1L(x(j)-x‾)2]]>求得，式中s‾=1LΣj=1Ls(j),x‾=1LΣj=1Lx(j),]]>其中：s为自主式水下机器人纵向速度原始信号，x为单稳系统输出信号，L为自主式水下机器人信号数据长度，j代表取L个数据中的第j个数据；(5)采用式Vi(t+1)＝w×Vi(t)+c1×R1×(Pi(t)‑Xi(t))+c2×R2×(G(t)‑Xi(t))所示的方式更新所有粒子的分形速度和位置，其中：t为当前迭代次数，i代表种群中第i个粒子，Xi(t)、Vi(t)分别为第t次迭代时第i个粒子的位置和飞行速度，Vi(t+1)为第t+1次迭代时第i个粒子的飞行速度，Vi(t)为第t次迭代时第i个粒子的飞行速度，c1＝2和c2＝2为学习因子，R1和R2为介于[0,1]之间相互独立均匀分布的随机数，Pi(t)为经过t次迭代第i个粒子迄今为止搜索到的最优解，G(t)为经过t次迭代整个种群所有粒子迄今为止搜索到的最优解，w为惯性权重、由式w(t)＝wmax‑(wmax‑wmin)×t/Tmax求得，式中wmax＝25为惯性权重的上限，wmin＝1为惯性权重的下限，t为当前迭代次数，Tmax为最大迭代次数，得到更新后的粒子飞行速度Vi(t+1)后，带入式Xi(t+1)＝Xi(t)+Vi(t+1)，式中Xi(t+1)为第t+1次迭代时第i个粒子的位置，Xi(t)为第t次迭代时第i个粒子的位置。经上述计算后，得到更新后粒子的位置和飞行速度；(6)判断是否符合终止条件，当满足下列条件中的任何一个时，迭代终止：1.进化代数达到设定的最大阈值Tmax；2.两代之间的最优适应度值之差小于预设精度，满足上述两个条件任何一个则迭代优化过程终止，进入下一步骤；两个条件均不满足则回到步骤(3)继续进行迭代优化；(7)经粒子迭代更新优化后，得到最优结构参数a和b的值并输出给单稳随机共振系统；(8)采用最优结构参数a和b的单稳态随机共振：根据步骤(7)得到的最优结构参数a和b的值，固化步骤(3)所述单稳态随机共振系统的结构参数值，对自主式水下机器人控制信号和状态信号进行随机共振处理。...

【技术特征摘要】
1.一种基于粒子群优化单稳态自适应随机共振的自主式水下机器人信号处理方法，其特征是：(1)首先，对自主式水下机器人采集到的数据进行滑动窗处理，当采集到数据长度为L＝250～500的传感器和控制器信号后开始检测，当再次采集到新的数据后，舍弃原数组第一个数据并将新采集回来的数据放在原数组的末尾，始终保持数据长度为L；(2)采用粒子群优化算法对单稳态随机共振系统结构参数进行优化，首先进行种群初始化，具体方法为：设置单稳系统结构参数a、b的搜索范围在[0,7]区间，Tmax＝100，随机初始化N＝50个粒子组成一个种群；(3)采用式所示的单稳态随机共振模型处理自主式水下机器人信号，式中x(t)为单稳随机共振系统输出信号，为中x(t)的导数，s(t)为步骤(1)得到的长度为L的数据，a、b为单稳随机共振系统结构参数；(4)采用互相关系数R作为适应度函数来评价粒子的个体适应度，R值越大则粒子的适应性越好；R由式求得，式中其中：x为单稳系统输出信号，L为自主式水下机器人信号数据长度，s(j)代表L个数据中的第j个数据；(5)采用式Vi(t+1)＝w×Vi(t)+c1×R1×(Pi(t)-Xi(t))+c2×R2×(G(t)-Xi(t))所示的方式更新所有粒子的分形速度和位置，其中：t为当前迭代次数，i代表种群中第i个粒子，Xi(t)、Vi(t)分别为第t次迭代时第i个粒子的位置和飞行速度，Vi(t+1)为第t+1次迭代时第i个粒子的飞行速度，Vi(t)为第t次迭代时第i个粒子的...

【专利技术属性】
技术研发人员：张铭钧，刘维新，刘星，殷宝吉，王玉甲，赵文德，姚峰，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：黑龙江;23