CN110263762B - 一种基于输出的半潜式海洋平台能量传递路径分析方法 - Google Patents

Abstract

本发明公开了一种基于输出的半潜式海洋平台能量传递路径分析方法，包括如下步骤：在振源处和目标点处布置加速度传感器，组成阵列，利用基于状态空间模型的短时Prony方法对振源处信号进行特征成分分析，确定特征频率；离散振源处和目标点信号

，利用复指数分解技术，引入一阶矩阵微分方程，提取特征频率处信号；求解特征频率附近重构信号的相关系数，提取特征分量；计算特征分量的幅值，并代入公式求解传递率，通过在空间上绘制传递率，绘制传递率空间分布图，既可以得到整个空间范围内特征信号成分的传递率的变化规律，从而揭示特征信号传递情况。该方法仅依赖于测试得到的输出响应信号，避免了传统方法中激励力以及传递函数的测量。

Description

一种基于输出的半潜式海洋平台能量传递路径分析方法

技术领域

本发明涉及一种基于输出的半潜式海洋平台能量传递路径分析方法。

背景技术

目前，结构的能量传递路径分析(TPA)主要分为两个方向，一种是根据振源和目标元的时域、频域传递函数TPA方法，通过分解多输入源传递到目标元的分量，确定能量的主要传递路径；另一类是基于传递权重的TPA方法，用相邻结构之间的耦合特征表征路径的传递系数，依据子系统之间的连接关系搜索从振源到目标元的路径。理论上，传统TPA方法可以用于海洋工程结构振源能量传递路径分析，但是海洋工程结构尺度庞大，结构复杂，激励力、传递函数等输入通常难以测量。另外，特征成分的激振能量经过长距离的传播和耗散以后，达到目标元的能量通常极其微弱，导致常用的激振-响应方法在目标元处测量得到的响应没有足够的信噪比，振声响应等级很可能低于背景噪声。因此，针对振源能量在海洋工程结构的传递特点，研究复杂系统的振源能量传递路径评估方法，以及相应的控制技术，对海洋工程结构噪声预报和控制具有重要意义。

随着海洋资源开发利用的进一步深入，海洋工程结构对动力装置功率要求也在不断提高，由主动力装置和辅机产生的振动和噪声日益成为焦点问题。开展振声能量传递路径研究可以为海洋工程结构设计提供有益的参考，提高振声控制措施的针对性，节约建造过程中的振声控制成本。能量传递路径分析可以帮助工程师检测噪声问题的根本原因，为每个关键结构部位设定性能目标，并快速评估设计改进，广泛适用于各种工程领域，如汽车、航空航天、海洋工程装备、建筑声学等领域。2005年，Vecchio将传递路径分析技术应用于汽车结构振动声荷载的识别，测量了振源到目标元的传递函数，并估计了结构内部的舒适性。2009年，Guasch通过建立能量传递网络，利用线性算法和路径算法，找到最主要的传递路径，并再图形中修改损耗因子和耦合损耗因子从而达到控制振声能量的目的。2012年，Sievi将能量传递路径分析技术应用到铁路行业中，通过应用串扰消除技术，评估了参考信号之间的相关性，较精确的对传递函数进行了估计。2014年吴彬利用VA ONE软件对海洋平台典型舱室结构在分体空调室外机空气噪声作用下的舱室噪声进行预报，对海洋工程舱室噪声预报及设备布置方面提供一定参考价值。2017年，张文春以统计能量分析为基础，将噪声传递路径问题转化为求解赋权图的最大权重路径问题，并将其用在船舶噪声模型中，为船舶降噪措施具体实施方案提供明确参考。2018年，Ping Zhu提出了一种使用人工激励的TPA方法，通过在振源处施加随机荷载，并在目标元处施加确定荷载，通过建立响应和荷载之间的关系实现能量传递路径的分析。但是在传统的能量路径分析方法中，频率响应函数的测量无法避免。虽然该技术已经较为成熟，但是在测量频率响应函数时，必须移除系统的有源部分，这种测量极为耗时，并且测量部分的边界条件也不再准确。同时，在测量力的过程中，需要较大的空间来安装复杂的装置，并且还需要专门的支撑臂来安装测量传感器，在应用到海洋工程中时，极为不便。

发明内容

本发明的目的是提供一种基于输出的半潜式海洋平台能量传递路径分析方法，该方法仅依赖于测试得到的输出响应信号，避免了传统方法中激励力以及传递函数的测量。

