CN113114128A

CN113114128A - 基于广义Bouc-Wen逆模型的压电前馈补偿方法

Info

Publication number: CN113114128A
Application number: CN202110514171.6A
Authority: CN
Inventors: 张泉; 高源蓬; 李清灵; 尹达一
Original assignee: Shanghai Institute of Technical Physics of CAS
Current assignee: Shanghai Institute of Technical Physics of CAS
Priority date: 2021-05-12
Filing date: 2021-05-12
Publication date: 2021-07-13
Anticipated expiration: 2041-05-12
Also published as: CN113114128B

Abstract

本发明公开了一种基于广义Bouc‑Wen逆模型的压电前馈补偿方法。针对压电执行器迟滞的率相关性，采用Hammerstein模型来描述压电率相关逆迟滞特性，以广义Bouc‑Wen非对称逆模型表征逆迟滞非线性，以ARX模型表征逆迟滞率相关特性；采用自适应天牛群智能优化算法对模型参数进行辨识；将建立Hammerstein率相关逆迟滞模型与压电执行器级联构成整体线性化系统来对迟滞非线性进行补偿。本发明通过直接建立广义Bouc‑Wen非对称逆迟滞模型和采用自适应天牛群智能优化算法，降低了建模难度和提高了模型准确度，从而使压电执行器能够进行高精度的定位。

Description

基于广义Bouc-Wen逆模型的压电前馈补偿方法

技术领域：

本发明涉及基于广义Bouc-Wen逆模型的压电前馈补偿方法。

背景技术：

在卫星激光通信和空间天文观测等领域中，为了达到秒级甚至是毫秒级的精度，大型光电跟踪系统广泛采用复合轴控制技术，其是实现大范围、高精度跟踪指标的最有效的控制结构。复合轴控制(Compound axis control)是二维关联控制系统的一种实现形式，以快摆镜(Fast Steering Mirror,FSM)为执行机构的控制系统是复合轴系统的精级跟踪系统，系统能够达到的精度主要由快摆镜执行机构决定的。快摆镜机构通常采用音圈电机(Voice Coil Motor,VCM)和压电执行器(piezoelectric actuator，PZT)作为驱动器。相比于前者，后者具有谐振频率高、位移分辨率大的优点，然而其固有的迟滞非线性会对定位精度产生不利影响。学界普遍采用建立前馈补偿器的方法进行压电迟滞补偿，前馈补偿器与被控对象级联构成整体线性化系统。压电执行器的前馈补偿主要有电荷控制补偿和电压控制补偿两种方式。电荷控制补偿方式是通过专门的电荷放大器来对压电执行器迟滞非线性进行抑制。电压控制补偿方式是引入建立好的压电逆迟滞数学模型来补偿迟滞非线性，相比于电荷控制，电压控制不需要额外的电路，且从控制理论的角度来看，电压控制很有意义。

电压控制补偿方式的重点就是压电迟滞模型的建立。从数学特点上来说，压电迟滞曲线具有四点性质：1)多值映射性：相同的输入电压在电压值升高和电压值下降阶段分别对应不同的位移输出；2)非局部记忆性：压电执行器的输出位移既和输入电压的瞬时值相关，又和输出位移的历史极值有关；3)率相关特性：随着输入电压频率的不断增大,压电迟滞曲线呈现频率相关性；4)非对称性：在输入电压升高和输入电压下降两个阶段的电压-位移曲线是不对称的。因此，对于压电陶瓷迟滞曲线的建模要从这四点性质出发才能保证拟合精度。目前，压电迟滞模型的建立和改进基本都是围绕这四点性质进行的。压电迟滞模型主要有微分方程模型和数学算子模型两类。数学算子迟滞模型无法表示压电执行器的动态特征，而微分方程迟滞模型既能表征迟滞的数学特性，还能描述压电执行器的动态特性。有代表性的微分方程迟滞模型是Bouc-Wen模型，其最早是由Bouc于1967年提出的，并在1976年由Wen完善，其将压电执行器等效成一个单自由度的质量-弹簧-阻尼系统。Bouc-Wen模型只利用了一个辅助的微分方程来描述各种迟滞形状，因而受到各类学者的青睐。

