具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”、“第三”、“第四”等(如果存在)是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。

下面以具体地实施例对本发明的技术方案进行详细说明。下面这几个具体的实施例可以相互结合，对于相同或相似的概念或过程可能在某些实施例不再赘述。

请参考图1，基于反电动势的最大转矩电流比控制方法，包括：

S101：构建目标关系；

所述目标关系表征了电机处于最大转矩电流比时，所述电机的ABC三相的反电动势与电流的关系；

其中的电流可以理解为在电机的定子中注入的电流，其中的反电动势可以理解为电机的转子上永磁体旋转过程中，最大转矩电流比可以理解为对平均输出转矩与定子电流的比值；

由于磁链在定子绕组上运动引起磁场变化而生成的感应电动势，进而通过实验数据，分析电机处于最大转矩电流比时，反电动势与电路的关系，得到目标关系。

一种实施方式中，所述电机包括以下至少之一：

表贴式永磁同步电机、内嵌式永磁同步电机、开关磁阻电机。

请参考图2，一种实施方式中，步骤S101，包括：

S1011：确定第一关系与第二关系；

所述第一关系表征了所述ABC三相的反电动势的k次谐波的幅值、电机转子电角度以及所述反电动势之间的关系；

所述第二关系表征了所述ABC三相的电流的k次谐波的幅值、电机转子电角度以及所述电流之间的关系；

一种举例中，第一关系可以通过对反电动势进行傅里叶变换，得到反电动势的k次谐波的幅值、电机转子电角度以及所述反电动势之间的关系，即反电动势可以由不用频率和幅值的正弦信号叠加表示；

进一步举例中，第一关系可以表示为：

同理，所述第二关系可以表示为：

其中，

E_a(θ)、E_b(θ)、E_c(θ)表示所述ABC三相的反电动势；

E_k表示所述反向电动势的k次谐波的幅值；

I_a(θ)、I_b(θ)、I_c(θ)表示所述ABC三相的电流；

I_k表示所述电流的k次谐波的幅值；

θp表示所述电流与所述反电动势的相角，即以度的形式表征反电动势与电流的相位差；

θ表示所述电机转子电角度。

S1012：根据所述第一关系和所述第二关系，确定所述电机的平均输出转矩与所述电流的k次谐波的幅值以及所述反电动势的k次谐波的幅值之间的第三关系；

所述平均输出转矩可以理解为平均输出功率与电机的转速的比值，平均输出功率可以通过反电动势与电流得到；

进一步举例中，当所述相角为零时，电机可以保证以最大的平均功率输出，此时，电机的平均输出转矩与所述电流的k次谐波的幅值以及所述反电动势的k次谐波的幅值之间的关系可以表示为：

对上述公式(3)进行简化，可以得到下述第三关系：

其中，

T_ave表征了所述平均输出转矩；

w_m表征了所述电机的转速；

通过上述公式(4)表征的第三关系可以看出，本发明一实施例中的得到是平均输出转矩与反电动势谐波的幅值、电流谐波的幅值之间的关系，而非部分方案中平均输出转矩与反电动势、电流之间的关系，进而在反电动势中包含大量谐波，而非正弦波时，也能够为电机提供最大的输出转矩电流比，保证电机的输出功率达到最大。

S1013：根据所述第二关系，确定第四关系；

所述第四关系表征了所述电流的有效值与所述电流的k次谐波的幅值之间的关系；

进一步举例中，当所述ABC三相的电流的有效值相同时，所述第四关系表征为：

其中，I_rms表征了所述电流的有效值。

S1014：根据所述第三关系以及所述第四关系，确定所述目标关系。

进一步举例中，步骤S1014，包括：

S10141：利用拉格朗日乘数法，以所述电机的总电流值为常数为约束条件，根据所述第三关系以及所述第四关系，求解所述转矩电流比的最大值，得到所述目标关系。

一种举例中，最大转矩电流比可以通过以下公式表征：

其中，MTPA表征了最大转矩电流比；

电机的总电流值为常数可以通过以下公式表示：

