基于最小二乘支持向量机的机电伺服系统摩擦补偿和动态面控制方法技术方案

一种基于最小二乘支持向量机的机电伺服系统摩擦补偿和动态面控制方法，包括：建立机电伺服系统模型和LuGre摩擦模型，初始化系统状态以及相关控制参数；通过最小二乘支持向量机逼近摩擦力，并补偿到系统中。根据逼近的摩擦力设计动态面控制器，保证系统跟踪误差快速稳定地收敛至零点，实现机电伺服系统的快速稳定控制。本发明专利技术提出一种基于最小二乘支持向量机的机电伺服系统摩擦补偿和动态面控制方法，解决摩擦不确定性的问题，利用最小二乘支持向量机的函数逼近功能对摩擦进行准确建模，进而根据建立的摩擦模型，结合动态面控制法设计自适应控制器。由于系统中的摩擦已经得到了补偿，系统状态能更好的收敛。

【技术实现步骤摘要】
基于最小二乘支持向量机的机电伺服系统摩擦补偿和动态面控制方法
本专利技术涉及一种基于最小二乘支持向量机的机电伺服系统摩擦补偿和动态面控制方法。特别是带有系统部分状态不可测、参数不确定以及外部扰动的机电伺服系统动态面控制方法。
技术介绍
机电伺服系统是以电动机作为动力驱动元件的伺服系统，广泛应用于飞行控制、火力控制等各种领域。但是，系统中的摩擦会影响伺服系统的控制精度，甚至严重降低机电伺服系统的性能，并且摩擦力的表现形式较为复杂，不易建模。因此，如何有效地控制和消除摩擦的不利影响，已成为机电控制中亟待解决的关键问题之一。动态面控制方法(DynamicSurfaceControl)在实现不确定非线性系统(特别是当干扰或不确定性不满足匹配条件时)的鲁棒控制或自适应控制方面有着明显的优越性，受到国内外学者的极大关注。利用系统的结构特性递推地构造整个系统的Lyapunov函数，使得控制器的设计结构化、系统化。针对反步控制法(BacksteppingControl)中可能导致微分项膨胀的问题，动态面控制方法在每一步设计中引入一阶积分滤波器，使得每一步设计的虚拟控制输入通过该滤波器，从而避免了系统中一些非线性函数的微分计算和控制器微分项的膨胀问题，具有较好的控制性能。
技术实现思路
为了克服现有技术的系统部分状态参数不确定的缺点，并消除摩擦对机电伺服系统性能的影响，本专利技术提出一种基于最小二乘支持向量机的机电伺服系统摩擦补偿和动态面控制方法，解决摩擦不确定性的问题，利用最小二乘支持向量机的函数逼近功能对摩擦进行准确建模，进而根据建立的摩擦模型，结合动态面控制法设计自适应控制器。由于系统中的摩擦已经得到了补偿，系统状态能更好的收敛。为了解决上述技术问题提出的技术方案为：一种基于最小二乘支持向量机的机电伺服系统摩擦补偿和动态面控制方法，其特征在于：包括如下步骤：步骤1，建立如式(1)所示的机电伺服系统模型，初始化系统状态以及控制参数；其中，θm，ωm为状态变量，分别表示电机输出轴位置和转速；J和D是折算到电机轴上的等效转动惯量和等效阻尼系数；Kt是电机扭矩常数；u是控制量；Tl是折算到电机轴上的负载扭矩；TN是折算到电机轴上的摩擦力；步骤2，建立非线性摩擦的LuGre模型，过程如下：2.1，对于摩擦采用LuGre模型：其中，σ0为鬃毛刚度系数，σ1为鬃毛阻尼系数，σ2为粘滞摩擦系数，z为接触表面鬃毛的平均变形量；2.2，将式(2)做如下分析：其中是非线性函数，i是电机转速ωm与负载转速ωl的传动比，FS表示最大静摩擦力矩，FC表示Coulomb摩擦力，ωs表示Stribeck速度；步骤3，应用最小二乘支持向量机逼近摩擦力TN，过程如下：3.1，通过试验得到TN的样本数据；根据伺服系统的实际工作条件，确定正转和反转的转动速度范围是[0,ωmax]和[ωmin,0]，ωmax表示正转最大速度，ωmin表示反转最大速度；令控制器输入为0，在速度范围内确定摩擦力矩的数据样本：其中，ω1...ωP为正转的转速采样样本数据，为反转的转速的采样样本数据，T1...TP和为所对应的转速下的摩擦力大小，P表示正转采样样本数据容量，Q表示反转采样样本数据容量；3.2，利用样本数据和最小二乘支持向量机建立模型：fsvrm(ωm)＝αφ(ωm)+b(4)其中，A＝φ(ωm)+1v/γ，δ＝P+Q，1v＝[1；1；...；1]，c和γ是大于0的常数，α，b为模型参数；式(1)中的TN用式(4)得到的fsvrm(ωm)来代替，则式(1)改写为：步骤4，用动态面的方法来设计控制器u，过程如下；4.1，定义跟踪误差：e1＝θref-θm(6)其中，e1为跟踪误差，θref为跟踪参考信号，θm为电机输出轴位置；对式(6)求导，得：其中，为跟踪误差的导数，为跟踪参考信号的导数，为电机输出轴位置的导数；4.2，定义虚拟控制器S及它的导数有如下关系：其中，k1，k为大于0的常数，是位置跟踪误差的积分作用，τ＞0表示时间常数；定义中间误差e2：e2＝S-ωm(9)将式(9)代入式(7)得到：对式(9)求导得：将式(4)，式(8)代入式(11)得：4.3，设计控制器u：其中k2为正的常数；4.4，设计李雅普诺夫函数：将式(6)，(10)，(12)，(13)代入式(14)，如果判定系统是稳定的。本专利技术结合最小二乘支持向量机和动态面控制技术，设计基于最小二乘支持向量机的动态面控制器，实现机电伺服系统的摩擦补偿和位置精确跟踪控制。本专利技术的技术构思为：机电伺服系统中由于存在摩擦力而导致控制精度不高。针对部分状态不可测、参数不确定以及存在外部扰动的机电伺服系统，运用最小二乘支持向量机，结合动态面控制理论，设计一种基于最小二乘支持向量机的机电伺服系统动态面控制方法，尽可能地消除了摩擦力对系统控制的影响。通过最小二乘支持向量机估计摩擦力，并补偿到系统中。设计动态面控制器保证系统跟踪误差快速稳定地收敛至零点，实现机电伺服系统的快速稳定控制。本专利技术的优点为：能精确估计并且补偿摩擦力，提高系统的跟踪精度和鲁棒性，改善系统的跟踪性能。附图说明图1为本专利技术的摩擦非线性模型输出曲线；图2为本专利技术的算法的基本流程；图3为本专利技术的控制系统响应曲线；图4为本专利技术的控制系统跟踪误差；图5为本专利技术的控制信号输出；图6为本专利技术的摩擦力的估计情况。具体实施方式下面结合附图对本专利技术做进一步说明。参照图1-图6，一种基于最小二乘支持向量机的机电伺服系统摩擦补偿和动态面控制方法，包括以下步骤：步骤1，建立如式(1)所示的机电伺服系统模型，初始化系统状态以及相关控制参数；其中，θm，ωm为状态变量，分别表示电机输出轴位置和转速；J和D是折算到电机轴上的等效转动惯量和等效阻尼系数；Kt是电机扭矩常数；u是控制量；Tl是折算到电机轴上的负载扭矩；TN是折算到电机轴上的摩擦力；步骤2，建立非线性摩擦的LuGre模型；2.1，对于摩擦采用LuGre模型：其中，σ0为鬃毛刚度系数，σ1为鬃毛阻尼系数，σ2为粘滞摩擦系数，z为接触表面鬃毛的平均变形量；2.2，将式(2)做如下分析：其中是非线性函数，i是电机转速ωm与负载转速ωl的传动比，FS表示最大静摩擦力矩，FC表示Coulomb摩擦力，ωs表示Stribeck速度；步骤3，应用最小二乘支持向量机逼近摩擦力TN；3.1，通过试验得到TN的样本数据；根据伺服系统的实际工作条件，可以确定正转和反转的转动速度范围是[0,ωmax]和[ωmin,0]，ωmax＞0表示正转最大速度，ωmin＜0表示反转最大速度；令控制器输入为0，我们可以在速度范围内确定摩擦力矩的数据样本：其中，ω1...ωP为正转的转速采样样本数据，为反转的转速的采样样本数据，T1...TP和为所对应的转速下的摩擦力大小，P表示正转采样样本数据容量，Q表示反转采样样本数据容量；3.2，利用样本数据和最小二乘支持向量机建立模型：fsvrm(ωm)＝αφ(ωm)+b(4)其中，A＝φ(ωm)+1v/γ，δ＝P+Q，1v＝[1；1；...；1]，c和γ是大于0的常数，α，b为模型参数；这时式(1)中的TN可以用式(4)得到的fsvrm(ωm)来代替，则式(1)改写为：步骤4，用动态面的方法来设计控制器u；4.1本文档来自技高网...

