基于谐波分解的永磁无刷电机高性能运行控制方法及系统

1.本发明涉及永磁无刷电机控制领域，特别是涉及一种基于谐波分解的永磁无刷电机高性能运行控制方法及系统。


背景技术：

2.永磁无刷电机因其结构简单、响应速度快、调速性能好等优点，在电力传动、工业控制、家用电器以及航空航天等各种场合得到广泛的研究与应用。永磁无刷电机分为永磁同步电机与无刷直流电机，对于理想状态下，永磁同步电机与无刷直流电机转子磁场在气隙中的分布分别为理想正弦波与理想梯形波。然而在实际生产制造过程中，电机本体设计并不能达到理想状态，转子永磁体磁链在气隙圆内的分布是由基波和一系列谐波叠加而成。
3.尽管理想正弦波形驱动方式相较于含谐波分量的非理想波形驱动方式简单，但是在含谐波分量的波形的驱动方式下产生的功率密度、输出转矩要比理想正弦波驱动的高。然而对于含谐波分量的非理想反电势而言，由于谐波分量存在会使得永磁无刷电机运行过程中产生相应的转矩脉动，从而影响电机的使用寿命和运行效率。此外许多高性能的驱动场合要求永磁无刷电机转矩具有较高的平滑度，在这种情况下，转矩脉动的存在就会严重影响整个系统的可靠性、稳定性。


技术实现要素：

4.本发明的目的是提供一种基于谐波分解的永磁无刷电机高性能运行控制方法及系统，以降低谐波分量导致的转矩脉动对电机性能的影响，提高永磁无刷电机的性能。
5.为实现上述目的，本发明提供了如下方案：
6.一种基于谐波分解的永磁无刷电机高性能运行控制方法，包括：
7.根据永磁无刷电机的转子位置，确定所述永磁无刷电机的实际转速；
8.根据所述永磁无刷电机的给定转速和实际转速，通过转速调节器确定合成电流给定值；
9.获取所述永磁无刷电机的性能选择参数；所述性能选择参数为不大于1的非负数，所述性能选择参数根据用户期望注入的谐波电流设定；
10.根据所述合成电流给定值和所述性能选择参数，确定基波电流峰值给定值；
11.根据所述基波电流峰值给定值、所述合成电流给定值和所述性能选择参数，确定五次谐波电流峰值给定值和七次谐波电流峰值给定值；
12.将所述五次谐波电流峰值给定值和所述七次谐波电流峰值给定值注入所述永磁无刷电机的定子电流中，以抑制五次谐波和七次谐波产生的转矩脉动及抑制超出七次的非零序谐波反电势与基波电流产生的转矩脉动。
13.可选的，所述性能选择参数为注入七次谐波的电流与注入总电流之间的比值；所述注入总电流为注入五次谐波的电流与注入七次谐波的电流之和；所述注入五次谐波的电
流和所述注入七次谐波的电流均为用户期望注入的谐波电流。
14.可选的，所述根据所述合成电流给定值和所述性能选择参数，确定基波电流峰值给定值，具体包括：
15.获取基波反电势常数、五次谐波反电势常数和七次谐波反电势常数；
16.根据所述基波反电势常数、所述五次谐波反电势常数和所述七次谐波反电势常数，确定五次谐波反电势倍数、七次谐波反电势倍数和最佳谐波电流倍数；
17.利用公式确定所述基波电流峰值给定值；其中，为基波电流峰值给定值，i
*
为合成电流给定值，α为性能选择参数，k5为五次谐波反电势倍数，k7为七次谐波反电势倍数，k
c
为最佳谐波电流倍数。
18.可选的，所述获取基波反电势常数、五次谐波反电势常数和七次谐波反电势常数，具体包括：
19.通过离线试验或有限元仿真，获得不同转子位置下的任一相反电势常数，所述反电势常数为反电势与永磁无刷电机转速的比值；
20.对由0开始增加的一个电周期的反电势常数与转子位置之间的数据，进行傅里叶分解，得到基波反电势常数、五次谐波反电势常数和七次谐波反电势常数。
21.可选的，所述根据所述基波反电势常数、所述五次谐波反电势常数和所述七次谐波反电势常数，确定五次谐波反电势倍数、七次谐波反电势倍数和最佳谐波电流倍数，具体包括：
22.利用公式确定所述五次谐波反电势倍数；
23.利用公式确定所述七次谐波反电势倍数；
24.式中，k
e1
为基波反电势常数，k
e5
为五次谐波反电势常数，k
e7
为七次谐波反电势常数；s5为与五次谐波反电势常数初始相位相关的系数，当五次谐波反电势常数与基波反电势常数的初相位之差为0时s5＝1，当五次谐波反电势常数与基波反电势常数的初相位之差为π时s5＝0；s7为与七次谐波反电势常数初始相位相关的系数，当七次谐波反电势常数与基波反电势常数的初相位之差为0时s7＝1，初相位之差为π时s7＝0；
25.根据由0开始增加的一个电周期的反电势常数与转子位置之间的数据，获得三相反电势常数k
ea
、k
eb
、k
ec
与转子位置θ之间的关系数据库；
26.根据三相反电势常数k
ea
、k
eb
、k
ec
与转子位置θ之间的关系数据库，确定q轴反电势常数k
eq
与转子位置θ之间的关系数据库，公式为：
27.k
eq
(θ)＝
‑
sinθ
·
k
ea
(θ)
‑
sin(θ
‑
2π/3)
·
k
eb
(θ)
‑
sin(θ+2π/3)
·
k
ec
(θ)
28.式中，k
eq
(θ)表示转子位置为θ时的q轴反电势常数；
29.确定超出七次的非零序q轴谐波反电势常数k
ehq
与转子位置θ之间的关系数据库，公式为：
[0030][0031]
式中，k
ehq
(θ)表示转子位置为θ时超出七次的非零序q轴谐波反电势常数；
[0032]
对于所有转子位置θ，通过仿真确定使|k
ehq
(θ)+k
c
k
e1
cos(6θ)|达到最小值时的k
c
的取值为最佳谐波电流倍数。
[0033]
可选的，所述根据所述基波电流峰值给定值、所述合成电流给定值和所述性能选择参数，确定五次谐波电流峰值给定值和七次谐波电流峰值给定值，具体包括：
[0034]
利用公式确定五次谐波电流峰值给定值；
[0035]
利用公式确定七次谐波电流峰值给定值；
[0036]
式中，表示五次谐波电流峰值给定值，表示七次谐波电流峰值给定值。
[0037]
可选的，其特征在于，所述将所述五次谐波电流峰值给定值和所述七次谐波电流峰值给定值注入所述永磁无刷电机的定子电流中，具体包括：
[0038]
根据所述基波电流峰值给定值、所述五次谐波电流峰值给定值、所述七次谐波电流峰值和所述转子位置，确定三相电流瞬时值的给定值；公式为：
[0039][0040]
式中，和为所述三相电流瞬时值的给定值；θ为转子位置；为基波电流峰值给定值，为五次谐波电流峰值给定值，为七次谐波电流峰值给定值；
[0041]
将三相电流瞬时值的给定值分别与三相电流的实际值作差，经电流调节器确定三相电压瞬时值的给定值；
[0042]
将所述三相电压瞬时值的给定值经过pwm生成单元输入到驱动所述永磁无刷电机运行的功率变换电路中。
[0043]
本发明还提供一种基于谐波分解的永磁无刷电机高性能运行控制系统，包括：
[0044]
转速确定模块，用于根据永磁无刷电机的转子位置，确定所述永磁无刷电机的实际转速；
[0045]
合成电流确定模块，用于根据所述永磁无刷电机的给定转速和实际转速，通过转速调节器确定合成电流给定值；
[0046]
性能选择参数获取模块，用于获取所述永磁无刷电机的性能选择参数；所述性能选择参数为不大于1的非负数，所述性能选择参数根据用户期望注入的谐波电流设定；
[0047]
基波电流峰值给定值确定模块，用于根据所述合成电流给定值和所述性能选择参数，确定基波电流峰值给定值；
[0048]
谐波电流峰值给定值确定模块，用于根据所述基波电流峰值给定值、所述合成电流给定值和所述性能选择参数，确定五次谐波电流峰值给定值和七次谐波电流峰值给定值；
[0049]
谐波电流注入模块，用于将所述五次谐波电流峰值给定值和所述七次谐波电流峰值给定值注入所述永磁无刷电机的定子电流中，以抑制五次谐波和七次谐波产生的转矩脉动及抑制超出七次的非零序谐波反电势与基波电流产生的转矩脉动。
[0050]
可选的，所述基波电流峰值给定值确定模块具体包括：
[0051]
反电势常数获取单元，用于获取基波反电势常数、五次谐波反电势常数和七次谐波反电势常数；
[0052]
谐波反电势倍数确定单元，用于根据所述基波反电势常数、所述五次谐波反电势常数和所述七次谐波反电势常数，确定五次谐波反电势倍数、七次谐波反电势倍数和最佳谐波电流倍数；
[0053]
基波电流峰值确定单元，用于利用公式确定所述基波电流峰值给定值；其中，为基波电流峰值给定值，i
*
为合成电流给定值，α为性能选择参数，k5为五次谐波反电势倍数，k7为七次谐波反电势倍数，k
c
为最佳谐波电流倍数；
[0054]
所述反电势常数获取单元，具体包括：
[0055]
反电势常数获取子单元，用于通过离线试验或有限元仿真，获得不同转子位置下的任一相反电势常数，所述反电势常数为反电势与永磁无刷电机转速的比值；
[0056]
傅里叶分解子单元，用于对由0开始增加的一个电周期的反电势常数与转子位置之间的数据，进行傅里叶分解，得到基波反电势常数、五次谐波反电势常数和七次谐波反电势常数；
[0057]
所述谐波反电势倍数确定单元，具体包括：
[0058]
五次谐波反电势倍数计算子单元，用于利用公式确定所述五次谐波反电势倍数；
[0059]
七次谐波反电势倍数计算子单元，用于利用公式确定所述七次谐波反电势倍数；
[0060]
式中，k
e1
为基波反电势常数，k
e5
为五次谐波反电势常数，k
e7
为七次谐波反电势常数；s5为与五次谐波反电势常数初始相位相关的系数，当五次谐波反电势常数与基波反电势常数的初相位之差为0时s5＝1，当五次谐波反电势常数与基波反电势常数的初相位之差为π时s5＝0；s7为与七次谐波反电势常数初始相位相关的系数，当七次谐波反电势常数与基波反电势常数的初相位之差为0时s7＝1，初相位之差为π时s7＝0；
[0061]
三相反电势常数与转子位置关系数据库获取子单元，用于根据由0开始增加的一个电周期的反电势常数与转子位置之间的数据，获得三相反电势常数k
ea
、k
eb
、k
ec
与转子位置θ之间的关系数据库；
[0062]
q轴反电势常数与转子位置关系数据库确定子单元，用于根据三相反电势常数k
ea
、k
eb
、k
ec
与转子位置θ之间的关系数据库，确定q轴反电势常数k
eq
与转子位置θ之间的关系数据库，公式为：
[0063]
k
eq
(θ)＝
‑
sinθ
·
k
ea
(θ)
‑
sin(θ
‑
2π/3)
·
k
eb
(θ)
‑
sin(θ+2π/3)
·
k
ec
(θ)
[0064]
式中，k
eq
(θ)表示转子位置为θ时的q轴反电势常数；
[0065]
超出七次的非零序q轴谐波反电势常数与转子位置关系数据库确定子单元，用于确定超出七次的非零序q轴谐波反电势常数k
ehq
与转子位置θ之间的关系数据库，公式为：
[0066][0067]
式中，k
ehq
(θ)表示转子位置为θ时超出七次的非零序q轴谐波反电势常数；
[0068]
最佳谐波电流倍数确定子单元，用于对于所有转子位置θ，通过仿真确定使|k
ehq
(θ)+k
c
k
e1
cos(6θ)|达到最小值时的k
c
的取值为最佳谐波电流倍数。
[0069]
可选的，所述谐波电流峰值给定值确定模块，具体包括：
[0070]
五次谐波电流峰值给定值计算单元，用于利用公式确定五次谐波电流峰值给定值；
[0071]
七次谐波电流峰值给定值计算单元，用于利用公式确定七次谐波电流峰值给定值；
[0072]
式中，表示五次谐波电流峰值给定值，表示七次谐波电流峰值给定值。
[0073]
可选的，所述谐波电流注入模块，具体包括：
[0074]
三相电流瞬时值确定单元，用于根据所述基波电流峰值给定值、所述五次谐波电流峰值给定值、所述七次谐波电流峰值和所述转子位置，确定三相电流瞬时值的给定值；公式为：
[0075][0076]
式中，和为所述三相电流瞬时值的给定值；θ为转子位置；为基波电流峰值给定值，为五次谐波电流峰值给定值，为七次谐波电流峰值给定值；
[0077]
三相电压瞬时值确定单元，用于将三相电流瞬时值的给定值分别与三相电流的实际值作差，经电流调节器确定三相电压瞬时值的给定值；
[0078]
永磁无刷电机驱动单元，用于将所述三相电压瞬时值的给定值经过pwm生成单元输入到驱动所述永磁无刷电机运行的功率变换电路中。
[0079]
根据本发明提供的具体实施例，本发明公开了以下技术效果：
[0080]
本发明针对非理想反电势作用下产生的转矩脉动进行抑制，研究了要使转矩脉动最小时反电势与相电流谐波分量间应满足的条件，并有效的利用谐波分量来抑制永磁无刷电机运行过程中产生的谐波转矩脉动，此外能够在抑制转矩脉动的同时根据性能选择参数α实现铜耗最小、无直轴电枢反应等高性能运行。本发明为非理想反电势情况下的电磁转矩特性提供了一种研究手段，为永磁无刷电机的转矩脉动抑制提供了一种新的方法，控制方法简单，谐波抑制效果良好，适用性强且无需增设外围电路。
附图说明
[0081]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0082]
图1为本发明基于谐波分解的永磁无刷电机高性能运行控制方法的流程示意图；
[0083]
图2为本发明具体实施案例中控制系统的结构框图；
[0084]
图3为本发明具体实施案例中基于谐波分解的永磁无刷电机高性能运行控制的框图；
[0085]
图4为本发明具体实施案例中三相反电势波形为理想梯形波时的反电势图；其中，(a)部分为a相反电势波形图，(b)部分为经傅里叶分解后的反电势频谱图，(c)部分为同步旋转d
‑
q坐标系下的反电势空间矢量图；
[0086]
图5为本发明具体实施案例中同步旋转d
‑
q坐标系下反电势和电流的空间矢量图；其中，(a)部分为五次、七次谐波与基波反电势反向空间矢量图，(b)部分为五次谐波与基波反电势同向空间矢量图，(c)部分为七次谐波与基波反电势同向空间矢量图，(d)部分为五次、七次谐波与基波反电势同向空间矢量图；
[0087]
图6为本发明具体实施案例中设定α＝0.5时的实验图；其中，(a)部分为a相电流及其基波、五次谐波、七次谐波分量波形，(b)部分为电磁转矩波形，(c)部分为转子永磁体产生的主磁通及其与定子电流产生的磁通进行合成后的气隙合成磁通轨迹。
具体实施方式
[0088]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0089]
本发明的目的是提供一种基于谐波分解的永磁无刷电机高性能运行控制方法及系统，通过利用谐波分量，将含有特定谐波分量的电流注入到永磁无刷电机的定子电流中来抑制谐波分量产生的转矩脉动，与此同时用户可自行通过对性能选择参数α的选择来实现永磁无刷电机的无直轴电枢反应、铜耗最小等高性能运行。
[0090]
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
[0091]
图1为本发明基于谐波分解的永磁无刷电机高性能运行控制方法的流程示意图，如图1所示，本发明基于谐波分解的永磁无刷电机高性能运行控制方法包括以下步骤：
[0092]
步骤100：根据永磁无刷电机的转子位置，确定永磁无刷电机的实际转速。
[0093]
步骤200：根据永磁无刷电机的给定转速和实际转速，通过转速调节器确定合成电流给定值。将给定转速n
*
与转速实际值n作差，经过转速调节器asr可以得到合成电流给定值i
*
。合成电流与给定电磁转矩之间的关系为：
[0094][0095]
式中，t
*
为电机给定电磁转矩；i
*
为合成电流；k
e1
为基波反电势常数，即基波反电势峰值e1与转速n的比值。
[0096]
步骤300：获取永磁无刷电机的性能选择参数。性能选择参数为不大于1的非负数，由用户根据对诸如效率、转矩电流比、铜耗及电枢反应等性能方面的需求而设定。具体的，性能选择参数为注入七次谐波的电流与注入总电流之间的比值，即i7为
[0109]
式中，i1＞0时，表示基波电流与基波反电势同相；反之，i1＜0时，表示基波电流与基波反电势反相；时，表示五次谐波电流与五次谐波反电势同相；反之，时，表示五次谐波电流与五次谐波反电势反相；时，表示七次谐波电流与七次谐波反电势同相；反之，时，表示七次谐波电流与七次谐波反电势反相。
[0110]
而且，注入的五次谐波电流峰值i5在注入过程中，分为i
50
与i
51
两个分量，七次谐波电流峰值i7在注入过程中，分为i
70
与i
71
两个分量。其中，i
50
与i
70
用于完全抵消五、七次谐波反电势与基波电流产生的6次转矩脉动，而i
51
和i
71
用于抑制超出七次的非零序谐波反电势与基波电流产生的转矩脉动，即：
[0111][0112]
此时，最佳谐波电流倍数k
c
等于i
51
+i
71
与基波电流峰值i1之比，具体推导过程在后续具体实施案例中进一步说明。
[0113]
基于上述方案，本发明还提供一种基于谐波分解的永磁无刷电机高性能运行控制系统，包括：
[0114]
转速确定模块，用于根据永磁无刷电机的转子位置，确定所述永磁无刷电机的实际转速。
[0115]
合成电流确定模块，用于根据所述永磁无刷电机的给定转速和实际转速，通过转速调节器确定合成电流给定值。
[0116]
性能选择参数获取模块，用于获取所述永磁无刷电机的性能选择参数；所述性能选择参数为不大于1的非负数。
[0117]
基波电流峰值给定值确定模块，用于根据所述合成电流给定值和所述性能选择参数，确定基波电流峰值给定值。
[0118]
谐波电流峰值给定值确定模块，用于根据所述基波电流峰值给定值、所述合成电流给定值和所述性能选择参数，确定五次谐波电流峰值给定值和七次谐波电流峰值给定值。
[0119]
谐波电流注入模块，用于将所述五次谐波电流峰值给定值和所述七次谐波电流峰值给定值注入所述永磁无刷电机的定子电流中，以抑制五次谐波和七次谐波产生的转矩脉动。
[0120]
作为具体实施例，本发明的基于谐波分解的永磁无刷电机高性能运行控制系统中，所述基波电流峰值给定值确定模块具体包括：
[0121]
反电势常数获取单元，用于获取基波反电势常数、五次谐波反电势常数和七次谐波反电势常数。
[0122]
谐波反电势倍数确定单元，用于根据所述基波反电势常数、所述五次谐波反电势常数和所述七次谐波反电势常数，确定五次谐波反电势倍数、七次谐波反电势倍数和最佳谐波电流倍数。
[0123]
基波电流峰值确定单元，用于利用公式确定所述基波电流峰值给定值；其中，为基波电流峰值给定值，i
*
为合成电流给定值，α为性能选择参数，k5为五次谐波反电势倍数，k7为七次谐波反电势倍数，k
c
为最佳谐波电流倍数。
[0124]
所述反电势常数获取单元，具体包括：
[0125]
反电势常数获取子单元，用于通过离线试验或有限元仿真，获得不同转子位置下的任一相反电势常数，所述反电势常数为反电势与永磁无刷电机转速的比值。
[0126]
傅里叶分解子单元，用于对由0开始增加的一个电周期的反电势常数与转子位置之间的数据，进行傅里叶分解，得到基波反电势常数、五次谐波反电势常数和七次谐波反电势常数。
[0127]
所述谐波反电势倍数确定单元，具体包括：
[0128]
五次谐波反电势倍数计算子单元，用于利用公式确定所述五次谐波反电势倍数。
[0129]
七次谐波反电势倍数计算子单元，用于利用公式确定所述七次谐波反电势倍数。
[0130]
式中，k
e1
为基波反电势常数，k
e5
为五次谐波反电势常数，k
e7
为七次谐波反电势常数；s5为与五次谐波反电势常数初始相位相关的系数，当五次谐波反电势常数与基波反电势常数的初相位之差为0时s5＝1，当五次谐波反电势常数与基波反电势常数的初相位之差为π时s5＝0；s7为与七次谐波反电势常数初始相位相关的系数，当七次谐波反电势常数与基波反电势常数的初相位之差为0时s7＝1，初相位之差为π时s7＝0。
[0131]
三相反电势常数与转子位置关系数据库获取子单元，用于根据由0开始增加的一个电周期的反电势常数与转子位置之间的数据，获得三相反电势常数k
ea
、k
eb
、k
ec
与转子位置θ之间的关系数据库；
[0132]
q轴反电势常数与转子位置关系数据库确定子单元，用于根据三相反电势常数k
ea
、k
eb
、k
ec
与转子位置θ之间的关系数据库，确定q轴反电势常数k
eq
与转子位置θ之间的关系数据库，公式为：
[0133]
k
eq
(θ)＝
‑
sinθ
·
k
ea
(θ)
‑
sin(θ
‑
2π/3)
·
k
eb
(θ)
‑
sin(θ+2π/3)
·
k
ec
(θ)
[0134]
式中，k
eq
(θ)表示转子位置为θ时的q轴反电势常数；
[0135]
超出七次的非零序q轴谐波反电势常数与转子位置关系数据库确定子单元，用于确定超出七次的非零序q轴谐波反电势常数k
ehq
与转子位置θ之间的关系数据库，公式为：
[0136][0137]
式中，k
ehq
(θ)表示转子位置为θ时超出七次的非零序q轴谐波反电势常数；
[0138]
最佳谐波电流倍数确定子单元，用于对于所有转子位置θ，通过仿真确定使|k
ehq
(θ)+k
c
k
e1
cos(6θ)|达到最小值时的k
c
的取值为最佳谐波电流倍数。
[0139]
作为具体实施例，本发明的基于谐波分解的永磁无刷电机高性能运行控制系统中，所述谐波电流峰值给定值确定模块，具体包括：
[0140]
五次谐波电流峰值给定值计算单元，用于利用公式确定五次谐波电流峰值给定值。
[0141]
七次谐波电流峰值给定值计算单元，用于利用公式确定七次谐波电流峰值给定值。
[0142]
式中，表示五次谐波电流峰值给定值，表示七次谐波电流峰值给定值。
[0143]
作为具体实施例，本发明的基于谐波分解的永磁无刷电机高性能运行控制系统中，所述谐波电流注入模块，具体包括：
[0144]
三相电流瞬时值确定单元，用于根据所述基波电流峰值给定值、所述五次谐波电流峰值给定值、所述七次谐波电流峰值和所述转子位置，确定三相电流瞬时值的给定值；公式为：
[0145][0146]
式中，和为所述三相电流瞬时值的给定值；为基波电流峰值给定值，为五次谐波电流峰值给定值，为七次谐波电流峰值给定值。
[0147]
三相电压瞬时值确定单元，用于将三相电流瞬时值的给定值分别与三相电流的实际值作差，经电流调节器确定三相电压瞬时值的给定值。
[0148]
永磁无刷电机驱动单元，用于将所述三相电压瞬时值的给定值经过pwm生成单元输入到驱动所述永磁无刷电机运行的功率变换电路中。
[0149]
下面提供一个具体实施案例来进一步说明本发明的上述方案。图2为本发明具体实施案例中控制系统的结构框图，图3为本发明具体实施案例中基于谐波分解的永磁无刷电机高性能运行控制控制框图，结合图2
‑
图3所示，本实施案例包括以下步骤：
[0150]
(1)通过位置传感器获得永磁无刷电机的转子位置θ，转子位置θ经过转速计算单元得到永磁无刷电机的实际转速n；
[0151]
(2)将转速给定值n
*
与转速实际值n作差，经过转速调节器asr得到合成电流给定值i
*
。合成电流与给定电磁转矩之间的关系为：
[0152][0153]
由此可知合成电流可作为转速调节器的输出。
[0154]
(3)根据合成电流给定值i
*
与性能选择参数α，计算出基波电流峰值给定值以及五、七次谐波电流峰值给定值此处公式与前述实施例中公式相同，不再赘述。
[0155]
(4)将和转子位置θ代入到三相电流计算公式中，得到三相电流瞬时值的给定
[0156]
(5)将三相电流瞬时值的给定分别与由电流传感器检测得到的三相电流的实际值i
a
、i
b
、i
c
作差，经过各自的电流调节器acar、acbr、accr，得到三相电压瞬时值的给定三相电压瞬时值的给定三相电压瞬时值的给定经过pwm生成单元输入到功率变换电路中，功率变换电路根据接收到的触发信号驱动永磁无刷电机运行。
[0157]
本实施例中涉及到的基波和五、七次谐波反电势常数k
e1
、k
e5
、k
e7
的获取方法如下：
[0158]
通过离线试验或有限元仿真，事先获得不同转子位置θ下的任一相反电势常数k
e
即任一相反电势e与转速n的比值，选取由0开始增加的一个电周期的反电势常数k
e
与转子位置θ数据，对其进行傅里叶分解得到基波和五、七次谐波反电势常数k
e1
、k
e5
、k
e7
。
[0159]
本实施例中涉及到的五、七次谐波反电势倍数k5、k7为带有符号的谐波反电势倍数，获取方法如下：
[0160]
根据永磁无刷电机绕组结构的特点，由前述经傅里叶分解得到的五、七次谐波反电势常数与基波反电势常数的初相位之差仅有0或π的可能，此时五次谐波反电势倍数k5的计算公式为
[0161][0162]
式中，当五次谐波反电势常数与基波反电势常数的初相位之差为0时s5代入1，初相位之差为π时s5代入0；
[0163]
同理，七次谐波反电势倍数k7的计算公式为
[0164][0165]
式中，当七次谐波反电势常数与基波反电势常数的初相位之差为0时s7代入1，初相位之差为π时s7代入0；
[0166]
最佳谐波电流倍数k
c
的获取方法如下：
[0167]
根据前述得到的由0开始增加的一个电周期的反电势常数k
e
与转子位置θ数据，获得三相反电势常数k
ea
、k
eb
、k
ec
与转子位置θ之间的关系数据库，把该数据代入下式
[0168]
k
eq
(θ)＝
‑
sinθ
·
k
ea
(θ)
‑
sin(θ
‑
2π/3)
·
k
eb
(θ)
‑
sin(θ+2π/3)
·
k
ec
(θ)
[0169]
计算得到q轴反电势常数k
eq
与转子位置θ之间的关系数据库，把该q轴反电势常数k
eq
代入下式
[0170][0171]
计算得到超出七次的非零序q轴谐波反电势常数k
ehq
与转子位置θ之间的关系数据库，对于所有转子位置θ，通过仿真确定使|k
ehq
(θ)+k
c
k
e1
cos(6θ)|达到最小值时的最佳谐波电流倍数k
c
。根据此处最佳谐波电流倍数k
c
的获取方法，可知最佳谐波电流倍数k
c
等于i
51
+i
71
与基波电流峰值i1之比。
[0172]
结合图4所示，本实施案例中以理想梯形波反电势为例，图4中(a)部分为a相反电势波形图，对其进行傅里叶分解，得到图4中(b)部分所示的经傅里叶分解后的反电势频谱图，经分解后的反电势基波峰值要远大于其他谐波峰值，且谐波的中主要含量为五次、七次分量。因此，由基波与五次、七次谐波分量相互作用产生的转矩脉动，要远大于其他分量间引起的转矩脉动。此外考虑到相电流基波分量的含量较高，其与超出七次的非零序谐波反电势作用产生转矩脉动亦不可忽略。因此本发明在着重考虑消除由基波分量与五、七次谐
波分量相互作用产生的六次转矩脉动的基础上，还考虑最大限度的抑制基波电流与超出七次的非零序谐波反电势作用产生的转矩脉动。
[0173]
以图4中(c)部分所示的同步旋转d
‑
q坐标系下反电势空间矢量图为例，转子位置处于零位置时，基波和五次、七次谐波的反电势空间矢量分别处于π/2、
‑
π/2、
‑
π/2；在该时刻以图5中(a)部分所示，基波和五次、七次谐波的电流空间矢量与基波反电势空间矢量同向，根据图示可以计算出基于谐波分解下的电磁转矩，其表达式如下：
[0174][0175]
式中，t为电磁转矩，e1、e5、e7分别为基波反电势与五次、七次谐波反电势峰值，i1、i5、i7分别为基波电流与五次、七次谐波电流峰值，e
h
为超出七次的非零序谐波反电势，ω为电机电角速度。
[0176]
由上式知，在非理想反电势情况下，电磁转矩除恒定分量外，还存在6次，12次以及e
h
引起的转矩脉动，考虑到12次转矩脉动以及e
h
与i5、i7相互作用产生转矩脉动对转矩平滑性的影响微乎其微，故可忽略；
[0177]
结合图6所示，图6(b)部分为本发明加入抑制策略且设定α＝0.5后输出的电磁转矩波形图，转矩输出平稳，转矩脉动得到了有效抑制；如图6(c)部分所示的当α设定为0.5时，转子永磁体产生的主磁通及其与定子电流产生的磁通进行合成后的气隙合成磁通轨迹，通过对比两者之间关系可以确定在设定α＝0.5时永磁电机无直轴电枢反应。
[0178]
因此，本发明针对非理想反电势作用下产生的转矩脉动进行抑制，有效的利用谐波分量来抑制永磁无刷电机运行过程中产生的谐波转矩脉动，此外能够在抑制转矩脉动的同时根据用户自行选择的性能选择参数α实现铜耗最小、无直轴电枢反应等高性能运行。
[0179]
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
[0180]
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。