CN113328668A - 一种基于离散全阶观测器的感应电机转速观测方法 - Google Patents
一种基于离散全阶观测器的感应电机转速观测方法 Download PDFInfo
- Publication number
- CN113328668A CN113328668A CN202110594045.6A CN202110594045A CN113328668A CN 113328668 A CN113328668 A CN 113328668A CN 202110594045 A CN202110594045 A CN 202110594045A CN 113328668 A CN113328668 A CN 113328668A
- Authority
- CN
- China
- Prior art keywords
- sampling period
- current sampling
- current
- axis component
- state variable
- 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.)
- Granted
Links
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
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
- H02P21/18—Estimation of position or speed
-
- 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
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/0003—Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
- H02P21/0017—Model reference adaptation, e.g. MRAS or MRAC, useful for control or parameter estimation
-
- 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
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/13—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
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
-
- 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
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/22—Current control, e.g. using a current control loop
-
- 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
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/24—Vector control not involving the use of rotor position or rotor speed sensors
- H02P21/26—Rotor flux based control
-
- 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
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/24—Vector control not involving the use of rotor position or rotor speed sensors
- H02P21/28—Stator flux based control
-
- 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)
- Control Of Ac Motors In General (AREA)
Abstract
一种基于离散全阶观测器的感应电机转速观测方法,它属于电机控制技术领域。本发明解决了在降低载波比的情况下,现有的无速度传感器感应电机离散全阶观测器的稳定性与观测精度差的问题。本发明首先建立全阶观测器连续域模型,其次基于全阶观测器连续域模型,分别预测一个采样周期区间内四个状态变量点的导数,并取四者的加权平均数;最后,根据当前状态变量值、加权平均数与采样周期长度Ts,估计得到下一个采样周期时状态变量值。采用本发明所设计的离散全阶观测器,即使在低开关频率时也能准确观测转速,在保持较低计算量的同时,实现了无速度传感器感应电机驱动系统在极低载波比下的高稳定高精度运行。本发明可以应用于感应电机转速观测。
Description
技术领域
本发明属于电机控制技术领域,具体涉及一种基于离散全阶观测器的感应电机转速观测方法。
背景技术
无速度传感器感应电机矢量控制技术具有高可靠性、结构简单、成本低廉等优点,广泛应用于石油钻机、矿车、数控切割机等各种工业场合。现如今,已有多种感应电机转速观测方法被提出,如:模型参考自适应系统、滑模观测器、自适应全阶观测器等。其中自适应全阶观测器通过合适的自适应率,利用测量的定子电流与其估计值之间的误差对估测变量进行校正。全阶观测器具有较强的参数鲁棒性,合理选择速度自适应律和观测器增益,可以使无速度传感器控制系统达到与配备速度传感器的控制系统类似的性能。
通常,全阶观测器在连续域中设计,但在基于工业微处理器的系统中,模型需进行离散化处理,而前向欧拉法因其方法简单而被作为最常用的离散方法。然而离散化必然会带来离散化误差,导致观测器性能恶化甚至不稳定。尤其在大功率或超高速驱动器中,为减小逆变器开关损耗,常常会让逆变器开关频率与电机运行频率之比减小,即降低载波比,此时基于前向欧拉法的全阶观测器模型的离散误差会随着载波比的减小急剧增大,导致观测器严重的震荡甚至不稳定。针对此问题,后向欧拉法和双线性法虽可提高观测器的估计精度和稳定性,但是其高计算量和隐式特性增加了系统的复杂性,所以寻求一种感应电机离散全阶观测器稳定性与精度提升策略对整个无速度传感器矢量控制系统的性能提升具有重要意义。
发明内容
本发明的目的是为解决在降低载波比的情况下,现有的无速度传感器感应电机离散全阶观测器的稳定性与观测精度差的问题,而提出了一种基于离散全阶观测器的感应电机转速观测方法。
本发明为解决上述技术问题所采取的技术方案是:一种基于离散全阶观测器的感应电机转速观测方法,所述方法具体包括以下步骤:
步骤一、将电机定子电流与转子磁链作为状态变量,建立感应电机数学模型;
基于建立的感应电机数学模型,推导出连续域全阶观测器的数学模型,并根据连续域全阶观测器数学模型计算当前采样周期起点时状态变量的一阶导数K(k);
步骤二、根据当前采样周期起点时的状态变量观测值和步骤一计算得到的K(k)估计当前采样周期中点时状态变量观测值,并根据当前采样周期中点时状态变量观测值和步骤一的连续域全阶观测器数学模型计算当前采样周期中点时状态变量的一阶导数K′(k);
步骤三、根据当前采样周期起点时的状态变量观测值和步骤二计算得到的K′(k)重新估计当前采样周期中点时状态变量观测值,并根据重新估计出的当前采样周期中点时状态变量观测值和步骤一的连续域全阶观测器数学模型重新计算当前采样周期中点时状态变量的一阶导数K″(k);
步骤四、根据当前采样周期起点时的状态变量观测值和步骤三计算得到的K″(k)估计当前采样周期终点时状态变量观测值,并根据当前采样周期终点时状态变量观测值和步骤一的连续域全阶观测器数学模型计算当前采样周期终点时状态变量的一阶导数K″′(k);
步骤五、根据步骤一至步骤四计算得到的K(k)、K′(k)、K″(k)以及K″′(k),计算K(k)、K′(k)、K″(k)和K″′(k)的加权平均数,根据计算出的加权平均数、当前采样周期起点时的状态变量观测值与采样周期长度,最终估计出下一个采样周期的状态变量观测值;
步骤六、根据估计出的下一个采样周期的状态变量观测值和采集的实际定子电流矢量,计算出在下一个采样周期感应电机的转速。
本发明的有益效果是:本发明提出了一种基于离散全阶观测器的感应电机转速观测方法,本发明首先建立全阶观测器连续域模型,其次基于全阶观测器连续域模型,分别预测一个采样周期区间内四个状态变量点的导数,并取四者的加权平均数;最后,根据当前状态变量值、加权平均数与采样周期长度Ts,估计得到下一个采样周期时状态变量值。本发明所设计的离散全阶观测器,即使在低开关频率(500Hz-1500Hz)时也能准确观测转速,当Ts=1/1200s时,基于本发明所提方法的最大离散误差小于2%,在保持较低计算量的同时,实现了无速度传感器感应电机驱动系统在极低载波比下的高稳定高精度运行。
附图说明
图1为本发明基于全阶观测器的无速度传感器感应电机矢量控制系统框图;
图中,ωr,ref为系统给定转速,为全阶观测器观测转速,为全阶观测器观测磁链,isd,rated为励磁电流给定值,isd,ref为弱磁分配的给定d轴电流分量,isq,ref为弱磁分配的给定q轴电流分量,isd为检测到的定子电流d轴分量,isq为检测到的定子电流q轴分量,为坐标变换旋转角,αβ为两相静止坐标系,dq为两相旋转坐标系,SVPWM为电压空间矢量控制,udc为直流电压,abc为三相坐标系,为定子电流矢量,为定子电流误差矢量,为定子电压矢量;
图2a为Ts=1/1200s时基于前向欧拉法的离散全阶观测器极点分布情况的示意图;
图2b为Ts=1/600s时基于前向欧拉法的离散全阶观测器极点分布情况的示意图;
图3a为Ts=1/1200s时基于预估校正法的离散全阶观测器极点分布情况的示意图;
图3b为Ts=1/600s时基于预估校正法的离散全阶观测器极点分布情况的示意图;
图4a为Ts=1/1200s时基于本发明所提出策略的离散全阶观测器极点分布情况的示意图;
图4b为Ts=1/600s时基于本发明所提出策略的离散全阶观测器极点分布情况的示意图;
图5为本发明所提出的离散全阶观测器与基于前向欧拉法和预估校正方法的离散全阶观测器的离散误差的对比分析图;
图6a为转速频率30Hz,空载下基于前向欧拉法的离散全阶观测器的定子电流α轴分量波形图;
图6b为转速频率30Hz,空载下基于预估校正法的离散全阶观测器的定子电流α轴分量波形图;
图6c为转速频率30Hz,空载下本发明所提出的离散全阶观测器的定子电流α轴分量波形图;
图7a为转速频率30Hz,额定负载下基于前向欧拉法的离散全阶观测器的定子电流α轴分量波形图;
图7b为转速频率30Hz,额定负载下基于预估校正法的离散全阶观测器的定子电流α轴分量波形图;
图7c为转速频率30Hz,额定负载下本发明所提出的离散全阶观测器的定子电流α轴分量波形图;
图8a为转速频率120Hz,空载下基于预估校正法的离散全阶观测器的定子电流α轴分量波形图;
图8b为转速频率120Hz,本发明所提出的离散全阶观测器的定子电流α轴分量波形图;
图9a为转速频率120Hz,轻载(4N·m)下基于预估校正法的离散全阶观测器的定子电流α轴分量波形图;
图9b为转速频率120Hz,轻载(4N·m)下本发明所提出的离散全阶观测器的定子电流α轴分量波形图。
具体实施方式
具体实施方式一、结合图1说明本实施方式。本实施方式所述的一种基于离散全阶观测器的感应电机转速观测方法,所述方法具体包括以下步骤:
步骤一、将电机定子电流与转子磁链作为状态变量,建立感应电机数学模型;
基于建立的感应电机数学模型,推导出连续域全阶观测器的数学模型,并根据连续域全阶观测器数学模型计算当前采样周期起点时状态变量的一阶导数K(k);
步骤二、根据当前采样周期起点时的状态变量观测值和步骤一计算得到的K(k)估计当前采样周期中点时状态变量观测值,并根据当前采样周期中点时状态变量观测值和步骤一的连续域全阶观测器数学模型计算当前采样周期中点时状态变量的一阶导数K′(k);
步骤三、根据当前采样周期起点时的状态变量观测值和步骤二计算得到的K′(k)重新估计当前采样周期中点时状态变量观测值,并根据重新估计出的当前采样周期中点时状态变量观测值和步骤一的连续域全阶观测器数学模型重新计算当前采样周期中点时状态变量的一阶导数K″(k);
步骤四、根据当前采样周期起点时的状态变量观测值和步骤三计算得到的K″(k)估计当前采样周期终点时状态变量观测值,并根据当前采样周期终点时状态变量观测值和步骤一的连续域全阶观测器数学模型计算当前采样周期终点时状态变量的一阶导数K″′(k);
步骤五、根据步骤一至步骤四计算得到的K(k)、K′(k)、K″(k)以及K″′(k),计算K(k)、K′(k)、K″(k)和K″′(k)的加权平均数,根据计算出的加权平均数、当前采样周期起点时的状态变量观测值与采样周期长度,最终估计出下一个采样周期的状态变量观测值;
步骤六、根据估计出的下一个采样周期的状态变量观测值和采集的实际定子电流矢量,观测出在下一个采样周期感应电机的转速。
具体实施方式二:本实施方式与具体实施方式一不同的是,所述步骤一中,将电机定子电流与转子磁链作为状态变量,建立感应电机数学模型;其具体过程为:
式中:定子电流矢量isα是在两相坐标系下的定子α轴分量,isβ是在两相坐标系下的定子β轴分量,转子磁链矢量λrα是在两相坐标系下的定子α轴分量,λrβ是在两相坐标系下的定子β轴分量,输入电压矢量u=[usα usβ]T,usα是u在两相坐标系下的定子α轴分量,usβ是u在两相坐标系下的定子β轴分量,A为系统矩阵,A11、A12、A21、A22为系统矩阵系数,x为状态变量矩阵, 为当前采样周期的定子电流矢量,为当前采样周期的转子磁链矢量,B为输入矩阵,b为电压项系数,t为时间变量。
具体实施方式三:本实施方式与具体实施方式二不同的是,所述系统矩阵系数A11、A12、A21、A22的表达式为:
A11=A′11I
A12=A′121I+A′122J
A21=A′21I
A22=A′221I+ωrJ
具体实施方式四:本实施方式与具体实施方式三不同的是,所述电压项系数b的表达式为:
b=b′I
具体实施方式五:本实施方式与具体实施方式四不同的是,所述基于建立的感应电机数学模型,推导出连续域全阶观测器的数学模型;其具体过程为:
由式(1)减去式(2),可以得状态误差方程如下:
根据李雅普诺夫稳定性理论,推导转速自适应律为:
其中,kp为比例系数,ki为积分系数,s为拉普拉斯算子。
具体实施方式六:本实施方式与具体实施方式五不同的是,所述根据连续域全阶观测器数学模型(公式(2))计算当前采样周期起点时状态变量的一阶导数K(k);其具体过程为:
将公式(5)展开为:
其中:“Ksα(k)”代表当前采样周期起点时两相坐标系下定子电流α轴分量的一阶导数,“Ksβ(k)”代表当前采样周期起点时两相坐标系下定子电流β轴分量的一阶导数,“Krα(k)”代表当前采样周期起点时两相坐标系下转子磁链α轴分量的一阶导数,“Krβ(k)”代表当前采样周期起点时两相坐标系下转子磁链β轴分量的一阶导数,代表当前采样周期定子电流α轴分量观测值,代表当前采样周期定子电流β轴分量观测值,代表当前采样周期转子磁链α轴分量观测值,代表当前采样周期转子磁链β轴分量观测值,uα(k)代表当前采样周期定子电压α轴分量,uβ(k)代表当前采样周期定子电压β轴分量。
具体实施方式七:本实施方式与具体实施方式六不同的是,所述步骤二的具体过程为:
将公式(7)和公式(8)联立得:
将公式(9)展开为:
其中:“K′sα(k)”代表当前采样周期中点时两相坐标系下定子电流α轴分量的一阶导数,“K′sβ(k)”代表当前采样周期中点时两相坐标系下定子电流β轴分量的一阶导数,“K′rα(k)”代表当前采样周期中点时两相坐标系下转子磁链α轴分量的一阶导数,“K′rβ(k)”代表当前采样周期中点时两相坐标系下转子磁链β轴分量的一阶导数。
具体实施方式八:本实施方式与具体实施方式七不同的是,所述步骤三的具体过程为:
将公式(11)和公式(12)联立得:
将公式(13)展开为:
其中:“K″sα(k)”代表重新估计出的当前采样周期中点时两相坐标系下定子电流α轴分量一阶导数,“K″sβ(k)”代表重新估计出的当前采样周期中点时两相坐标系下定子电流β轴分量一阶导数,“K″rα(k)”代表重新估计出的当前采样周期中点时两相坐标系下转子磁链α轴分量一阶导数,“K″rβ(k)”代表重新估计出的当前采样周期中点时两相坐标系下转子磁链β轴分量一阶导数。
具体实施方式九:本实施方式与具体实施方式八不同的是,所述步骤四的具体过程为:
将公式(15)和公式(16)联立得:
将公式(17)展开为:
其中:uα(k+1)代表下一采样周期定子电压α轴分量,uβ(k+1)代表下一采样周期定子电压β轴分量,“K″′sα(k)”代表当前采样周期终点时两相坐标系下定子电流α轴分量的一阶导数,“K″′sβ(k)”代表当前采样周期终点时两相坐标系下定子电流β轴分量的一阶导数,“K″′rα(k)”代表当前采样周期终点时两相坐标系下转子磁链α轴分量的一阶导数,“K″′rβ(k)”代表当前采样周期终点时两相坐标系下转子磁链β轴分量的一阶导数。
具体实施方式十:本实施方式与具体实施方式九不同的是,所述步骤五的具体过程为:
将公式(19)和公式(20)联立得:
将公式(21)展开为:
采用以下实施例验证本发明的有益效果:
为进一步提升低载波比时离散全阶观测器的稳定性与精度,本发明首先分析离散全阶观测器极点分布情况对离散全阶观测器的稳定性进行分析。为保证离散全阶观测器是稳定的,观测器的所有极点必须位于单位圆内。图2a和图2b为基于前向欧拉法的离散全阶观测器极点分布情况;图3a和图3b为基于预估校正法的离散全阶观测器极点分布情况;图4为本发明的离散全阶观测器极点分布情况。
由图2a可知,对于基于前向欧拉法的离散全阶观测器,在Ts=1/1200s时,离散全阶观测器的临界稳定点为ωr=70Hz,此时载波比为17。由图2b可知,在Ts=1/600s时,离散全阶观测器的临界稳定点为ωr=50Hz,此时载波比为12。因此在低载波比时,采用前向欧拉法的离散全阶观测器不能在全速范围内(0-200Hz)稳定运行。由图3a可知,对于预估校正法的离散全阶观测器,在Ts=1/1200s时,离散全阶观测器的临界稳定点为ωr=180Hz,此时载波比为6.7。由图3b可知,在Ts=1/600s时,离散全阶观测器的临界稳定点为ωr=110Hz,此时载波比为5.1。因此在低载波比时,采用预估校正方法的离散全阶观测器不能在全速范围内(0-200Hz)稳定运行。
本发明的离散全阶观测器极点分布情况如图4a和图4b所示,由图中可以看出,所有极点均在单位圆内,即使在载波比为3时,本发明所提的离散全阶观测器也能在全速范围内稳定运行。
本发明采用矩阵的2-范数对全阶观测器近似离散模型与精确离散模型的误差进行定量分析。矩阵的2-范数定义为:
式中:A∈Rn×n,λmax(ATA)为ATA的最大特征值。
根据公式(2),可以看到全阶观测器模型为一阶常微分方程,计算机通常用数值方法对微分方程求解,即通过某种离散方法,将微分方程转换差分方程后对其求解。定义Ts为离散步长,基于非齐次矩阵的矩阵方程定理,可得离散全阶观测器的精确离散形式为:
其中:φ为离散全阶观测器的状态转移矩阵。
根据离散观测器数学模型,经过数学推导,可得基于前向欧拉法的离散全阶观测器的状态转移矩阵为:
φ1=1+ATs (27)
基于预估校正法的离散全阶观测器的状态转移矩阵为:
本发明所提出的离散全阶观测器的状态转移矩阵为:
分别将公式(27)、(28)、(29)的状态转移矩阵φ1、φ2、φ3带入到公式(24)的φ中,得到如图5所示的基于本发明所提出策略的离散全阶观测器与基于前向欧拉法和预估校正方法的离散全阶观测器的离散误差的对比分析。水平轴代表着转子转速频率,单位为Hz,垂直轴用百分比表示离散误差。从图5中可以看出,随着转速的增大,三种离散方法的离散误差随着转子转速的上升和离散步长的增加而增大。但是,可以明显看出,相较于基于前向欧拉法和预估校正法的离散全阶观测器,本发明所提出的离散全阶观测器的离散误差要低的多。当Ts=1/1200s时,基于本发明所提方法近似的离散全阶观测器的最大误差小于2%,而基于前向欧拉法和预估校正法的离散全阶观测器的最大误差分别在20%和60%以下。即使在Ts=1/600s的情况下,基于本发明所提出策略的离散全阶观测器在150hz时的误差也在10%内,远小于基于前向欧拉法和预估校正法的离散观测器。
实验部分:
在开关频率为600Hz,不同电机转速、负载情况下,分析在采用不同离散全阶观测器时相应的定子电流α轴分量的实验波形及其估计值。图6a为转速频率30Hz,空载下基于前向欧拉法的离散全阶观测器的定子电流α轴分量波形;图6b为转速频率30Hz,空载下基于预估校正法的离散全阶观测器的定子电流α轴分量波形;图6c为转速频率30Hz,空载下本发明所提出的离散全阶观测器的定子电流α轴分量波形。图7a为转速频率30Hz,额定负载下基于前向欧拉法的离散全阶观测器的定子电流α轴分量波形;图7b为转速频率30Hz,额定负载下基于预估校正法的离散全阶观测器的定子电流α轴分量波形;图7c为转速频率30Hz,额定负载下本发明所提出的离散全阶观测器的定子电流α轴分量波形。
在开关频率600Hz,给定转速频率30Hz时,由图6a、图6b、图6c、图7a、图7b和图7c可以看到,基于前向欧拉法的离散全阶观测器对定子电流估计值的幅值与相位均有很大误差,且空载时估计误差较负载时大,离散观测器精度差,此时感应电机在无速度传感器模式下已无法高精度运行。但是基于预估校正法和本发明所提出的离散全阶观测器对定子电流仍估计精确,离散观测器精度良好。
图8a为转速频率120Hz,空载下基于预估校正法的离散全阶观测器的定子电流α轴分量波形;图8b为转速频率120Hz,本发明所提出的离散全阶观测器的定子电流α轴分量波形。图9a为转速频率120Hz,轻载(4N·m)下基于预估校正法的离散全阶观测器的定子电流α轴分量波形;图9b为转速频率120Hz,轻载(4N·m)下本发明所提出的离散全阶观测器的定子电流α轴分量波形。
在开关频率600Hz,给定转速频率120Hz时,由图8a、图8b、图9a和图9b可以看到,无论是否带载,随着转速的升高,基于预估校正法法的离散全阶观测器对定子电流估计值的幅值误差与相位误差逐渐增大,空载时估计误差较大,离散观测器精度差,此时感应电机在无速度传感器模式下已无法高精度运行。但是本发明所提出的离散全阶观测器对定子电流仍估计精确,实现了无传感器感应电机驱动系统在低载波比下的高精度运行。
以上对所提出的一种感应电机离散全阶观测器稳定性与精度提升策略进行了详细介绍,以上实例介绍仅用帮助理解本发明的具体方法及核心思想;对于相关领域的一般技术人员,依照本发明的思想,在上述说明的基础上还可做出其他改变,这里无法对所有方法逐一例举,凡是依照本发明所进行的明显改变和变形都属于本发明所附的权利要求的保护范围。
Claims (10)
1.一种基于离散全阶观测器的感应电机转速观测方法,其特征在于,所述方法具体包括以下步骤:
步骤一、将电机定子电流与转子磁链作为状态变量,建立感应电机数学模型;
基于建立的感应电机数学模型,推导出连续域全阶观测器的数学模型,并根据连续域全阶观测器数学模型计算当前采样周期起点时状态变量的一阶导数K(k);
步骤二、根据当前采样周期起点时的状态变量观测值和步骤一计算得到的K(k)估计当前采样周期中点时状态变量观测值,并根据当前采样周期中点时状态变量观测值和步骤一的连续域全阶观测器数学模型计算当前采样周期中点时状态变量的一阶导数K′(k);
步骤三、根据当前采样周期起点时的状态变量观测值和步骤二计算得到的K′(k)重新估计当前采样周期中点时状态变量观测值,并根据重新估计出的当前采样周期中点时状态变量观测值和步骤一的连续域全阶观测器数学模型重新计算当前采样周期中点时状态变量的一阶导数K″(k);
步骤四、根据当前采样周期起点时的状态变量观测值和步骤三计算得到的K″(k)估计当前采样周期终点时状态变量观测值,并根据当前采样周期终点时状态变量观测值和步骤一的连续域全阶观测器数学模型计算当前采样周期终点时状态变量的一阶导数K″′(k);
步骤五、根据步骤一至步骤四计算得到的K(k)、K′(k)、K″(k)以及K″′(k),计算K(k)、K′(k)、K″(k)和K″′(k)的加权平均数,根据计算出的加权平均数、当前采样周期起点时的状态变量观测值与采样周期长度,最终估计出下一个采样周期的状态变量观测值;
步骤六、根据估计出的下一个采样周期的状态变量观测值和采集的实际定子电流矢量,观测出在下一个采样周期感应电机的转速。
2.根据权利要求1所述的一种基于离散全阶观测器的感应电机转速观测方法,其特征在于,所述步骤一中,将电机定子电流与转子磁链作为状态变量,建立感应电机数学模型;其具体过程为:
6.根据权利要求5所述的一种基于离散全阶观测器的感应电机转速观测方法,其特征在于,所述根据连续域全阶观测器数学模型计算当前采样周期起点时状态变量的一阶导数K(k);其具体过程为:
将公式(5)展开为:
9.根据权利要求8所述的一种基于离散全阶观测器的感应电机转速观测方法,其特征在于,所述步骤四的具体过程为:
将公式(15)和公式(16)联立得:
将公式(17)展开为:
其中:uα(k+1)代表下一采样周期定子电压α轴分量,uβ(k+1)代表下一采样周期定子电压β轴分量,“K″′sα(k)”代表当前采样周期终点时两相坐标系下定子电流α轴分量的一阶导数,“K″′sβ(k)”代表当前采样周期终点时两相坐标系下定子电流β轴分量的一阶导数,“K″′rα(k)”代表当前采样周期终点时两相坐标系下转子磁链α轴分量的一阶导数,“K″′rβ(k)”代表当前采样周期终点时两相坐标系下转子磁链β轴分量的一阶导数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110594045.6A CN113328668B (zh) | 2021-05-28 | 2021-05-28 | 一种基于离散全阶观测器的感应电机转速观测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110594045.6A CN113328668B (zh) | 2021-05-28 | 2021-05-28 | 一种基于离散全阶观测器的感应电机转速观测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113328668A true CN113328668A (zh) | 2021-08-31 |
CN113328668B CN113328668B (zh) | 2022-01-14 |
Family
ID=77422266
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110594045.6A Active CN113328668B (zh) | 2021-05-28 | 2021-05-28 | 一种基于离散全阶观测器的感应电机转速观测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113328668B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103701386A (zh) * | 2014-01-03 | 2014-04-02 | 哈尔滨工业大学 | 基于观测磁链误差的异步电机无速度传感器的全阶磁链观测器的获取方法 |
CN106330043A (zh) * | 2015-07-03 | 2017-01-11 | 日立(中国)研究开发有限公司 | 永磁同步电动机的控制方法及控制装置 |
CN107994826A (zh) * | 2017-12-12 | 2018-05-04 | 合肥工业大学 | 一种基于误差加权的全阶观测器无速度传感器控制系统 |
CN109274303A (zh) * | 2018-11-29 | 2019-01-25 | 合肥工业大学 | 一种基于扩展控制集的异步电机模型预测控制方法 |
US20190199258A1 (en) * | 2017-12-22 | 2019-06-27 | Inalfa Roof Systems Group B.V. | Controlling a motor of a closure and/or blind in a vehicle body based on a disturbance observer signal |
CN110492815A (zh) * | 2019-08-23 | 2019-11-22 | 哈尔滨工业大学 | 无传感器感应电机极低速稳定性与动态性同步优化方法 |
CN111082726A (zh) * | 2019-12-19 | 2020-04-28 | 南京航空航天大学 | 一种永磁电机伺服系统的电流控制方法 |
-
2021
- 2021-05-28 CN CN202110594045.6A patent/CN113328668B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103701386A (zh) * | 2014-01-03 | 2014-04-02 | 哈尔滨工业大学 | 基于观测磁链误差的异步电机无速度传感器的全阶磁链观测器的获取方法 |
CN106330043A (zh) * | 2015-07-03 | 2017-01-11 | 日立(中国)研究开发有限公司 | 永磁同步电动机的控制方法及控制装置 |
CN107994826A (zh) * | 2017-12-12 | 2018-05-04 | 合肥工业大学 | 一种基于误差加权的全阶观测器无速度传感器控制系统 |
US20190199258A1 (en) * | 2017-12-22 | 2019-06-27 | Inalfa Roof Systems Group B.V. | Controlling a motor of a closure and/or blind in a vehicle body based on a disturbance observer signal |
CN109274303A (zh) * | 2018-11-29 | 2019-01-25 | 合肥工业大学 | 一种基于扩展控制集的异步电机模型预测控制方法 |
CN110492815A (zh) * | 2019-08-23 | 2019-11-22 | 哈尔滨工业大学 | 无传感器感应电机极低速稳定性与动态性同步优化方法 |
CN111082726A (zh) * | 2019-12-19 | 2020-04-28 | 南京航空航天大学 | 一种永磁电机伺服系统的电流控制方法 |
Non-Patent Citations (2)
Title |
---|
CHENG LUO ET AL.: "Operating-Point Tracking Method for Sensorless Induction Motor Stability Enhancement in Low-Speed Regenerating Mode", 《 IEEE TRANSACTIONS ON INDUSTRIAL ELECTRONICS》 * |
马继先 等: "基于全阶状态观测的无速度传感器矢量控制仿真研究", 《江苏科技大学学报(自然科学版)》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113328668B (zh) | 2022-01-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110460281B (zh) | 一种三电平永磁同步电机双矢量模型预测磁链控制方法 | |
US6194865B1 (en) | Control method and system for electric rotary machine | |
CN209844868U (zh) | 永磁同步电机无差拍电流预测控制系统 | |
US9525377B2 (en) | System and method of rotor time constant online identification in an AC induction machine | |
CN110995076A (zh) | 一种永磁同步电机模型预测电流控制方法 | |
CN112422004B (zh) | 一种永磁同步电机弱磁控制模式下的扰动抑制方法 | |
CN112713834B (zh) | 一种永磁同步电机无位置传感器控制方法及系统 | |
US20030128009A1 (en) | Motor control apparatus and method | |
CN113691186B (zh) | 一种永磁同步电机无位置传感器控制转子位置角补偿方法 | |
CN110995072B (zh) | 一种电机转子位置估测方法 | |
CN110912485B (zh) | 一种考虑铁损电阻的永磁同步电机无位置传感器控制方法 | |
CN112422002B (zh) | 一种鲁棒性的永磁同步电机单电流传感器预测控制方法 | |
CN111711392B (zh) | 永磁同步电机单电流传感器预测控制和参数扰动抑制方法 | |
CN110535389B (zh) | 一种永磁同步电机系统过采样预测电流控制方法 | |
CN111478637A (zh) | 电机控制方法及电机控制系统 | |
CN111478632B (zh) | 一种提高永磁同步电机参数鲁棒性的无观测器控制方法 | |
CN113794419A (zh) | 一种基于增量式电机模型的无差拍电流预测鲁棒控制方法 | |
CN113328668B (zh) | 一种基于离散全阶观测器的感应电机转速观测方法 | |
CN106712629A (zh) | 一种永磁同步电机的电流控制方法 | |
CN114301361B (zh) | 一种基于母线电流控制的无电解电容永磁同步电机驱动系统控制方法 | |
CN111181462B (zh) | 基于变步长神经网络的表贴式永磁同步电机参数辨识方法 | |
JP2005160199A (ja) | 3相交流電動機の制御装置および制御方法 | |
CN109842337B (zh) | 一种磁链差异量处理方法及电机控制装置 | |
CN115133832B (zh) | 一种表贴式永磁同步电机的参数实时修正方法 | |
CN110460279B (zh) | 一种基于辅助变量的扩展转速自适应观测器低速控制方法 |
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 |