【技术保护点】
一种基于最小二乘支持向量机的机电伺服系统摩擦补偿和动态面控制方法，其特征在于：包括如下步骤：步骤1，建立如式(1)所示的机电伺服系统模型，初始化系统状态以及控制参数；dθmdt=ωmJdωmdt=Ktu-Dωm-TN-Tl---(1)]]>其中，θm，ωm为状态变量，分别表示电机输出轴位置和转速；J和D是折算到电机轴上的等效转动惯量和等效阻尼系数；Kt是电机扭矩常数；u是控制量；Tl是折算到电机轴上的负载扭矩；TN是折算到电机轴上的摩擦力；步骤2，建立非线性摩擦的LuGre模型，过程如下：2.1，对于摩擦采用LuGre模型：TN=σ0z+σ1z.+σ2ωl---(2)]]>其中，σ0为鬃毛刚度系数，σ1为鬃毛阻尼系数，σ2为粘滞摩擦系数，z为接触表面鬃毛的平均变形量；2.2，将式(2)做如下分析：z.=ωl-|ωl|g(ωl)z---(3)]]>其中g(ωl)=FC+(FS-FC)e-(ωl/ωs)2σ0]]>是非线性函数，ωl=ωmi,]]>i是电机转速ωm与负载转速ωl的传动比，FS表示最大静摩擦力矩，FC表示Coulomb摩擦力，ωs表示Stribeck速度；步骤3，应用最小二乘支持向量机逼近摩擦力TN，过程如下：3.1，通过试验得到TN的样本数据；根据伺服系统的实际工作条件，确定正转和反转的转动速度范围是[0,ωmax]和[ωmin,0]，ωmax＞0表示正转最大速度，ωmin＜0表示反转最大速度；令控制器输入为0，在速度范围内确定摩擦力矩的数据样本：(ω1T1),...,(ωi,Ti),...,(ωP,TP)(ω‾1,T‾1),...,(ω‾j,T‾j),...,(ω‾Q,T‾Q)]]>其中，ω1...ωP为正转的转速采样样本数据，为反转的转速的采样样本数据，T1...TP和为所对应的转速下的摩擦力大小，P表示正转采样样本数据容量，Q表示反转采样样本数据容量；3.2，利用样本数据和最小二乘支持向量机建立模型：fsvrm(ωm)＝αφ(ωm)+b   (4)其中，φ(ωm)=e||ω-c||2δ2,b=(1vTAyT)/(1vTA1v),α=yT-bA,A=φ(ωm)+1v/γ,]]>ω=[ω1,ω2,...,ωP,ω‾1,ω‾2,...,ω‾N],]]>δ＝P+Q，1v＝[1；1；...；1]，c和γ是大于0的常数，α，b为模型参数；式(1)中的TN用式(4)得到的fsvrm(ωm)来代替，则式(1)改写为：dθmdt=ωmJdωmdt=Ktu-Dωm-fsvrm(ωm)-Tl---(5)]]>步骤4，用动态面的方法来设计控制器u，过程如下；4.1，定义跟踪误差：e1＝θref‑θm   (6)其中，e1为跟踪误差，θref为跟踪参考信号，θm为电机输出轴位置；对式(6)求导，得：e.1=θ.ref-θ.m=θ.ref-ωm---(7)]]>其中，为跟踪误差的导数，为跟踪参考信号的导数，为电机输出轴位置的导数；4.2，定义虚拟控制器S及它的导数有如下关系：S.τ+S=k1e1+kχ+θ.ref---(8)]]>其中，k1，k为大于0的常数，是位置跟踪误差的积分作用，τ＞0表示时间常数；定义中间误差e2：e2＝S‑ωm   (9)将式(9)代入式(7)得到：e.1=e2-k1e1-kχ---(10)]]>对式(9)求导得：e.2=S.-ω.m---(11)]]>将式(4)，式(8)代入式(11)得：e.2=k1e1+kχ+θ.ref-Sτ2-1J(Ktu-Dωm-fsvrm(ωm)+Tl)---(12)]]>4.3，设计控制器u：u=1a0(e1+S.+k2e2+a1ωm+a2fsvrm(ωm)+a2ΔT)---(13)]]>其中a0=KtJ,a1=DJ,a2=1J,]]>k2为正的常数；4.4，设...