传统Bouc-Wen模型只能描述关于原点对称的迟滞曲线，且针对的都是单一频率迟滞曲线，无法满足压电执行器变频控制的要求。虽然可以通过改变模型中的迟滞函数来实现非对称性表征，但是微分方程中过多的未知参数会导致数值解的发散，不利于控制，且非对称模型求逆过程复杂。而且，传统Bouc-Wen模型描述的是固定频率下的压电迟滞曲线，无法表征压电迟滞曲线的率相关性。因此，需要针对上述问题对传统Bouc-Wen模型进行改进，并与压电执行器级联构成整体线性化系统来对迟滞非线性进行补偿。

发明内容：

针对上述应用背景，本发明提出了基于广义Bouc-Wen逆模型的压电动态迟滞前馈补偿方法，包括以下步骤：

1)、根据逆函数定理，建立直接Bouc-Wen逆迟滞模型，并在模型中引入了多项式来表征压电迟滞曲线非对称性；

2)、采用Hammerstein模型来描述压电率相关逆迟滞特性，以广义Bouc-Wen非对称逆模型表征逆迟滞非线性，以ARX模型表征逆迟滞率相关特性；

3)、以模型输出数据与实际采集数据的均方根误差为目标函数，采用自适应天牛群智能优化算法对广义Bouc-Wen逆迟滞模型和ARX率相关模型参数进行辨识；

4)、将建立Hammerstein率相关逆迟滞模型与压电执行器级联构成整体线性化系统来对迟滞非线性进行补偿。

具体地，所述步骤1中：

1)、根据逆函数定理，得到以位移为自变量和以驱动电压为因变量的Bouc-Wen逆迟滞模型，表达式为：

u＝d₁y+d₂h (1)

其中，u为驱动电压，h为迟滞状态变量，y为位移，α，β，γ是迟滞曲线形状系数，d₁，d₂为系数；

2)、为了反映压电迟滞曲线的非对称性质，在逆模型中引入多项式，表达式为：

其中，f’和g’是关于y的多项式，α’，β’，γ’是逆迟滞曲线形状系数，p_i’，q_i’表示多项式f’和g’的系数；

3)、离散化的广义Bouc-Wen逆迟滞模型表达式为：

h(t)-h(t-1)＝α[y(t)-y(t-1)]...

-β|[y(t)-y(t-1)]|h(t-1)...

-γ[y(t)-y(t-1)]|h(t-1)|

具体地，所述步骤2中：通过串联非线性静态模块和线性动态模块的Hammerstein模型来描述压电逆迟滞非线性，以广义Bouc-wen逆迟滞模型表征逆迟滞静态非线性，以ARX模型表征逆迟滞率相关特性。

ARX模型是离散系统传递函数模型，如式(4)所示：

式(4)对应的差分方程如式(5)所示：

out(t)+a₁out(t-1)+a₂out(t-2)+…+a_nout(t-n)

＝b₁in(t-d)+b₂in(t-d-1)+…+b_min(t-d-m+1)+ε(t) (5)

其中，in(t)是输入信号，in(t-1)是前一采样时刻输入信号，out(t)是输出信号，out(t-1)是前一采样时刻输出信号，ε(t)为残差信号。

ARX模型的建立是由低频的压电执行器输入输出信号辨识出广义Bouc-Wen逆模型来表征逆迟滞率无关非线性部分，再将扫频压电执行器输出位移信号代入Bouc-Wen逆迟滞模型得到中间量v(t)，然后根据v(t)和输入电压u(t)得到ARX模型。

具体地，所述步骤3中：

1)：在群体优化算法的启发下，通过对天牛群体觅食行为的模拟来对广义Bouc-Wen逆迟滞模型参数和ARX模型参数进行辨识。假设在D维度待优化问题中，第i只天牛的位置和速度分别表示为X_i＝[x_i，1，x_i，2，…，x_i，D]和V_i＝[v_i，1，v_i，2，…，v_i，D]。确定在第t次迭代中每只天牛所经过的最佳位置pbest以及群体所发现的最佳位置gbest，通过跟踪这两个最佳位置按照式(6)更新每只天牛的速度和位置。

