基于Duhem逆模型的压电动态迟滞前馈补偿方法

技术领域：

本发明涉及一种基于duhem逆模型的压电动态迟滞前馈补偿方法。

背景技术：
：

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

电压控制补偿方式的重点就是压电迟滞模型的建立。从数学特点上来说，压电迟滞曲线具有四点性质：1)多值映射性：相同的输入电压在电压值升高和电压值下降阶段分别对应不同的位移输出；2)非局部记忆性：压电执行器的输出位移既和输入电压的瞬时值相关，又和输出位移的历史极值有关；3)率相关特性：随着输入电压频率的不断增大,压电迟滞曲线呈现频率相关性；4)非对称性：在输入电压升高和输入电压下降两个阶段的电压-位移曲线是不对称的。因此，对于压电陶瓷迟滞曲线的建模要从这四点性质出发才能保证拟合精度。目前，压电迟滞模型的建立和改进基本都是围绕这四点性质进行的。压电迟滞模型主要有微分方程模型和数学算子模型两类。数学算子迟滞模型无法表示压电执行器的动态特征，而微分方程迟滞模型既能表征迟滞的数学特性，还能描述压电执行器的动态特性。有代表性的微分方程迟滞模型是duhem模型，其最早是由物理学家p.duhem和stefanini在十九世纪末提出的。之后，coleman和hodgdon将对duhem模型进行深入研究，得到简化的duhem模型微分方程表达式。duhem模型最大的优点在于其函数表达式简单，且输出与输入之间是速率相关的，通过调整模型参数可以准确反映不同复杂情况下压电陶瓷驱动器的迟滞非线性，符合实际应用的要求。

传统duhem模型根据weierstrass第一逼近定理，采用多项式对duhem模型中的分段连续函数进行近似，但是多项式的引入导致逆模型求解过程非常复杂，不利于实时对压电迟滞非线性进行补偿。而且，传统duhem模型描述的是固定频率下的压电迟滞曲线，无法表征压电迟滞曲线的率相关性。因此，需要针对上述问题对传统duhem模型进行改进，并与压电执行器级联构成整体线性化系统来对迟滞非线性进行补偿。

技术实现要素：
：

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

1)、根据逆函数定理，建立直接duhem逆迟滞模型，并采用多项式函数对模型中的分段连续函数进行逼近；

2)、采用hammerstein模型来描述压电率相关逆迟滞特性，以duhem逆模型表征逆迟滞非线性，以arx模型表征逆迟滞率相关特性；

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

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

具体地，所述步骤1中：

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

其中，u为驱动电压，y为位移，α是常数，f,g为分段连续函数；

2)、由weierstrass第一逼近定理可知，对于duhem逆模型中的连续函数f(·)和g(·)，在满足一定精度条件下可以应用多项式函数进行逼近，表达式为：

其中，f和g是关于y的多项式，pi,qi表示多项式f’和g的系数；

3)、对式(1)进行离散化处理，duhem逆模型离散化表达式为：

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

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

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

out(t)+a1out(t-1)+a2out(t-2)+…+anout(t-n)

＝b1in(t-d)+b2in(t-d-1)+…+bmin(t-d-m+1)+ε(t)(5)

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

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

具体地，所述步骤3中：

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

vi(t)＝ωvi(t-1)+c1r1[pbesti-xi(t-1)]+c1r1[gbest-xi(t-1)],

xi(t)＝xi(t-1)+λvi(t-1)+(1-k)δ(t-1)(6)

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

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

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

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

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

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

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

其中，xr(t)为第t次迭代的天牛右天线空间坐标，xl(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逆模型参数α,pi和qi，以及arx模型参数进行辨识。

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

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

附图说明：

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

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

图3是0.5hz频率下duhem逆迟滞曲线与实际逆迟滞曲线。

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

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

具体实施方式：

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

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

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

借助自适应天牛群优化算法辨识模型参数α,p0,p1,q1,q0。设自适应天牛群优化算法的迭代次数为300，种群数np＝120，加速常数c1＝2.8,c2＝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的单频变幅值正弦电压信号驱动下的duhem逆迟滞模型参数值如表1所示。

表1duhem逆迟滞模型参数

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

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

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

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

表2迟滞曲线拟合误差

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