本发明是一种对海洋结构物振动能量传递评估的新思路，该方法基于阵列的加速度传感器测得的输出响应信号，利用STPT方法对振源点振动信号进行特征成分分析，确定特征频率。利用复指数形式的极值、留数分解，提取振源处和目标点处特征成分信号的复指数表达式，通过计算相关系数并重构特征成分信号的时间序列，去除伪分量，之后量化提取误差，并进行重复提取，最终计算振源处到目标点处的传递率，从而揭示特征信号成分振动能量在空间中的传递情况。本发明通过基于复指数分解的分析技术而非传统的能量传递路径分析方法，克服了传统方法难以分析微弱能量特征成分的不足，并避免了傅里叶变换导致的周期性假设与实际不符、固定频率导致的漏频等问题，最终重建得到海洋结构振源特征成分信号的真实振动能量传递路径。

为解决上述技术问题，本发明采用如下具体技术方案：

一种基于输出的半潜式海洋平台能量传递路径分析方法，包括如下步骤：

1)在振源处和目标点处布置加速度传感器，组成阵列，并利用基于状态空间模型的短时Prony方法对振源处信号进行特征成分分析，确定特征频率；

2)离散振源处和目标点信号x_i(t)，利用复指数分解技术，引入一阶矩阵微分方程，提取特征频率处信号；

3)求解特征频率附近重构信号的相关系数，提取特征分量；

4)计算特征分量的幅值，并代入公式求解传递率，通过在空间上绘制传递率，绘制传递率空间分布图，既可以得到整个空间范围内特征信号成分的传递率的变化规律，从而揭示特征信号传递情况。

进一步的，所述步骤1)具体包括：

S1、在时间轴上设置一个长度为w的窗函数，并设置滑动距离h，且h≤w，通过滑动窗口可以把离散信号y(t)分解成若干段,每段数据长度都为w，且认为每段数据内的频率不随时间变换，则第q段信号表示为：

y_q(t)＝y_q(aΔt)(1)

式中，a＝(q-1)d，(q-1)d+1，...，(q-1)d+w-1，w为窗口长度，d为窗口每次滑动的距离，

N为原数据总长度，Δt表示采样信号的离散时间间隔；

S2、将每段数据y_q(t)分解成复指数形式，并求解其频率、阻尼等特征参数，之后画出时间-频率-幅值的三维图形，即可确定信号的时频特性，从而确定振源处的特征成分信号。

进一步的，所述步骤2)具体包括：

S3、为了计算特征信号能量有效值，必须获得特征信号

的表达式，将实测加速度信号分解为复指数形式，即

对于等间隔离散信号，式(4)可以进一步表示为：

其中，k＝1，2，3，…，N。

S4、引入

构造矩阵

式中，ξ和η分别为矩阵C的行数和列数。并分别对H_m(0)和H_m(1)进行奇异值分解，可以得到系统状态矩阵的一个实现为：

矩阵A_m的特征根为矩阵A_m的特征根为

其中l＝1、2、3n，由公式/>

可计算复指数序列中的复指数/>

进而由/>

可求得/>

对应的频率f_m，i，其中Im表示虚数部分。

通过最小二乘法求解

通过求得的

和/>

即可重构对应于特征频率的信号，

其中x_m，l表示m点处，频率在特征频率附近的重构信号，其中l＝1,2,3n。

进一步的，所述步骤3)具体包括：

S5、经过分解，我们可以得到特征频率附近的若干个分量。通过进行相关性分析进行筛选。计算所有特征频率附近的重构信号x_m,c和测量信号x_m之间的相关系数

式中,cov(x_m,c,x_m)＝E(x_m,cx_m)-E(x_m,c)E(x_m)(9)

其中，

和/>

代表x_m,c和x_m的方差其中c＝1,2,3…，cov(x_m,c,x_m)代表x_m,c和x_m的协方差，E代表期望。

通过上面的计算，可以得到所有分量的相关系数，一般情况下，相关系数

数值越大，说明该分量与原信号的相关性越大，通过比较相关系数，并重构该特征成份信号，若该重构的分量在时域上不衰减且相关系数较大，则认为该信号为目标点m处的特征成分信号x_m，c。

更进一步的，还包括：

S6、通过一次提取并不能保证特征分量被完全提取，为了量化提取误差，保证提取准确性，在L2-范数基础上，引入变量

式中，B_m,c代表重构信号的幅值，