其中，C表示所述电机的总电流值。

对于公式(6)、(7)，利用拉格朗日乘数法求的电流与反电动势的关系，以及拉格朗日乘数与反电动势谐波的幅值的关系，进而，在实际应用中，通过对反电动势谐波的分析，得到拉格朗日乘数，再根据反电动势得到电流，实现电机的最大转矩电流比。

一种举例中，所述目标关系表征为：

其中，

γ表示拉格朗日乘数。

以上实施方式中，得到的目标关系可以看出，本发明一实施例中，电流与反电动势以及反电动势的谐波有关，进而，在控制电机的反电动势为梯形波的同时，能够控制电机的输出效率达到最高，且控制方法简单，对MCU的要求也不高。

S102：获取目标电机的反电动势；

S103：根据所述目标电机的反电动势以及所述目标关系，确定所述目标电机的电流。

进一步举例中，步骤S103中，在确定目标电机的电流时，需要先对获取到的目标电机的反电动势经过快速傅里叶变换，得到目标电机的反电动势的多次谐波分量，进而根据得到的谐波分量以及反电动势确定目标电机的电流。

请参考图4，一种实施方式中，所述方法还包括：

S104：根据所述目标电机的电流，控制所述目标电机。

本发明一实施例中还提供了一种电机控制系统，包括控制模块、电机，所述控制模块连接所述电机，所述控制模块用于实现本发明第一方面及其可选方案所述的基于反电动势的最大转矩电流比控制方法。

下面结合图5至图12，详细阐述本发明一实施例中所述方法的积极效果：

图5为现有技术中通过矢量控制法控制电机的原理图，具体控制过程如下：

通过在直(d)轴和交(q)轴给定参考电流(例如图5中的I_{d_Ref}与I_{q_Ref})，将给定参考电流与电机定子中电流采样值进行比较(其中，电流采样值经过CLARK变换得到I_α和I_β，将I_α和I_β经过PARK变换得到I_d和I_q)，产生误差信号，误差信号经过PI控制器产生d轴和q轴参考电压(例如图5中的U_{d_Ref}和U_{q_Ref})，所产生的参考电压经过PARK逆变换生成α-β坐标系下的参考电压(U_{α_Ref}和U_{β_Ref})，最后根据α-β坐标系下的参考电压，空间矢量调制算法可以产生相应的PWM信号(例如图5中的PWM 1-6)驱动逆变器以控制电机，通过选择合理的PI控制器的相关参数，最终电机定子中流过的电流为正弦波。

图5中，电机的平均输出转矩的表达式为：

其中，θ为电机转子电角度，E_x(θ)和I_x(θ)为电机任意一相的反电动势和相电流，w_m为电机机械转速。

由傅里叶级数展开可知，E_x(θ)和I_x(θ)是不同频率和幅值的正弦信号的叠加。

基于正弦函数的正交性，对于反电动势为正弦波的电机(例如图6所示，其中A、B、C对应的正弦波分别为电机三相的反电动势的波形)，此时，在电机的定子中注入正弦电流可以提供最大的输出转矩电流比，使电机效率达到最优；

对于反电动势为非正弦波的电机(例如图7所示，其中A、B、C对应的正弦波分别为电机三相的反电动势的波形)，反电动势中含有大量谐波，若此时在电机中注入频率为基波的正弦电流，并不能提供最大的输出转矩电流比，进而无法使得电机达到最大的输出效率。

本发明一实施例中，通过对反电动势与电流的分析，构建了表征电流与反电动势之间关系的目标关系，即公式(8)所示，进而在对电机注入电流时，通过分析电机的反电动势的谐波，计算出能够达到最大转矩电流比的最优的定子电流，进而实现对电机的控制。

请参考图8，一种举例中，电机转速为1200rpm时，获取到电机的反电动势，对反电动势进行快速傅里叶变换，得到反电动势的多次谐波的分量，可例如图9中所示，

由图9可以看出，反电动势的基波分量的幅值为11.257V，三次谐波分量的幅值为1.334V，五次谐波的分量的幅值为-0.36V，根据公式(8)中表征的目标关系，可以计算得到电机处于最大转矩电流比时的电流，电流的波形可例如图10中I₂对应的波形，I₁对应的波形为正弦波；

由图10可知，电流的基波分量的幅值为0.991A，三次谐波分量的幅值为0.127A，五次谐波的分量的幅值为-0.032A，

图11中为电机运行中反电动势的波形，根据电流与反电动势，可以得到电机的平均输出转矩，平均输出转矩可例如图12中所示，其中，灰度值偏低的波形为在定子中输入正弦波时平均输出转矩对应的波形，此时，平均输出转矩为1.065N.M，灰度值偏高的波形为在定子中输入通过公式(8)得到的电流时，平均输出转矩对应的波形，此时，平均输出转矩为1.086N.M，可见，通过本发明得到的电流注入定子中时，在电流有效值相同的情况下，相比于在定子中注入正弦波，平均输出转矩提升了2％。

需要说明的是，本发明是实施例提供的方法，不仅可以用于对电机的控制，还可以用于电机转矩脉动分析中。

请参考图13，本发明一实施例还提供了一种基于反电动势的最大转矩电流比控制装置2，包括：

关系构建模块201，用于构建目标关系，所述目标关系表征了电机处于最大转矩电流比时，所述电机的ABC三相的反电动势与电流的关系；

数据获取模块202，用于获取目标电机的反电动势；

电流确定模块203，用于根据所述目标电机的反电动势以及所述目标关系，确定所述目标电机的电流。

可选的，关系构建模块201，具体用于：

确定第一关系与第二关系，所述第一关系表征了所述ABC三相的反电动势的k次谐波的幅值、电机转子电角度以及所述反电动势之间的关系；所述第二关系表征了所述ABC三相的电流的k次谐波的幅值、电机转子电角度以及所述电流之间的关系；

根据所述第一关系和所述第二关系，确定所述电机的平均输出转矩与所述电流的k次谐波的幅值以及所述反电动势的k次谐波的幅值之间的第三关系；

根据所述第二关系，确定第四关系，所述第四关系表征了所述电流的有效值与所述电流的k次谐波的幅值之间的关系；

根据所述第三关系以及所述第四关系，确定所述目标关系。

可选的，所述第一关系表征为：

所述第二关系表征为：

其中，

E_a(θ)、E_b(θ)、E_c(θ)表示所述ABC三相的反电动势；

E_k表示所述反向电动势的k次谐波的幅值；

I_a(θ)、I_b(θ)、I_c(θ)表示所述ABC三相的电流；

I_k表示所述电流的k次谐波的幅值；

θp表示所述电流与所述反电动势的相角；

θ表示所述电机转子电角度。

可选的，当所述相角为零时，所述第三关系表征为：

当所述ABC三相的电流的有效值相同时，所述第四关系表征为：

其中，

T_ave表征了所述平均输出转矩；

w_m表征了所述电机的转速；

I_rms表征了所述电流的有效值。

可选的，根据所述第三关系以及所述第四关系，确定所述目标关系，包括：

利用拉格朗日乘数法，以所述电机的总电流值为常数为约束条件，根据所述第三关系以及所述第四关系，求解所述转矩电流比的最大值，得到所述目标关系。

可选的，所述目标关系表征为：

其中，

γ表示拉格朗日乘数；

C表示所述电机的总电流值。

请参考图14，一种实施方式中，所述装置还包括：

电机控制模块204，用于根据所述目标电机的电流，控制所述目标电机。

可选的，所述电机包括以下至少之一：

表贴式永磁同步电机、内嵌式永磁同步电机、开关磁阻电机。

请参考图15，本发明还提供了一种电子设备30，包括

处理器31；以及，

存储器32，用于存储处理器的可执行指令；

其中，处理器31配置为经由执行可执行指令来执行以上所涉及的方法。

处理器31能够通过总线33与存储器32通讯。

本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现以上所涉及的方法。

最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。