CN103362741A - 基于adaline技术的风力发电机组系统辨识方法 - Google Patents

Abstract

本发明公开了一种基于ADALINE技术的风力发电机组系统辨识方法。本发明针对风力发电机组运行机理，确定系统输入输出；然后基于LPV方法，依照有效风速划分稳定工作点；选取合适激励信号，以桨叶节距角或发电机转矩为测试输入信号，输出为发电机转速，在稳定工作点和介于工作点间过渡阶段进行激励，获取全域输入输出数据；提取测试数据中所有稳定工作点输入输出数据，采用基于ADALINE技术的辨识方法，获得各稳定工作点下的简单模型；再通过模型插值方法建立全局系统的LPV模型。本发明算法简单，容易实施，具有良好的运行速度和较低的计算量，模型性能及稳定性能够得到保证。

Description

基于ADALINE技术的风力发电机组系统辨识方法

技术领域

本发明属于风力发电、智能控制和系统辨识技术领域，具体涉及一种基于神经网络中ADALINE技术的系统辨识方法，并结合LPV方法用于复杂、非线性的风力发电机组系统辨识。

背景技术

随着风力发电技术不断发展，以往小型的定桨风力发电机组逐渐被大型变桨、变速风力发电机组取代。这种变桨、变速风力发电机组结构更加复杂、功率输出曲线更加优化，对控制算法提出了更高要求，尤其在工作点和工作点之间切换的过渡阶段。而控制算法的精准性强烈依赖于对象模型的精准性。针对于风力发电机组系统，其主要包含：空气动力子系统，描述自然界的风如何转化为桨叶转动的驱动力；机械子系统，可划分为传动结构和支撑结构两个功能块，其中传动结构将桨叶上的空气动力转矩传递到发电机轴；电力子系统，描述发电机轴上如何将机械能转换为电能；执行子系统，对桨叶伺服行为进行建模。要建立风力发电机组系统各子系统乃至整机系统模型，严格理论意义上的解决方案是采用机理建模方法或者非线性系统辨识方法，但这两种方法并不易于实现。一方面，采用机理建模方法，即利用微分方程和代数方程来理解和描述风力发电机组系统的物理、化学过程需要耗费大量的人力物力，而且对于复杂非线性过程而言几乎不切实际。另一方面，采用非线性系统辨识方法，经验地依赖输入输出数据辨识非线性模型强烈依赖于激励信号，但对于风力发电机系统而言，激励测试可能会产生大量过激的扰动和产品损耗代价，而且难以在全局工作域上进行持续激励测试。为此，需要寻找更好的系统辨识技术，建立既能准确表述系统结构、特征，又能容易辨识参数、实施控制算法的模型。

发明内容

本发明的目的是针对现有技术的不足，提出一种基于ADALINE(ADAptiveLINear Element)技术的风力发电机组系统辨识方法。

本发明主要针对风力发电机组系统由风能捕获至传动系统再至发电机转速这一过程，其中涉及空气动力学子系统、机械子系统、电力子系统部分模块，对不同风速状态下系统输入（桨距或发电机转矩）和系统输出（发电机转速）的模型进行辨识，提出一种结构较简单的系统辨识方法。本发明的对象为变桨或变速风力发电机组，针对实际工况具有调节桨叶节距角β和发电机转矩T_g以控制发电机转速Ω两种工作机制，因此能将辨识分为桨距环辨识和转矩环辨识两个环节。若为桨距环辨识时，系统输入为桨叶节距角β，输出为发电机转速Ω，风速V为外部输入；若为转矩环辨识时，系统输入为发电机转矩T_g，输出为发电机转速Ω，风速V为外部输入。根据LPV(Linear parameter varying)方法，需要找寻一个合适调度参数p，以表征系统工作状态和划分系统稳定工作点。针对风力发电机组，将外部输入风速V表达为两个部分

其中为表征风力发电机组运行状态的有效风速（以下简称有效风速），v表征在该有效风速上的干扰叠加，且V的不同取值决定了系统的不同工作点，故选取有效风速作为调度参数变量。针对风力发电机组运行的某个有效风速段

根据系统特性划分出m个有效风速点，分别标记为

对应系统稳定工作点为WP_k(k=1,2,…,m)，m为自然数。

