CN114172436B - 基于观测器的异步电动机指令滤波离散控制方法 - Google Patents

基于观测器的异步电动机指令滤波离散控制方法 Download PDF

Info

Publication number
CN114172436B
CN114172436B CN202111506688.7A CN202111506688A CN114172436B CN 114172436 B CN114172436 B CN 114172436B CN 202111506688 A CN202111506688 A CN 202111506688A CN 114172436 B CN114172436 B CN 114172436B
Authority
CN
China
Prior art keywords
asynchronous motor
function
formula
neural network
discrete
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202111506688.7A
Other languages
English (en)
Other versions
CN114172436A (zh
Inventor
于金鹏
张加伟
徐雨梦
刘加朋
于慧慧
王栋诚
林文娟
朱一平
司晨怡
朱钰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qingdao University
Original Assignee
Qingdao University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Qingdao University filed Critical Qingdao University
Priority to CN202111506688.7A priority Critical patent/CN114172436B/zh
Publication of CN114172436A publication Critical patent/CN114172436A/zh
Application granted granted Critical
Publication of CN114172436B publication Critical patent/CN114172436B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P23/00Arrangements or methods for the control of AC motors characterised by a control method other than vector control
    • H02P23/0004Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
    • H02P23/0018Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control using neural networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P23/00Arrangements or methods for the control of AC motors characterised by a control method other than vector control
    • H02P23/12Observer control, e.g. using Luenberger observers or Kalman filters
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P25/00Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
    • H02P25/02Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details characterised by the kind of motor
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P2207/00Indexing scheme relating to controlling arrangements characterised by the type of motor
    • H02P2207/01Asynchronous machines

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Control Of Electric Motors In General (AREA)

Abstract

本发明属于异步电动机位置跟踪控制技术领域,具体公开了一种基于观测器的异步电动机指令滤波离散控制方法。该异步电动机指令滤波离散控制方法针对异步电动机容易出现的输入饱和问题,同时结合欧拉方法建立异步电动机离散系统模型;构造降维观测器估计异步电动机的转子位置、角速度;在传统的反步法中引入指令滤波控制技术,克服了反步控制算法中存在的“计算复杂性”和因果矛盾问题;利用RBF神经网络处理异步电动机离散系统中的高阶非线性项,构造了基于观测器的异步电动机指令滤波离散控制器。本发明方法能够克服输入饱和的影响,利于提高系统控制的准确性,保证系统能够快速稳定运行。

Description

基于观测器的异步电动机指令滤波离散控制方法
技术领域
本发明属于异步电动机位置跟踪控制技术领域,尤其涉及一种基于观测器的异步电动机指令滤波离散控制方法。
背景技术
异步电动机(Induction Motors,简称IMs),也叫感应电动机,是由气隙旋转磁场与转子绕组感应电流相互作用产生电磁转矩,从而实现机电能量转换为机械能量的一种交流电机。近几年来,随着信息技术,自动化技术以及微电子技术的迅速发展,异步电动机普遍应用于工农业生产中。然而,由于异步电动机系统是高阶、强耦合、非线性的多变量系统,且容易受到负载扰动、参数变化等不确定因素的影响。
因此,如何对异步电动机复杂系统进行有效控制成为一项具有挑战性的工作。
近些年,非线性系统控制方法取得了很大发展,如滑模控制、自适应控制、鲁棒控制和反步法控制等。然而,这些技术大多针对异步电机连续系统提出的,针对离散系统的控制方法相对较少。在计算机技术高速发展的时代,实际工程应用过程中大多采用离散控制方法,且离散控制在系统可实现性及稳定性上均优于连续系统。因此,针对异步电动机离散系统构建控制方法有着十分重要的现实意义。此外,一些异步电动机驱动系统的变流器母线电压的幅值受限会导致输入电压饱和问题,输入电压不稳定会使控制系统性能下降,影响电动机的调速性能。因此,需要对电压进行约束,在对异步电动机进行位置跟踪控制过程中考虑输入饱和的影响非常具有现实意义。
反步法作为非线性系统的有效控制方法之一,已经被广泛应用到异步电动机控制系统中,并取得了良好的控制效果。然而,反步法应用到离散系统中时,对虚拟控制函数进行连续求差分的过程,会引起“计算爆炸”和“因果矛盾”问题。对于一些具有高度非线性和参数不确定的高阶系统,非线性项会导致控制器设计变得十分复杂。在上述背景下,动态面技术和指令滤波技术被提出用来弥补反步法的不足。其中,指令滤波技术通过引入误差补偿机制,既能消除动态面技术产生的滤波误差,也能解决反步法产生的“计算复杂性”问题。此外,在系统实际运行过程中,诸如成本限制和传感器使用过程中产生振动引起性能下降等原因导致的系统状态不可测或测量不准确是很常见的。降维观测器以其结构简单、维数低等优点,用来估计电机系统中难以测量的状态量。
发明内容
本发明的目的在于提出一种基于观测器的异步电动机指令滤波离散控制方法,考虑异步电动机离散系统实际运行中存在的输入饱和问题,以实现电动机快速稳定的位置跟踪控制。
本发明为了实现上述目的,采用如下技术方案:
本发明提出了一种基于观测器的异步电动机指令滤波离散控制方法,该方法针对异步电动机容易出现的输入饱和问题,结合欧拉方法建立异步电动机离散系统模型;构造降维观测器估计异步电动机的转子位置、角速度;根据指令滤波技术和反步法设计基于观测器的异步电动机指令滤波离散控制器,在传统的反步法中引入指令滤波控制技术,以克服反步控制算法中存在的“计算复杂性”和因果矛盾问题,同时引入补偿信号,消除滤波误差;利用RBF神经网络处理异步电动机离散系统中的高阶非线性项,结合自适应控制方法解决系统中存在的参数未知和输入饱和问题,构造基于观测器的异步电动机指令滤波离散控制器。
本发明具有如下优点:
(1)本发明方法是针对离散系统提出的控制方法,具有较高的稳定性和可实现性。
(2)本发明方法综合考虑异步电动机离散系统的输入饱和问题,能够保证设计的控制系统在电压不稳定的情况下仍能稳定运行,有利于保护设备及人身安全,有效地解决了输入饱和情况下异步电动机离散系统的位置跟踪控制问题。
(3)本发明方法采用指令滤波技术,有效地避免了在传统反步法中存在的“计算复杂性”和因果矛盾问题,同时引入补偿信号,消除了滤波误差。
(4)本发明方法使用径向基函数(RBF)神经网络来处理电动机系统中未知的非线性项,有效地解决了异步电动机的高度非线性控制问题,最终达到更加准确的控制精度。
(5)本发明方法采用降维观测器估计转子位置、角速度,提高了系统的可靠性。
(6)本发明方法不需要根据异步电动机离散系统的不同而修改控制器的参数,原理上可以实现对所有型号和功率的异步电动机离散系统的稳定控制,在控制过程中减少了对异步电动机离散系统参数的测量,利于实现异步电动机离散系统的位置跟踪控制。
附图说明
图1为本发明实施例中基于观测器的异步电动机指令滤波离散控制器、坐标变换单元和SVPWM逆变器组成的复合被控对象的示意图;
图2是采用本发明控制方法后转子角度和转子角度设定值跟踪仿真图;
图3是采用本发明控制方法后转子角度和转子角度设定值跟踪误差仿真图;
图4是采用本发明控制方法后d轴对称饱和非线性输入ud(k)和d轴定子电压仿真图;
图5是采用本发明控制方法后q轴对称饱和非线性输入uq(k)和q轴定子电压仿真图;
图6是采用本发明控制方法后转子位置x1(k)的实际值和观测值仿真图;
图7是采用本发明控制方法后转子角速度x2(k)的实际值和观测值仿真图。
具体实施方式
下面结合附图以及具体实施方式对本发明作进一步详细说明:
图1为本发明实施例中基于观测器的异步电动机指令滤波离散控制器、坐标变换单元和SVPWM逆变器组成的复合被控对象的示意图,图中涉及到的部件有基于观测器的异步电动机指令滤波离散控制器1、坐标变换单元2、SVPWM逆变器3、转速检测单元4和电流检测单元5,ω表示转子角速度,Uα和Uβ表示两相旋转坐标系下的电压,U、V和W表示三相交流电压。其中,转速检测单元4和电流检测单元5主要用于检测异步电动机离散系统转速相关变量的电流值,通过实际测量的电流和转速变量作为输入,通过基于观测器的异步电动机指令滤波离散控制器1进行电压控制,最终转换为三相电控制异步电动机的位置跟踪控制。为了设计一个更加有效的控制器,建立异步电动机离散系统模型是十分必要的。
基于观测器的异步电动机指令滤波离散控制方法,包括如下步骤:
步骤1.建立异步电动机离散系统动态数学模型。
在同步旋转坐标系下,异步电动机离散系统动态数学模型如公式(1)所示。
Figure BDA0003404677270000031
其中,k为异步电动机离散系统的步数,
Figure BDA0003404677270000032
θ(k)、θ(k+1)分别表示异步电动机离散系统第k步、k+1步的转子角度;ω(k)、ω(k+1)分别表示异步电动机离散系统第k步、k+1步的角速度;id(k)、id(k+1)分别表示异步电动机离散系统第k步、k+1步的d轴电流;iq(k)、iq(k+1)分别表示异步电动机离散系统第k步、k+1步的q轴电流;ud(k)表示异步电动机离散系统第k步d轴的对称饱和非线性输入,uq(k)表示异步电动机离散系统第k步q轴对称饱和非线性输入;ψd(k)表示转子磁链;Δt表示异步电动机离散系统的采样周期;np表示磁极对数;Lm表示互感系数;Ls和Lr分别为定子和转子侧等效电感;TL表示异步电动机离散系统的负载转矩;J表示转动惯量;Rs和Rr分别表示定子和转子的等效电阻。
为简化以上异步电动机离散系统动态数学模型,定义如下新的变量:
Figure BDA0003404677270000041
由式(2)定义的新变量,将式(1)通过欧拉方法,得到异步电动机离散系统模型如下:
Figure BDA0003404677270000042
由于uq(k)和ud(k)的基本特性相同,为表述方便,定义u(k)代指uq(k)和ud(k)。
根据u(k)的特性,u(k)描述为:
Figure BDA0003404677270000043
其中,umax>0、umin<0是未知的输入饱和常数。
vq(k)为q轴定子的输入电压,vd(k)为d轴定子的输入电压,由于vq(k)和vd(k)的基本特性相同,为表述方便,定义v(k)代指vq(k)和vd(k)。定义如下平滑函数g(v(k)):
Figure BDA0003404677270000044
由式(5),定义u(k)=g(v(k))+Y(v(k)) (6)
式中,Y(v(k))是有界函数;Y(v(k))的边界为:
|Y(v(k))|=|u(k)-g(v(k))|≤max{umax(1-tan(1),umin(tan(1)-1)}=D。
其中,D为正常数且D>0;根据中值定理得知,存在一个常数λ,0<λ<1,使得:
Figure BDA0003404677270000051
其中,v(0)和g(v(0))分别表示定子输入电压v(k)以及平滑函数g(v(k))的初始值,
Figure BDA0003404677270000052
vλ(k)=λv(k)+(1-λ)v(0)。
那么u(k)被重新定义为:
Figure BDA0003404677270000053
其中,v(k+1)表示k+1步定子的输入电压,g(v(k+1))表示k+1步的平滑函数。
设定v(0)=0,g(v(0))=0,则式(6)改写为:
Figure BDA0003404677270000054
其中,/>
Figure BDA0003404677270000055
则异步电动机离散系统第k步q轴的对称饱和非线性输入uq(k)为:
Figure BDA0003404677270000056
则异步电动机离散系统第k步d轴的对称饱和非线性输入ud(k)为:
Figure BDA0003404677270000057
步骤2.设计降维观测器,估计异步电动机的转子位置、角速度。
对于公式(3)示出的异步电动机驱动系统动态模型(1),其动态方程如下:
Figure BDA0003404677270000058
其中,未知非线性函数f2(k)=x2(k)-x3(k)+a1Δtx3(k)x4(k)+a2ΔtTL。y(k)表示k时刻降维观测器的输出;对于f2(k),考虑用RBF神经网络技术解决。假设f(Z)在紧集ΩZ中是一个连续的函数,存在径向基函数神经网络WTS(Z)使得f(Z)=WTS(Z)+τ;τ是逼近误差且满足|τ|≤ε,ε为任意小的正常数;
Figure BDA0003404677270000059
是输入向量,q是神经网络输入维数,Rq为实数向量集;W∈Rp是权重向量,p是神经网络节点数,p为正整数且p>1,Rp为实数向量集;S(Z)=[s1(Z),...,sp(Z)]T∈Rp为基函数向量。sc(Z)为高斯函数,sc(Z)的表达式为:
Figure BDA00034046772700000510
其中,c=1,...,p,μc是高斯函数sc(Z)接受域的中心,ηc为高斯函数sc(Z)的宽度。根据径向基函数神经网络的定义,存在径向基函数神经网络
Figure BDA0003404677270000061
使得
Figure BDA0003404677270000062
W2表示权重向量,S2(Z2(k))表示基函数向量,式(10)表示为:
Figure BDA0003404677270000063
其中,
Figure BDA0003404677270000064
分别表示第k步时异步电动机的转子角度和角速度的观测值;τ2为逼近误差且满足|τ2|≤ε2,ε2为任意小的正常数。
则降维观测器如公式(12)所示;
Figure BDA0003404677270000065
其中,
Figure BDA0003404677270000066
为第k步降维观测器输出的估计值,/>
Figure BDA0003404677270000067
分别表示第k+1步时异步电动机的转子角度和角速度的观测值,g1、g2为正常数;/>
Figure BDA0003404677270000068
Φ2>0且为未知常数,/>
Figure BDA0003404677270000069
为神经网络权重/>
Figure BDA00034046772700000610
的范数;自适应律/>
Figure BDA00034046772700000611
是Φ2的估计值,定义/>
Figure BDA00034046772700000612
定义h(k)=[h1(k),h2(k)]Τ,且
Figure BDA00034046772700000613
那么h(k+1)表示为:
Figure BDA00034046772700000614
其中,
Figure BDA00034046772700000615
选择Lyapunov函数V0(k)=hΤ(k)Ph(k),其中,P为正定矩阵,PΤ=P。
对V0(k)取一阶差分得:
Figure BDA00034046772700000616
其中,||P||表示正定矩阵P的范数,||S2(Z2(k))||表示神经网络径向基函数S2(Z2(k))的范数,Q为矩阵,且Q=P-3AΤPA。
步骤3.根据指令滤波和反步法设计基于观测器的异步电动机指令滤波离散控制方法。
定义指令滤波器的公式如下:
Figure BDA0003404677270000071
式中,ξ为指令滤波器的时间常数,ωn为正常数;d1(k)、d2(k)分别为第k步时指令滤波器的输出信号,d1(k+1)、d2(k+1)分别为指令滤波器在第k+1步时的输出信号。
α(k)为指令滤波器的输入信号,即虚拟控制函数;如果虚拟控制函数α(k)满足:
Figure BDA0003404677270000072
且/>
Figure BDA0003404677270000073
对于任意的k≥1成立。
其中,
Figure BDA0003404677270000074
且d1(0)=α(0),d2(0)=0,α(0)为α(k)的初始值,d1(0)为d1(k)的初始值,d2(0)为d2(k)的初始值。
那么对于任意的Q>0,有0<ξ≤1和ωn>0,使得|d1(k)-αi(k)|≤Q,Δd1(k)=|d1(k+1)-d1(k)|有界;其中,αi(k)表示虚拟控制函数,i=1,2,3。
虚拟控制函数α(k)的表达式将在下面的控制方法设计中给出。
定义误差信号e1(k)、e2(k)、e3(k)、e4(k)、e5(k)为:
Figure BDA0003404677270000075
定义误差补偿信号ξ1(k)、ξ2(k)、ξ3(k)、ξ4(k)、ξ5(k)为:
Figure BDA0003404677270000076
其中,x1d(k)、x4d(k)为给定的期望信号;α1d(k)、α2d(k)、α3d(k)为指令滤波器的输出信号;vj(k)表示补偿误差,j=1,2,3,4,5。
对离散系统每一步都选取一个Lyapunov函数来构建虚拟控制函数,具体步骤如下:
步骤3.1.选择Lyapunov函数
Figure BDA0003404677270000077
那么V1(k)的一阶差分为:
Figure BDA0003404677270000078
其中,x1d(k+1)、ξ1(k+1)分别为k+1步的期望信号和误差补偿信号。
选取虚拟控制函数α1(k)和误差补偿信号ξ1(k)为:
Figure BDA0003404677270000081
其中,t1为常数,且|t1|≤1;根据虚拟控制函数α1(k)、误差补偿信号ξ1(k)、误差信号
Figure BDA0003404677270000082
以及误差补偿信号ξ2(k)=e2(k)-v2(k)和式(16)得:
Figure BDA0003404677270000083
步骤3.2.选择Lyapunov函数
Figure BDA0003404677270000084
那么V2(k)的一阶差分为:
Figure BDA0003404677270000085
选取虚拟控制函数α2(k)、误差补偿信号ξ2(k)和自适应律
Figure BDA0003404677270000086
为:
Figure BDA0003404677270000087
Figure BDA0003404677270000088
其中,t2为常数,且|t2|≤1,γ2和δ2为正参数,将式(19)和(20)代入式(18)中,得:
Figure BDA0003404677270000089
步骤3.3.选择Lyapunov函数
Figure BDA00034046772700000810
那么V3(k)的一阶差分为:
Figure BDA00034046772700000811
其中,f3(k)为一未知的非线性函数,且
Figure BDA00034046772700000812
由径向基函数神经网络逼近原理得知,对于给定的任意ε3>0,存在径向基函数神经网络
Figure BDA0003404677270000091
使得/>
Figure BDA0003404677270000092
W3表示权重向量,S3(Z3(k))表示基函数向量;Z3(k)=[x1(k),x2(k),x3(k),x4(k),x5(k)]Τ;τ3表示逼近误差,并满足不等式|τ3|≤ε3,||W3||是向量W3的范数。取误差补偿信号ξ3(k)=0,则异步电动机q轴定子输入电压vq(k),对称饱和非线性输入信号uq(k)和自适应律/>
Figure BDA0003404677270000093
分别为:
Figure BDA0003404677270000094
Figure BDA0003404677270000095
其中,γ3和δ3为正参数,
Figure BDA0003404677270000096
表示神经网络权重向量W3(k)的估计值,||S3(Z3(k))||表示神经网络径向基函数S3(Z3(k))的范数;定义/>
Figure BDA0003404677270000097
Φ3为未知常数,且Φ3>0,/>
Figure BDA0003404677270000098
为神经网络权重向量/>
Figure BDA0003404677270000099
的范数,自适应律/>
Figure BDA00034046772700000910
是Φ3的估计值,定义估计误差/>
Figure BDA00034046772700000911
将式(23)、(24)代入式(22)中得将式(23)、(24)代入式(22)中并结合杨氏不等式得:
Figure BDA00034046772700000912
步骤3.4.选择Lyapunov函数
Figure BDA00034046772700000913
取V4(k)的一阶差分,得到:
Figure BDA00034046772700000914
选取虚拟控制函数α3(k)和误差补偿信号ξ4(k)为:
Figure BDA00034046772700000915
其中,t4为常数,且|t4|≤1,将式(27)代入式(26)中,得:
Figure BDA00034046772700000916
步骤3.5.选择Lyapunov函数
Figure BDA00034046772700000917
其中,M为正常数。
取V5(k)的一阶差分,得到:
Figure BDA0003404677270000101
其中,f5(k)为一未知的非线性函数,且
Figure BDA0003404677270000102
由径向基函数神经网络逼近原理得知,对于给定的任意ε5>0,存在径向基函数神经网络
Figure BDA0003404677270000103
使得:/>
Figure BDA0003404677270000104
Z5(k)=[x1(k),x2(k),x3(k),x4(k),x5(k)]Τ,τ5表示逼近误差,并满足不等式|τ5|≤ε5。取误差补偿信号ξ5(k)=0则异步电动机d轴定子输入电压vd(k),对称饱和非线性输入信号ud(k)和自适应律/>
Figure BDA0003404677270000105
分别为:
Figure BDA0003404677270000106
Figure BDA0003404677270000107
其中,γ5和δ5为正参数,
Figure BDA0003404677270000108
表示神经网络权重向量W5(k)的估计值,||S5(Z5(k))||表示神经网络径向基函数S5(Z5(k))的范数,/>
Figure BDA0003404677270000109
定义/>
Figure BDA00034046772700001010
为神经网络权重向量/>
Figure BDA00034046772700001011
的范数Φ5为未知常数且Φ5>0,估计误差/>
Figure BDA00034046772700001012
是Φ5的估计值;将式(30)、(31)代入式(29)中得将式(30)、(31)代入式(29)中并结合杨氏不等式得:
Figure BDA00034046772700001013
步骤4.对构建的基于观测器的异步电动机指令滤波离散控制方法进行稳定性分析。
选择Lyapunov函数:
Figure BDA00034046772700001014
取V(k)的一阶差分,得到:
Figure BDA0003404677270000111
根据
Figure BDA0003404677270000112
和式(20)、(24)、(31)得到:
Figure BDA0003404677270000113
其中,m=2,3,5,δm为正参数,由基函数向量S(Z)的定义得知,||Sm(Zm(k))||2≤lm,lm表示径向基函数神经网络Sm(Zm(k))的节点数;由杨氏不等式,得:
Figure BDA0003404677270000114
由式(9)、(23)、(30)及杨氏不等式得:
Figure BDA0003404677270000115
将式(36)、(37)代入式(35)中得到:
Figure BDA0003404677270000116
Figure BDA0003404677270000117
Figure BDA0003404677270000118
将式(38)、(39)、(40)代入式(33)中得到:
Figure BDA0003404677270000121
其中:
Figure BDA0003404677270000122
选择合适的采样周期Δt,控制参数ξ,ωn235235,g1和g2,使下列不等式成立:
Figure BDA0003404677270000123
κ=2,3,5,且矩阵Q=P-3AΤPA正定,只需满足下列条件:
Figure BDA0003404677270000124
则有ΔV(k)≤0成立,这意味着对任意σ>0,有
Figure BDA0003404677270000125
成立。
假设
Figure BDA0003404677270000126
有界且/>
Figure BDA0003404677270000127
i3=1,2,3,i4=1,2,4,那么补偿信号/>
Figure BDA0003404677270000128
最终一致有界,i1=1,2,3,4,5,由v1(k)=e1(k)-ξ1(k)和ξ1(k)是有界的,则跟踪误差e1(k)亦是有界的。
由以上分析得到,在定子的输入电压vq(k)和vd(k)的作用下,异步电动机离散系统的误差变量e1(k)能够收敛到原点的一个充分小的邻域内,并保证其他信号有界。
下面在Matlab环境下对所建立的基于观测器的异步电动机指令滤波离散控制方法进行仿真,以验证所提出控制方法的可行性。
电动机及负载参数为:
J=0.0586Kg·m2,Rs=0.1Ω,Rr=0.15Ω,Ls=Lr=0.0699H,Lm=0.068H。
跟踪参考信号为:x1d(k)=cos(Δtkπ/2);期望转子磁链信号为:x4d(k)=1。
采样周期:Δt=0.0025s,负载转矩为:
Figure BDA0003404677270000131
选择的控制参数为:ξ=0.98,ωn=500,γ2=2.88,γ3=0.1,γ5=0.1,δ2=1.0,δ3=0.1,δ5=0.1,t1=t2=t4=0.1,g1=0.0025,g2=300。
选择RBF神经网络如下:神经网络
Figure BDA0003404677270000132
和/>
Figure BDA0003404677270000133
包含11个中心平均分布在在[-10,10]内节点,宽度均为2。
相应的仿真结果如图2至图7所示。其中:
图2是采用本发明控制方法后转子角度和转子角度设定值跟踪仿真图。
图3是采用本发明控制方法后转子角度和转子角度设定值跟踪误差仿真图。
由图2和图3的仿真结果表明,由图2和图3的仿真结果表明,本发明方法能够准确跟踪异步电动机转子位置,具有很小的跟踪误差,并且具有较强的抗扰动能力。
图4是采用本发明控制方法后d轴对称饱和非线性输入ud(k)和d轴定子电压仿真图。
图5是采用本发明控制方法后q轴对称饱和非线性输入uq(k)和q轴定子电压仿真图。
由图4和图5的仿真结果表明,本发明方法能够有效的减少输入饱和带来的不利影响,位置跟踪误差小、跟踪效果好。
图6是采用本发明控制方法后转子位置x1(k)的实际值和观测值仿真图。
图7是采用本发明控制方法后转子角速度x2(k)的实际值和观测值仿真图。
由图6和图7的仿真结果表明,采用本发明方法控制后降维观测器能够有效地估计转子位置及角速度。
综上所述,本发明提出的基于观测器的异步电动机指令滤波离散控制方法,能够在考虑输入饱和的情况下,快速稳定地跟踪参考信号。
当然,以上说明仅仅为本发明的较佳实施例,本发明并不限于列举上述实施例,应当说明的是,任何熟悉本领域的技术人员在本说明书的教导下,所做出的所有等同替代、明显变形形式,均落在本说明书的实质范围之内,理应受到本发明的保护。

Claims (1)

1.基于观测器的异步电动机指令滤波离散控制方法,其特征在于,包括如下步骤:
步骤1.建立异步电动机离散系统动态数学模型;
在同步旋转坐标系下,异步电动机离散系统动态数学模型如公式(1)所示;
Figure FDA0003404677260000011
其中,k为异步电动机离散系统的步数,
Figure FDA0003404677260000012
θ(k)、θ(k+1)分别表示异步电动机离散系统第k步、k+1步的转子角度;ω(k)、ω(k+1)分别表示异步电动机离散系统第k步、k+1步的角速度;id(k)、id(k+1)分别表示异步电动机离散系统第k步、k+1步的d轴电流;iq(k)、iq(k+1)分别表示异步电动机离散系统第k步、k+1步的q轴电流;ud(k)表示异步电动机离散系统第k步d轴的对称饱和非线性输入,uq(k)表示异步电动机离散系统第k步q轴对称饱和非线性输入;ψd(k)表示转子磁链;Δt表示异步电动机离散系统的采样周期;np表示磁极对数;Lm表示互感系数;Ls和Lr分别为定子和转子侧等效电感;TL表示异步电动机离散系统的负载转矩;J表示转动惯量;Rs和Rr分别表示定子和转子的等效电阻;
为简化以上异步电动机离散系统动态数学模型,定义如下新的变量:
Figure FDA0003404677260000013
由式(2)定义的新变量,将式(1)通过欧拉方法,得到异步电动机离散系统模型如下:
Figure FDA0003404677260000021
由于uq(k)和ud(k)的基本特性相同,为表述方便,定义u(k)代指uq(k)和ud(k);
根据u(k)的特性,u(k)描述为:
Figure FDA0003404677260000022
其中,umax>0、umin<0是未知的输入饱和常数;vq(k)和vd(k)均为定子的输入电压,由于vq(k)和vd(k)的基本特性相同,为表述方便,定义v(k)代指vq(k)和vd(k);
定义如下平滑函数g(v(k)):
Figure FDA0003404677260000023
由式(5),定义u(k)=g(v(k))+Y(v(k)) (6)
式中,Y(v(k))是有界函数;Y(v(k))的边界为:
|Y(v(k))|=|u(k)-g(v(k))|≤max{umax(1-tan(1),umin(tan(1)-1)}=D;
其中,D为正常数且D>0;根据中值定理得知,存在一个常数λ,0<λ<1,使得:
Figure FDA0003404677260000024
其中,v(0)和g(v(0))分别表示定子输入电压v(k)以及平滑函数g(v(k))的初始值,
Figure FDA0003404677260000025
vλ(k)=λv(k)+(1-λ)v(0);
那么u(k)被重新定义为:
Figure FDA0003404677260000026
其中,v(k+1)表示k+1步定子的输入电压,g(v(k+1))表示k+1步的平滑函数;设定v(0)=0,g(v(0))=0,则式(6)改写为:
Figure FDA0003404677260000027
其中,/>
Figure FDA0003404677260000028
则异步电动机离散系统第k步q轴的对称饱和非线性输入uq(k)为:
Figure FDA0003404677260000029
则异步电动机离散系统第k步d轴的对称饱和非线性输入ud(k)为:
Figure FDA0003404677260000031
步骤2.设计降维观测器,估计异步电动机的转子位置和角速度;
对于公式(3)给出的异步电动机离散系统模型,其动态方程如下:
Figure FDA0003404677260000032
其中,未知非线性函数f2(k)=x2(k)-x3(k)+a1Δtx3(k)x4(k)+a2ΔtTL;y(k)表示k时刻降维观测器的输出;假设f(Z)在紧集ΩZ中是一个连续的函数,存在径向基函数神经网络WTS(Z)使得f(Z)=WTS(Z)+τ;τ是逼近误差且满足|τ|≤ε,ε为任意小的正常数;
Figure FDA0003404677260000033
是输入向量,q是神经网络输入维数,Rq为实数向量集;W∈Rp是权重向量,p是神经网络节点数,p为正整数且p>1,Rp为实数向量集;S(Z)=[s1(Z),...,sp(Z)]T∈Rp为基函数向量;sc(Z)为高斯函数,sc(Z)的表达式为:
Figure FDA0003404677260000034
其中,c=1,...,p,μc是高斯函数sc(Z)接受域的中心,ηc为高斯函数sc(Z)的宽度;
根据径向基函数神经网络的定义,存在径向基函数神经网络
Figure FDA0003404677260000035
使得
Figure FDA0003404677260000036
W2表示权重向量,S2(Z2(k))表示基函数向量,式(10)表示为:
Figure FDA0003404677260000037
其中,
Figure FDA0003404677260000038
Figure FDA0003404677260000039
和/>
Figure FDA00034046772600000310
分别表示第k步时异步电动机的转子角度和角速度的观测值;τ2为逼近误差且满足|τ2|≤ε2,ε2为任意小的正常数;
则降维观测器如公式(12)所示;
Figure FDA00034046772600000311
其中,
Figure FDA0003404677260000041
为第k步降维观测器输出的估计值,/>
Figure FDA0003404677260000042
分别表示第k+1步时异步电动机的转子角度和角速度的观测值,g1、g2为正常数;/>
Figure FDA0003404677260000043
Φ2>0且为未知常数,/>
Figure FDA0003404677260000044
为神经网络权重/>
Figure FDA0003404677260000045
的范数;
自适应律
Figure FDA0003404677260000046
是Φ2的估计值,定义/>
Figure FDA0003404677260000047
定义h(k)=[h1(k),h2(k)]T,且
Figure FDA0003404677260000048
那么h(k+1)表示为:
Figure FDA0003404677260000049
其中,
Figure FDA00034046772600000417
选择Lyapunov函数V0(k)=hT(k)Ph(k),其中,P为正定矩阵,PT=P;
对V0(k)取一阶差分得:
Figure FDA00034046772600000411
其中,||P||表示正定矩阵P的范数,||S2(Z2(k))||表示神经网络径向基函数S2(Z2(k))的范数,Q为矩阵,且Q=P-3ATPA;
步骤3.根据指令滤波和反步法设计基于观测器的异步电动机指令滤波离散控制方法;
定义指令滤波器的公式如下:
Figure FDA00034046772600000412
式中,ξ为指令滤波器的时间常数,ωn为正常数;d1(k)、d2(k)分别为第k步时指令滤波器的输出信号,d1(k+1)、d2(k+1)分别为指令滤波器在第k+1步时的输出信号;
α(k)为指令滤波器的输入信号,即虚拟控制函数;如果虚拟控制函数α(k)满足:
Figure FDA00034046772600000413
且/>
Figure FDA00034046772600000414
对于任意的k≥1成立;其中,/>
Figure FDA00034046772600000415
Figure FDA00034046772600000416
且d1(0)=α(0),d2(0)=0,α(0)为α(k)的初始值,d1(0)为d1(k)的初始值,d2(0)为d2(k)的初始值;那么对于任意的Q>0,有0<ξ≤1和ωn>0,使得|d1(k)-αi(k)|≤Q,Δd1(k)=|d1(k+1)-d1(k)|有界;其中,αi(k)表示虚拟控制函数,i=1,2,3;
定义误差信号e1(k)、e2(k)、e3(k)、e4(k)、e5(k)为:
Figure FDA0003404677260000051
定义误差补偿信号ξ1(k)、ξ2(k)、ξ3(k)、ξ4(k)、ξ5(k)为:
Figure FDA0003404677260000052
其中,x1d(k)、x4d(k)为给定的期望信号;α1d(k)、α2d(k)、α3d(k)为指令滤波器的输出信号;vj(k)表示补偿误差,j=1,2,3,4,5;
对离散系统每一步都选取一个Lyapunov函数来构建虚拟控制函数,具体步骤如下:
步骤3.1.选择Lyapunov函数
Figure FDA0003404677260000053
那么V1(k)的一阶差分为:
Figure FDA0003404677260000054
其中,x1d(k+1)、ξ1(k+1)分别为k+1步的期望信号和误差补偿信号;
选取虚拟控制函数α1(k)和误差补偿信号ξ1(k)为:
Figure FDA0003404677260000055
其中,t1为常数,且|t1|≤1;根据虚拟控制函数α1(k)、误差补偿信号ξ1(k)、误差信号
Figure FDA0003404677260000056
以及误差补偿信号ξ2(k)=e2(k)-v2(k)和式(16)得:
Figure FDA0003404677260000057
步骤3.2.选择Lyapunov函数
Figure FDA0003404677260000058
那么V2(k)的一阶差分为:
Figure FDA0003404677260000061
选取虚拟控制函数α2(k)、误差补偿信号ξ2(k)和自适应律
Figure FDA0003404677260000062
为:
Figure FDA0003404677260000063
Figure FDA0003404677260000064
其中,t2为常数,且|t2|≤1,γ2和δ2为正参数,将式(19)和(20)代入式(18)中,得:
Figure FDA0003404677260000065
步骤3.3.选择Lyapunov函数
Figure FDA0003404677260000066
那么V3(k)的一阶差分为:
Figure FDA0003404677260000067
其中,f3(k)为一未知的非线性函数,且
Figure FDA0003404677260000068
由径向基函数神经网络逼近原理得知,对于给定的任意ε3>0,存在径向基函数神经网络W3 TS3(Z3(k)),使得f3(k)=W3 TS3(Z3(k))+τ3;W3表示权重向量,S3(Z3(k))表示基函数向量;Z3(k)=[x1(k),x2(k),x3(k),x4(k),x5(k)]T;τ3表示逼近误差,并满足不等式|τ3|≤ε3,||W3||是向量W3的范数;取误差补偿信号ξ3(k)=0,则异步电动机q轴定子输入电压vq(k),对称饱和非线性输入uq(k)和自适应律
Figure FDA0003404677260000069
分别为:
Figure FDA00034046772600000610
Figure FDA00034046772600000611
其中,γ3和δ3为正参数,
Figure FDA00034046772600000612
Figure FDA00034046772600000613
表示神经网络权重向量W3(k)的估计值,||S3(Z3(k))||表示神经网络径向基函数S3(Z3(k))的范数;
定义Φ3=||W3 T||,Φ3为未知常数,且Φ3>0,||W3 T||为神经网络权重向量W3 T的范数,自适应律
Figure FDA0003404677260000071
是Φ3的估计值,定义估计误差/>
Figure FDA0003404677260000072
将式(23)、(24)代入式(22)中得将式(23)、(24)代入式(22)中并结合杨氏不等式得:
Figure FDA0003404677260000073
步骤3.4.选择Lyapunov函数
Figure FDA0003404677260000074
取V4(k)的一阶差分,得到:
Figure FDA0003404677260000075
选取虚拟控制函数α3(k)和误差补偿信号ξ4(k)为:
Figure FDA0003404677260000076
其中,t4为常数,且|t4|≤1,将式(27)代入式(26)中,得:
Figure FDA0003404677260000077
步骤3.5.选择Lyapunov函数
Figure FDA0003404677260000078
其中,M为正常数;
取V5(k)的一阶差分,得到:
Figure FDA0003404677260000079
其中,f5(k)为一未知的非线性函数,且
Figure FDA00034046772600000710
由径向基函数神经网络逼近原理得知,对于给定的任意ε5>0,存在径向基函数神经网络W5 TS5(Z5(k)),使得:f5(k)=W5 T||S5(Z5(k))||+τ5;Z5(k)=[x1(k),x2(k),x3(k),x4(k),x5(k)]T,τ5表示逼近误差,并满足不等式|τ5|≤ε5;取误差补偿信号ξ5(k)=0,则异步电动机d轴定子输入电压vd(k),对称饱和非线性输入ud(k)和自适应律
Figure FDA00034046772600000711
分别为:
Figure FDA0003404677260000081
Figure FDA0003404677260000082
其中,γ5和δ5为正参数,
Figure FDA0003404677260000083
表示神经网络权重向量W5(k)的估计值,||S5(Z5(k))||表示神经网络径向基函数S5(Z5(k))的范数,/>
Figure FDA0003404677260000084
定义Φ5=||W5 T||,||W5 T||为神经网络权重向量W5 T的范数Φ5为未知常数且Φ5>0,估计误差/>
Figure FDA0003404677260000085
Figure FDA0003404677260000086
是Φ5的估计值;将式(30)、(31)代入式(29)中得将式(30)、(31)代入式(29)中并结合杨氏不等式得:
Figure FDA0003404677260000087
步骤4.对构建的基于观测器的异步电动机指令滤波离散控制方法进行稳定性分析;
选择Lyapunov函数:
Figure FDA0003404677260000088
取V(k)的一阶差分,得到:
Figure FDA0003404677260000089
根据
Figure FDA00034046772600000810
和式(20)、(24)、(31)得到:
Figure FDA00034046772600000811
其中,m=2,3,5,δm为正参数,由基函数向量S(Z)的定义得知,||Sm(Zm(k))||2≤lm,lm表示径向基函数神经网络Sm(Zm(k))的节点数;由杨氏不等式,得:
Figure FDA0003404677260000091
由式(9)、(23)、(30)及杨氏不等式得:
Figure FDA0003404677260000092
将式(36)、(37)代入式(35)中得到:
Figure FDA0003404677260000093
Figure FDA0003404677260000094
Figure FDA0003404677260000095
将式(38)、(39)、(40)代入式(33)中得到:
Figure FDA0003404677260000096
其中:
Figure FDA0003404677260000097
选择合适的采样周期Δt,控制参数ξ,ωn235235,g1和g2,使下列不等式成立:
Figure FDA0003404677260000101
其中,κ=2,3,5,且矩阵Q=P-3ATPA正定,只需满足下列条件:
Figure FDA0003404677260000102
则有ΔV(k)≤0成立,这意味着对任意σ>0,有
Figure FDA0003404677260000103
成立;
假设
Figure FDA0003404677260000104
有界且/>
Figure FDA0003404677260000105
i3=1,2,3,i4=1,2,4,那么补偿信号/>
Figure FDA0003404677260000106
最终一致有界,i1=1,2,3,4,5,由v1(k)=e1(k)-ξ1(k)和ξ1(k)是有界的,则跟踪误差e1(k)亦是有界的;
由以上分析得到,在定子的输入电压vq(k)和vd(k)的作用下,异步电动机离散系统的误差变量e1(k)能够收敛到原点的一个充分小的邻域内,并保证其他信号有界。
CN202111506688.7A 2021-12-10 2021-12-10 基于观测器的异步电动机指令滤波离散控制方法 Active CN114172436B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111506688.7A CN114172436B (zh) 2021-12-10 2021-12-10 基于观测器的异步电动机指令滤波离散控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111506688.7A CN114172436B (zh) 2021-12-10 2021-12-10 基于观测器的异步电动机指令滤波离散控制方法

Publications (2)

Publication Number Publication Date
CN114172436A CN114172436A (zh) 2022-03-11
CN114172436B true CN114172436B (zh) 2023-06-23

Family

ID=80485357

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111506688.7A Active CN114172436B (zh) 2021-12-10 2021-12-10 基于观测器的异步电动机指令滤波离散控制方法

Country Status (1)

Country Link
CN (1) CN114172436B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004166408A (ja) * 2002-11-13 2004-06-10 Yoichi Hayashi 永久磁石同期電動機制御方法
EP3160037A1 (en) * 2015-10-20 2017-04-26 ABB Schweiz AG Method for identifying the discrete instantaneous angular speed of an electromechanical system
CN106788052A (zh) * 2016-11-30 2017-05-31 青岛大学 基于观测器的异步电机命令滤波误差补偿模糊控制方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004166408A (ja) * 2002-11-13 2004-06-10 Yoichi Hayashi 永久磁石同期電動機制御方法
EP3160037A1 (en) * 2015-10-20 2017-04-26 ABB Schweiz AG Method for identifying the discrete instantaneous angular speed of an electromechanical system
CN106788052A (zh) * 2016-11-30 2017-05-31 青岛大学 基于观测器的异步电机命令滤波误差补偿模糊控制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于观测器的异步电动机命令滤波反步技术;牛浩;马玉梅;于金鹏;史发涛;孙凤龙;;青岛大学学报(工程技术版)(第03期);全文 *

Also Published As

Publication number Publication date
CN114172436A (zh) 2022-03-11

Similar Documents

Publication Publication Date Title
Saghafinia et al. Adaptive fuzzy sliding-mode control into chattering-free IM drive
Tan et al. Neural network observers and sensorless robust optimal control for partially unknown PMSM with disturbances and saturating voltages
CN108092567B (zh) 一种永磁同步电动机转速控制系统及方法
CN106788046B (zh) 永磁同步电机命令滤波有限时间模糊控制方法
Ammar Performance improvement of direct torque control for induction motor drive via fuzzy logic-feedback linearization: Simulation and experimental assessment
CN106788052B (zh) 基于观测器的异步电机命令滤波误差补偿模糊控制方法
Lin et al. Adaptive fuzzy-neural-network control for a DSP-based permanent magnet linear synchronous motor servo drive
CN108390597A (zh) 带有扰动观测器的永磁同步电机非线性预测控制器设计
CN109921698B (zh) 考虑铁损的永磁同步电机随机命令滤波神经网络控制方法
CN109728755A (zh) 一种pmsm反演终端滑模控制方法
Zerikat et al. High-performance sensorless vector control of induction motor drives using artificial intelligent technique
CN110401391A (zh) 异步电动机随机系统模糊自适应动态面控制方法
Salim et al. Sensorless passivity based control for induction motor via an adaptive observer
Boukhlouf et al. Hardware-in-the-loop implementation of an unknown input observer for synchronous reluctance motor
CN113659895B (zh) 基于指令滤波的永磁同步电动机全状态约束有限时间控制方法
Adigintla et al. Design and analysis of a speed controller for fractional‐order‐modeled voltage‐source‐inverter‐fed induction motor drive
Barambones et al. Adaptive sensorless robust control of AC drives based on sliding mode control theory
Tsai et al. Implementation of a DSP‐based speed‐sensorless adaptive control for permanent‐magnet synchronous motor drives with uncertain parameters using linear matrix inequality approach
CN114172436B (zh) 基于观测器的异步电动机指令滤波离散控制方法
Tiwari et al. ANN based RF-MRAS speed estimation of induction motor drive at low speed
CN111293941B (zh) 一种考虑铁损的永磁同步电动机有限时间动态面控制方法
CN110492809B (zh) 基于神经网络逼近的异步电动机动态面离散容错控制方法
Ganjewar Novel approaches in sensorless induction motor drive for industrial applications
Wei et al. Model-Free Predictive Control Using Sinusoidal Generalized Universal Model for PMSM Drives
Wang et al. A high performance permanent magnet synchronous motor servo system using predictive functional control and Kalman filter

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant