一种用于电力系统多输入多输出线性化模型辨识的低频段多正弦信号设计方法技术方案

本发明专利技术公开了一种用于电力系统多输入多输出线性化模型辨识的低频段多正弦信号设计方法，通过指定多正弦信号的ak、ωk及采样长度N，运用算法求解信号的一组参数

A Low Frequency Band Multi-Sinusoidal Signal Design Method for Multi-Input and Multi-Output Linearized Model Identification of Power System

The invention discloses a low-frequency multi-sinusoidal signal design method for multi-input and multi-output linearized model identification of power system. By specifying ak, _k and sampling length N of multi-sinusoidal signals, a set of parameters of signals are solved by an algorithm.

【技术实现步骤摘要】
一种用于电力系统多输入多输出线性化模型辨识的低频段多正弦信号设计方法
本专利技术涉及电力系统辨识和信号处理的交叉
，具体涉及一种用于电力系统多输入多输出线性化模型辨识的低频段多正弦信号设计方法。
技术介绍
随着互联电网的发展，电力系统在0.1至2.5赫兹频段范围内的低频振荡问题突出，需要通过安装电力系统低频振荡控制器加以解决。电力系统低频振荡控制器设计依赖电力系统多输入多输出线性化模型，在实际电力系统工程设计中需要首先解决多输入多输出线性化模型辨识问题。为了保障辨识过程中电力系统能安全稳定运行，实际工程通常采用小幅扰动信号对电力系统进行激励，然后采集电力系统激励输入与响应输出信号，通过系统辨识算法进行相应的线性化模型辨识。此时，作为辨识工作的源头，采用的小幅扰动信号就成为决定电力系统多输入多输出线性化模型辨识成功与否的关键因素。目前，实际工程采用的小幅扰动信号有经过低通滤波器的白噪声信号和有限频带的伪随机信号。两类信号在低频振荡关心的0.1至2.5赫兹频段内能量较大，但是在该关心频段以外存在一定的能量，导致对电力系统进行多输入多输出线性化模型辨识时，输入信号能量不够集中，系统输出响应信号的信噪比降低，从而影响了电力系统多输入多输出线性化模型辨识的精度。
技术实现思路
本专利技术的目的是为了解决现有技术中的上述缺陷，提供一种用于电力系统多输入多输出线性化模型辨识的低频段多正弦信号设计方法。本专利技术的目的可以通过采取如下技术方案达到：一种用于电力系统多输入多输出线性化模型辨识的低频段多正弦信号设计方法，所述的的低频段多正弦信号设计方法包括以下步骤：初始启动阶段：S1、确定低频段多正弦信号的时域长度为N，预期的波峰因子为Cr,set，对应的傅里叶系数记录变量为FS，其幅值记录变量为Mag、相位记录变量为Pha，初始FS、Mag、Pha均为长度为N的零向量；S2、确定低频段多正弦信号的谐波分量频率ωk以及频域幅值ak，其中k＝0,1,…,Nk-1，Nk≤N/2，然后将幅值ak按照谐波分量频率ωk对应于幅值记录变量Mag的位置填入幅值记录变量Mag，从而得到初始低频段多正弦信号的幅频特性；S3、设定低频段多正弦信号的谐波分量初始相位为一组范围在-π到π之间均匀分布的随机数，其中k＝0,1,…,Nk-1，然后将相位按照谐波分量频率ωk对应于相位记录变量Pha的位置填入相位记录变量Pha，从而得到初始低频段多正弦信号的相位特性；S4、将幅值记录变量Mag与相位记录变量Pha进行复数合成，得到初始低频段多正弦信号的频域特性FS，其中关心频段ωk对应的傅里叶系数为S5、对初始低频段多正弦信号的频域特性FS进行傅立叶逆变换，得到低频段多正弦信号的初始时域波形u(t)，并将初始时域波形u(t)存储在变量Signal中，变量Signal的长度为N；S6、计算变量Signal的波峰因子Cr，并将其存储在变量CF中，若变量CF小于预期的波峰因子Cr,set，计Signal_min＝Signal，然后转入步骤S13；否则，设定迭代循环次数Num为0，设置最大迭代次数为Num_max，Num_max>1，进入如下迭代循环阶段；迭代循环阶段：S7、若迭代循环次数Num小于最大迭代次数Num_max，则将迭代循环次数Num加1，然后转入步骤S8，否则，直接进入步骤S13；S8、将时域波形中绝对值超过最大值max(abs(u(t)))的90％的低频段多正弦信号波形值设定为最大值max(abs(u(t)))的90％，并保持符号不变，得到更新后的时域波形u(t)*，并对变量Signal进行替换更新，其中abs()表示取绝对值操作，max()表示取最大值操作；S9、对更新后的变量Signal进行傅立叶变换，得到傅立叶系数FS*，其相位特性为Pha*，用幅值记录变量Mag和相位特性Pha*合成作为更新后的频域特性FS；S10、对更新后的频域特性FS进行傅立叶逆变换，得到更新后的时域波形u(t)，并对变量Signal进行替换更新；S11、计算更新的变量Signal的波峰因子Cr，若波峰因子Cr小于预期的波峰因子Cr,set，计Signal_min＝Signal，然后进入步骤S13；否则进入步骤S12；S12、若波峰因子Cr小于上一次迭代过程的变量CF，则计变量CF＝Cr，Signal_min＝Signal，返回步骤S7；否则直接返回步骤S7；迭代循环退出；S13、将变量Signal_min作为最终得到的时域波形。进一步地，所述的低频段多正弦信号的时域波形表达式为其中t是采样时间，Ak、ωk和分别为第k个正弦谐波分量的时域幅值、频率和相位，Nk为正弦谐波频率分量的个数；所述的变量Signal的波峰因子按照波峰因子Cr的定义式计算，其中，波峰因子Cr的定义式如下：其中N是低频段多正弦信号的总采样次数(即信号的总长度)，max()代表取括号内数组的最大值。显然，Cr的大小反映了信号在时域的波动情况：给定信号u(t)的频域特性，信号u(t)的时域波动越小时，则其Cr越小。本专利技术相对于现有技术具有如下的优点及效果：本专利技术提出的一种用于电力系统多输入多输出线性化模型辨识的低频段多正弦信号设计方法，通过指定信号的ak、ωk及采样长度N，运用算法求解信号的一组参数其中k＝0,1,…,Nk-1，Nk≤N/2，使得u(t)的最大值尽可能小，从而使输入信号既能满足时域幅值限制要求，使其对系统的扰动较小，又能将输入信号的能量集中在关心的频段，而不关心频段的能量几乎为零。与传统的小幅激励信号设计方法相比，该方法解决了输入信号能量不集中，输出响应信号的信噪比低的问题，该设计好的低频段多正弦信号相比其他的低频段激励信号更适用于电力系统多输入多输出线性化模型辨识。附图说明图1是本专利技术的一个实施例中三种低频段输入信号的时域波形；图2是本专利技术的一个实施例中三种低频段输入信号的频域波形；图3是本专利技术的一个实施例所依据的电力系统接线图；图4是本专利技术的一个实施例中三种输入信号对应的输出信号的频域波形；图5是本专利技术公开的一种用于电力系统多输入多输出线性化模型辨识的低频段多正弦信号设计方法的流程步骤图。具体实施方式为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。实施例本专利技术提出的一种用于电力系统多输入多输出线性化模型辨识的低频段多正弦信号设计方法，如图5所示，包括以下步骤：初始启动(1)确定低频段多正弦信号的时域长度为N，预期的波峰因子为Cr,set，对应的傅里叶系数记录变量为FS，其幅值记录变量为Mag、相位记录变量为Pha，初始FS、Mag、Pha均为长度为N的零向量。(2)确定低频段多正弦信号的谐波分量频率ωk以及频域幅值ak，其中k＝0,1,…,Nk-1，Nk≤N/2，然后将幅值ak按照谐波分量频率ωk对应于幅值记录变量Mag的位置填入幅值记录变量Mag，从而得到初始低频段多正弦信号的幅频特性。(3)设定低频段多正本文档来自技高网...