本发明解决其技术问题所采用的技术方案包括如下步骤：

步骤1.系统测试数据的获取：

根据风力发电机组系统特性，选取激励信号,在稳定工作点WP_k和介于两个稳定工作点WP_k之间的过渡阶段均进行激励，以获取该工作域的全域输入数据、输出数据；采样时间选取系统内部采样时间；桨距环辨识时激励信号为桨叶节距角β，即输入数据为桨叶节距角β（转矩环则为发电机转矩T_g），输出数据为发电机转速Ω，同时包含外部输入风速V；

所述的激励信号为正弦信号与白噪声信号的叠加信号；

步骤2.基于ADALINE技术对稳定工作点进行辨识：

风力发电机组系统描述如下（实际风力发电机系统应当不限于此表达式）：

y(t)=G(p,q^-1)u(t)+v(t)    （1）

其中， $G(p,q^{-1})=\frac{B(p,q^{-1})}{A(p,q^{-1})}=\frac{b_1\left(p\right)q^{-1}+\·\·\·+b_{n_b}\left(p\right)q^{{-n}_b}}{1+a_1\left(p\right)q^{-1}+\·\·\·+a_{n_a}\left(p\right)q^{{-n}_a}},$ 为系统的传递函数，u(t),y(t)分别为系统输入、系统输出；v(t)为噪声信号，其均值为零，方差为有限值；q^-1为单位延迟控制器；p为调度参数，表征系统工作状态，且设定其区间跨度为[p_min,p_max]；a_i,b_j(1≤i≤n_a,1≤j≤n_b,n_a≤n_b)为传递函数分母多项式A(p,q^-1)、分子多项式B(p,q^-1)的系数，i,j,n_a,n_b均为自然数。

在稳定工作点WP_k(k=1,2,…,m)上，其调度参数为一常数p_k(k=1,2,…,m)，此时系统方程描述为：

$y\left(t\right)=\frac{b_1\left(p_k\right)q^{-1}+\·\·\·+{b_n}_b\left(p_k\right)q^{{-n}_b}}{1+a_1\left(p_k\right)q^{-1}+\·\·\·{+a_n}_a\left(p_k\right)q^{{-n}_a}}u\left(t\right)-v\left(t\right)---\left(2\right)$

在ADALINE单元的输入端信号加入TDL，即单位延迟控制器，连接到系统输入和系统输出，即有：

X=[x₁ x₂ … x_H]^T=[u(t-1) …u(t-n_b) y(t-1) … y(t-n_a)]^T    （3）

其中，X为ADALINE的输入，t为迭代指数，对应的ADALINE权值W为：

$W={\left[\begin{array}{cccc}{w}_1& w_2& \·\·\·& w_H\end{array}\right]}^T={\left[\begin{array}{cccccc}{\hat{b}}_1& \·\·\·& {\hat{b}}_{n_b}& -{\hat{a}}_1& \·\·\·& -{\hat{a}}_{n_a}\end{array}\right]}^T---\left(4\right)$

(1≤i≤n_a,1≤j≤n_b,n_a≤n_b)为式（2）中传递函数分母多项式、分子多项式各系数a_i,b_j(1≤i≤n_a,1≤j≤n_b,n_a≤n_b)的估计值，H=n_a+n_b为ADALINE输入个数。

则ADALINE的输出为：

$\hat{y}={\hat{y}}^k=\underset{i=1}{\overset{H}{\mathrm{\Σ}}}{w}_i{x}_i=X^{T}W,p=p_k---\left(5\right)$

通过学习算法，使ADALINE输出和系统输出不断逼近，在一定精度条件下能够认为权值的估计值即为系统传递函数分母多项式、分子多项式各系数，由此确定各稳态工作点模型为：

${\hat{y}}^k\left(t\right)={\hat{G}}_k\left(q^{-1}\right)u\left(t\right),p=p_k---\left(6\right)$

其中

即为各稳态工作点的传递函数估计，不同于式（1），式（6）中的传递函数为定结构和定参数，表达式中不含调度参数变量p。

步骤3.基于ADALINE技术对全局LPV模型进行辨识：

结合LPV方法，将系统全局模型表达为（6）中模型的插值：