s代表迭代次数。

进一步的，所述步骤4)具体包括：

S7、从振源m处到目标点n处的传递率T_mn等于两者响应的幅值的比值，在传统的方法中，幅值是通过傅里叶变换得到的，但是傅里叶变换存在频率泄露的问题，可能导致求得的幅值不准确，不同于基于傅里叶变换求解传递率的方法，幅值由公式(13)求得：

S8、重复利用步骤2)和3)，可以振源m处和目标点n处的特征分量的极值和留数，然后利用公式(13)，即可获得振源点和所有目标点处提取出的特征成份信号幅值；

S9、通过传递率的定义，从振源点到目标点的传递率可以表示为：

B_m,c＝T_mnB_n,c (14)

即

写成矩阵形式为

通过在空间上绘制传递率，既可以得到整个空间范围内特征信号成分的传递率的变化规律，从而揭示特征信号传递情况。

本发明提供了一种对半潜式海洋平台振源能量传递路径评估的新思路——基于阵列的加速度传感器测试得到的输出响应信号，利用STPT方法进行特征成分分析，确定特征成份信号频率，引入一阶状态空间模型，提取特征成份信号的复指数表达式，通过计算相关系数并重构特征信号成分的时间序列，去除伪分量，之后量化提取误差，重复提取特征成分信号，最终计算振源处到目标点处的传递率，从而揭示特征信号成分振动能量在空间中的传递情况；在工程上，为包括半潜式海洋平台在内的浮式海洋结构物的振源能量传递衰减评估提供了一种新的分析方法，能够为相关结构的噪声预报和控制等工作提供新的技术手段，具有一定的工程应用前景。

与现有技术相比，本发明的有益技术效果：

1)本发明基于阵列的加速度传感器测得的输出响应信号，具有测量速度快，空间分辨率好等优势，避免了传统方法中激励力及频率响应函数等参数的测量。通过将加速度传感器按照一定排列规则布置在目标点处，具有较强的指向性，可以用来重建真实的能量分布，从而揭示真实的能量振动传递细节信息。

2)本发明中，利用STPT方法对振源点信号进行特征成分分析，克服传统STFT及小波变换等时频分辨率比较低的缺点，并且避免希尔伯特-黄变换在对信号EMD分解时存在模态混叠的现象，该方法在分析存在大量噪音的大型复杂海上结构时振动信号时，其时频分析结果相较于传统方法更加清晰，工程应用前景较好。

3)本发明通过复指数分解技术提取特征信号成分，并重构特征信号成分的时间序列，同时求得能量幅值，进而求解振源点到目标点的传递率，提供一种新的海洋结构振源能量传递路径评估方法，增加了工程上进行能量传递路径评估的方法。

4)由于海洋工程结构尺度庞大，结构复杂，特征信号成分的激振能量经过长距离的传播和耗散以后，达到目标点的能量通常极其微弱，甚至可能低于背景噪声，难以识别。通过复指数分解技术，针对振源能量在海洋工程结构的传递特点，实现对复杂系统的振源能量传递路径的评估。

附图说明

图1为半潜式海洋平台DP舱平面示意图及传感器安装示意图；

图2为推进器处不同转速下测试数据时域图，其中，(a)为在推进器转速为620rpm时的时域图，(b)为在推进器转速为680rpm时的时域图；

图3为推进器在680rmp转速下振动信号时频图；

图4为1号—9号传感器146Hz特征信号成分时域和频域图，其中，a对应图1中上边缘处146Hz特征信号的时域图，b对应图1中左边缘处146Hz特征信号的时域图，c对应图1中下边缘处146Hz特征信号的时域图，d对应图1中右边缘处146Hz特征信号的时域图；

图5为本发明方法得到的振动能量传递路径空间分布结果；

附图标记说明：1-传感器。

具体实施方式

下面结合附图说明对本发明作进一步说明。

一种基于输出的半潜式海洋平台能量传递路径分析方法，包括如下步骤：

1)在振源处和目标点处布置加速度传感器，组成阵列，并利用基于状态空间模型的短时Prony方法对振源处信号进行特征成分分析，确定特征频率；

具体包括：

S1、在时间轴上设置一个长度为w的窗函数，并设置滑动距离h，且h≤w，通过滑动窗口可以把离散信号y(t)分解成若干段,每段数据长度都为w，且认为每段数据内的频率不随时间变换，则第q段信号表示为：

