CN114172436B - 基于观测器的异步电动机指令滤波离散控制方法 - Google Patents
基于观测器的异步电动机指令滤波离散控制方法 Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 90
- 238000001914 filtration Methods 0.000 title claims abstract description 33
- 238000013528 artificial neural network Methods 0.000 claims abstract description 39
- 229920006395 saturated elastomer Polymers 0.000 claims description 15
- 230000003044 adaptive effect Effects 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000009499 grossing Methods 0.000 claims description 6
- 238000013178 mathematical model Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000006467 substitution reaction Methods 0.000 claims description 4
- 230000004907 flux Effects 0.000 claims description 3
- 230000001537 neural effect Effects 0.000 claims description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 230000009471 action Effects 0.000 claims description 2
- 230000001360 synchronised effect Effects 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 11
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 238000004422 calculation algorithm Methods 0.000 abstract description 2
- 238000012545 processing Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 16
- 230000008569 process Effects 0.000 description 5
- 238000001514 detection method Methods 0.000 description 4
- 239000002131 composite material Substances 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000006698 induction Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 229920002430 Fibre-reinforced plastic Polymers 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000012271 agricultural production Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 239000011151 fibre-reinforced plastic Substances 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000004377 microelectronic Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P23/00—Arrangements or methods for the control of AC motors characterised by a control method other than vector control
- H02P23/0004—Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
- H02P23/0018—Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control using neural networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P23/00—Arrangements or methods for the control of AC motors characterised by a control method other than vector control
- H02P23/12—Observer control, e.g. using Luenberger observers or Kalman filters
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P25/00—Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
- H02P25/02—Arrangements 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P2207/00—Indexing scheme relating to controlling arrangements characterised by the type of motor
- H02P2207/01—Asynchronous 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)所示。
其中,k为异步电动机离散系统的步数,θ(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分别表示定子和转子的等效电阻。
为简化以上异步电动机离散系统动态数学模型,定义如下新的变量:
由式(2)定义的新变量,将式(1)通过欧拉方法,得到异步电动机离散系统模型如下:
由于uq(k)和ud(k)的基本特性相同,为表述方便,定义u(k)代指uq(k)和ud(k)。
其中,umax>0、umin<0是未知的输入饱和常数。
vq(k)为q轴定子的输入电压,vd(k)为d轴定子的输入电压,由于vq(k)和vd(k)的基本特性相同,为表述方便,定义v(k)代指vq(k)和vd(k)。定义如下平滑函数g(v(k)):
由式(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,使得:
其中,v(k+1)表示k+1步定子的输入电压,g(v(k+1))表示k+1步的平滑函数。
则异步电动机离散系统第k步q轴的对称饱和非线性输入uq(k)为:
则异步电动机离散系统第k步d轴的对称饱和非线性输入ud(k)为:
步骤2.设计降维观测器,估计异步电动机的转子位置、角速度。
对于公式(3)示出的异步电动机驱动系统动态模型(1),其动态方程如下:
其中,未知非线性函数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)+τ;τ是逼近误差且满足|τ|≤ε,ε为任意小的正常数;是输入向量,q是神经网络输入维数,Rq为实数向量集;W∈Rp是权重向量,p是神经网络节点数,p为正整数且p>1,Rp为实数向量集;S(Z)=[s1(Z),...,sp(Z)]T∈Rp为基函数向量。sc(Z)为高斯函数,sc(Z)的表达式为:
其中,c=1,...,p,μc是高斯函数sc(Z)接受域的中心,ηc为高斯函数sc(Z)的宽度。根据径向基函数神经网络的定义,存在径向基函数神经网络使得W2表示权重向量,S2(Z2(k))表示基函数向量,式(10)表示为:
则降维观测器如公式(12)所示;
其中,为第k步降维观测器输出的估计值,/>分别表示第k+1步时异步电动机的转子角度和角速度的观测值,g1、g2为正常数;/>Φ2>0且为未知常数,/>为神经网络权重/>的范数;自适应律/>是Φ2的估计值,定义/>
选择Lyapunov函数V0(k)=hΤ(k)Ph(k),其中,P为正定矩阵,PΤ=P。
对V0(k)取一阶差分得:
其中,||P||表示正定矩阵P的范数,||S2(Z2(k))||表示神经网络径向基函数S2(Z2(k))的范数,Q为矩阵,且Q=P-3AΤPA。
步骤3.根据指令滤波和反步法设计基于观测器的异步电动机指令滤波离散控制方法。
定义指令滤波器的公式如下:
式中,ξ为指令滤波器的时间常数,ωn为正常数;d1(k)、d2(k)分别为第k步时指令滤波器的输出信号,d1(k+1)、d2(k+1)分别为指令滤波器在第k+1步时的输出信号。
α(k)为指令滤波器的输入信号,即虚拟控制函数;如果虚拟控制函数α(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)的表达式将在下面的控制方法设计中给出。
其中,x1d(k)、x4d(k)为给定的期望信号;α1d(k)、α2d(k)、α3d(k)为指令滤波器的输出信号;vj(k)表示补偿误差,j=1,2,3,4,5。
对离散系统每一步都选取一个Lyapunov函数来构建虚拟控制函数,具体步骤如下:
其中,x1d(k+1)、ξ1(k+1)分别为k+1步的期望信号和误差补偿信号。
选取虚拟控制函数α1(k)和误差补偿信号ξ1(k)为:
其中,t2为常数,且|t2|≤1,γ2和δ2为正参数,将式(19)和(20)代入式(18)中,得:
其中,f3(k)为一未知的非线性函数,且
由径向基函数神经网络逼近原理得知,对于给定的任意ε3>0,存在径向基函数神经网络使得/>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)和自适应律/>分别为:
其中,γ3和δ3为正参数,表示神经网络权重向量W3(k)的估计值,||S3(Z3(k))||表示神经网络径向基函数S3(Z3(k))的范数;定义/>Φ3为未知常数,且Φ3>0,/>为神经网络权重向量/>的范数,自适应律/>是Φ3的估计值,定义估计误差/>将式(23)、(24)代入式(22)中得将式(23)、(24)代入式(22)中并结合杨氏不等式得:
选取虚拟控制函数α3(k)和误差补偿信号ξ4(k)为:
其中,t4为常数,且|t4|≤1,将式(27)代入式(26)中,得:
取V5(k)的一阶差分,得到:
其中,f5(k)为一未知的非线性函数,且
由径向基函数神经网络逼近原理得知,对于给定的任意ε5>0,存在径向基函数神经网络使得:/>Z5(k)=[x1(k),x2(k),x3(k),x4(k),x5(k)]Τ,τ5表示逼近误差,并满足不等式|τ5|≤ε5。取误差补偿信号ξ5(k)=0则异步电动机d轴定子输入电压vd(k),对称饱和非线性输入信号ud(k)和自适应律/>分别为:
其中,γ5和δ5为正参数,表示神经网络权重向量W5(k)的估计值,||S5(Z5(k))||表示神经网络径向基函数S5(Z5(k))的范数,/>定义/>为神经网络权重向量/>的范数Φ5为未知常数且Φ5>0,估计误差/>是Φ5的估计值;将式(30)、(31)代入式(29)中得将式(30)、(31)代入式(29)中并结合杨氏不等式得:
步骤4.对构建的基于观测器的异步电动机指令滤波离散控制方法进行稳定性分析。
取V(k)的一阶差分,得到:
其中,m=2,3,5,δm为正参数,由基函数向量S(Z)的定义得知,||Sm(Zm(k))||2≤lm,lm表示径向基函数神经网络Sm(Zm(k))的节点数;由杨氏不等式,得:
将式(36)、(37)代入式(35)中得到:
将式(38)、(39)、(40)代入式(33)中得到:
选择合适的采样周期Δt,控制参数ξ,ωn,γ2,γ3,γ5,δ2,δ3,δ5,g1和g2,使下列不等式成立:
κ=2,3,5,且矩阵Q=P-3AΤPA正定,只需满足下列条件:
由以上分析得到,在定子的输入电压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。
选择的控制参数为:ξ=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。
相应的仿真结果如图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)所示;
其中,k为异步电动机离散系统的步数,θ(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分别表示定子和转子的等效电阻;
为简化以上异步电动机离散系统动态数学模型,定义如下新的变量:
由式(2)定义的新变量,将式(1)通过欧拉方法,得到异步电动机离散系统模型如下:
由于uq(k)和ud(k)的基本特性相同,为表述方便,定义u(k)代指uq(k)和ud(k);
其中,umax>0、umin<0是未知的输入饱和常数;vq(k)和vd(k)均为定子的输入电压,由于vq(k)和vd(k)的基本特性相同,为表述方便,定义v(k)代指vq(k)和vd(k);
由式(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,使得:
则异步电动机离散系统第k步q轴的对称饱和非线性输入uq(k)为:
则异步电动机离散系统第k步d轴的对称饱和非线性输入ud(k)为:
步骤2.设计降维观测器,估计异步电动机的转子位置和角速度;
对于公式(3)给出的异步电动机离散系统模型,其动态方程如下:
其中,未知非线性函数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)+τ;τ是逼近误差且满足|τ|≤ε,ε为任意小的正常数;是输入向量,q是神经网络输入维数,Rq为实数向量集;W∈Rp是权重向量,p是神经网络节点数,p为正整数且p>1,Rp为实数向量集;S(Z)=[s1(Z),...,sp(Z)]T∈Rp为基函数向量;sc(Z)为高斯函数,sc(Z)的表达式为:
其中,c=1,...,p,μc是高斯函数sc(Z)接受域的中心,ηc为高斯函数sc(Z)的宽度;
则降维观测器如公式(12)所示;
选择Lyapunov函数V0(k)=hT(k)Ph(k),其中,P为正定矩阵,PT=P;
对V0(k)取一阶差分得:
其中,||P||表示正定矩阵P的范数,||S2(Z2(k))||表示神经网络径向基函数S2(Z2(k))的范数,Q为矩阵,且Q=P-3ATPA;
步骤3.根据指令滤波和反步法设计基于观测器的异步电动机指令滤波离散控制方法;
定义指令滤波器的公式如下:
式中,ξ为指令滤波器的时间常数,ωn为正常数;d1(k)、d2(k)分别为第k步时指令滤波器的输出信号,d1(k+1)、d2(k+1)分别为指令滤波器在第k+1步时的输出信号;
α(k)为指令滤波器的输入信号,即虚拟控制函数;如果虚拟控制函数α(k)满足:且/>对于任意的k≥1成立;其中,/> 且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;
其中,x1d(k)、x4d(k)为给定的期望信号;α1d(k)、α2d(k)、α3d(k)为指令滤波器的输出信号;vj(k)表示补偿误差,j=1,2,3,4,5;
对离散系统每一步都选取一个Lyapunov函数来构建虚拟控制函数,具体步骤如下:
其中,x1d(k+1)、ξ1(k+1)分别为k+1步的期望信号和误差补偿信号;
选取虚拟控制函数α1(k)和误差补偿信号ξ1(k)为:
其中,t2为常数,且|t2|≤1,γ2和δ2为正参数,将式(19)和(20)代入式(18)中,得:
其中,f3(k)为一未知的非线性函数,且
由径向基函数神经网络逼近原理得知,对于给定的任意ε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)和自适应律分别为:
将式(23)、(24)代入式(22)中得将式(23)、(24)代入式(22)中并结合杨氏不等式得:
选取虚拟控制函数α3(k)和误差补偿信号ξ4(k)为:
其中,t4为常数,且|t4|≤1,将式(27)代入式(26)中,得:
取V5(k)的一阶差分,得到:
其中,f5(k)为一未知的非线性函数,且
由径向基函数神经网络逼近原理得知,对于给定的任意ε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)和自适应律分别为:
其中,γ5和δ5为正参数,表示神经网络权重向量W5(k)的估计值,||S5(Z5(k))||表示神经网络径向基函数S5(Z5(k))的范数,/>定义Φ5=||W5 T||,||W5 T||为神经网络权重向量W5 T的范数Φ5为未知常数且Φ5>0,估计误差/> 是Φ5的估计值;将式(30)、(31)代入式(29)中得将式(30)、(31)代入式(29)中并结合杨氏不等式得:
步骤4.对构建的基于观测器的异步电动机指令滤波离散控制方法进行稳定性分析;
取V(k)的一阶差分,得到:
其中,m=2,3,5,δm为正参数,由基函数向量S(Z)的定义得知,||Sm(Zm(k))||2≤lm,lm表示径向基函数神经网络Sm(Zm(k))的节点数;由杨氏不等式,得:
将式(36)、(37)代入式(35)中得到:
将式(38)、(39)、(40)代入式(33)中得到:
选择合适的采样周期Δt,控制参数ξ,ωn,γ2,γ3,γ5,δ2,δ3,δ5,g1和g2,使下列不等式成立:
其中,κ=2,3,5,且矩阵Q=P-3ATPA正定,只需满足下列条件:
由以上分析得到,在定子的输入电压vq(k)和vd(k)的作用下,异步电动机离散系统的误差变量e1(k)能够收敛到原点的一个充分小的邻域内,并保证其他信号有界。
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)
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 | 青岛大学 | 基于观测器的异步电机命令滤波误差补偿模糊控制方法 |
-
2021
- 2021-12-10 CN CN202111506688.7A patent/CN114172436B/zh active Active
Patent Citations (3)
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)
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 |