$\hat{y}\left(t\right)=\underset{k=1}{\overset{m}{\mathrm{\Σ}}}{\mathrm{\α}}_k\left(p\right){\hat{y}}^k\left(t\right)=\underset{k=1}{\overset{m}{\mathrm{\Σ}}}{\mathrm{\α}}_k\left(p\right){\hat{G}}_k\left(q^{-1}\right)u\left(t\right)---\left(7\right)$

其中α_k(p)(1≤k≤m)为基于调度参数p的插值函数，插值函数可为多项式函数、分段线性函数、三次样条函数或高斯函数。

若插值函数为三次样条函数，则表达式如下：

${\mathrm{\α}}_k\left(p\right)={\mathrm{\γ}}_1^k+{\mathrm{\γ}}_2^{k}p+\underset{i=2}{\overset{M-1}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i+1}^k{|p-{\tilde{p}}_i^k|}^3---\left(8\right)$

其中M为三次样条函数的阶数，M为自然数；

为调度参数取值范围[p_min,p_max]内给定的某一序列，针对任意工作点WP_k(k=1,2,…,m)，

能够取不同序列，这里简化为取相同序列

采用最小二乘法进行辨识，则式（8）中需要估计的系数构成的向量θ为：

$\mathrm{\θ}={\left[\begin{array}{cccccccccc}{\mathrm{\γ}}_1^1& \·\·\·& {\mathrm{\γ}}_M^1& {\mathrm{\γ}}_1^2& \·\·\·& {\mathrm{\γ}}_M^2& \·\·\·& {\mathrm{\γ}}_1^m& \·\·\·& {\mathrm{\γ}}_M^m\end{array}\right]}^T---\left(9\right)$

向量θ对应的调度参数向量φ^k(t)为：

${\mathrm{\φ}}^k\left(t\right)=\left[\begin{array}{ccccc}1& p\left(t\right)& {|p\left(t\right)-{\tilde{p}}_2^{*}|}^3& \·\·\·& {|p\left(t\right)-{\tilde{p}}_{M-1}^{*}|}^3\end{array}\right]---\left(10\right)$

将（9）和（10）代入（7）中可得：

$\hat{y}\left(t\right)=\left[\begin{array}{ccc}{\mathrm{\φ}}^1\left(t\right){\hat{y}}^1\left(t\right)& \·\·\·& {\mathrm{\φ}}^m\left(t\right){\hat{y}}^m\left(t\right)\end{array}\right]\mathrm{\θ}---\left(11\right)$

则θ通过最小二乘准则获得的估计值

如下：

$\widehat{\mathrm{\θ}}=\underset{\mathrm{\θ}}{\min }\underset{t=1}{\overset{N}{\mathrm{\Σ}}}{(y\left(t\right)-\hat{y}\left(t\right))}^2---\left(12\right)$

直接对式（12）进行运算，若数据量较为庞大、三次样条的阶数较高，则会产生病态或者难以求解的问题；所以采用以下方法结构替换：

将全局系统输入u通过式（6）中传递函数获得的估计值

与对应的序列Ψ(p)运算所得的结果

作为ADALINE输入。这里取Ψ(p)为（10）中调度参数向量，

为Ψ(p)与各个稳态工作点模型输出

乘积，即：

$\mathrm{\Ψ}\left(p\right)=\left[\begin{array}{ccccc}1& p& {|p-{\tilde{p}}_2^{*}|}^3& \·\·\·& {|p-{\tilde{p}}_{M-1}^{*}|}^3\end{array}\right]---\left(13\right)$

${\mathrm{\ψ}}^k(p,{\hat{y}}^k)=\mathrm{\Ψ}\left(p\right){\hat{y}}^k---\left(14\right)$

则全局LPV模型的输出

的表达式为：

$\hat{y}=\left[\begin{array}{cccc}\mathrm{\Ψ}\left(p\right){\hat{y}}^1& \mathrm{\Ψ}\left(p\right){\hat{y}}^2& \·\·\·& \mathrm{\Ψ}\left(p\right){\hat{y}}^m\end{array}\right]W---\left(15\right)$

将（15）与（11）相比较，能够得出W的估计即为θ的估计。如此，则建立好（7）中全局系统的LPV模型。模型的精确性能够通过迭代停止准则和基于ADALINE的网络规模改进来调节，以达到要求的模型精度。

