CN115347842A - 基于有约束模型预测控制的永磁同步电机电流的控制方法 - Google Patents

基于有约束模型预测控制的永磁同步电机电流的控制方法 Download PDF

Info

Publication number
CN115347842A
CN115347842A CN202211249061.2A CN202211249061A CN115347842A CN 115347842 A CN115347842 A CN 115347842A CN 202211249061 A CN202211249061 A CN 202211249061A CN 115347842 A CN115347842 A CN 115347842A
Authority
CN
China
Prior art keywords
current
permanent magnet
magnet synchronous
region
error signal
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
Application number
CN202211249061.2A
Other languages
English (en)
Other versions
CN115347842B (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.)
Xuanzhi Electronic Technology Shanghai Co ltd
Original Assignee
Xuanzhi Electronic Technology Shanghai Co ltd
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 Xuanzhi Electronic Technology Shanghai Co ltd filed Critical Xuanzhi Electronic Technology Shanghai Co ltd
Priority to CN202211249061.2A priority Critical patent/CN115347842B/zh
Priority to PCT/CN2022/129516 priority patent/WO2024077682A1/zh
Publication of CN115347842A publication Critical patent/CN115347842A/zh
Application granted granted Critical
Publication of CN115347842B publication Critical patent/CN115347842B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/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/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
    • 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
    • H02P25/024—Synchronous motors controlled by supply frequency
    • H02P25/026—Synchronous motors controlled by supply frequency thereby detecting the rotor position
    • 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
    • H02P6/00—Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/34—Modelling or simulation for control purposes
    • 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
    • H02P2205/00—Indexing scheme relating to controlling arrangements characterised by the control loops
    • H02P2205/01—Current loop, i.e. comparison of the motor current with a current reference
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Ac Motors In General (AREA)
  • Control Of Motors That Do Not Use Commutators (AREA)

Abstract

本公开的实施例涉及基于有约束模型预测控制的永磁同步电机电流的控制方法。在该方法中，获取采样时刻下的永磁同步电机的定子电流测量数据；基于定子电流测量数据和定子电流设定值，生成电流跟随误差信号；确定电流跟随误差信号处于电流平面状态分区的六边形内部区域、第一区域还是第二区域。响应于确定电流跟随误差信号在电流平面状态分区中所处的区域为第一区域，将电流跟随误差信号向第一区域的边进行映射以便得到更新后的电流跟随误差信号；以及生成电流环控制电压。该方法根据系统状态给出了当前控制时刻预测控制率的解析表达式，克服了传统有约束预测控制算法在线计算量大的缺陷，从而降低对芯片计算能力的要求。

Description

基于有约束模型预测控制的永磁同步电机电流的控制方法
技术领域
本公开的实施例总体涉及永磁同步电机的电流环控制领域，并且更具体地涉及一种基于有约束模型预测控制的永磁同步电机电流的控制方法。
背景技术
电机将电能转化为机械能，是现代工业中最重要的原动力。电机的耗能占据了工业总耗能的60~70%，其控制效率的提高具有重大意义。
传统上，电机控制包括内环控制和外环控制。其中，内环控制为电流环控制，是对电信号的控制；外环控制是机械控制，是对机械信号的控制。对电信号和机械信号的控制在时间尺度上差别很大。通常情况下，对电信号的控制对应的时间尺度为毫秒级，而对机械信号的控制对应的时间尺度为秒级。因此，将机电控制分解成内外环的递阶控制可以实现对于电机的有效控制。
对于电流环控制问题，传统上主要采用解耦的PI（比例和积分）控制方法。但PI控制方法存在动态响应慢等问题。尤其是当需要实现快速控制的时候，往往会导致PI控制器上产生积分饱和，引起系统超调，导致系统进入稳态很慢，相应的控制带宽难以达到较高的水平。
在电流环控制中采用有约束模型预测控制（MPC）可以实现快速控制。但是，传统MPC需要实时在线求解复杂优化问题，涉及非常大的计算量，为了在很短的采样间隔内完成计算，需要采用高性能芯片，成本较高。工业上，受限于硬件成本，传统MPC通常难以实现高开关频率下的控制。
综上，目前，在电流环控制中采用MPC技术，往往涉及非常大的计算量，导致实现难度大、成本高。
发明内容
针对上述问题，本公开提供了一种基于有约束模型预测控制的永磁同步电机电流的控制方法、计算设备，将预测控制理论应用于永磁同步电机的电流环控制，从而实现了对电机电流的高带宽控制。本公开根据系统状态给出了当前控制时刻预测控制率的解析表达式从而显著减少了对运算资源的消耗，克服了传统有约束预测控制算法在线计算量大的缺陷。因此本公开中设计的算法可以在低成本芯片上实现。
根据本公开的第一方面，提供一种基于有约束模型预测控制的永磁同步电机电流的控制方法。该基于有约束模型预测控制的永磁同步电机电流的控制方法包括：获取采样时刻下的永磁同步电机的定子电流测量数据；基于定子电流测量数据和定子电流设定值，生成电流跟随误差信号；确定电流跟随误差信号在电流平面状态分区中所处的区域为六边形内部区域、第一区域还是第二区域，电流平面状态分区包括六边形内部区域和包围六边形内部区域的第一区域和第二区域响应于确定电流跟随误差信号在电流平面状态分区中所处的区域为第一区域，将电流跟随误差信号向第一区域的边进行映射以便得到更新后的电流跟随误差信号；以及基于更新后的电流跟随误差信号，生成用于控制永磁同步电机的定子电流的电流环控制电压，以使得永磁同步电机的定子电流跟踪参考电流。
在一些实施例中，该基于有约束模型预测控制的永磁同步电机电流的控制方法还包括：响应于确定电流跟随误差信号在电流平面状态分区中所处的区域为第二区域，基于六边形内部区域的对应两条边相交的顶点，得到更新后的电流跟随误差信号，以便用于生成用于控制永磁同步电机的定子电流的电流环控制电压。
在一些实施例中，该基于有约束模型预测控制的永磁同步电机电流的控制方法还包括：响应于确定电流跟随误差信号在电流平面状态分区中所处的区域为六边形内部区域，基于电流跟随误差信号生成用于控制永磁同步电机的定子电流的电流环控制电压，以使得永磁同步电机的定子电流跟踪参考电流。
在一些实施例中，生成用于控制永磁同步电机的定子电流的电流环控制电压包括：基于更新后的电流跟随误差信号，经由永磁同步电机的离散模型对应的电流环控制的解析解，生成用于控制永磁同步电机的定子电流的电流环控制电压。
在一些实施例中，永磁同步电机的离散模型对应的电流环控制的解析解为经由以下方式确定：构建永磁同步电机的离散模型；以定子电流跟随误差作为状态变量，以定子电压跟随误差作为控制输入，对离散模型进行变量代换以将离散模型转换为关于状态变量和控制输入的线性时不变模型；以及根据预定约束条件、目标函数和线性时不变模型，确定永磁同步电机的离散模型对应的电流环控制的解析解。
在一些实施例中，预定约束条件包括关于永磁同步电机的原始电压、永磁同步电机的稳态电压、永磁同步电机的母线电压以及控制输入的惩罚矩阵的约束条件。
在一些实施例中，第一区域包括分别垂直于六边形内部区域的边并向六边形内部区域的外部延伸的六个第一子区域，第二区域包括分别被六个第一子区域分隔的六个第二子区域。
在一些实施例中，定子电流设定值包括永磁同步电机的稳态电流数据。
在一些实施例中，永磁同步电机的离散模型的系统状态矩阵以及输入矩阵均为伸缩旋转矩阵。
根据本公开的第二方面，提供一种计算设备。该计算设备包括：至少一个处理器；以及与至少一个处理器通信连接的存储器；存储器存储有可被至少一个处理器执行的指令，指令被至少一个处理器执行，以使至少一个处理器能够执行根据本公开的第一方面的方法。
根据本公开的第三方面，提供一种计算机可读存储介质。该计算机可读存储介质上存储有计算机程序，计算机程序被机器执行时实现根据本公开的第一方面的方法。
应当理解，本部分所描述的内容并非旨在标识本公开的实施例的关键或重要特征，也不用于限制本公开的范围。本公开的其它特征将通过以下的说明书而变得容易理解。
附图说明
结合附图并参考以下详细说明，本公开各实施例的上述和其他特征、优点及方面将变得更加明显。在附图中，相同或相似的附图标注表示相同或相似的元素。
图1示出了用于实现根据本公开的实施例的基于有约束模型预测控制的永磁同步电机电流的控制方法的计算设备的示意图。
图2示出了本公开的实施例的基于有约束模型预测控制的永磁同步电机电流的控制方法的流程图。
图3示出了本公开的实施例的电流平面状态分区的示意图。
图4示出了本公开的实施例的电流跟随误差信号处于第一区域的示意图。
图5示出了本公开的实施例的电流跟随误差信号处于第二区域的示意图。
图6示出了本公开的实施例的用于确定永磁同步电机的离散模型对应的电流环控制的解析解的方法的流程图。
图7示出了本公开的实施例的电流平面状态分区的一种示意图。
图8示出了本公开的实施例的电流平面状态分区的另一种示意图。
图9示出了本公开的实施例的电流平面状态分区的又一种示意图。
图10示出了可以用来实施本公开内容的实施例的基于有约束模型预测控制的永磁同步电机电流的控制方法的示例电子设备的示意性框图。
具体实施方式
以下结合附图对本公开的示范性实施例做出说明，其中包括本公开实施例的各种细节以助于理解，应当将它们认为仅仅是示范性的。因此，本领域普通技术人员应当认识到，可以对这里描述的实施例做出各种改变和修改，而不会背离本公开的范围和精神。同样，为了清楚和简明，以下的描述中省略了对公知功能和结构的描述。
在本文中使用的术语“包括”及其变形表示开放性包括，即“包括但不限于”。除非特别申明，术语“或”表示“和/或”。术语“基于”表示“至少部分地基于”。术语“一个示例实施例”和“一个实施例”表示“至少一个示例实施例”。术语“另一实施例”表示“至少一个另外的实施例”。术语“第一”、“第二”等等可以指代不同的或相同的对象。下文还可能包括其他明确的和隐含的定义。
如前文所描述，在电流环控制中采用MPC技术，往往涉及非常大的计算量。为了便于说明，以下以表贴式永磁同步电机为例进行说明。应当理解，本公开的技术方案同样适用于内嵌式永磁同步电机等多种类型的永磁同步电机。表贴式永磁同步电机的数学模型如公式（1）所示：
Figure 845226DEST_PATH_IMAGE002
（1）
其中，
Figure 90263DEST_PATH_IMAGE004
，
Figure 666738DEST_PATH_IMAGE006
表征对变量t（表征时间）求导数,i s 表征永磁同步电机的定子电流，u s 表征永磁同步电机的定子电压，R s 表征永磁同步电机的电阻，L s 表征永磁同步电机的电感，
Figure 124264DEST_PATH_IMAGE007
表征永磁同步电机的永磁链系数，
Figure 673057DEST_PATH_IMAGE009
表征永磁同步电机的转速，T s 表征永磁同步电机的采样时间。
根据公式（1）所示的连续动力学方程，在两相旋转坐标系下生成永磁同步电机的离散模型，该离散模型如以下公式（2）所示：
Figure 557836DEST_PATH_IMAGE011
（2）
其中，
Figure 887187DEST_PATH_IMAGE012
，
Figure 682492DEST_PATH_IMAGE013
，
Figure 148108DEST_PATH_IMAGE015
，F表征系统矩阵，B表征输入矩阵，I表征单位矩阵，k表征第k个采样时刻，k为大于或者等于0 的整数。
为了实现电流环控制，需要设计一个电流控制器，即控制定子电流
Figure 203789DEST_PATH_IMAGE016
跟踪参考电流
Figure 426960DEST_PATH_IMAGE017
。在模型预测控制理论的框架下，该跟踪问题可以被描述为如公式（3）所示的优化问题：
Figure 23026DEST_PATH_IMAGE019
Figure 343149DEST_PATH_IMAGE021
（3）
其中，min表征取最小值，
Figure 835310DEST_PATH_IMAGE022
表征受限制于，
Figure 545777DEST_PATH_IMAGE023
表征目标函数，
Figure 679955DEST_PATH_IMAGE024
表征永磁同步电机的稳态电压，
Figure 792268DEST_PATH_IMAGE025
表征采样时刻k的永磁同步电机的定子电流，
Figure 720910DEST_PATH_IMAGE026
表征永磁同步电机的稳态电流，
Figure 509219DEST_PATH_IMAGE028
表征采样时刻k的永磁同步电机的电流跟随误差，Q表征电流跟随误差的权重矩阵，
Figure 650350DEST_PATH_IMAGE030
表征采样时刻k的永磁同步电机的定子电压，
Figure 945065DEST_PATH_IMAGE024
表征永磁同步电机的稳态电压，
Figure 185553DEST_PATH_IMAGE031
表征永磁同步电机的定子电流的初始值，
Figure 995246DEST_PATH_IMAGE033
表征实数域。
Figure 674489DEST_PATH_IMAGE035
为预测步长，
Figure 89290DEST_PATH_IMAGE037
为控制步长。
对公式（3）求解涉及非常大的运算量，为了在很短的采样间隔内完成计算，需要采用高性能芯片，因此实现难度较大，成本较高。
为了至少部分地解决上述问题以及其他潜在问题中的一个或者多个，本公开的示例实施例提出了一种用于控制存储器访问权限的方案。在本公开方案中，首先获取采样时刻下的永磁同步电机的定子电流测量数据，然后，基于该定子电流测量数据和定子电流设定值，生成电流跟随误差信号。通过将定子电流测量数据转换为电流跟随误差信号，可以便于简化运算，快速确定电流跟随误差信号所处的区域。然后，确定电流跟随误差信号在电流平面状态分区中所处的区域为六边形内部区域、第一区域还是第二区域。其中，电流平面状态分区包括六边形内部区域和包围该六边形内部区域的第一区域和第二区域。如果确定该电流跟随误差信号在电流平面状态分区中所处的区域为第一区域，则将该电流跟随误差信号向第一区域的边进行映射以便得到更新后的电流跟随误差信号。通过将电流跟随误差信号向第一区域的边进行映射的方式得到更新后的电流跟随误差信号，可以有效简化获取更新后的电流跟随误差信号的算法，减少运算量，提高效率。然后，基于更新后的电流跟随误差信号，生成用于控制永磁同步电机的定子电流的电流环控制电压，以使得永磁同步电机的定子电流跟踪参考电流。由此，不但可以显著减少获取电流环控制电压过程的运算量，在很短的采样间隔内完成相关计算，而且能够显著减少针对运算资源的占用。
在下文中，将结合附图更详细地描述本方案的具体示例。
图1示出了用于实现根据本公开的实施例的基于有约束模型预测控制的永磁同步电机电流的控制方法的控制系统100的示意图。控制系统100例如包括计算设备120、电流测量装置130、电流环140以及永磁同步电机150。计算设备120例如用于执行根据本公开的实施例的基于有约束模型预测控制的永磁同步电机电流的控制方法，以便生成用于控制永磁同步电机的定子电流的电流环控制电压，以使得永磁同步电机的定子电流跟踪参考电流；电流测量装置130例如用于对永磁同步电机150进行测量，以得到永磁同步电机的定子电流测量数据；电流环140例如用于接收计算设备120输出的电流环控制电压，以便控制永磁同步电机150。
计算设备120可以具有一个或多个处理单元，包括诸如GPU（Graphics ProcessingUnit，图形处理器）、FPGA（Field Programmable Gate Array，现场可编程逻辑门阵列）和ASIC（Application Specific Integrated Circuit，专用集成电路）等的专用处理单元以及诸如CPU（Central Processing Unit，中央处理器）的通用处理单元。另外，在每个计算设备120上也可以运行着一个或多个虚拟机。在一些实施例中，计算设备120例如包括测量数据获取单元102、电流跟随误差信号生成单元104、区域确定单元106、更新信号获取单元108、电流环控制电压生成单元110。应当理解，图1示出的计算设备120仅仅是示例性的，而不应当构成对本公开所描述的实现的功能和范围的任何限制。
关于测量数据获取单元102，其用于获取采样时刻下的永磁同步电机的定子电流测量数据。测量数据获取单元102例如从电流测量装置130处接收测量得到的永磁同步电机的定子电流测量数据。
关于电流跟随误差信号生成单元104，其用于基于定子电流测量数据和定子电流设定值，生成电流跟随误差信号。
关于区域确定单元106，其用于确定电流跟随误差信号在电流平面状态分区中所处的区域为六边形内部区域、第一区域还是第二区域，电流平面状态分区包括六边形内部区域和包围六边形内部区域的第一区域和第二区域。
关于更新信号获取单元108，其用于响应于确定电流跟随误差信号在电流平面状态分区中所处的区域为第一区域，将电流跟随误差信号向第一区域的边进行映射以便得到更新后的电流跟随误差信号。
关于电流环控制电压生成单元110，其用于基于更新后的电流跟随误差信号，生成用于控制永磁同步电机的定子电流的电流环控制电压，以使得永磁同步电机的定子电流跟踪参考电流。
以下结合图2~5，详细说明本公开实施例的基于有约束模型预测控制的永磁同步电机电流的控制方法。图2示出了本公开的实施例的基于有约束模型预测控制的永磁同步电机电流的控制方法200的流程图。方法200可由如图1所示的计算设备120执行，也可以在图10所示的电子设备1000处执行。应当理解的是，方法200还可以包括未示出的附加步骤和/或可以省略所示出的步骤，本公开的范围在此方面不受限制。
在步骤202处，计算设备获取采样时刻下的永磁同步电机的定子电流测量数据。其中，定子电流测量数据例如以
Figure 766259DEST_PATH_IMAGE038
表征。
在步骤204处，计算设备基于定子电流测量数据和定子电流设定值，生成电流跟随误差信号。
在一些实施例中，定子电流设定值例如为永磁同步电机的稳态电流，以
Figure 932DEST_PATH_IMAGE039
表征；电流跟随误差信号
Figure 280603DEST_PATH_IMAGE040
例如为定子电流跟随误差
Figure 222014DEST_PATH_IMAGE042
，即，
Figure 400710DEST_PATH_IMAGE044
。
在步骤206处，计算设备确定电流跟随误差信号在电流平面状态分区中所处的区域为六边形内部区域、第一区域还是第二区域。电流平面状态分区包括六边形内部区域和包围六边形内部区域的第一区域和第二区域。图3示出了本公开的实施例的电流平面状态分区的示意图。应当理解，参照图3所示，电流平面状态分区基于两相旋转坐标系构建。电流平面状态分区包括六边形内部区域H1和包围六边形内部区域H1的第一区域和第二区域。其中，六边形内部区域H1参照以下公式（4）定义：
Figure 388258DEST_PATH_IMAGE046
（4）
其中，
Figure 674883DEST_PATH_IMAGE048
，并且，
Figure 798697DEST_PATH_IMAGE050
以及，
Figure 879785DEST_PATH_IMAGE052
，
Figure 354629DEST_PATH_IMAGE054
。
应当理解，b表征六边形内部区域H1的外接圆的半径，1表征元素全为1的列向量，
Figure 241682DEST_PATH_IMAGE055
表征永磁同步电机的母线电压，
Figure 892106DEST_PATH_IMAGE056
表征输入矩阵B对应的缩放系数，
Figure 144096DEST_PATH_IMAGE057
表征输入矩阵B的特征旋转角度，
Figure 634465DEST_PATH_IMAGE058
表征永磁同步电机的采样时间，R表征惩罚矩阵，
Figure 997313DEST_PATH_IMAGE059
表征永磁同步电机的转子角度，
Figure 33402DEST_PATH_IMAGE061
，r表征关于控制输入的惩罚权重，
Figure 253031DEST_PATH_IMAGE062
表征系统矩阵F的缩放系数，
Figure 233625DEST_PATH_IMAGE063
表征矩阵的逆矩阵
Figure 337847DEST_PATH_IMAGE065
的行列式的平方根，
Figure 290760DEST_PATH_IMAGE066
表征永磁同步电机的稳态电压。可以理解，
Figure 884552DEST_PATH_IMAGE068
的每一行，也即
Figure 352443DEST_PATH_IMAGE070
分别定义了六边形内部区域H1的每一条边。应当理解，在一些实施例中，该六边形为正六边形。
六边形内部区域H1的中心HC由以下公式（5）定义：
Figure 591182DEST_PATH_IMAGE072
（5）
其中，
Figure 664180DEST_PATH_IMAGE073
表征系统矩阵F的逆矩阵。
在六边形内部区域H1中，中心HC到每条边的距离D1如以下公式（6）所示：
Figure 100978DEST_PATH_IMAGE075
（6）
第一区域包括分别垂直于六边形内部区域的边并向六边形内部区域的外部延伸的六个第一子区域，例如，六个第一子区域分别为子区域I、子区域II、子区域III、子区域IV、子区域V以及子区域VI；第二区域包括分别被六个第一子区域分隔的六个第二子区域，例如，六个第二子区域分别为子区域Z1、子区域Z2、子区域Z3、子区域Z4、子区域Z5以及子区域Z6。
定子电流测量数据
Figure 259426DEST_PATH_IMAGE076
具有D轴分量
Figure 564506DEST_PATH_IMAGE077
以及Q轴分量
Figure 226431DEST_PATH_IMAGE078
。相应地，电流跟随误差信号
Figure 834130DEST_PATH_IMAGE079
具有D轴分量
Figure 745454DEST_PATH_IMAGE080
以及Q轴分量
Figure 588645DEST_PATH_IMAGE081
，其中，
Figure 370657DEST_PATH_IMAGE083
，
Figure 414836DEST_PATH_IMAGE085
。
Figure 568385DEST_PATH_IMAGE087
为稳态电流
Figure 684108DEST_PATH_IMAGE089
的D轴分量，
Figure 320626DEST_PATH_IMAGE091
为稳态电流
Figure 660341DEST_PATH_IMAGE088
的Q轴分量。
根据电流跟随误差信号
Figure 749519DEST_PATH_IMAGE092
的D轴分量
Figure 75459DEST_PATH_IMAGE093
以及Q轴分量
Figure 832062DEST_PATH_IMAGE081
，计算设备确定电流跟随误差信号在电流平面状态分区中所处的区域为六边形内部区域、第一区域还是第二区域。
在步骤208处，如果计算设备确定电流跟随误差信号在电流平面状态分区中所处的区域为第一区域，计算设备将电流跟随误差信号向第一区域的边进行映射以便得到更新后的电流跟随误差信号。
图4示出了本公开的实施例的电流跟随误差信号处于第一区域的示意图。参照图4，电流跟随误差信号x处于其中一个子区域II，计算设备确定电流跟随误差信号x在电流平面状态分区中所处的区域为第一区域。因此，计算设备将电流跟随误差信号x向第一区域的边进行映射以便得到更新后的电流跟随误差信号x*。应当理解，将电流跟随误差信号x向第一区域的边进行映射，例如包括：在两相旋转坐标系内，经过电流跟随误差信号x构建第一区域的边s1的垂线，以便获取垂足（例如 ，图中x*对应的点）对应的数据x*作为更新后的电流跟随误差信号。也即，更新后的电流跟随误差信号x*是电流跟随误差信号x沿着子区域II对应的边s1的法向量的垂直投影（即正交投影）。应当理解，通过将电流跟随误差信号向第一区域的边进行映射的方式得到更新后的电流跟随误差信号，可以有效简化获取更新后的电流跟随误差信号的算法，减少运算量，提高效率。
在步骤210处，如果计算设备确定电流跟随误差信号在电流平面状态分区中所处的区域为第二区域，计算设备基于六边形内部区域的对应两条边相交的顶点，得到更新后的电流跟随误差信号。图5示出了本公开的实施例的电流跟随误差信号处于第二区域的示意图。参照图5，电流跟随误差信号x处于其中一个子区域Z2，计算设备确定电流跟随误差信号x在电流平面状态分区中所处的区域为第二区域。因此，计算设备基于六边形内部区域H1的对应两条边相交的顶点，得到更新后的电流跟随误差信号x*。应当理解，更新后的电流跟随误差信号x*对应的顶点p1，为六边形内部区域H1的对应两条边相交的顶点，其中对应两条边分别为六边形内部区域H1的边s1和s2。可以理解，更新后的电流跟随误差信号x*对应的顶点p1为电流跟随误差信号x所处的子区域Z2对应的顶点p1，也即，六边形内部区域H1的与电流跟随误差信号x的距离最近的顶点。应当理解，顶点p1与子区域Z2具有确定的对应关系。而且，在电流平面状态分区构建完成后，顶点p1对应的数值即可确定。因此，计算设备一旦确定电流跟随误差信号x处于子区域Z2，即可根据顶点p1与子区域Z2的对应关系确定更新后的电流跟随误差信号x*的数值，无需进行额外的数据运算，效率极高。
在步骤212处，计算设备基于更新后的电流跟随误差信号，生成用于控制永磁同步电机的定子电流的电流环控制电压，以使得永磁同步电机的定子电流跟踪参考电流。
在一些实施例中，基于更新后的电流跟随误差信号，生成用于控制永磁同步电机的定子电流的电流环控制电压，例如包括：将更新后的电流跟随误差信号x*代入以下公式（7），以得到电流环控制电压：
Figure 77098DEST_PATH_IMAGE095
（7）
其中，
Figure 184732DEST_PATH_IMAGE097
表征电流环控制电压，
Figure 314362DEST_PATH_IMAGE099
表征系统矩阵的缩放系数，
Figure 928401DEST_PATH_IMAGE100
表征矩阵H的逆矩阵
Figure 813181DEST_PATH_IMAGE102
的行列式的平方根，
Figure 939269DEST_PATH_IMAGE104
表征旋转角度，r表征关于控制输入
Figure 200486DEST_PATH_IMAGE105
的惩罚权重，R表征关于控制输入
Figure 400523DEST_PATH_IMAGE105
的惩罚矩阵，
Figure 128307DEST_PATH_IMAGE009
表征永磁同步电机的转速，
Figure 741691DEST_PATH_IMAGE106
表征永磁同步电机的采样时间。关于控制输入
Figure 947545DEST_PATH_IMAGE105
，其例如为定子电压跟随误差，即采样时刻下的永磁同步电机的定子电压测量数据
Figure 64405DEST_PATH_IMAGE108
与永磁同步电机的稳态电压
Figure 556567DEST_PATH_IMAGE109
的差值，即
Figure 326421DEST_PATH_IMAGE111
。
其中，公式（7）为基于模型预测控制的永磁同步电机的离散模型对应的电流环控制的解析解。关于永磁同步电机的离散模型对应的电流环控制的解析解，例如可以经由以下方式确定：构建永磁同步电机的离散模型；以定子电流跟随误差作为状态变量，以定子电压跟随误差作为控制输入，对离散模型进行变量代换以将离散模型转换为关于状态变量和控制输入的线性时不变模型；以及根据预定约束条件、目标函数和线性时不变模型，确定永磁同步电机的离散模型对应的电流环控制的解析解。关于确定永磁同步电机的离散模型对应的电流环控制的解析解的方法，下文将结合图6详细说明，此处不再赘述。
在步骤214处，如果计算设备确定电流跟随误差信号在电流平面状态分区中所处的区域为六边形内部区域，计算设备基于电流跟随误差信号生成用于控制永磁同步电机的定子电流的电流环控制电压，以使得永磁同步电机的定子电流跟踪参考电流。也即，如果电流跟随误差信号所处的区域为六边形内部区域，则计算设备直接基于电流跟随误差信号生成用于控制永磁同步电机的定子电流的电流环控制电压。
在一些实施例中，基于电流跟随误差信号，生成用于控制永磁同步电机的定子电流的电流环控制电压，例如包括：将电流跟随误差信号代入公式（7），以得到电流环控制电压，其中，公式（7）中变量x*的取值为计算设备获取得到的电流跟随误差信号。
根据上述方案，根据定子电流测量数据确定对应的电流跟随误差信号，以便基于电流跟随误差信号确定更新后的电流跟随误差信号。进一步地，在上述方案中，通过构建表征电流矢量轨迹边界的六边形内部区域，以及划分包围六边形内部区域的第一区域和第二区域，可以便于简便、快速地确定电流跟随误差信号所属的区域，以便快速确定更新后的电流跟随误差信号。具体的，当电流跟随误差信号在电流平面状态分区中所处的区域为六边形内部区域，则直接基于电流跟随误差信号通过公式（7）计算得到电流环控制电压；当电流跟随误差信号处于第一区域，则将电流跟随误差信号向第一区域的边进行映射以便得到更新后的电流跟随误差信号，从而基于更新后的电流跟随误差信号，生成用于控制永磁同步电机的定子电流的电流环控制电压，以使得永磁同步电机的定子电流跟踪参考电流；当电流跟随误差信号处于第二区域，则基于六边形内部区域的对应两条边相交的顶点，得到更新后的电流跟随误差信号，以便用于生成用于控制永磁同步电机的定子电流的电流环控制电压。尤其是，当电流跟随误差信号在电流平面状态分区中所处的区域为第一区域时，通过将电流跟随误差信号向第一区域的边进行映射的方式得到更新后的电流跟随误差信号，可以显著降低运算复杂度、减小运算量。综上，上述方案可以避免电流环控制过程中，对电流跟踪问题对应的优化问题求解的复杂过程，以及所涉及的非常大的运算量，可以快速获得用于控制永磁同步电机的定子电流的电流环控制电压。在基于具备同等算力的计算设备的条件下，传统方案需要耗时大约100微秒完成电流环控制所需的运算，而基于本公开的上述方案，耗时降低到大约5微秒，速度提高非常显著。而且，鉴于运算复杂度降低、运算量明显减小，针对运算资源的占用显著减少。
图6示出了本公开的实施例的用于确定永磁同步电机的离散模型对应的电流环控制的解析解的方法600的流程图。方法600可由如图1所示的计算设备120执行，也可以在图10所示的电子设备1000处执行。应当理解的是，方法600还可以包括未示出的附加步骤和/或可以省略所示出的步骤，本公开的范围在此方面不受限制。
在步骤602处，计算设备构建永磁同步电机的离散模型。
在一些实施例中，计算设备构建的永磁同步电机的离散模型如前文提及的公式（2）所示：
Figure 929440DEST_PATH_IMAGE113
（2）
在步骤604处，计算设备以定子电流跟随误差作为状态变量，以定子电压跟随误差作为控制输入，对离散模型进行变量代换以将离散模型转换为关于状态变量和控制输入的线性时不变模型。
在一些实施例中，系统矩阵F可以表示为以下公式（8）：
Figure 307332DEST_PATH_IMAGE115
（8）
其中，
Figure 704816DEST_PATH_IMAGE116
表征系统矩阵F对应的缩放系数，R表征惩罚矩阵。类似地，输入矩阵B可以表示为以下公式（9）：
Figure 230475DEST_PATH_IMAGE118
（9）
其中，
Figure 168344DEST_PATH_IMAGE119
表征输入矩阵B对应的缩放系数，
Figure 931900DEST_PATH_IMAGE120
表征输入矩阵B的特征旋转角度。
从而，系统矩阵F与输入矩阵B均为伸缩旋转矩阵，系统矩阵F与输入矩阵B是可交换的，即，
Figure 31443DEST_PATH_IMAGE121
。
对于公式（2）所示的离散模型，稳态时，该离散模型满足
Figure 44399DEST_PATH_IMAGE122
，即
Figure 726572DEST_PATH_IMAGE124
，从而可以由公式（2）得到如下公式（10）：
Figure 141372DEST_PATH_IMAGE125
（10）
针对公式（10），计算设备以定子电流跟随误差作为状态变量，并以定子电压跟随误差作为控制输入，对公式（10）进行变量代换，从而将该离散模型转换为关于状态变量和控制输入的线性时不变模型。
例如，设置状态变量
Figure 146237DEST_PATH_IMAGE127
，设置控制输入
Figure 646489DEST_PATH_IMAGE129
。经过变量代换后，公式（10）转换成如下公式（11）：
Figure 129423DEST_PATH_IMAGE131
（11）
应当理解，公式（11）表征关于状态变量x k 和控制输入v k 的线性时不变模型。
在步骤606处，计算设备根据预定约束条件、目标函数和线性时不变模型，确定永磁同步电机的离散模型对应的电流环控制的解析解。
关于预定约束条件，其例如包括关于控制输入（即定子电压跟随误差）的约束条件，以及关于状态变量（即定子电流跟随误差）的约束条件。
预定约束条件例如包括关于永磁同步电机的原始电压、永磁同步电机的稳态电压、永磁同步电机的母线电压以及控制输入的惩罚矩阵的约束条件。
关于目标函数，例如包括目标函数
Figure 929889DEST_PATH_IMAGE132
，目标函数
Figure 574497DEST_PATH_IMAGE132
可以参照公式（15）所示的函数，此处不再赘述。
在该永磁同步电机的驱动系统中，原始控制输入
Figure DEST_PATH_IMAGE133
（即永磁同步电机的定子电压测量数据）受到电压六边形的约束。例如，对于预测步长为N p 的模型预测控制器来说，其控制输入
Figure 827623DEST_PATH_IMAGE134
受到如以下公式（12）所示的约束：
Figure 845739DEST_PATH_IMAGE136
（12）
其中，
Figure 235132DEST_PATH_IMAGE138
，
Figure 378538DEST_PATH_IMAGE140
，
Figure 853381DEST_PATH_IMAGE141
；
Figure 412539DEST_PATH_IMAGE143
表征永磁同步电机的稳态电压，R表征控制输入的惩罚矩阵，
Figure 656438DEST_PATH_IMAGE145
表征永磁同步电机的转子角度，
Figure 439587DEST_PATH_IMAGE146
表征控制输入的初始值，
Figure 136147DEST_PATH_IMAGE147
表征永磁同步电机的转速。
在一些实施例中，原始电流
Figure 767504DEST_PATH_IMAGE148
（即永磁同步电机的定子电流测量数据）受到电流圆（用于约束电流的圆形区域）的约束，即
Figure 538014DEST_PATH_IMAGE150
，其中，
Figure 960905DEST_PATH_IMAGE152
表征用于约束原始电流的最大电流值定子电流测量数据。
在一些实施例中，使用
Figure 410341DEST_PATH_IMAGE154
表征的
Figure 576880DEST_PATH_IMAGE155
多边形对电流圆进行近似以构成对原始电流
Figure 592110DEST_PATH_IMAGE156
的约束。其中，
Figure 185902DEST_PATH_IMAGE157
Figure 857055DEST_PATH_IMAGE158
为大于或者等于6的整数，
Figure 499389DEST_PATH_IMAGE158
多边形为正多边形。
Figure 103545DEST_PATH_IMAGE159
的取值越大，则该
Figure 805922DEST_PATH_IMAGE160
多边形越逼近电流圆。
相应地，状态变量
Figure 695862DEST_PATH_IMAGE161
受到的约束如以下公式（13）所示：
Figure 942DEST_PATH_IMAGE162
（13）
其中，
Figure 928446DEST_PATH_IMAGE163
表征永磁同步电机的稳态电流，
Figure DEST_PATH_IMAGE165
，
Figure DEST_PATH_IMAGE167
表征用于约束原始电流
Figure 926358DEST_PATH_IMAGE168
的最大电流值，1表征元素全为1的列向量。
当永磁同步电机的状态变量的初始状态为
Figure DEST_PATH_IMAGE169
，控制输入序列为
Figure DEST_PATH_IMAGE171
时，状态变量
Figure 962316DEST_PATH_IMAGE172
将按照以下公式（14）进行演变：
Figure 946453DEST_PATH_IMAGE174
（14）
其中，
Figure 731394DEST_PATH_IMAGE176
，
Figure 103469DEST_PATH_IMAGE178
，
Figure 236510DEST_PATH_IMAGE180
，
Figure 352234DEST_PATH_IMAGE182
，
即，
Figure DEST_PATH_IMAGE183
。
于是，公式（5）中的目标函数J 0 可以表示为以下公式（15）：
Figure DEST_PATH_IMAGE185
（15）
其中，
Figure 175702DEST_PATH_IMAGE186
，
Figure DEST_PATH_IMAGE187
，r为控制输入
Figure 187520DEST_PATH_IMAGE134
的惩罚权重，
Figure 93945DEST_PATH_IMAGE188
表征输入矩阵B的缩放系数，I表征单位矩阵。其中，惩罚权重r表征实际控制输入和目标值之差在整个目标函数或者代价函数中的占比。惩罚权重r越小，允许的控制输入的变化越大，相应的系统响应动态越快，反之亦然。
相应地，公式（5）所示的优化问题被重新构造为以下公式（16）所示的更新后的优化问题：
Figure 747780DEST_PATH_IMAGE190
（16）
其中，
Figure 176487DEST_PATH_IMAGE192
，
Figure 218261DEST_PATH_IMAGE194
，
Figure 591474DEST_PATH_IMAGE196
。
基于公式（16），可以得到公式（16）所示的优化问题的解析解。
为了便于说明，以下将公式（16）所示的更新后的优化问题称为“优化问题（16）”。
当优化问题（16）无约束时，其解析解如以下公式（17）所示：
Figure DEST_PATH_IMAGE197
（17）
当优化问题（16）无约束时，根据模型预测控制理论，只选取
Figure 517842DEST_PATH_IMAGE198
第一步作为预测控制率用于进行电流环控制，该预测控制率即为基于公式（7）所示的公式得到的控制永磁同步电机的定子电流的电流环控制电压；
Figure 863372DEST_PATH_IMAGE200
（7）
其中，
Figure DEST_PATH_IMAGE201
表征系统矩阵F的缩放系数，
Figure 78977DEST_PATH_IMAGE202
表征矩阵H的逆矩阵
Figure DEST_PATH_IMAGE203
的行列式的平方根，
Figure 673907DEST_PATH_IMAGE204
表征输入矩阵B的特征旋转角度，x*表征更新后的电流跟随误差信号。
可以理解，当优化问题（16）有约束时，本公开根据约束被满足的数量，将电流平面进行了划分，具体如图7所示。图7示出了本公开的实施例的电流平面状态分区的一种示意图。其中，花括号{}中的数字代表优化问题（16）中不等式约束等号成立的约束序号。应当理解，优化问题（16）中的不等式组
Figure 403965DEST_PATH_IMAGE206
包括多个不等式约束，当其中例如第k个不等式约束中等号对应的条件成立时，则称“第m个不等式约束等号成立”，其中对应的“约束序号”为m。例如，如图7所示，区域0代表无约束的情况；区域1代表第1个不等式约束等号成立，因此，区域1中花括号{}中标记的约束序号为1，表征当电流跟随误差信号处于区域1时，第1个不等式约束等号成立；类似地，区域13对应的约束序号包括1和7；楔形区域7对应的约束序号包括1和2。关于其他区域的约束序号，可参照图7所示，此处不再赘述。
可以证明，优化问题（16）在无约束情况下对应的电流平面状态分区（criticalregion）是0号区域对应的六边形内部区域（即六边形内部区域H1）。
当优化问题（16）仅存在一个控制约束时（即仅有一个不等式约束等号成立时），可以证明，总是第一步的控制约束率先被满足。在这种情况下，优化问题（16）的电流平面状态分区是无约束情况下所确定的六边形内部区域沿着其六条边向外的垂直延伸所形成的区域，即图7中的区域1-区域6六个区域。
当优化问题（16）存在两个控制约束时（即具有两个不等式约束等号成立时），会出现两种情况。
在第一种情况中，以区域13所示的矩形区域为例，电流跟随误差信号处于区域13中时，对于由电压六边形的其中一条边所对应的控制约束，在连续两步（即
Figure 135161DEST_PATH_IMAGE207
和
Figure 987579DEST_PATH_IMAGE208
）都可以满足不等式约束等号成立，这条边即为区域1对应的边s3。可以证明，区域13为沿着区域1向电压六边形外部垂直延伸的区域。区域14-区域18五个区域与区域13类似，均属于前述第一种情况，此处不再赘述。区域14-区域18分别为沿着区域2-区域6向电压六边形外部垂直延伸的区域。也即，在第一种情况中，优化问题（16）的电流平面状态分区包括区域13-区域18六个区域。
在第二种情况中，以区域7所示的楔形区域为例，电流跟随误差信号处于区域7中时，对于由电压六边形的其中两条边所对应的控制约束在同一步
Figure 210750DEST_PATH_IMAGE207
可以满足不等式约束等号成立，这两条边为区域1对应的边s3以及区域2对应的边s4。区域8-区域12五个区域与区域7类似，均属于前述第二种情况，此处不再赘述。也即，在第二种情况中，优化问题（16）的电流平面状态分区包括区域7-区域12六个区域。
图8示出了本公开的实施例的电流平面状态分区的另一种示意图。可以证明，当优化问题（16）仅存在一个控制约束时，例如，电流跟随误差信号位于图8中区域2的位置p2时，其对应的受约束情况下的预测控制率，等于区域0的位置p3所对应的无约束情况下的预测控制率。其中位置p3是位置p2沿着六边形内部区域的边s4（即区域2对应的边s4）的法向量的垂直投影。图8中所示区域1、区域3-区域6在这方面与区域2相同，此处不再赘述。
图9示出了本公开的实施例的电流平面状态分区的又一种示意图。可以证明，当优化问题（16）存在两个控制约束时，例如，电流跟随误差信号位于图9中区域14的位置p4时，其对应的受约束情况下的预测控制率，等于区域0的位置p5所对应的无约束情况下的预测控制率，其中位置p5是位置p4沿着六边形内部区域的边s4（即区域14对应的边s4）的法向量的垂直投影。图9中所示区域13、区域15-区域18在这方面与区域14相同，此处不再赘述。
以此类推，当优化问题（16）中的有任意个控制约束时，电流平面状态分区可以被划分成如图3所示的区域分布：内部的电流六边形内部区域区域、第一区域（子区域I、子区域II、子区域III、子区域IV、子区域V以及子区域VI共同构成的区域），以及第二区域（子区域Z1、子区域Z2、子区域Z3、子区域Z4、子区域Z5以及子区域Z6共同构成的区域）。
其中，对于电流跟随误差信号处于第一区域的情况，计算设备将电流跟随误差信号向第一区域的边进行映射以便得到更新后的电流跟随误差信号，并将更新后的电流跟随误差信号代入公式（7）以得到预测控制率，即可得到对应的用于控制永磁同步电机的定子电流的电流环控制电压；对于电流跟随误差信号处于第二区域的情况，计算设备基于六边形内部区域的对应两条边相交的顶点，得到更新后的电流跟随误差信号，并将更新后的电流跟随误差信号代入公式（7）以得到预测控制率，即可得到对应的用于控制永磁同步电机的定子电流的电流环控制电压; 对于电流跟随误差信号处于六边形内部区域的情况，计算设备将电流跟随误差信号直接代入公式（7）（其中变量的取值为计算设备获取得到的电流跟随误差信号），即可得到对应的用于控制永磁同步电机的定子电流的电流环控制电压。
本质上，优化问题（16）是一个复杂的有约束优化问题，对其求解通常需要耗费大量的计算资源。本公开通过分析得到了优化问题（16）的解析解（即公式（7））。基于此，只要根据采样获得的电流数据，即可快速得到有约束优化问题的预测控制率以作为控制永磁同步电机的定子电流的电流环控制电压，用于电流环控制，省略了通过复杂的在线优化算法求解的过程。基于此本公开的技术方案，例如可以在中低端控制芯片上实现有约束预测控制算法，极大改善电机的动态响应性能，达到快速高效的控制性能。
图10示出了可以用来实施本公开内容的实施例的基于有约束模型预测控制的永磁同步电机电流的控制方法的示例电子设备1000的示意性框图。如图所示，电子设备1000包括中央处理单元（即，CPU 1001），其可以根据存储在只读存储器（即，ROM 1002）中的计算机程序指令或者从存储单元1008加载到随机存取存储器（即，RAM 1003）中的计算机程序指令，来执行各种适当的动作和处理。在RAM 1003中，还可存储电子设备1000操作所需的各种程序和数据。CPU 1001、ROM 1002以及RAM 1003通过总线1004彼此相连。输入/输出接口（即，I/O接口1005）也连接至总线1004。
电子设备1000中的多个部件连接至I/O接口1005，包括：输入单元1006，例如键盘、鼠标、麦克风等；输出单元1007，例如各种类型的显示器、扬声器等；存储单元1008，例如磁盘、光盘等；以及通信单元1009，例如网卡、调制解调器、无线通信收发机等。通信单元1009允许电子设备1000通过诸如因特网的计算机网络和/或各种电信网络与其他设备交换信息/数据。
上文所描述的各个过程和处理，例如方法200和600，可由CPU 1001执行。例如，在一些实施例中，方法200和600可被实现为计算机软件程序，其被有形地包含于机器可读介质，例如存储单元1008。在一些实施例中，计算机程序的部分或者全部可以经由ROM 1002和/或通信单元1009而被载入和/或安装到电子设备1000上。当计算机程序被加载到RAM1003并由CPU 1001执行时，可以执行上文描述的方法200和600的一个或多个动作。
本公开涉及方法、装置、系统、电子设备、计算机可读存储介质和/或计算机程序产品。计算机程序产品可以包括用于执行本公开的各个方面的计算机可读程序指令。
计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是――但不限于――电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合适的组合。计算机可读存储介质的更具体的例子（非穷举的列表）包括：便携式计算机盘、硬盘、随机存取存储器（RAM）、只读存储器（ROM）、可擦式可编程只读存储器（EPROM或闪存）、静态随机存取存储器（SRAM）、便携式压缩盘只读存储器（CD-ROM）、数字多功能盘（DVD）、记忆棒、软盘、机械编码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合适的组合。这里所使用的计算机可读存储介质不被解释为瞬时信号本身，诸如无线电波或者其他自由传播的电磁波、通过波导或其他传输媒介传播的电磁波（例如，通过光纤电缆的光脉冲）、或者通过电线传输的电信号。
这里所描述的计算机可读程序指令可以从计算机可读存储介质下载到各个计算/处理设备，或者通过网络、例如因特网、局域网、广域网和/或无线网下载到外部计算机或外部存储设备。网络可以包括铜传输电缆、光纤传输、无线传输、路由器、防火墙、交换机、网关计算机和/或边缘计算设备。每个计算/处理设备中的网络适配卡或者网络接口从网络接收计算机可读程序指令，并转发该计算机可读程序指令，以供存储在各个计算/处理设备中的计算机可读存储介质中。
用于执行本公开操作的计算机程序指令可以是汇编指令、指令集架构（ISA）指令、机器指令、机器相关指令、微代码、固件指令、状态设置数据、或者以一种或多种编程语言的任意组合编写的源代码或目标代码，所述编程语言包括面向对象的编程语言—诸如Smalltalk、C++等，以及常规的过程式编程语言—诸如“C”语言或类似的编程语言。计算机可读程序指令可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络—包括局域网(LAN)或广域网(WAN)—连接到用户计算机，或者，可以连接到外部计算机（例如利用因特网服务提供商来通过因特网连接）。在一些实施例中，通过利用计算机可读程序指令的状态信息来个性化定制电子电路，例如可编程逻辑电路、现场可编程门阵列（FPGA）或可编程逻辑阵列（PLA），该电子电路可以执行计算机可读程序指令，从而实现本公开的各个方面。
这里参照根据本公开实施例的方法、装置（系统）和计算机程序产品的流程图和/或框图描述了本公开的各个方面。应当理解，流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合，都可以由计算机可读程序指令实现。
这些计算机可读程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理单元，从而生产出一种机器，使得这些指令在通过计算机或其它可编程数据处理装置的处理单元执行时，产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。也可以把这些计算机可读程序指令存储在计算机可读存储介质中，这些指令使得计算机、可编程数据处理装置和/或其他设备以特定方式工作，从而，存储有指令的计算机可读介质则包括一个制造品，其包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的各个方面的指令。
也可以把计算机可读程序指令加载到计算机、其它可编程数据处理装置、或其它设备上，使得在计算机、其它可编程数据处理装置或其它设备上执行一系列操作步骤，以产生计算机实现的过程，从而使得在计算机、其它可编程数据处理装置、或其它设备上执行的指令实现流程图和/或框图中的一个或多个方框中规定的功能/动作。
附图中的流程图和框图显示了根据本公开的多个实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或指令的一部分，所述模块、程序段或指令的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。
以上已经描述了本公开的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。

