一种基于龙格库塔法的压电陶瓷驱动器控制方法技术

本发明专利技术属于压电驱动器控制技术领域，公开一种基于龙格库塔法的压电陶瓷驱动器控制方法。该方法包括S1.采用机电耦合非线性动力系统模型对压电陶瓷驱动器进行数学建模；S2.根据所得模型确定初值，即令t＝0，计算出X(0)为初值；S3.采用四阶龙格库塔法求出λ；S4.判断四阶龙格库塔法满足绝对稳定条件，求出步长h；S5.采用四阶龙格库塔法对压电陶瓷驱动器机电耦合非线性动力系统模型求解，计算出最优控制状态变量x(t)，求出性能指标函数的最小值，得到全局最优解。通过判断选定或者自适应更换步长，使得四阶龙格库塔法绝对稳定，应用到压电陶瓷驱动器机电耦合非线性动力系统模型上，简单且控制精度高。

A piezoelectric ceramic actuator control method based on Runge Kutta method

The invention belongs to the field of piezoelectric actuator control technology, and discloses a control method of piezoelectric ceramic actuator based on Runge Kutta method. The method includes S1. Mathematical modeling of piezoelectric ceramic actuator using electromechanical coupled non-linear dynamic system model; S2. Determining the initial value according to the obtained model, i.e. t=0, calculating X(0) as the initial value; S3. Using the fourth-order Runge-Kutta method to calculate lambda; S4. Judging that the fourth-order Runge-Kutta method satisfies the absolute stability condition and calculates the step length h; The fourth-order Runge-Kutta method is used to solve the electromechanical coupling nonlinear dynamic system model of piezoelectric ceramic actuator, and the optimal control state variable x(t) is calculated. The minimum value of performance index function is obtained, and the global optimal solution is obtained. The fourth-order Runge-Kutta method is absolutely stable by judging the selected or self-adapting step size. It is applied to the electromechanical coupling nonlinear dynamic system model of piezoelectric ceramic actuator with simple and high control accuracy.