本发明有益效果如下：

本发明算法简单，容易实施，具有良好的运行速度和较低的计算量，模型性能及稳定性能够得到保证。

附图说明

图1变桨变速风力发电机组系统结构示意图

图2基于ADALINE技术辨识稳定工作点的方法结构示意图

图3基于ADALINE技术辨识全局LPV模型的方法结构示意图

具体实施方式

下面结合附图内容对本发明做进一步说明。

本发明主要针对风力发电机组系统由风能捕获至传动系统再至发电机转速这一过程，其中涉及空气动力学子系统、机械子系统、电力子系统部分模块，对不同风速状态下系统输入（桨距或发电机转矩）和系统输出（发电机转速）的模型进行辨识，提出一种结构较简单的系统辨识方法。本发明的对象为变桨或变速风力发电机组，针对实际工况具有调节桨叶节距角β和发电机转矩T_g以控制发电机转速Ω两种工作机制，因此能将辨识分为桨距环辨识和转矩环辨识两个环节。若为桨距环辨识时，系统输入为桨叶节距角β，输出为发电机转速Ω，风速V为外部输入；若为转矩环辨识时，系统输入为发电机转矩T_g，输出为发电机转速Ω，风速V为外部输入。根据LPV(Linear parameter varying)方法，需要找寻一个合适调度参数p，以表征系统工作状态和划分系统稳定工作点。针对风力发电机组，将外部输入风速V表达为两个部分

其中

为表征风力发电机组运行状态的有效风速（以下简称有效风速），v表征在该有效风速上的干扰叠加，且的不同取值决定了系统的不同工作点，故选取有效风速作为调度参数变量。针对风力发电机组运行的某个有效风速段

根据系统特性划分出m个有效风速点，分别标记为

对应系统稳定工作点为WP_k(k=1,2,…,m)，m为自然数。

实施例

如图1所示，系统输入为桨叶节距角β和发电机转矩T_g，系统输出为发电机转速Ω；V为风速，λ为叶尖速比；R为风轮半径，T_a,F_a为风轮捕获风能产生的转矩和推力，f(λ,β,V²)为某一非线性函数。根据LPV方法，需要找寻一个合适调度参数p，以表征系统工作状态和划分系统稳定工作点。针对风力发电机组，将外部输入风速V表达为两个部分

其中为表征风力发电机组运行状态的有效风速（以下简称有效风速），v表征在该有效风速上的干扰叠加，且

的不同取值决定了系统的不同工作点，故选取有效风速

作为调度参数变量。针对风力发电机组运行的某个有效风速段

根据系统特性划分出m个有效风速点，分别标记为

对应系统稳定工作点为WP_k(k=1,2,…,m)，m为自然数。

步骤1.系统测试数据的获取：

本实施例针对桨距环辨识进行说明，则选取激励信号为桨叶节距角β，信号类型为正弦信号与白噪声信号的叠加，在稳定工作点WP_k和介于两个稳定工作点WP_k之间的过渡阶段均进行激励，以获取该工作域的全域输入数据、输出数据；采样时间选取系统内部采样时间。获取的全域数据，包括输入数据为桨叶节距角β、输出数据为发电机转速Ω和外部输入风速V数据。

步骤2.基于ADALINE技术对稳定工作点进行辨识：

本实施例针对桨距环辨识，对应风力发电机组高风速工作段的变桨工作状态。此时系统输入为桨叶节距角β，输出为发电机转速Ω，调度参数为有效风速

。为简化说明，将给定的系统方程描述替换输入、输出及调度参数后如下，实际风力发电机系统应当不限于此表达式：

$\mathrm{\Ω}\left(t\right)=G(\stackrel{\‾}{V},q^{-1})\mathrm{\β}\left(t\right)+v\left(t\right)---\left(1\right)$

其中 $G(\stackrel{\‾}{V},q^{-1})=\frac{B(\stackrel{\‾}{V},q^{-1})}{A(\stackrel{\‾}{V},q^{-1})}=\frac{b_1\left(\stackrel{\‾}{V}\right)q^{-1}+\·\·\·+{b_n}_b\left(\stackrel{\‾}{V}\right)q^{{-n}_b}}{1+a_1\left(\stackrel{\‾}{V}\right)q^{-1}+\·\·\·+{a_n}_a\left(\stackrel{\‾}{V}\right)q^{{-n}_a}}$ 为系统的传递函数，β(t),Ω(t)分别为系统输入、系统输出。有效风速

