CN111181458A - 基于扩展卡尔曼滤波器的表贴式永磁同步电机转子磁链观测方法 - Google Patents
基于扩展卡尔曼滤波器的表贴式永磁同步电机转子磁链观测方法 Download PDFInfo
- Publication number
- CN111181458A CN111181458A CN202010062837.4A CN202010062837A CN111181458A CN 111181458 A CN111181458 A CN 111181458A CN 202010062837 A CN202010062837 A CN 202010062837A CN 111181458 A CN111181458 A CN 111181458A
- Authority
- CN
- China
- Prior art keywords
- permanent magnet
- flux linkage
- current
- state
- magnet synchronous
- 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.)
- Pending
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/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
- H02P21/141—Flux 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
- 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
- H02P25/022—Synchronous motors
-
- 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/05—Synchronous machines, e.g. with permanent magnets or DC excitation
- H02P2207/055—Surface mounted magnet motors
Abstract
本发明提供基于扩展卡尔曼滤波器的表贴式永磁同步电机转子磁链观测方法,以矢量控制为基础,依据表贴式永磁同步电机两相静止坐标系下的数学模型,得出用于转子磁链观测的数学模型,将其代入扩展卡尔曼滤波算法中,实现对表贴式永磁同步电机转子磁链的观测。本发明的基于扩展卡尔曼滤波器的表贴式永磁同步电机转子磁链观测方法,通过算法形式的观测器代替传统物理上实体形式的参数测量方法,不仅可以实现无传感器监测永磁同步电机的转子磁链,而且结构简单,易操作,成本低,具有较强的自适应能力和抗干扰能力。
Description
技术领域
本发明涉及永磁同步电机技术领域,具体涉及基于扩展卡尔曼滤波器的表贴式永磁同步电机转子磁链观测方法。
背景技术
永磁同步电机(PMSM:Permanent Magnet Synchronous Motor)具有高功率、高能量密度、结构简单、在众多工业领域得到越来越广泛的应用。在高温或大电流条件下永磁同步电机的永磁体将会发生不可逆失磁,这将会导致电机的空载电势下降,性能退化,甚至可能造成电机的损坏,严重影响到电机的效率、性能和可靠性,因此对永磁同步电机的永磁体磁链的实时在线观测显得尤为重要。
而传统的离线参数辨识方法,需要用到另一台反拖被测电机,测量不同转速下的反电势大小,通过磁链与转速的关系计算永磁体磁链,这种方法需要拖动电机、测量设备和专用的连接器,操作复杂,工程代价高,时间消耗多,而且不能实现转子永磁体磁链的实时在线观测。
发明内容
为解决上述现有技术中存在的问题,本发明提供基于扩展卡尔曼滤波器的表贴式永磁同步电机转子磁链观测方法,其具有预测性、自适应能力、抗干扰性、结构简单、易于软件实现等优点。
为了达到上述目的,本发明采用如下技术方案:
基于扩展卡尔曼滤波器的表贴式永磁同步电机转子磁链观测方法,包含以下步骤:
步骤1、获取表贴式永磁同步电机的实际转速与电机三相电流iabc,三相电流iabc经过Clarke变换变换得到静止两相坐标系下电流iα和iβ,通过Park变换得到旋转两相坐标系下id和iq;
步骤2、对电机的给定转速与电机的实际转速的差值进行转速调节器(ASR)调节确定q轴给定电流iq *,d轴电流为给定电流id *;
步骤3、将q轴给定电流iq *与反馈电流iq的差值和d轴给定电流id *与反馈电流id的差值经过电流调节器(ACR)调节后,得到d和q轴实时电压ud、uq;
步骤4、d和q轴实时电压ud、uq经过反PARK变换得到α和β轴实时电压uα、uβ;
步骤5、对α和β轴实时电压uα、uβ进行SVPWM调制得到脉宽调制波形,将脉宽度制波形送入逆变器对永磁同步电机进行控制,得到实时的三相电流iabc和三相电压uabc;
步骤6、将实时电流iα、iβ和实时电压uα、uβ送入基于扩展卡尔曼滤波器(EKF)的转子磁链观测器中,实现永磁同步电机转子磁链的实时观测。
进一步地,所述表贴式永磁同步电机中,转子磁路结构对称,又由于永磁材料的磁导率与气隙的磁导率相接近,于是定子直轴电感Ld等于交轴电感Lq,在两相静止坐标系α-β下,永磁同步电机的电压方程为:
将式(1)变换为电流方程可得:
又由于在实际情况下,电机转子磁链波动时间远远大于电机控制系统动态过程时间,所以转子磁链可以看作稳态值来计算,对其求导可以得到以下方程:
由式(2)、(3)可以写成如下形式:
可以将式(4)写成如下状态方程:
由式(4)、(5)可得:
为构建卡尔曼滤波器状态观测器数字化系统,将式(5)进行离散化处理可得:
整理得:
式中,A'=I+AT为离散化的系数矩阵,B'=BT为离散化的系统输入矩阵,C'=C为离散化的系统输出矩阵;
可以得到离散化的系统矩阵A',B',C'如下:
离散化状态方程(8)是确定性方程,但是在实际系统中,模型参数存在不确定性和可变性,定子电压和电流中不可避免的存在测量噪声,将这些不确定因素纳入到系统噪声矢量V和测量噪声W中,于是式(8)可以改成:
式中,V(k)为系统噪声,W(k)为测量噪声,这两个噪声都是零均值的白噪声。
进一步地,步骤6的具体过程如下:
EKF算法离散化的数学模型为:
式中,x(k+1)表示k+1时刻状态估计值,x(k)表示k时刻状态估计值,V(k)为系统噪声矢量,W(k)为测量噪声,A'=I+AT为离散化的系数矩阵,B'=BT为离散化的系统输入矩阵,C'=C为离散化的系统输出矩阵,T为采样周期;
扩展卡尔曼滤波观测器的状态估计主要分为两个状态:预测和校正
(1)状态预测:
(2)计算协方差矩阵:
(3)计算增益矩阵:
式中,K(k+1)是增益矩阵,主要用来完成对状态矢量估计的校正,R是测量噪声矢量W的协方差矩阵;
(4)状态矢量估计:
(5)计算估计误差协方差矩阵:
与现有技术相比,本发明的有益技术效果为:本发明的基于扩展卡尔曼滤波器的表贴式永磁同步电机转子磁链观测方法,通过算法形式的观测器代替传统物理上实体形式的参数测量方法,通过EKF算法对永磁同步电机的转子磁链实时监测,不仅结构简单,易于操作,成本低;而且具有较强的自适应能力和抗干扰能力。
附图说明
图1为本发明实施例的矢量控制系统示意图;
图2为本发明实施例的基于EKF算法的转子磁链观测流程示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例的附图,对本发明实施例的技术方案进行清楚、完整地描述。显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于所描述的本发明的实施例,本领域普通技术人员所获得的所有其他实施例,都属于本发明保护的范围。
以下结合说明书附图对本发明的技术方案做进一步的详细说明。
如图1所示,基于扩展卡尔曼滤波器的表贴式永磁同步电机转子磁链观测方法,包含以下步骤:
步骤1、获取表贴式永磁同步电机的实际转速与电机三相电流iabc,三相电流iabc经过Clarke变换变换得到静止两相坐标系下电流iα和iβ,通过Park变换得到旋转两相坐标系下id和iq;
步骤2、对电机的给定转速与电机的实际转速的差值进行转速调节器(ASR)调节确定q轴给定电流iq *,d轴电流为给定电流id *;
步骤3、将q轴给定电流iq *与反馈电流iq的差值和d轴给定电流id *与反馈电流id的差值经过电流调节器(ACR)调节后,得到d和q轴实时电压ud、uq;
步骤4、d和q轴实时电压ud、uq经过反PARK变换得到α和β轴实时电压uα、uβ;
步骤5、对α和β轴实时电压uα、uβ进行SVPWM调制得到脉宽调制波形,将脉宽度制波形送入逆变器对永磁同步电机进行控制,得到实时的三相电流iabc和三相电压uabc;
步骤6、将实时电流iα、iβ和实时电压uα、uβ送入基于扩展卡尔曼滤波器(EKF)的转子磁链观测器中,实现永磁同步电机转子磁链的实时观测。
其中,表贴式永磁同步电机中,转子磁路结构对称,又由于永磁材料的磁导率与气隙的磁导率相接近,于是定子直轴电感Ld等于交轴电感Lq,在两相静止坐标系α-β下,永磁同步电机的电压方程为:
将式(1)变换为电流方程可得:
又由于在实际情况下,电机转子磁链波动时间远远大于电机控制系统动态过程时间,所以转子磁链可以看作稳态值来计算,对其求导可以得到以下方程:
由式(2)、(3)可以写成如下形式:
可以将式(4)写成如下状态方程:
由式(4)、(5)可得:
为构建卡尔曼滤波器状态观测器数字化系统,将式(5)进行离散化处理可得:
整理得:
式中,A'=I+AT为离散化的系数矩阵,B'=BT为离散化的系统输入矩阵,C'=C为离散化的系统输出矩阵;
可以得到离散化的系统矩阵A',B',C'如下:
离散化状态方程(8)是确定性方程,但是在实际系统中,模型参数存在不确定性和可变性,定子电压和电流中不可避免的存在测量噪声,将这些不确定因素纳入到系统噪声矢量V和测量噪声W中,于是式(8)可以改成:
式中,V(k)为系统噪声,W(k)为测量噪声,这两个噪声都是零均值的白噪声。
步骤6的具体过程如下:
EKF算法离散化的数学模型为:
式中,x(k+1)表示k+1时刻状态估计值,x(k)表示k时刻状态估计值,V(k)为系统噪声矢量,W(k)为测量噪声,A'=I+AT为离散化的系数矩阵,B'=BT为离散化的系统输入矩阵,C'=C为离散化的系统输出矩阵,T为采样周期;
扩展卡尔曼滤波观测器的状态估计主要分为两个状态:预测和校正
(1)状态预测:
(2)计算协方差矩阵:
(3)计算增益矩阵:
式中,K(k+1)是增益矩阵,主要用来完成对状态矢量估计的校正,R是测量噪声矢量W的协方差矩阵;
(4)状态矢量估计:
(5)计算估计误差协方差矩阵:
本发明的基于扩展卡尔曼滤波器的表贴式永磁同步电机转子磁链观测方法,通过算法形式的观测器代替传统物理上实体形式的参数测量方法,通过EKF算法对永磁同步电机的转子磁链实时监测,不仅结构简单,操作简单易行,成本低,且具有较强的自适应能力和抗干扰能力。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明所述原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (3)
1.基于扩展卡尔曼滤波器的表贴式永磁同步电机转子磁链观测方法,其特征在于,包含以下步骤:
步骤1、获取表贴式永磁同步电机的实际转速与电机三相电流iabc,三相电流iabc经过Clarke变换变换得到静止两相坐标系下电流iα和iβ,通过Park变换得到旋转两相坐标系下id和iq;
步骤2、对电机的给定转速与电机的实际转速的差值进行转速调节器(ASR)调节确定q轴给定电流iq *,d轴电流为给定电流id *;
步骤3、将q轴给定电流iq *与反馈电流iq的差值和d轴给定电流id *与反馈电流id的差值经过电流调节器(ACR)调节后,得到d和q轴实时电压ud、uq;
步骤4、d和q轴实时电压ud、uq经过反PARK变换得到α和β轴实时电压uα、uβ;
步骤5、对α和β轴实时电压uα、uβ进行SVPWM调制得到脉宽调制波形,将脉宽度制波形送入逆变器对永磁同步电机进行控制,得到实时的三相电流iabc和三相电压uabc;
步骤6、将实时电流iα、iβ和实时电压uα、uβ送入基于扩展卡尔曼滤波器(EKF)的转子磁链观测器中,实现对表贴式永磁同步电机转子磁链的实时观测。
2.根据权利要求1所述的基于扩展卡尔曼滤波器的表贴式永磁同步电机转子磁链观测方法,其特征在于,所述表贴式永磁同步电机中,转子磁路结构对称,又由于永磁材料的磁导率与气隙的磁导率相接近,于是定子直轴电感Ld等于交轴电感Lq,在两相静止坐标系α-β下,永磁同步电机的电压方程为:
将式(1)变换为电流方程可得:
又由于在实际情况下,电机转子磁链波动时间远远大于电机控制系统动态过程时间,所以转子磁链可以看作稳态值来计算,对其求导可以得到以下方程:
由式(2)、(3)可以写成如下形式:
可以将式(4)写成如下状态方程:
由式(4)、(5)可得:
为构建卡尔曼滤波器状态观测器数字化系统,将式(5)进行离散化处理可得:
整理得:
x(k+1)=A'x(k)+B'u(k) (8)
y(k)=C'x(k)
式中,A'=I+AT为离散化的系数矩阵,B'=BT为离散化的系统输入矩阵,C'=C为离散化的系统输出矩阵;
可以得到离散化的系统矩阵A',B',C'如下:
离散化状态方程(8)是确定性方程,但是在实际系统中,模型参数存在不确定性和可变性,定子电压和电流中不可避免的存在测量噪声,将这些不确定因素纳入到系统噪声矢量V和测量噪声W中,于是式(8)可以改成:
式中,V(k)为系统噪声,W(k)为测量噪声,这两个噪声都是零均值的白噪声。
3.根据权利要求1所述的基于扩展卡尔曼滤波器的表贴式永磁同步电机转子磁链观测方法,其特征在于,步骤6的具体过程如下:
EKF算法离散化的数学模型为:
式中,x(k+1)表示k+1时刻状态估计值,x(k)表示k时刻状态估计值,V(k)为系统噪声矢量,W(k)为测量噪声,A'=I+AT为离散化的系数矩阵,B'=BT为离散化的系统输入矩阵,C'=C为离散化的系统输出矩阵,T为采样周期;
扩展卡尔曼滤波观测器的状态估计主要分为两个状态:预测和校正
(1)状态预测:
(2)计算协方差矩阵:
(3)计算增益矩阵:
式中,K(k+1)是增益矩阵,主要用来完成对状态矢量估计的校正,R是测量噪声矢量W的协方差矩阵;
(4)状态矢量估计:
(5)计算估计误差协方差矩阵:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010062837.4A CN111181458A (zh) | 2020-01-20 | 2020-01-20 | 基于扩展卡尔曼滤波器的表贴式永磁同步电机转子磁链观测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010062837.4A CN111181458A (zh) | 2020-01-20 | 2020-01-20 | 基于扩展卡尔曼滤波器的表贴式永磁同步电机转子磁链观测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111181458A true CN111181458A (zh) | 2020-05-19 |
Family
ID=70651172
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010062837.4A Pending CN111181458A (zh) | 2020-01-20 | 2020-01-20 | 基于扩展卡尔曼滤波器的表贴式永磁同步电机转子磁链观测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111181458A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112671284A (zh) * | 2020-10-28 | 2021-04-16 | 珠海格力电器股份有限公司 | 一种永磁同步电机、控制方法、控制器及存储介质 |
CN113037165A (zh) * | 2021-03-12 | 2021-06-25 | 上海金脉电子科技有限公司 | 一种永磁同步电机磁链系数的修正方法及装置 |
CN113078852A (zh) * | 2021-04-13 | 2021-07-06 | 东南大学 | 一种永磁同步电机参数的实时辨识方法 |
CN114722528A (zh) * | 2022-03-31 | 2022-07-08 | 中国矿业大学 | 一种基于自适应扩展卡尔曼滤波的感应电机无传感器矢量控制方法 |
CN114744941A (zh) * | 2022-06-09 | 2022-07-12 | 浙江大学 | 一种永磁同步电机永磁体退磁在线监测方法及系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110138301A (zh) * | 2019-04-29 | 2019-08-16 | 江苏大学 | 一种基于扩展卡尔曼滤器的永磁同步电机无位置控制方法 |
-
2020
- 2020-01-20 CN CN202010062837.4A patent/CN111181458A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110138301A (zh) * | 2019-04-29 | 2019-08-16 | 江苏大学 | 一种基于扩展卡尔曼滤器的永磁同步电机无位置控制方法 |
Non-Patent Citations (1)
Title |
---|
胡跃伟: "《基于转子磁链观测器的可控磁通永磁同步电机矢量控制研究》", 《中国优秀硕士论文全文数据库(电子期刊)》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112671284A (zh) * | 2020-10-28 | 2021-04-16 | 珠海格力电器股份有限公司 | 一种永磁同步电机、控制方法、控制器及存储介质 |
CN112671284B (zh) * | 2020-10-28 | 2024-04-16 | 珠海格力电器股份有限公司 | 一种永磁同步电机、控制方法及控制器 |
CN113037165A (zh) * | 2021-03-12 | 2021-06-25 | 上海金脉电子科技有限公司 | 一种永磁同步电机磁链系数的修正方法及装置 |
CN113037165B (zh) * | 2021-03-12 | 2023-01-03 | 上海金脉电子科技有限公司 | 一种永磁同步电机磁链系数的修正方法及装置 |
CN113078852A (zh) * | 2021-04-13 | 2021-07-06 | 东南大学 | 一种永磁同步电机参数的实时辨识方法 |
CN114722528A (zh) * | 2022-03-31 | 2022-07-08 | 中国矿业大学 | 一种基于自适应扩展卡尔曼滤波的感应电机无传感器矢量控制方法 |
CN114744941A (zh) * | 2022-06-09 | 2022-07-12 | 浙江大学 | 一种永磁同步电机永磁体退磁在线监测方法及系统 |
CN114744941B (zh) * | 2022-06-09 | 2022-09-23 | 浙江大学 | 一种永磁同步电机永磁体退磁在线监测方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110198150B (zh) | 一种永磁同步电机多参数在线辨识方法 | |
Tsuji et al. | A sensorless vector control system for induction motors using q-axis flux with stator resistance identification | |
CN110350835B (zh) | 一种永磁同步电机无位置传感器控制方法 | |
CN111181458A (zh) | 基于扩展卡尔曼滤波器的表贴式永磁同步电机转子磁链观测方法 | |
Feng et al. | Hybrid terminal sliding-mode observer design method for a permanent-magnet synchronous motor control system | |
CN111193448B (zh) | 基于扩展卡尔曼滤波器的永磁同步电机负载转矩观测方法 | |
CN110022106B (zh) | 一种基于高频信号注入的永磁同步电机无位置传感器控制方法 | |
CN101149423A (zh) | 永磁同步电机永磁磁场畸变实时检测与分析方法及其装置 | |
CN107994826A (zh) | 一种基于误差加权的全阶观测器无速度传感器控制系统 | |
Sun et al. | Sensorless control technique of PMSM based on RLS on-line parameter identification | |
CN112217428B (zh) | 水下机器人推进系统无位置传感器控制方法 | |
CN114598206B (zh) | 永磁同步电机宽速域转子位置观测器设计方法 | |
CN113691186B (zh) | 一种永磁同步电机无位置传感器控制转子位置角补偿方法 | |
CN111987961A (zh) | 一种永磁同步电机无位置传感器直接转矩控制方法 | |
CN113659904A (zh) | 一种基于非奇异快速终端滑模观测器的spmsm无传感器矢量控制方法 | |
Zhang et al. | An improved sensorless control strategy of ship IPMSM at full speed range | |
Lam et al. | Torque ripple minimization in PM synchronous motors an iterative learning control approach | |
CN114465543B (zh) | 一种永磁同步电机无位置传感器控制方法 | |
Comanescu | Implementation of time-varying observers used in direct field orientation of motor drives by trapezoidal integration | |
Gan et al. | An adaptive nonlinear extended state observer for the sensorless speed control of a PMSM | |
CN113708673B (zh) | 高速开关电机驱动控制方法 | |
CN114465530A (zh) | 永磁同步电机的速度控制方法及系统 | |
CN113890438A (zh) | 一种基于内置式永磁同步电机的无速度传感器控制方法 | |
CN111654223A (zh) | 一种混合励磁磁通切换电机无位置传感器复合控制方法 | |
Shi et al. | Speed sensorless vector control oflnduction motor based on adaptive full-order flux observer |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200519 |
|
RJ01 | Rejection of invention patent application after publication |