式中，a＝(q-1)d，(q-1)d+1，...，(q-1)d+w-1，w为窗口长度，d为窗口每次滑动的距离，N为原数据总长度，Δt表示采样信号的离散时间间隔；

S2、将每段数据y_q(t)分解成复指数形式，并求解其频率、阻尼等特征参数，之后画出时间-频率-幅值的三维图形，即可确定信号的时频特性，从而确定振源处的特征成分信号。

2)离散振源处和目标点信号x_i(t)，利用复指数分解技术，引入一阶矩阵微分方程，提取特征频率处信号；

具体包括：

S3、为了计算特征信号能量有效值，必须获得特征信号

的表达式，将实测加速度信号分解为复指数形式，即

对于等间隔离散信号，式(4)可以进一步表示为：

其中，k＝1，2，3，…，N。

S4、引入

构造矩阵

式中，ξ和η分别为矩阵C的行数和列数。并分别对H_m(0)和H_m(1)进行奇异值分解，可以得到系统状态矩阵的一个实现为：

矩阵A_m的特征根为矩阵A_m的特征根为

其中l＝1、2、3n。由公式/>

可计算复指数序列中的复指数/>

进而由/>

可求得/>

对应的频率f_m，i，其中Im表示虚数部分。

通过最小二乘法求解

通过求得的

和/>

即可重构对应于特征频率的信号，

其中x_m，l表示m点处，频率在特征频率附近的重构信号，其中l＝1,2,3n。

3)求解特征频率附近重构信号的相关系数，提取特征分量；

具体包括：

S5、经过分解，我们可以得到特征频率附近的若干个分量。通过进行相关性分析进行筛选。计算所有特征频率附近的重构信号x_m,c和测量信号x_m之间的相关系数

式中,cov(x_m,c,x_m)＝E(x_m,cx_m)-E(x_m,c)E(x_m)(9)

其中，

和/>

代表x_m,c和x_m的方差其中c＝1,2,3…，cov(x_m,c,x_m)代表x_m,c和x_m的协方差，E代表期望；

通过上面的计算，可以得到所有分量的相关系数，一般情况下，相关系数

数值越大，说明该分量与原信号的相关性越大，通过比较相关系数，并重构该特征成份信号，若该重构的分量在时域上不衰减且相关系数较大，则认为该信号为目标点m处的特征成分信号x_m，c；

还包括：

S6、通过一次提取并不能保证特征分量被完全提取，为了量化提取误差，保证提取准确性，在L2-范数基础上，引入变量

式中，B_m,c代表重构信号的幅值，

s代表迭代次数。

4)计算特征分量的幅值，并代入公式求解传递率，通过在空间上绘制传递率，绘制传递率空间分布图，既可以得到整个空间范围内特征信号成分的传递率的变化规律，从而揭示特征信号传递情况；

具体包括：

S7、从振源m处到目标点n处的传递率T_mn等于两者响应的幅值的比值，在传统的方法中，幅值是通过傅里叶变换得到的，但是傅里叶变换存在频率泄露的问题，可能导致求得的幅值不准确，不同于基于傅里叶变换求解传递率的方法，幅值由公式(13)求得：

S8、重复利用步骤2)和3)，可以振源m处和目标点n处的特征分量的极值和留数，然后利用公式(13)，即可获得振源点和所有目标点处提取出的特征成份信号幅值；

S9、通过传递率的定义，从振源点到目标点的传递率可以表示为：

B_m,c＝T_mnB_n,c (14)

即

写成矩阵形式为

通过在空间上绘制传递率，既可以得到整个空间范围内特征信号成分的传递率的变化规律，从而揭示特征信号传递情况。

其中：

步骤1)的创新性体现在：传统的传递路径分析方法中，需要测量频率响应函数，测试时需要移除系统的有源部分，这种测量极为耗时，并且测量部分的边界条件也不再准确。同时，在测量力的过程中，需要较大的空间来安装复杂的装置，并且还需要专门的支撑臂来安装测量传感器。加速度传感器较其他类型的传感器而言，具有磁性基底，可以实现快速的安装和拆除，更适用于海洋工程结构的振动测试。通过初步主观评估，将加速度传感器按照一定排列规则组成阵列，这种阵列具有较强的指向性，可以较好的测定能量的空间分布。利用STPT方法对振源处信号进行特征成分分析，克服传统STFT及小波变换不能同时获得高的时间和频率分辨率的缺点，同时也不存希尔伯特-黄变换存在模态混叠的现象，时频特性更加准确且分辨率较高。