为调度参数，表征不同工作状态。v(t)为噪声信号，其均值为零，方差为有限值。q^-1为单位延迟控制器。设定有效风速段为

a_i,b_j(1≤i≤n_a,1≤j≤n_b,n_a≤n_b)为传递函数分母多项式

分子多项式

的系数，i,j,n_a,n_b均为自然数。

针对稳定工作点WP_k(k=1,2,…,m)，对应调度参数为

则系统方程描述为：

$\mathrm{\Ω}\left(t\right)=\frac{b_1\left({\stackrel{\‾}{V}}_k\right)q^{-1}+{\·\·\·b_n}_b\left({\stackrel{\‾}{V}}_k\right)q^{{-n}_b}}{1+a_1\left({\stackrel{\‾}{V}}_k\right)q^{-1}+\·\·\·+{a_n}_a\left({\stackrel{\‾}{V}}_k\right)q^{{-n}_a}}\mathrm{\β}\left(t\right)+v\left(t\right)---\left(2\right)$

这里由于

已知，则分子分母中各系数为常数。采用如图2所示方法结构进行辨识，即有：

X=[x₁ x₂ … x_H]^T=[β(t-1) … β(t-n_b) Ω(t-1) … Ω(t-n_a)]^T   （3）

其中，X为ADALINE的输入，t为迭代指数，对应的ADALINE权值W为：

$W={\left[\begin{array}{cccc}{w}_1& w_2& \·\·\·& w_H\end{array}\right]}^T={\left[\begin{array}{cccccc}{\hat{b}}_1& \·\·\·& {\hat{b}}_{n_b}& -{\hat{a}}_1& \·\·\·& -{\hat{a}}_{n_a}\end{array}\right]}^T---\left(4\right)$

(1≤i≤n_a,1≤j≤n_b,n_a≤n_b)为（2）中传递函数分母多项式、分子多项式各系数a_i,b_j(1≤i≤n_a,1≤j≤n_b,n_a≤n_b)的估计,H=n_a+n_b为ADALINE输入个数。

则ADALINE的输出

为：

$\widehat{\mathrm{\Ω}}={\widehat{\mathrm{\Ω}}}^k=\underset{i=1}{\overset{H}{\mathrm{\Σ}}}{w}_i{x}_i=X^{T}W,\stackrel{\‾}{V}={\stackrel{\‾}{V}}_k---\left(5\right)$

通过学习算法，使ADALINE输出和系统输出不断逼近，在一定精度条件下能够认为权值的估计值即为系统传递函数分母多项式、分子多项式各系数，由此确定各稳态工作点模型为：

${\widehat{\mathrm{\Ω}}}^k\left(t\right)={\hat{G}}_k\left(q^{-1}\right)\mathrm{\β}\left(t\right),\stackrel{\‾}{V}={\stackrel{\‾}{V}}_k---\left(6\right)$

其中即为各稳态工作点的传递函数估计，不同于式（1），式（6）中的传递函数为定结构和定参数，表达式中不含调度参数变量。

步骤3.基于ADALINE技术对全局LPV模型进行辨识：

结合LPV方法，将系统全局模型表达为（6）中模型的插值：

$\widehat{\mathrm{\Ω}}\left(t\right)=\underset{k=1}{\overset{m}{\mathrm{\Σ}}}{\mathrm{\α}}_k\left(\stackrel{\‾}{V}\right){\widehat{\mathrm{\Ω}}}^k\left(t\right)=\underset{k=1}{\overset{m}{\mathrm{\Σ}}}{\mathrm{\α}}_k\left(\stackrel{\‾}{V}\right)G_k\left(q^{-1}\right)\mathrm{\β}\left(t\right)---\left(7\right)$

其中

为基于调度参数

的插值函数，插值函数可为多项式函数、分段线性函数、三次样条函数或高斯函数等。

若插值函数为三次样条函数，则表达式如下：