Claims (11)

1.一种基于有约束模型预测控制的永磁同步电机电流的控制方法，其特征在于，包括：
获取采样时刻下的永磁同步电机的定子电流测量数据；
基于所述定子电流测量数据和定子电流设定值，生成电流跟随误差信号；
确定所述电流跟随误差信号在电流平面状态分区中所处的区域为六边形内部区域、第一区域还是第二区域，电流平面状态分区包括六边形内部区域和包围所述六边形内部区域的第一区域和第二区域；
响应于确定所述电流跟随误差信号在电流平面状态分区中所处的区域为第一区域，将所述电流跟随误差信号向第一区域的边进行映射以便得到更新后的电流跟随误差信号；以及
基于更新后的电流跟随误差信号，生成用于控制永磁同步电机的定子电流的电流环控制电压，以使得永磁同步电机的定子电流跟踪参考电流。
2.根据权利要求1所述的方法，其特征在于，还包括：
响应于确定所述电流跟随误差信号在电流平面状态分区中所处的区域为第二区域，基于六边形内部区域的对应两条边相交的顶点，得到更新后的电流跟随误差信号，以便用于生成用于控制永磁同步电机的定子电流的电流环控制电压。
3.根据权利要求1所述的方法，其特征在于，还包括：
响应于确定所述电流跟随误差信号在电流平面状态分区中所处的区域为六边形内部区域，基于所述电流跟随误差信号生成用于控制永磁同步电机的定子电流的电流环控制电压，以使得永磁同步电机的定子电流跟踪参考电流。
4.根据权利要求1所述的方法，其特征在于，其中生成用于控制永磁同步电机的定子电流的电流环控制电压包括：
基于更新后的电流跟随误差信号，经由永磁同步电机的离散模型对应的电流环控制的解析解，生成用于控制永磁同步电机的定子电流的电流环控制电压。
5.根据权利要求4所述的方法，其特征在于，永磁同步电机的离散模型对应的电流环控制的解析解为经由以下方式确定：
构建永磁同步电机的离散模型；
以定子电流跟随误差作为状态变量，以定子电压跟随误差作为控制输入，对离散模型进行变量代换以将离散模型转换为关于所述状态变量和控制输入的线性时不变模型；以及
根据预定约束条件、目标函数和线性时不变模型，确定永磁同步电机的离散模型对应的电流环控制的解析解。
6.根据权利要求5所述的方法，其特征在于，预定约束条件包括关于永磁同步电机的原始电压、永磁同步电机的稳态电压、永磁同步电机的母线电压以及控制输入的惩罚矩阵的约束条件。
7.根据权利要求5所述的方法，其特征在于，第一区域包括分别垂直于所述六边形内部区域的边并向所述六边形内部区域的外部延伸的六个第一子区域，第二区域包括分别被六个第一子区域分隔的六个第二子区域。
8.根据权利要求1所述的方法，其特征在于，定子电流设定值包括永磁同步电机的稳态电流数据。
9.根据权利要求5所述的方法，其特征在于，永磁同步电机的离散模型的系统状态矩阵以及输入矩阵均为伸缩旋转矩阵。
10.一种计算设备，其特征在于，包括：
至少一个处理器；以及
与所述至少一个处理器通信连接的存储器；
所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行根据权利要求1-9中任一项所述的方法。
11.一种计算机可读存储介质，其特征在于，计算机可读存储介质上存储有计算机程序，所述计算机程序被机器执行时实现根据权利要求1-9中任一项所述的方法。
CN202211249061.2A 2022-10-12 2022-10-12 基于有约束模型预测控制的永磁同步电机电流的控制方法 Active CN115347842B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202211249061.2A CN115347842B (zh) 2022-10-12 2022-10-12 基于有约束模型预测控制的永磁同步电机电流的控制方法
PCT/CN2022/129516 WO2024077682A1 (zh) 2022-10-12 2022-11-03 基于有约束模型预测控制的永磁同步电机电流的控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211249061.2A CN115347842B (zh) 2022-10-12 2022-10-12 基于有约束模型预测控制的永磁同步电机电流的控制方法