【技术特征摘要】
1.一种基于最小二乘支持向量机的机电伺服系统摩擦补偿和动态面控制方法，其特征在于：包括如下步骤：步骤1，建立如式(1)所示的机电伺服系统模型，初始化系统状态以及控制参数；其中，θm，ωm为状态变量，分别表示电机输出轴位置和转速；J和D是折算到电机轴上的等效转动惯量和等效阻尼系数；Kt是电机扭矩常数；u是控制量；Tl是折算到电机轴上的负载扭矩；TN是折算到电机轴上的摩擦力；步骤2，建立非线性摩擦的LuGre模型，过程如下：2.1，对于摩擦采用LuGre模型：其中，σ0为鬃毛刚度系数，σ1为鬃毛阻尼系数，σ2为粘滞摩擦系数，z为接触表面鬃毛的平均变形量；2.2，将式(2)做如下分析：其中是非线性函数，i是电机转速ωm与负载转速ωl的传动比，FS表示最大静摩擦力矩，FC表示Coulomb摩擦力，ωs表示Stribeck速度；步骤3，应用最小二乘支持向量机逼近摩擦力TN，过程如下：3.1，通过试验得到TN的样本数据；根据伺服系统的实际工作条件，确定正转和反转的转动速度范围是[0,ωmax]和[ωmin,0]，ωmax＞0表示正转最大速度，ωmin＜0表示反转最大速度；令控制器输入为0，在正转和反转的转动速度范围内确定摩擦力矩的数据样本：其中，ω1...ωP为正转的转速采样样本数据，为反转的转速的采样样本数据，T1...TP和为所对应的转速下的摩擦力大小，P表示正转采样样本数据容量，Q表示反转采样样本数据容量；3.2，利用样本数据和最小二乘支持向量机建立模型：fsvrm(ωm)＝αφ(ωm)+b(4)其中，A...

【专利技术属性】
技术研发人员：陈强，楼成林，
申请(专利权)人：浙江工业大学，
类型：发明
国别省市：浙江;33