${\mathrm{\α}}_k\left(\stackrel{\‾}{V}\right){=\mathrm{\γ}}_1^k+{\mathrm{\γ}}_2^k\stackrel{\‾}{V}+\underset{i=2}{\overset{M-1}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i+1}^k{|\stackrel{\‾}{V}-{\tilde{V}}_i^k|}^3---\left(8\right)$

其中M为三次样条的阶数，M为自然数；为调度参数取值范围

内给定序列，针对任意工作点WP_k(k=1,2,…,m)可以取不同序列，这里简化为取同样序列

采用最小二乘法进行辨识，则式（8）中需要估计的系数构成的向量θ为：

$\mathrm{\θ}={\left[\begin{array}{cccccccccc}{\mathrm{\γ}}_1^1& \·\·\·& {\mathrm{\γ}}_M^1& {\mathrm{\γ}}_1^2& \·\·\·& {\mathrm{\γ}}_M^2& \·\·\·& {\mathrm{\γ}}_1^m& \·\·\·& {\mathrm{\γ}}_M^m\end{array}\right]}^T---\left(9\right)$

向量θ对应的调度参数向量φ^k(t)为：

${\mathrm{\φ}}^k\left(t\right)=\left[\begin{array}{ccccc}1& \stackrel{\‾}{V}\left(t\right)& {|\stackrel{\‾}{V}\left(T\right)-{\tilde{V}}_2^{*}|}^3& \·\·\·& {|\stackrel{\‾}{V}\left(T\right)-{\tilde{V}}_{M-1}^{*}|}^3\end{array}\right]---\left(10\right)$

将（9）和（10）代入（7）中可得：

$\widehat{\mathrm{\Ω}}\left(t\right)=\left[\begin{array}{ccc}{\mathrm{\φ}}^1\left(t\right){\widehat{\mathrm{\Ω}}}^1\left(t\right)& \·\·\·& {\mathrm{\φ}}^m\left(t\right){\widehat{\mathrm{\Ω}}}^m\left(t\right)\end{array}\right]\mathrm{\θ}---\left(11\right)$

则θ通过最小二乘准则获得的估计值

如下：

$\widehat{\mathrm{\θ}}=\underset{\mathrm{\θ}}{\min }\underset{t=1}{\overset{N}{\mathrm{\Σ}}}{(\mathrm{\Ω}\left(t\right)-\widehat{\mathrm{\Ω}}\left(t\right))}^2---\left(12\right)$

直接进行（12）的运算，若数据量较为庞大、三次样条的阶数较高，则会产生病态或者难以求解的问题。因此采用图3中的结构进行估计。

如图3所示，全局数据包含系统测试过程中的全局输入、全局输出及全局调度参数。左上角

为式（6）中所得到的稳态工作点的模型，则将全局系统输入u通过式（6）中传递函数

获得的估计值

与对应的序列Ψ(p)运算所得的结果

作为ADALINE输入。

在此方法结构中，关键的设计为和

针对三次样条，取

为调度参数向量，为

与各个稳态工作点模型输出

乘积，即：

$\mathrm{\Ψ}\left(\stackrel{\‾}{V}\right)=\left[\begin{array}{ccccc}1& \stackrel{\‾}{V}& {|\stackrel{\‾}{V}-{\stackrel{\‾}{V}}_2^{*}|}^3& \·\·\·& {|\stackrel{\‾}{V}-{\tilde{V}}_{M-1}^{*}|}^3\end{array}\right]---\left(13\right)$

${\mathrm{\ψ}}^k(\stackrel{\‾}{V},{\widehat{\mathrm{\Ω}}}^k)=\mathrm{\Ψ}\left(\stackrel{\‾}{V}\right){\widehat{\mathrm{\Ω}}}^k---\left(14\right)$

则全局LPV模型的输出的表达式为：

$\stackrel{\‾}{\mathrm{\Ω}}=\left[\begin{array}{cccc}\mathrm{\Ψ}\left(\stackrel{\‾}{V}\right){\widehat{\mathrm{\Ω}}}^1& \mathrm{\Ψ}\left(\stackrel{\‾}{V}\right){\widehat{\mathrm{\Ω}}}^2& \·\·\·& \mathrm{\Ψ}\left(\stackrel{\‾}{V}\right){\widehat{\mathrm{\Ω}}}^m\end{array}\right]W---\left(15\right)$