Publications (2)

Publication Number Publication Date
CN115347842A true CN115347842A (zh) 2022-11-15
CN115347842B CN115347842B (zh) 2023-01-24

Family

ID=83957755

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211249061.2A Active CN115347842B (zh) 2022-10-12 2022-10-12 基于有约束模型预测控制的永磁同步电机电流的控制方法

Country Status (2)

Country Link
CN (1) CN115347842B (zh)
WO (1) WO2024077682A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024077682A1 (zh) * 2022-10-12 2024-04-18 旋智电子科技（上海）有限公司 基于有约束模型预测控制的永磁同步电机电流的控制方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103326598A (zh) * 2013-05-20 2013-09-25 河南师范大学 三相整流器快速模型预测控制方法
CN104242768A (zh) * 2014-09-11 2014-12-24 天津大学 一种应用于多电机控制系统的有限状态模型预测控制方法
CN110557074A (zh) * 2019-08-27 2019-12-10 天津大学 用于级联h桥逆变器永磁电机系统的预测电流控制方法
CN111969914A (zh) * 2020-07-21 2020-11-20 北方工业大学 永磁同步电机无差拍电流预测控制方法及设备、存储介质
WO2022127167A1 (zh) * 2020-12-14 2022-06-23 泛仕达机电股份有限公司 一种抑制单电阻采样永磁同步电机低速噪音的方法
WO2022151609A1 (zh) * 2021-01-12 2022-07-21 东南大学 交替执行采样与控制程序的双三相永磁同步电机控制方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001159521A (ja) * 1999-12-01 2001-06-12 Canon Inc 角速度検出装置
JP3636098B2 (ja) * 2001-06-06 2005-04-06 東芝三菱電機産業システム株式会社 電力変換装置の制御回路
SE536142C2 (sv) * 2011-11-25 2013-05-28 Comsys Ab Aktivt filter för resonansreduktion
CN103532408B (zh) * 2013-09-23 2016-08-17 辽宁工程技术大学 基于电压矢量法的中压cs-apf双滞环容错控制方法
CN115347842B (zh) * 2022-10-12 2023-01-24 旋智电子科技（上海）有限公司 基于有约束模型预测控制的永磁同步电机电流的控制方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103326598A (zh) * 2013-05-20 2013-09-25 河南师范大学 三相整流器快速模型预测控制方法
CN104242768A (zh) * 2014-09-11 2014-12-24 天津大学 一种应用于多电机控制系统的有限状态模型预测控制方法
CN110557074A (zh) * 2019-08-27 2019-12-10 天津大学 用于级联h桥逆变器永磁电机系统的预测电流控制方法
CN111969914A (zh) * 2020-07-21 2020-11-20 北方工业大学 永磁同步电机无差拍电流预测控制方法及设备、存储介质
WO2022127167A1 (zh) * 2020-12-14 2022-06-23 泛仕达机电股份有限公司 一种抑制单电阻采样永磁同步电机低速噪音的方法
WO2022151609A1 (zh) * 2021-01-12 2022-07-21 东南大学 交替执行采样与控制程序的双三相永磁同步电机控制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张志鑫: "永磁同步电机系统转矩脉动抑制的研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024077682A1 (zh) * 2022-10-12 2024-04-18 旋智电子科技（上海）有限公司 基于有约束模型预测控制的永磁同步电机电流的控制方法