【技术保护点】
1.一种用于电力系统多输入多输出线性化模型辨识的低频段多正弦信号设计方法，其特征在于，所述的的低频段多正弦信号设计方法包括以下步骤：初始启动阶段：S1、确定低频段多正弦信号的时域长度为N，预期的波峰因子为Cr,set，对应的傅里叶系数记录变量为FS，其幅值记录变量为Mag、相位记录变量为Pha，初始FS、Mag、Pha均为长度为N的零向量；S2、确定低频段多正弦信号的谐波分量频率ωk以及频域幅值ak，其中k＝0,1,…,Nk‑1，Nk≤N/2，然后将幅值ak按照谐波分量频率ωk对应于幅值记录变量Mag的位置填入幅值记录变量Mag，从而得到初始低频段多正弦信号的幅频特性；S3、设定低频段多正弦信号的谐波分量初始相位

【技术特征摘要】
1.一种用于电力系统多输入多输出线性化模型辨识的低频段多正弦信号设计方法，其特征在于，所述的的低频段多正弦信号设计方法包括以下步骤：初始启动阶段：S1、确定低频段多正弦信号的时域长度为N，预期的波峰因子为Cr,set，对应的傅里叶系数记录变量为FS，其幅值记录变量为Mag、相位记录变量为Pha，初始FS、Mag、Pha均为长度为N的零向量；S2、确定低频段多正弦信号的谐波分量频率ωk以及频域幅值ak，其中k＝0,1,…,Nk-1，Nk≤N/2，然后将幅值ak按照谐波分量频率ωk对应于幅值记录变量Mag的位置填入幅值记录变量Mag，从而得到初始低频段多正弦信号的幅频特性；S3、设定低频段多正弦信号的谐波分量初始相位为一组范围在-π到π之间均匀分布的随机数，其中k＝0,1,…,Nk-1，然后将相位按照谐波分量频率ωk对应于相位记录变量Pha的位置填入相位记录变量Pha，从而得到初始低频段多正弦信号的相位特性；S4、将幅值记录变量Mag与相位记录变量Pha进行复数合成，得到初始低频段多正弦信号的频域特性FS，其中关心频段ωk对应的傅里叶系数为S5、对初始低频段多正弦信号的频域特性FS进行傅立叶逆变换，得到低频段多正弦信号的初始时域波形u(t)，并将初始时域波形u(t)存储在变量Signal中，变量Signal的长度为N；S6、计算变量Signal的波峰因子Cr，并将其存储在变量CF中，若变量CF小于预期的波峰因子Cr,set，计Signal_min＝Signal，然后转入步骤S13；否则，设定迭代循环次数Num为0，设置最大迭代次数为Num_max，Num_max>1，进入如下迭代循环阶段；迭代循环阶段：S7、若迭代循环次数Num小于最大迭代次...

【专利技术属性】
技术研发人员：张俊勃，曾繁宏，
申请(专利权)人：华南理工大学，
类型：发明
国别省市：广东,44