将（15）与（11）相比较，可以看到W的估计即为θ的估计。如此，则建立好（7）中全局系统的LPV模型。模型的精确性可以通过迭代停止准则和基于ADALINE的网络规模改进来调节，以达到要求的模型精度。

Claims (1)

1.基于ADALINE技术的风力发电机组系统辨识方法，其特征在于包括如下步骤：

步骤1.系统测试数据的获取：

根据风力发电机组系统特性，选取激励信号,在稳定工作点WP_k和介于两个稳定工作点之间的过渡阶段均进行激励，以获取该工作域的全域输入数据、输出数据；采样时间选取系统内部采样时间；桨距环辨识时激励信号为桨叶节距角β，即输入数据为桨叶节距角β，输出为发电机转速Ω，风速V为外部输入；转矩环辨识时，系统输入为发电机转矩T_g，输出数据为发电机转速Ω，风速V为外部输入；

所述的激励信号为正弦信号与白噪声信号的叠加信号；

步骤2.基于ADALINE技术对稳定工作点进行辨识：

风力发电机组系统描述如下：

y(t)=G(p,q^-1)u(t)+v(t)    （1）

其中， $G(p,q^{-1})=\frac{B(p,q^{-1})}{A(p,q^{-1})}=\frac{b_1\left(p\right)q^{-1}+\·\·\·+{b_n}_b\left(p\right)q^{{-n}_b}}{1+a_1\left(p\right)q^{-1}+\·\·\·+{a_n}_a\left(p\right)q^{{-n}_a}},$ 为系统的传递函数，u(t),y(t)分别为系统输入、系统输出；v(t)为噪声信号，其均值为零，方差为有限值；q^-1为单位延迟控制器；p为调度参数，表征系统工作状态，且设定其区间跨度为[p_min,p_max]；a_i,b_j(1≤i≤n_a,1≤j≤n_b,n_a≤n_b)为传递函数分母多项式A(p,q^-1)、分子多项式B(p,q^-1)的系数，i,j,n_a,n_b均为自然数；在稳定工作点WP_k(k=1,2,…,m)上，其调度参数为一常数p_k(k=1,2,…,m)，m为自然数，此时系统方程描述为：

$y\left(t\right)=\frac{b_1\left(p_k\right)q^{-1}+\·\·\·+{b_n}_b\left(p_k\right)q^{{-n}_b}}{1+a_1\left(p_k\right)q^{-1}+\·\·\·+{a_n}_a\left(p_k\right)q^{{-n}_a}}u\left(t\right)+v\left(t\right)---\left(2\right)$

在ADALINE单元的输入端信号加入TDL，即单位延迟控制器，连接到系统输入和系统输出，即有：

X=[x₁ x₂ … x_H]^T=[u(t-1) … u(t-n_b) y(t-1) … y(t-n_a)]^T   （3）

其中，X为ADALINE的输入，t为迭代指数，对应的ADALINE权值W为：

$W={\left[\begin{array}{cccc}{w}_1& w_2& \·\·\·& w_H\end{array}\right]}^T={\left[\begin{array}{cccccc}{\hat{b}}_1& \·\·\·& {\hat{b}}_{n_b}& -{\hat{a}}_1& \·\·\·& -{\hat{a}}_{n_a}\end{array}\right]}^T---\left(4\right)$

(1≤i≤n_a,1≤j≤n_b,n_a≤n_b)为式（2）中传递函数分母多项式、分子多项式各系数a_i,b_j(1≤i≤n_a,1≤j≤n_b,n_a≤n_b)的估计值，H=n_a+n_b为ADALINE输入个数；

则ADALINE的输出

为：

$\hat{y}={\hat{y}}^k=\underset{i=1}{\overset{H}{\mathrm{\Σ}}}{w}_i{x}_i=X^{T}W,p=p_k---\left(5\right)$

通过学习算法，在一定精度条件下能够认为权值的估计值即为系统传递函数分母多项式、分子多项式各系数，由此确定各稳态工作点模型为：

${\hat{y}}^k\left(t\right)={\hat{G}}_k\left(q^{-1}\right)u\left(t\right),p=p_k---\left(6\right)$

其中

即为各稳态工作点的传递函数估计，且式（6）中的传递函数为定结构和定参数，表达式中不含调度参数变量p；