【技术实现步骤摘要】
一种基于龙格库塔法的压电陶瓷驱动器控制方法
本专利技术属于压电驱动器控制
，更具体地，涉及一种基于龙格库塔法的压电陶瓷驱动器控制方法。
技术介绍
压电陶瓷驱动器是一种利用压电逆压电效应制作的微/纳米级的超精密定位装置，具有位移分辨力高、响应快、无噪音、体积小、刚性大和质量轻等优点，所以被广泛用于精密定位技术、超声马达技术和纳米技术等领域，但是压电陶瓷驱动器的迟滞性是影响其精度的关键因素，其迟滞性不仅会降低其系统的控制精度，还会导致其系统的不稳定。对压电陶瓷驱动器迟滞性进行建模有很多学者，典型的有Wiener模型，Preisach模型和Prandtl-Ishlinslii算子迟滞蠕变模型等，因为对压电陶瓷驱动器迟滞性模型和控制算法的研究具有重要的实践意义，本专利在压电陶瓷驱动器机电耦合非线性动力系统建模基础下，利用龙格库塔法最优控制方法对该模型进行控制求解。对于在压电陶瓷驱动器机电耦合非线性动力系统模型下，目前还没有一种稳定的最优控制方法对其控制求解。
技术实现思路
为了解决上述现有技术存在的不足和缺点，提供一种基于龙格库塔法的压电陶瓷驱动器控制方法。在压电陶瓷驱动器机电耦合非线性动力系统模型下，用四阶龙格库塔控制方法对其求解，通过自适应选定步长，使得四阶龙格库塔法在压电陶瓷驱动器机电耦合非线性动力系统模型计算中绝对稳定，从而提高压电陶瓷驱动器的定位精度。本专利技术的目的通过下述技术方案来实现：一种基于龙格库塔法的压电陶瓷驱动器控制方法，包括如下具体步骤：S1.采用机电耦合非线性动力系统模型对压电陶瓷驱动器进行数学建模，得到压电陶瓷驱动器的模型如下：性能指标函数为：约束条件为：其中，Q、R、X、F、G、U表示矩阵，为X的导数，x(t)为状态变量，u(t)为控制变量，xdes为期望位移，x为实际输出位移，Q为压电陶瓷驱动器的输出总电荷，U为输入电压，R，G为权矩阵，F为与状态变量与控制变量有关的函数矩阵。S2.根据步骤S1所得模型确定初值，即令t＝0，计算出初始状态数值X(0)；S3.采用四阶龙格库塔法将X(0)代入式(3)中求出λ，所述四阶龙格库塔法计算公式如下：其中，h为状态变量离散化后的步长，ym为第m次迭代由式(3)得到的函数值，xm表示第m次迭代状态变量分量的数值，由K1＝f(x0，y0)＝λX(0)，得λ为初始状态数值X(0)与K1的比值；S4.将步骤S3求出的λ代入式(4)中判断四阶龙格库塔法满足绝对稳定条件，求出步长h；所述四阶龙格库塔法满足绝对稳定的条件为：其中，h为状态变量离散化后的步长，如果首次选定的步长h满足式(4)，则四阶龙格库塔法绝对稳定；由式(4)求出四阶龙格库塔法的绝对稳定区间为(-2.785，0)，如果首次选定的步长h不满足式(4)，则自适应的更换h，使得-2.785＜λh＜0，以使四阶龙格库塔法满足绝对稳定，通过以下方式更换步长h：若λ为实数，则步长h满足若λ不为实数，则步长h满足-2.785＜λh＜0；S5.采用四阶龙格库塔法对压电陶瓷驱动器机电耦合非线性动力系统模型求解，计算出最优控制状态变量x(t)，求出性能指标函数的最小值，得到全局最优解。进一步地，步骤S1中所述的性能指标函数的值的具体计算步骤如下：S11.把性能指标函数式(1)转化为S12.使M趋于正无穷，根据即当M趋于正无穷，求性能指标函数J；S13.令M→+∞，求出极限即为性能指标函数J的最小值。进一步地，步骤S5中所述的最优控制状态变量x(t)的具体计算步骤如下：S21.选定M＝nh，其中n≥2；S22.用四阶龙格库塔法式(3)计算状态变量x(t)，记t0＝0，更新自变量时间t，t的方程为tn+1＝tn+h；S23.判断自变量时间tn+1是否小于M，若不满足tn+1＜M，则计算最优控制状态变量x(t)，若满足tn+1＜M，则回到步骤S22更新自变量时间t，重复步骤S22使得tn+1＝M；S24.计算出状态变量x(t)关于M的值，通过x(t)计算出控制变量u(t)的值；S25.把状态变量x(t)和控制变量u(t)代入性能指标函数式(1)中，求出性能指标函数的最小值。与现有技术相比，本专利技术具有以下有益效果：本专利技术基于龙格库塔法的压电陶瓷驱动器控制方法，通过判定四阶龙格库塔法是否绝对稳定，从而自适应地确定龙格库塔法离散化连续状态变量的步长大小，并且把该离散技术应用到基于机电耦合非线性动力系统模型的压电陶瓷驱动器控制方法上，该方法计算精度较高，提高压电驱动系统的位移跟踪精度，可有效抑制建模误差带来的大幅值超现象。附图说明图1为本专利技术基于龙格库塔法的压电陶瓷驱动器控制方法的流程图。具体实施方式下面结合具体实施例进一步说明本专利技术的内容，但不应理解为对本专利技术的限制。图1为本专利技术基于龙格库塔法的压电陶瓷驱动器控制方法的流程图。如图1所示，一种基于龙格库塔法的压电陶瓷驱动器控制方法，包括如下具体步骤：S1.采用机电耦合非线性动力系统模型对压电陶瓷驱动器进行数学建模，得到压电陶瓷驱动器的模型如下：性能指标函数为：约束条件为：其中，Q、R、X、F、G、U表示矩阵，为X的导数，x(t)为状态变量，u(t)为控制变量，xdes为期望位移，x为实际输出位移，Q为压电陶瓷驱动器的输出总电荷，U为输入电压，R，G为权矩阵，F为与状态变量与控制变量有关的函数矩阵。S2.根据步骤S1所得模型确定初值，即令t＝0，计算出初始状态数值X(0)；S3.采用四阶龙格库塔法将X(0)代入式(3)中求出λ，所述四阶龙格库塔法计算公式如下：K1＝f(xm，ym)K4＝f(xM+h，ym+hK3)其中，h为状态变量离散化后的步长，ym为第m次迭代由式(3)得到的函数值，xm表示第m次迭代状态变量分量的数值，由K1＝f(x0，y0)＝λX(0)，得λ为初始状态数值X(0)与K1的比值；S4.将步骤S3求出的λ代入式(4)中判断四阶龙格库塔法满足绝对稳定条件，求出步长h；所述四阶龙格库塔法满足绝对稳定的条件为：其中，h为状态变量离散化后的步长，如果首次选定的步长h满足式(4)，则四阶龙格库塔法绝对稳定；由式(4)求出四阶龙格库塔法的绝对稳定区间为(-2.785，0)，如果首次选定的步长h不满足式(4)，则自适应的更换h，使得-2.785＜λh＜0，以使四阶龙格库塔法满足绝对稳定，通过以下方式更换步长h：若λ为实数，则步长h满足若λ不为实数，则步长h满足-2.785＜λh＜0；S5.采用四阶龙格库塔法对压电陶瓷驱动器机电耦合非线性动力系统模型求解，计算出最优控制状态变量x(t)，求出性能指标函数的最小值，得到全局最优解。具体地，步骤S1中所述的性能指标函数的值的具体计算步骤如下：S11.把性能指标函数式(1)转化为S12.使M趋于正无穷，根据即当M趋于正无穷，求性能指标函数J；S13.令M→+∞，求出极限即为性能指标函数J的最小值。具体地，步骤S5中所述的最优控制状态变量x(t)的具体计算步骤如下：S21.选定M＝nh，其中n≥2；S22.用四阶龙格库塔法式(3)计算状态变量x(t)，记t0＝0，更新自变量时间t，t的方程为tn+1＝tn+h；S23.判断自变量时间tn+1是否小于M，若不满足tn+1＜M，则计算最优控制状态变量x(t)，若满足tn+1本文档来自技高网...