v_i(t)＝ωv_i(t-1)+c₁r₁[pbest_i-x_i(t-1)]+c₁r₁[gbest-x_i(t-1)]，

x_i(t)＝x_i(t-1)+λv_i(t-1)+(1-k)δ(t-1) (6)

其中，ω为惯性权重系数，c₁和c₂为正的加速常数，r₁和r₂为0到1之间均匀分布的随机数，λ为正常数。在天牛群优化中，以每只天牛运动速度的方向表示搜索方向。通过设置每只天牛的位置范围[x_min，x_max]，则可以对每只天牛的移动进行适当的限制。

采用线性变化的权重，让惯性权重从最大值ω_max线性减小到最小值ω_min，随算法迭代次数的变化公式为：

其中，ω_max，ω_min分别表示ω的最大值和最小值，t表示当前迭代步数，T表示最大迭代步数；

2)、每只天牛的搜索行为是通过迭代机制来模拟的，如式(8)所示：

x(t)＝x(t-1)+δ(t)，

其中，x(t)为第t次迭代天牛的位置，δ(t)表示天牛运动位置的增量，step(t)为第t次迭代天牛的步长，sign(.)为符号函数，f(.)为表示天牛天线处气味强度的目标函数。

天牛左边和右边天线的空间坐标为：

其中，x_r(t)为第t次迭代的天牛右天线空间坐标，x_l(t)为第t次迭代天牛左天线的空间坐标，d(t)表示第t次迭代左右天线之间的距离。

天牛天线距离d(t)和步长step(t)需要随着迭代次数的增加而减小，如式(10)和式(11)所示：

step(t)＝eta·step(t-1) (10)

d(t)＝step(t)/c (11)

其中，c为常数，需要根据实际问题设置。

为了避免算法早熟，采用自适应的衰减因子，如式(12)所示：

其中，t表示当前迭代步数，T表示最大迭代步数；

3)、初始化自适应天牛群优化算法的迭代次数，种群数NP，加速常数c1和c2，惯性权重最大值ω_max，惯性权重最小值ω_min，常数λ，初始步长step(1)，天线距离衰减因子c，衰减因子初值eta(1)，以及每只天牛位置范围。以模型输出电压数据与实际采集数据的均方根误差为目标函数，对Duhem逆模型参数α,β,γ,p_i和q_i，以及ARX模型参数进行辨识。

具体地，所述步骤4中：基于Hammerstein率相关逆迟滞模型建立了压电执行器动态迟滞前馈补偿器，目标位移通过该前馈补偿器的解算得到压电执行器驱动电压，经过数模转化模块(DA converter,DAC)与压电驱动模块对压电执行器进行控制。

本发明可以避免复杂的正迟滞模型求逆运算，只需要采用智能优化算法辨识出模型参数就可以得到非对称、率相关逆迟滞模型，该模型可以直接作为压电执行器控制系统的前馈补偿器。自适应天牛群优化算法在进行Hammerstein率相关逆迟滞模型参数时不易陷入局部最优，拟合精度高，收敛速度快。通过本发明所提出动态迟滞前馈补偿方法可以实现压电执行器的高精度定位和控制。

附图说明：

图1是基于广义Bouc-Wen逆模型的压电执行器动态迟滞前馈补偿系统结构图。

图2是压电执行器控制系统框图。

图3是0.5Hz频率下广义Bouc-wen逆迟滞曲线与实际逆迟滞曲线。

图4是Hammerstein率相关逆迟滞曲线与实际逆迟滞曲线；其中，图(a)是1Hz频率下曲线；图(b)是10Hz频率下曲线；图(c)是50Hz频率下曲线；图(d)是100Hz频率下曲线。

图5是基于动态迟滞前馈补偿器的实时跟踪及误差曲线。

具体实施方式：

下面结合具体实施例对本发明作进一步描述，在此发明的示意性实施例以及说明用来解释本发明，但并不作为对本发明的限定。