步骤3.基于ADALINE技术对全局LPV模型进行辨识：

结合LPV方法，将系统全局模型表达为（6）中模型的插值：

$\hat{y}\left(t\right)=\underset{k=1}{\overset{m}{\mathrm{\Σ}}}{\mathrm{\α}}_k\left(p\right){\hat{y}}^k\left(t\right)=\underset{k=1}{\overset{m}{\mathrm{\Σ}}}{\mathrm{\α}}_k\left(p\right){\hat{G}}_k\left(q^{-1}\right)u\left(t\right)---\left(7\right)$

其中α_k(p)(1≤k≤m)为基于调度参数p的插值函数，插值函数为多项式函数、分段线性函数、三次样条函数或高斯函数；

若插值函数为三次样条函数，则表达式如下：

${\mathrm{\α}}_k\left(p\right)={\mathrm{\γ}}_1^k+{\mathrm{\γ}}_2^{k}p+\underset{i=2}{\overset{M-1}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i+1}^k{|p-{\tilde{p}}_i^k|}^3---\left(8\right)$

其中M为三次样条函数的阶数，M为自然数；

为调度参数取值范围[p_min,p_max]内给定的某一序列，针对任意工作点WP_k(k=1,2,…,m)，

能够取不同序列，这里简化为取相同序列

采用最小二乘法进行辨识，则式（8）中需要估计的系数构成的向量θ为：

$\mathrm{\θ}={\left[\begin{array}{cccccccccc}{\mathrm{\γ}}_1^1& \·\·\·& {\mathrm{\γ}}_M^1& {\mathrm{\γ}}_1^2& \·\·\·& {\mathrm{\γ}}_M^2& \·\·\·& {\mathrm{\γ}}_1^m& \·\·\·& {\mathrm{\γ}}_M^m\end{array}\right]}^T---\left(9\right)$

向量θ对应的调度参数向量φ^k(t)为：

${\mathrm{\φ}}^k\left(t\right)=\left[\begin{array}{ccccc}1& p\left(t\right)& {|p\left(t\right)-{\tilde{p}}_2^{*}|}^3& \·\·\·& {|p\left(t\right)-{\tilde{p}}_{M-1}^{*}|}^3\end{array}\right]---\left(10\right)$

将（9）和（10）代入（7）中可得：

$\hat{y}\left(t\right)=\left[\begin{array}{ccc}{\mathrm{\φ}}^1\left(t\right){\hat{y}}^1\left(t\right)& \·\·\·& {\mathrm{\φ}}^m\left(t\right){\hat{y}}^m\left(t\right)\end{array}\right]\mathrm{\θ}---\left(11\right)$

则θ通过最小二乘准则获得的估计值

如下：

$\widehat{\mathrm{\θ}}=\underset{\mathrm{\θ}}{\min }\underset{t=1}{\overset{N}{\mathrm{\Σ}}}{(y\left(t\right)-\hat{y}\left(t\right))}^2---\left(12\right)$

采用以下方法结构替换式（12）：

将全局系统输入u(t)通过式（6）中传递函数

获得的估计值

与对应的序列Ψ(p)运算所得的结果

作为ADALINE输入，其中Ψ(p)为（10）中调度参数向量，为Ψ(p)与各个稳态工作点模型输出

乘积，即：

$\mathrm{\Ψ}\left(p\right)=\left[\begin{array}{ccccc}1& p& {|p-{\tilde{p}}_2^{*}|}^3& \·\·\·& {|p-{\tilde{p}}_{M-1}^{*}|}^3\end{array}\right]---\left(13\right)$

${\mathrm{\Ψ}}^k(p,{\hat{y}}^k)=\mathrm{\Ψ}\left(p\right){\hat{y}}^k---\left(14\right)$

则全局LPV模型的输出

的表达式为：

$\hat{y}=\left[\begin{array}{cccc}\mathrm{\Ψ}\left(p\right){\hat{y}}^1& \mathrm{\Ψ}\left(p\right){\hat{y}}^2& \·\·\·& \mathrm{\Ψ}\left(p\right){\hat{y}}^m\end{array}\right]W---\left(15\right)$

将（15）与（11）相比较，得出W的估计即为θ的估计，则（7）中全局系统的LPV模型建立完成。

Images