步骤2)的创新性体现在：在传统的振动信号分析方法中，由于傅里叶变换收到频率分辨率的影响，只能处理谐波信号，本方法通过将信号表示为复指数的形式，意味着特征信号成分可以是谐波或者是非谐波的形式，用其来处理海洋工程结构物实测信号，可以增加对工程实测信号分析的准确度。通过将高阶微分方程转化为一阶微分方程，克服了传统Prony方法对信号进行分解时需要求解多项式方程的根的病态问题，使得求得的结果更加精确稳定。同时在每次求解过程中，仅提取出我们所关心的特征信号成分，理论上每次只提取一个特征信号成分，因此在计算过程中，本方法具有更高的效率。

步骤3)中创新性体现在：通过计算相关系数，可以快速缩小特征成分范围，极大的减少人为筛选过程的工作量。同时，通过在时域进行特征成分信号的重构，避免了伪分量成份的干扰。其次，通过量化提取误差，进行多次提取，可以对特征成分信号进行充分提取，提高了计算精度，使得计算结果更加准确。

步骤4)中创新性体现在：通过极值、留数求解振源点和目标点幅值，避免了传统傅里叶变换在计算特征信号成分幅值时的能量泄露问题，使得计算结果更为准确，进而利用幅值来计算传递率。最后，通过整合空间中各个路径的传递率，形成传递率矩阵，并用其来表示特征成分信号在海洋工程结构中的振动传递情况，可以增加海洋工程结构中，能量传递路径评估的准确性。

具体实施案例：

本实施例采用某中深水半潜式钻井平台推进器舱振动测试数据进行计算分析，传感器1的具体布置位置如图1所示，1个传感器位于推进器基座处，4个传感器分别在推进器舱的4个舱壁上。推进器转速设为620rpm，等推进器运转平稳以后，开始测量振动加速度信号，测试时间48s，然后将推进器转速提高到680rmp，推进器运转平稳后，测量振动加速度信号，测量时间为30s，采样频率均设置为2560Hz。图2给出测量信号在不同转速下的时域图，其中(a)为在推进器转速为620rpm时的时域图，图(b)为转速为680rmp时的时域图。

结果对比：

首先用基于状态空间模型的短时Prony方法对推进器处的加速度信号进行时频分析，时频分析时的窗口长度设为2000，阶次设为200，窗口移动距离设为500，图3给出了推进器转速为680rmp时，推进器振动信号的时频图，从图中可以看出振动信号在146Hz和292Hz能量比较稳定。因此在本实施例中，对146Hz和292Hz特征信号进行分析，研究其能量传递路径。图4给出了经过特征信号提取后，四个舱壁传感器处146Hz特征信号的时域图，并利用极值、留数求解振源处到目标点处的传递率。同理，可以得到推进器转速为620rmp的相关数据，下表1给出了在不同转速下，不同特征频率的传递率计算结果。图5即为本发明方法所得到的推进器振动能量向四个舱壁的传递结果，从图中可以看出振动能量的主要传递方向和变化规律：较大的能量向结构的中心部分传递，同时由于结构右下部份对应平台的外部结构，主要用于抵抗波浪等外部荷载的作用，因此设计强度较大，振动传递相较于其他部分较弱，从结果可以看出，结果与平台的结构设计是相一致的。

表1在不同转速下，不同特征频率的传递率计算结果

	T21	T31	T41	T51
					620rmp-135Hz	0.3846	0.3670	0.1591	0.0812
620rmp-270Hz	0.4630	0.3722	0.1279	0.0236
					620rmp-146Hz	0.3444	0.2931	0.1182	0.0392
620rmp-292Hz	0.4404	0.3281	0.1156	0.0781

以上所述的实施例仅是对本发明的优选方式进行描述，并非对本发明的范围进行限定，在不脱离本发明设计精神的前提下，本领域普通技术人员对本发明的技术方案做出的各种变形和改进，均应落入本发明权利要求书确定的保护范围内。

Claims (5)

1.一种基于输出的半潜式海洋平台能量传递路径分析方法，其特征在于，包括如下步骤：

1)在振源处和目标点处布置加速度传感器，组成阵列，并利用基于状态空间模型的短时Prony方法对振源处信号进行特征成分分析，确定特征频率；