搭建实际控制系统用于数据采集和压电执行器定位控制。压电执行器控制系统由SGS微位移传感器，SGS信号调理模块，主控模块，压电执行器驱动模块，压电执行器组成，其中主控模块采用半实物实时仿真平台。主控模块产生压电执行器驱动信号，经过驱动模块来控制压电执行器，SGS信号调理模块将检测的实际位移反馈给主控模块，控制系统框图如图2所示。

以最大位移为30μm，最大驱动电压100V的压电执行器为被控对象。首先利用广义Bouc-wen逆模型对Hammerstein动态逆迟滞模型中的非线性部分进行建模。对压电执行器施加最大幅度为80V，频率为0.5Hz的单频变幅值正弦电压信号，采集对应的输出位移信号，采样率为10kHz。采用广义Bouc-wen逆迟滞模型进行该迟滞曲线拟合，模型中的多项式阶数取1。

借助自适应天牛群优化算法辨识模型参数α,β,γ,p₀,p₁,q₁,q₀,。设自适应天牛群优化算法的迭代次数为300，种群数NP＝120，加速常数c₁＝2.8,c₂＝1.3，权重最大值ω_max＝0.9，权重最小值ω_min＝0.4，λ＝0.95，初始步长step(1)＝2，c＝2，衰减因子初值eta(1)＝0.95，位置范围设为[-10,10]。寻优算法的核心问题是选取目标函数：

其中，F为模型驱动电压与实际驱动电压的均方根误差，N为数据样本数。自适应天牛群优化算法得到在0.5Hz的单频变幅值正弦电压信号驱动下的广义Bouc-Wen逆迟滞模型参数值如表1所示。

表1广义Bouc-Wen逆迟滞模型参数

将表1中的模型参数值代入式(3)，得到的逆迟滞拟合曲线，拟合精度为0.3172V，相对误差0.39％，如图3所示。

压电执行器控制系统生成幅值为80V，频率范围为1～100Hz的扫频驱动信号，得到1～100Hz扫频微位移信号。采集的输出位移y(k)代入前文建立的广义Bouc-wen逆迟滞模型得到ARX模型的输入v(k)，ARX模型的输出为施加给压电执行器的输入电压u(k)。

动态线性系统的阶次选择2阶，利用自适应天牛群优化算法辨识出式(5)的ARX逆迟滞率相关模型，如式(26)所示：

所建立的动态迟滞逆模型可以有效描述出1～100Hz频率范围内压电逆迟滞曲线，将实际采集的逆迟滞曲线与所建立Hammerstein模型拟合的逆迟滞曲线进行比较，如图4所示，拟合逆迟滞曲线与实际曲线的均方根误差和相对误差如表2所示。

表2迟滞曲线拟合误差

在搭建的验证平台进行目标位移跟踪试验，验证该动态迟滞前馈补偿器在压电执行器定位中有效性。实时跟踪频率范围为1～100Hz，最大位移为24μm的变频变幅值位移信号，跟踪误差值均方根值为0.2932μm，相对误差为1.22％，有效抑制了压电执行器的迟滞非线性，基于动态迟滞前馈补偿器跟踪和误差曲线如图5所示。

Claims

1.一种基于广义Bouc-Wen逆模型的压电前馈补偿方法，其特征在于包括以下步骤：

2.根据权利要求1所述的基于广义Bouc-Wen逆模型的压电前馈补偿方法，其特征在于，步骤1)中所述Bouc-Wen逆迟滞模型直接建立过程具体为：

u＝d₁y+d₂h (1)

2)、在逆模型中引入多项式，表达式为：

其中，

和g是关于y的多项式，α，β，γ是逆迟滞曲线形状系数，p_i，q_i表示多项式

和g的系数；

3)、离散化的广义Bouc-Wen逆迟滞模型表达式为：

h(t)-h(t-1)＝α[y(t)-y(t-1)]...-β|[y(t)-y(t-1)]|h(t-1)...-γ[y(t)-y(t-1)]|h(t-1)|

。

3.根据权利要求1所述的基于广义Bouc-Wen逆模型的压电前馈补偿方法，其特征在于，步骤2)中所述的ARX模型是离散系统传递函数模型，如式(4)所示：

