[发明专利]一种低频振荡快速衰减信号分量参数的提取方法

权利要求书

1.一种Prony法低频振荡快速衰减信号分量参数的提取方法，其特征在于：包括以下步骤：

A、读入滤波后的测量数据

按照采样间隔Δt和采样时长T读入滤波后的测量数据，所述的采样时长T>10s，采样间隔Δt＝0.005s；

B、划分时段

把采样时长T分成两个时段，分段时刻为t_s，第一时段时长为t₁，t₁的最佳取值范围为0.5s≤t₁≤1.4s，剩余时段划为第二时段t₂，则第二时段时长t₂＝T-t₁；

C、在第二时段内分析慢速衰减信号分量参数

第二时段测量数据用来分析低频振荡中的慢速衰减信号分量参数，所述的信号分量参数为信号分量的频率、衰减因子、幅值和初相参数；采样间隔为Δt₂，Δt₂＝2k₂Δt,k₂＝5,6,…,14，从t_s时刻开始每间隔Δt₂从读入的测量数据中取一个值，形成第二时段的测量数据，采用扩展Prony法分析低频振荡中的慢速衰减信号分量参数；在第二时段内，以t_s时刻为扩展Prony法分析的0时刻，所取测量数据的个数N₂为

式中，T为原始测量数据的采样时长，t_s为分段时刻，Δt₂为第二时段测量数据的采样间隔；为向上取整符号；

D、调整慢速衰减信号分量参数

步骤C分析得到的慢速衰减信号分量参数是以分段时刻t_s为计时起点的信号分量参数，需要转换成以0时刻为计时起点的信号分量参数，各信号分量的频率和衰减因子不变，幅值和初相为

式中，A_i和θ_i分别是以t_s时刻为计时起点的第i个信号分量的幅值和初相，A_i0和θ_i0分别是以0时刻为计时起点的第i个信号分量的幅值和初相，f_i和α_i分别为第i个信号分量的频率和衰减因子；p为低频振荡中慢速衰减信号分量的个数；

根据式(2)得到的θ_i0需要调整为反余弦函数的主值，即使得θ_i0在[-π,π]范围内；

E、第一时段测量数据处理

步骤D得到以0时刻为计时起点的各信号分量的频率、衰减因子、幅值和初相参数后，从第一时段测量数据减去已经分析出的慢速衰减信号分量值，只保留快速衰减信号分量；

F、在第一时段内分析快速衰减信号分量参数

利用第一时段已处理过的测量数据分析低频振荡中的快速衰减信号分量参数，采样间隔为Δt₁，Δt₁＝k₁Δt,k₁＝1,2，从0s开始每间隔Δt₁从已处理过的测量数据中取一个值，形成第一时段的测量数据，采用扩展Prony法分析低频振荡中的快速衰减信号分量参数；第一时段所取测量数据的个数N₁为

G、分析结果输出

将第二时段得到的慢速衰减信号分量参数与第一时段得到的快速衰减信号分量参数合在一起，作为分段Prony法得到的所有信号分量参数。