【技术保护点】
1.一种基于龙格库塔法的压电陶瓷驱动器控制方法，其特征在于，包括如下具体步骤：S1.采用机电耦合非线性动力系统模型对压电陶瓷驱动器进行数学建模，得到压电陶瓷驱动器的模型如下：性能指标函数为：

【技术特征摘要】
1.一种基于龙格库塔法的压电陶瓷驱动器控制方法，其特征在于，包括如下具体步骤：S1.采用机电耦合非线性动力系统模型对压电陶瓷驱动器进行数学建模，得到压电陶瓷驱动器的模型如下：性能指标函数为：约束条件为：其中，Q、R、X、F、G、U表示矩阵，为X的导数，x(t)为状态变量，u(t)为控制变量，xdes为期望位移，x为实际输出位移，Q为压电陶瓷驱动器的输出总电荷，U为输入电压，R和G为权矩阵，F为与状态变量与控制变量有关的函数矩阵；S2.根据步骤S1所得模型确定初值，即令t＝0，计算出初始状态数值X(0)；S3.采用四阶龙格库塔法将X(0)代入式(3)中求出λ，所述四阶龙格库塔法计算公式如下：其中，h为状态变量离散化后的步长，xm表示第m次迭代状态变量分量的数值，ym为第m次迭代由式(3)得到的函数值，由K1＝f(x0，y0)＝λX(0)，得λ为初始状态数值X(0)与K1的比值；S4.将步骤S3求出的λ代入式(4)中判断四阶龙格库塔法是否满足绝对稳定条件，求出步长h；所述四阶龙格库塔法满足绝对稳定的条件为：其中，h为状态变量离散化后的步长，如果首次选定的步长h满足式(4)，则四阶龙格库塔法绝对稳定；由式(4)求出四阶龙格库塔法的绝对稳定区间为(-2.785，0)，如果首次选定的步长h不满足式(4)，则自适应的更换h，使得-2.785＜λh＜0，以使四阶龙格库塔法满足绝...

【专利技术属性】
技术研发人员：麦恒，陈学松，张丽丽，蔡述庭，
申请(专利权)人：广东工业大学，
类型：发明
国别省市：广东,44