式(4)对应的差分方程如式(5)所示：

out(t)+a₁out(t-1)+a₂out(t-2)+…+a_nout(t-n)

＝b₁in(t-d)+b₂in(t-d-1)+…+b_min(t-d-m+1)+ε(t) (5)

其中，in(t)是输入信号，in(t-1)是前一采样时刻输入信号，out(t)是输出信号，out(t-1)是前一采样时刻输出信号，ε(t)为残差信号；

ARX模型的建立是由低频的压电执行器输入输出信号辨识出广义Bouc-Wen逆模型来表征逆迟滞率无关非线性部分，再将扫频压电执行器输出位移信号代入Bouc-Wen逆迟滞模型得到中间量v(t)，然后根据v(t)和输入电压乩(t)得到ARX模型。

4.根据权利要求1所述的基于广义Bouc-Wen逆模型的压电前馈补偿方法，其特征在于，步骤3)中所述自适应天牛群智能优化算法对广义Bouc-Wen逆迟滞模型和ARX率相关模型参数辨识的具体过程为：

1)：在群体优化算法的启发下，通过对天牛群体觅食行为的模拟来对广义Bouc-Wen逆迟滞模型参数和ARX模型参数进行辨识，假设在D维度待优化问题中，第i只天牛的位置和速度分别表示为X_i＝[x_i，1，x_i，2，…，x_i，D]和V_i＝[v_i，1，v_i，2，…，V_i，D]；确定在第t次迭代中每只天牛所经过的最佳位置pbest以及群体所发现的最佳位置gbest，通过跟踪这两个最佳位置按照式(6)更新每只天牛的速度和位置；

v_i(t)＝ωv_i(t-1)+c₁r₁[pbest_i-x_i(t-1)]+c₁r₁[gbest-x_i(t-1)]，

x_i(t)＝x_i(t-1)+λv_i(t-1)+(1-k)δ(t-1) (6)

其中，ω为惯性权重系数，c₁和c₂为正的加速常数，r₁和r₂为0到1之间均匀分布的随机数，λ为正常数。在天牛群优化中，以每只天牛运动速度的方向表示搜索方向，通过设置每只天牛的位置范围[x_min，x_max]，则对每只天牛的移动进行适当的限制；

x(t)＝x(t-1)+δ(t)，

其中，x(t)为第t次迭代天牛的位置，δ(t)表示天牛运动位置的增量，step(t)为第t次迭代天牛的步长，sign(.)为符号函数，f(.)为表示天牛天线处气味强度的目标函数；

天牛左边和右边天线的空间坐标为：

其中，x_r(t)为第t次迭代的天牛右天线空间坐标，x_l(t)为第t次迭代天牛左天线的空间坐标，d(t)表示第t次迭代左右天线之间的距离；

step(t)＝eta·step(t-1) (10)

d(t)＝step(t)/c (11)

其中，c为常数，需要根据实际问题设置；

采用自适应的衰减因子，如式(12)所示：

其中，t表示当前迭代步数，T表示最大迭代步数；

3)、初始化自适应天牛群优化算法的迭代次数，种群数NP，加速常数c1和c2，惯性权重最大值ω_max，惯性权重最小值ω_min，常数λ，初始步长step(1)，天线距离衰减因子c，衰减因子初值eta(1)，以及每只天牛位置范围；以模型输出电压数据与实际采集数据的均方根误差为目标函数，对Duhem逆模型参数α，β，γ，p_i和q_i，以及ARX模型参数进行辨识。

5.根据权利要求1所述的基于广义Bouc-Wen逆模型的压电前馈补偿方法，其特征在于，所述步骤4)中所述的将建立Hammerstein率相关逆迟滞模型与压电执行器级联构成整体线性化系统来对迟滞非线性进行补偿的方法为：基于Hammerstein率相关逆迟滞模型建立了压电执行器动态迟滞前馈补偿器，目标位移通过该前馈补偿器的解算得到压电执行器驱动电压，经过数模转化模块与压电驱动模块对压电执行器进行控制。