Also Published As

Publication number Publication date
WO2024077682A1 (zh) 2024-04-18
CN115347842B (zh) 2023-01-24

Similar Documents

Publication Publication Date Title
Bindal et al. Torque ripple reduction of induction motor using dynamic fuzzy prediction direct torque control
CN115347842B (zh) 基于有约束模型预测控制的永磁同步电机电流的控制方法
CN110112974B (zh) 电机控制方法、控制器、存储介质及电机驱动系统
Ahuja et al. Design of fractional order PID controller for DC motor using evolutionary optimization techniques
Mazaheri et al. Performance evaluation of nonlinear Kalman filtering techniques in low speed brushless DC motors driven sensor-less positioning systems
Dursun et al. Sensorless control application of PMSM with a novel adaptation mechanism
CN114301336A (zh) 一种永磁同步电机直接转矩预测控制方法
CN113364386B (zh) 永磁同步电机的基于强化学习的h∞电流控制方法及系统
CN116502373B (zh) 流道计算网格生成方法及装置
Leuer et al. Fast online model predictive control of IPMSM using parallel computing on FPGA
JP2014220938A (ja) モータ制御装置
CN104156604B (zh) 基于遗传算法的滤波器边界频带控制方法及其装置
Campos et al. Nonlinear speed sensorless control of a surface-mounted PMSM based on a Thau observer
Blanchini et al. Robust constrained model predictive control of fast electromechanical systems
CN115270363B (zh) 基于多代理模型的离心压气机叶片优化设计方法及装置
Thakar et al. Fractional-order PI controller for permanent magnet synchronous motor: A design-based comparative study
Nguyen et al. Asynchronous machine rotor speed estimation using a tabulated numerical approach
Allouche et al. Fuzzy tracking control for indirect field-oriented induction machine using integral action performance
CN114499312A (zh) 一种双三相永磁同步电机无模型电流预测控制方法及系统
CN114584026A (zh) 一种优化永磁同步电机参数鲁棒性的方法及相关设备
CN112398402A (zh) 基于快速有效集的电机模型预测控制方法、装置及控制器
WO2020227891A1 (zh) 电机控制方法、控制器、存储介质及电机驱动系统
CN109039168A (zh) 无刷直流电机的多尺度近似显式模型预测控制方法
Azadi et al. Wavenet based vector control of a permanent magnet synchronous motor drive
CN114337431B (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