2)离散振源处和目标点信号x_i(t)，利用复指数分解技术，引入一阶矩阵微分方程，提取特征频率处信号；

3)求解特征频率附近重构信号的相关系数，提取特征分量；

4)计算特征分量的幅值，并代入公式求解传递率，通过在空间上绘制传递率，绘制传递率空间分布图，既可以得到整个空间范围内特征信号成分的传递率的变化规律，从而揭示特征信号传递情况；

所述步骤4)具体包括：

S7、从振源m处到目标点n处的传递率T_mn等于两者响应的幅值的比值，在传统的方法中，幅值是通过傅里叶变换得到的，但是傅里叶变换存在频率泄露的问题，可能导致求得的幅值不准确，不同于基于傅里叶变换求解传递率的方法，幅值由公式(13)求得：

S8、重复利用步骤2)和3)，可以振源m处和目标点n处的特征分量的极值和留数，然后利用公式(13)，即可获得振源点和所有目标点处提取出的特征成份信号幅值；

S9、通过传递率的定义，从振源点到目标点的传递率可以表示为：

B_m,c＝T_mnB_n,c (14)

即

写成矩阵形式为

通过在空间上绘制传递率，既可以得到整个空间范围内特征信号成分的传递率的变化规律，从而揭示特征信号传递情况。

2.根据权利要求1所述的基于输出的半潜式海洋平台能量传递路径分析方法，其特征在于，所述步骤1)具体包括：

S1、在时间轴上设置一个长度为w的窗函数，并设置滑动距离h，且h≤w，通过滑动窗口可以把离散信号y(t)分解成若干段，每段数据长度都为w，且认为每段数据内的频率不随时间变换，则第q段信号表示为：

y_q(t)＝y_q(aΔt) (1)

式中，a＝(q-1)d，(q-1)d+1，...，(q-1)d+w-1，w为窗口长度，d为窗口每次滑动的距离，N为原数据总长度，Δt表示采样信号的离散时间间隔；

S2、将每段数据y_q(t)分解成复指数形式，并求解其频率、阻尼等特征参数，之后画出时间-频率-幅值的三维图形，即可确定信号的时频特性，从而确定振源处的特征成分信号。

3.根据权利要求1所述的基于输出的半潜式海洋平台能量传递路径分析方法，其特征在于，所述步骤2)具体包括：

S3、为了计算特征信号能量有效值，必须获得特征信号

的表达式，将实测加速度信号分解为复指数形式，即

对于等间隔离散信号，式(4)可以进一步表示为：

其中，k＝1，2，3，…，N；

S4、引入

构造矩阵

式中，ξ和η分别为矩阵C的行数和列数，并分别对H_m(0)和H_m(1)进行奇异值分解，可以得到系统状态矩阵的一个实现为：

矩阵A_m的特征根为矩阵A_m的特征根为

其中l＝1、2、3n，由公式/>

可计算复指数序列中的复指数/>

进而由/>

可求得/>

对应的频率f_m，i，其中Im表示虚数部分，通过最小二乘法求解/>

通过求得的

和/>

即可重构对应于特征频率的信号

其中x_m，l表示m点处，频率在特征频率附近的重构信号，其中l＝1,2,3n。

4.根据权利要求1所述的基于输出的半潜式海洋平台能量传递路径分析方法，其特征在于，所述步骤3)具体包括：

S5、经过分解，我们可以得到特征频率附近的若干个分量，通过进行相关性分析进行筛选，计算所有特征频率附近的重构信号x_m，c和测量信号x_m之间的相关系数

式中,cov(x_m,c,x_m)＝E(x_m,cx_m)-E(x_m,c)E(x_m)(9)

其中，

和/>

代表x_m,c和x_m的方差，其中c＝1,2,3…，cov(x_m,c,x_m)代表x_m,c和x_m的协方差，E代表期望；通过上面的计算，可以得到所有分量的相关系数，一般情况下，相关系数

数值越大，说明该分量与原信号的相关性越大，通过比较相关系数，并重构该特征成份信号，若该重构的分量在时域上不衰减且相关系数较大，则认为该信号为目标点m处的特征成分信号x_m，c。

5.根据权利要求4所述的基于输出的半潜式海洋平台能量传递路径分析方法，其特征在于，还包括：

S6、通过一次提取并不能保证特征分量被完全提取，为了量化提取误差，保证提取准确性，在L2-范数基础上，引入变量

式中，B_m,c代表重构信号的幅值，

s代表迭代次数。

Images