CN107612436A - 基于永磁电机谐波反电势的转子位置计算方法 - Google Patents
基于永磁电机谐波反电势的转子位置计算方法 Download PDFInfo
- Publication number
- CN107612436A CN107612436A CN201710851734.4A CN201710851734A CN107612436A CN 107612436 A CN107612436 A CN 107612436A CN 201710851734 A CN201710851734 A CN 201710851734A CN 107612436 A CN107612436 A CN 107612436A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mover
- theta
- mtr
- 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
- 238000000205 computational method Methods 0.000 title claims abstract description 14
- 238000010276 construction Methods 0.000 claims abstract description 4
- 238000000034 method Methods 0.000 claims description 14
- 238000000605 extraction Methods 0.000 abstract description 5
- 230000000694 effects Effects 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000005389 magnetism Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Landscapes
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Control Of Motors That Do Not Use Commutators (AREA)
Abstract
本发明公开一种基于永磁电机谐波反电势的转子位置计算方法。首先对5、7次谐波反电势进行标幺化处理;然后,通过对标幺化谐波反电势进行和差化积处理,得到含有转子位置正余弦作为因数的变量;接着,利用积化和差、三倍角公式构造解调参考波,对提取得到的谐波反电势进行变换,得到幅值恒定的转子位置的正弦与余弦值,从而解调出转子位置的正余弦值;最后,利用反正切运算,通过正弦与余弦值得到转子位置信息,从而得到转子位置。本发明可以规避基于基波反电势计算谐波转子位置计算方法中对电机参数敏感性,可以得到鲁棒性更强的无传感器算法,便于及时获悉转子实时所处位置并对转子位置进行调整,从而保证永磁电机处于高效的工作状态。
Description
技术领域
本发明公开一种转子位置计算方法,特别涉及一种基于永磁电机谐波反电势的转子位置的计算方法,属于永磁电机无传感器算法领域。
背景技术
近20年,永磁电机的无传感器算法得到了广泛的应用。在很多应用场合下,无传感器控制策略是必不可少的。有的场合受限于空间或重量的限制,不能安装旋转位置传感器。而诸如风电等大型永磁电机应用的场合,更是无法在其轴上安装旋转位置编码器。不过,无传感器控制技术应用的最大优势在于降低成本。在一般永磁电机的应用场合,旋转位置传感器往往是设备成本中比重很大的一部分,因此,省去此类传感器无疑可以大大提高永磁电机驱动系统的经济性。
通过电机电压、电流及反电势等变量计算永磁电机转子位置的方法是永磁电机无位置传感器算法的基础。传统的无传感器方法均是基于基波反电势或基波磁链来计算转子位置角度。一般使用在固定坐标系下或者观测旋转坐标系下得到的转子磁链或者反电势的正交分量,通过反正切函数直接得到转子位置角度。但是由于这种方法基于基波模型,而电机为了正常工作必须有基波电流,所以不可避免地会受到电机定子参数,如定子电阻与电枢电感的影响。这就对转子位置提取的精度产生了负面影响。
发明内容
为了克服现有技术的不足,本发明提供了一种基于永磁电机谐波反电势的转子位置计算方法,可以不必通过基波反电势或基波磁链来获得转子位置信息。
一种基于永磁电机谐波反电势的转子位置计算方法,步骤如下:
(1)对谐波反电势进行标幺化处理;
(2)通过对标幺化谐波反电势进行和差化积得到含有转子位置正余弦作为因数的变量;
(3)构造解调参考波;
(4)从步骤(2)所得变量中,使用步骤(3)构造的解调参考波解调出幅值恒定的转子位置正余弦信息;
(5)利用转子位置正余弦信息得到转子位置。
所述步骤(1)包括如下步骤:
(1A)对如式(1)所示的α-β坐标系下的5次谐波反电势7次谐波反电势进行标幺化处理,得到如式(2)、式(3)的标幺制5、7次谐波反电势:
其中,分别表示在α方向、β方向的分量,分别表示在α方向、β方向的分量,分别表示5、7次谐波反电势标幺化后的结果,ωr表示转子角速度,θr表示转子所处位置,λ5、λ7分别表示五次谐波、七次谐波产生的磁链大小。
所述步骤(2)包括如下步骤:
(2A)使用和差化积公式,对5、7次谐波反电势,即式(2)、式(3)做如下运算:
所述步骤(3)包括如下步骤:
(3A)需构造更多与sin(6θr)、cos(6θr)相关的变量作为解调参考波,以解调出sin(6θr)与cos(6θr),对式(2)、式(3)使用积化和差公式,得:
(3B)联立步骤(3A)中得到的式(8)至式(11),得:
(3C)由三倍角公式:
将式(12)、式(13)分别代入式(14)、式(15)可得:
cos(6θr)=-3cos(2θr)+4cos3(2θr) (16)
sin(6θr)=3sin(2θr)-4sin3(2θr) (17)。
所述步骤(4)包括如下步骤:
(4A)由式(16)、式(17)以及式(4)—(7)解调出转子位置的正余弦值
所述步骤(5)包括如下步骤:
(5A)使用反正切函数通过转子位置正余弦值得到转子位置:
本发明的有益效果:
本发明实现了一种基于谐波反电势的转子位置计算方法,可以规避基于基波反电势计算谐波转子位置计算方法中对电机参数敏感性。而对于谐波反电势的提取,则可以采用将对应次谐波抑制至0的方式,从而提高反电势提取的鲁棒性。此方法为永磁电机无传感器算法提供了新的转子位置计算方式,在与对参数鲁棒的谐波反电势提取方法相结合后可以得到鲁棒性更强的无传感器算法。
附图说明
图1为本发明基于永磁电机谐波反电势的转子位置计算方法的一种流程图。
图2为5、7次谐波反电势标幺化后在α轴方向的分量sin(5θr)、sin(7θr)的曲线。
图3为通过提出的计算方法计算求得的转子所处位置θr随时间t的变化曲线与设定的参考转子位置随时间t的变化曲线。
具体实施方式
无传感器算法中的转子位置信息可以不通过基波反电势来提取,而通过谐波反电势来提取。由于谐波反电势对于电机出力的作用很小,所以我们发现若通过将谐波电流抑制为0的方法提取谐波反电势,则获得的转子信息可以是对参数不敏感的。
下面结合附图和实例对本发明作进一步的阐述。
本发明基于永磁电机谐波反电势的转子位置计算方法的一种流程图如图1所示。
首先,为得到幅值相等的5次谐波反电势与7次谐波反电势,以进行后续的变换,需对谐波反电势进行标幺化处理。如图2所示为5、7次谐波反电势经过标幺化处理后在α轴方向的分量sin(5θr)、sin(7θr)。
其中,分别表示在α方向、β方向的分量,分别表示在α方向、β方向的分量,分别表示5、7次谐波反电势标幺化后的结果,ωr表示转子角速度,θr表示转子所处位置,λ5、λ7分别表示五次谐波、七次谐波产生的磁链大小。
考虑到5θr=6θr-θr、7θr=6θr+θr,因此使用和差化积公式能够得到与目标θr有关的条件,故对5、7次谐波反电势做如下运算:
此时得到的有关基波量的条件由于系数部分为变量sin(6θr)、cos(6θr)而非常数,所以无法直接得到θr。因此必须得到6θr的正弦与余弦的信息,从而将sinθr、cosθr的系数变换为常数。为得到6θr的正弦与余弦的信息,首先考虑得到2θr:
经过简单的消元计算可得:
为了得到sin(6θr)、cos(6θr)以确定前文已得的有关基波量的条件的系数部分将sin(2θr)、cos(2θr)代入三倍角公式:
可得:
cos(6θr)=-3cos(2θr)+4cos3(2θr)
sin(6θr)=3sin(2θr)-4sin3(2θr)。
此时,sin(6θr)、cos(6θr)的信息已经被构造出来,故可以将sin(6θr)、cos(6θr)作为参考波与上文得到的条件相乘。构造如下:
在已知了sinθr与cosθr后,可方便地通过反正切变换得到转子位置θr的值:
为了能够更直观、快速地计算求得的转子所处位置θr,在线计算时应选择调用atan2函数而非普通的的反正切函数atan。这是因为atan2函数具有解出四象限角度的能力,也即其输出范围是-π至π之间,而非像atan一样在-π/2至π/2之间。
另外,由于此种方法涉及很多三角函数运算,因此在在线计算的时候需调用CPU的快速计算函数库,如TI公司DSP所使用的fastRTS_fbu32.h库,其三角函数与反三角函数的计算才用高精度的查表法,而非一般的math.h库在计算三角函数时采用的数值计算方式,因而可以极大地提高本算法的计算速度。
实测转子所处位置θr与参考转子位置如图3所示,由于实测转子所处位置θr处于-π至π之间,而参考转子位置处于0至2π之间,因此印证了通过所提方法计算求得的转子所处位置θr与参考转子位置是吻合的,并由此证明了所提方法的有效性。
Claims (6)
1.一种基于永磁电机谐波反电势的转子位置计算方法,其特征在于,步骤如下:
(1)对谐波反电势进行标幺化处理;
(2)通过对标幺化谐波反电势进行和差化积得到含有转子位置正余弦作为因数的变量;
(3)构造解调参考波;
(4)从步骤(2)所得变量中,使用步骤(3)构造的解调参考波解调出幅值恒定的转子位置正余弦信息;
(5)利用转子位置正余弦信息得到转子位置。
2.如权利要求1所述的方法,其特征在于:所述步骤(1)包括如下步骤:
(1A)对如式(1)所示的α-β坐标系下的5次谐波反电势7次谐波反电势进行标幺化处理,得到如式(2)、式(3)的标幺制5、7次谐波反电势:
<mrow>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mn>5</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&omega;</mi>
<mi>r</mi>
</msub>
<msub>
<mi>&lambda;</mi>
<mn>5</mn>
</msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mo>(</mo>
<mn>5</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<mn>5</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mn>7</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&omega;</mi>
<mi>r</mi>
</msub>
<msub>
<mi>&lambda;</mi>
<mn>7</mn>
</msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mo>-</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mo>(</mo>
<mn>7</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>cos</mi>
<mo>(</mo>
<mn>7</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mn>5</mn>
</msub>
<msqrt>
<mrow>
<msubsup>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&alpha;</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&beta;</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mfrac>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&alpha;</mi>
<mi>n</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&beta;</mi>
<mi>n</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mo>(</mo>
<mn>5</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<mn>5</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mn>7</mn>
</msub>
<msqrt>
<mrow>
<msubsup>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&alpha;</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&beta;</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mfrac>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&alpha;</mi>
<mi>n</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&beta;</mi>
<mi>n</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mo>-</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mo>(</mo>
<mn>7</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>cos</mi>
<mo>(</mo>
<mn>7</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,分别表示在α方向、β方向的分量,分别表示在α方向、β方向的分量,分别表示5、7次谐波反电势标幺化后的结果,ωr表示转子角速度,θr表示转子所处位置,λ5、λ7分别表示五次谐波、七次谐波产生的磁链大小。
3.如权利要求1所述的方法,其特征在于:所述步骤(2)包括如下步骤:
(2A)使用和差化积公式,对5、7次谐波反电势,即式(2)、式(3)做如下运算:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&alpha;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&alpha;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>5</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>7</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mo>-</mo>
<mn>2</mn>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mn>6</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>sin&theta;</mi>
<mi>r</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&beta;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&beta;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>5</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>7</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>=</mo>
<mn>2</mn>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>6</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>cos&theta;</mi>
<mi>r</mi>
</msub>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&alpha;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&alpha;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>5</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>7</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mn>2</mn>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>6</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>cos&theta;</mi>
<mi>r</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&beta;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&beta;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>5</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mn>7</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mn>2</mn>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>6</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>sin&theta;</mi>
<mi>r</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
4.如权利要求1所述的方法,其特征在于:所述步骤(3)包括如下步骤:
(3A)需构造更多与sin(6θr)、cos(6θr)相关的变量作为解调参考波,以解调出sin(6θr)与cos(6θr),对式(2)、式(3)使用积化和差公式,得:
<mrow>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&alpha;</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&alpha;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>&lsqb;</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>12</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&beta;</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&beta;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>&lsqb;</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>12</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&alpha;</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&beta;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>&lsqb;</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>12</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&beta;</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&alpha;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>&lsqb;</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>12</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
(3B)联立步骤(3A)中得到的式(8)至式(11),得:
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&beta;</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&beta;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&alpha;</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&alpha;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&alpha;</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&beta;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&beta;</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&alpha;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
(3C)由三倍角公式:
将式(12)、式(13)分别代入式(14)、式(15)可得:
cos(6θr)=-3cos(2θr)+4cos3(2θr) (16)
sin(6θr)=3sin(2θr)-4sin3(2θr) (17)。
5.如权利要求1所述的方法,其特征在于:所述步骤(4)包括如下步骤:
(4A)由式(16)、式(17)以及式(4)—(7)解调出转子位置的正余弦值
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mi>&alpha;</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&alpha;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&alpha;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mn>6</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&beta;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&beta;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mn>6</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mn>2</mn>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mn>6</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>sin&theta;</mi>
<mi>r</mi>
</msub>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mn>6</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mn>6</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>sin&theta;</mi>
<mi>r</mi>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mn>6</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mn>2</mn>
<msub>
<mi>sin&theta;</mi>
<mi>r</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mi>&beta;</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&beta;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&beta;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mn>6</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>5</mn>
<mi>&alpha;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>7</mn>
<mi>&alpha;</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mn>6</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mn>2</mn>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mn>6</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>cos&theta;</mi>
<mi>r</mi>
</msub>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mn>6</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mn>6</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>cos&theta;</mi>
<mi>r</mi>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mn>6</mn>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mn>2</mn>
<msub>
<mi>cos&theta;</mi>
<mi>r</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>19</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
6.如权利要求1所述的方法,其特征在于:所述步骤(5)包括如下步骤:
(5A)使用反正切函数通过转子位置正余弦值得到转子位置:
<mrow>
<msub>
<mi>&theta;</mi>
<mi>r</mi>
</msub>
<mo>=</mo>
<mi>arctan</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mi>&alpha;</mi>
</msub>
<msub>
<mover>
<mi>E</mi>
<mo>^</mo>
</mover>
<mi>&beta;</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>arctan</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>sin&theta;</mi>
<mi>r</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>cos&theta;</mi>
<mi>r</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>20</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710851734.4A CN107612436B (zh) | 2017-09-19 | 2017-09-19 | 基于永磁电机谐波反电势的转子位置计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710851734.4A CN107612436B (zh) | 2017-09-19 | 2017-09-19 | 基于永磁电机谐波反电势的转子位置计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107612436A true CN107612436A (zh) | 2018-01-19 |
CN107612436B CN107612436B (zh) | 2019-11-22 |
Family
ID=61061201
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710851734.4A Expired - Fee Related CN107612436B (zh) | 2017-09-19 | 2017-09-19 | 基于永磁电机谐波反电势的转子位置计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107612436B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110380733A (zh) * | 2019-05-21 | 2019-10-25 | 日本电产凯宇汽车电器(江苏)有限公司 | 一种汽车刹车系统控制器用旋转变压器解码方法 |
CN112095218A (zh) * | 2019-06-18 | 2020-12-18 | 美名格-艾罗有限公司 | 纱线输送装置和用于将纱线输送到纺织机的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103199779A (zh) * | 2013-04-22 | 2013-07-10 | 哈尔滨工业大学 | 基于自适应滤波的内置式永磁同步电机转子位置观测装置及观测方法 |
US20140159632A1 (en) * | 2012-12-10 | 2014-06-12 | Electronics And Telecommunications Research Institute | Motor driving module and brushless dc motor system |
CN106841901A (zh) * | 2017-03-09 | 2017-06-13 | 浙江大学 | 一种变频器驱动内嵌式永磁同步电机定子匝间短路故障诊断方法 |
CN106877768A (zh) * | 2017-02-20 | 2017-06-20 | 杭州富生电器有限公司 | 多相永磁电机转子位置辨识方法 |
CN107070318A (zh) * | 2017-05-09 | 2017-08-18 | 天津大学 | 一种永磁无刷直流电机换相转矩脉动的谐波分析法 |
US9739815B2 (en) * | 2012-07-20 | 2017-08-22 | Brose Fahrzeugteile Gmbh & Co. Kommanditsellschaft, Wurzburg | Method for determining the rotor position of an electrically-commuted multi-phase direct current motor |
-
2017
- 2017-09-19 CN CN201710851734.4A patent/CN107612436B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9739815B2 (en) * | 2012-07-20 | 2017-08-22 | Brose Fahrzeugteile Gmbh & Co. Kommanditsellschaft, Wurzburg | Method for determining the rotor position of an electrically-commuted multi-phase direct current motor |
US20140159632A1 (en) * | 2012-12-10 | 2014-06-12 | Electronics And Telecommunications Research Institute | Motor driving module and brushless dc motor system |
CN103199779A (zh) * | 2013-04-22 | 2013-07-10 | 哈尔滨工业大学 | 基于自适应滤波的内置式永磁同步电机转子位置观测装置及观测方法 |
CN106877768A (zh) * | 2017-02-20 | 2017-06-20 | 杭州富生电器有限公司 | 多相永磁电机转子位置辨识方法 |
CN106841901A (zh) * | 2017-03-09 | 2017-06-13 | 浙江大学 | 一种变频器驱动内嵌式永磁同步电机定子匝间短路故障诊断方法 |
CN107070318A (zh) * | 2017-05-09 | 2017-08-18 | 天津大学 | 一种永磁无刷直流电机换相转矩脉动的谐波分析法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110380733A (zh) * | 2019-05-21 | 2019-10-25 | 日本电产凯宇汽车电器(江苏)有限公司 | 一种汽车刹车系统控制器用旋转变压器解码方法 |
CN112095218A (zh) * | 2019-06-18 | 2020-12-18 | 美名格-艾罗有限公司 | 纱线输送装置和用于将纱线输送到纺织机的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107612436B (zh) | 2019-11-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101355337B (zh) | 基于磁场正交控制的永磁同步电动机的驱动控制方法 | |
CN103199779B (zh) | 基于自适应滤波的内置式永磁同步电机转子位置观测装置及观测方法 | |
CN103199788B (zh) | 永磁同步电机转子位置观测器 | |
CN106788127B (zh) | 基于二维查表与插值法的逆变器非线性谐波补偿方法 | |
CN105743414A (zh) | 电力转换装置、控制装置及载波频率的改变方法 | |
CN107612436A (zh) | 基于永磁电机谐波反电势的转子位置计算方法 | |
CN109039204A (zh) | 车用永磁同步电机无位置传感器模型预测控制系统及方法 | |
CN108574444A (zh) | 一种用于永磁同步电机转子初始位置检测方法 | |
CN108900127A (zh) | 考虑交叉耦合效应的ipmsm低速段无位置传感器控制方法 | |
CN107769636A (zh) | 一种无位置传感器永磁同步电机转子位置检测方法 | |
CN103078572A (zh) | 一种永磁同步电机高精度转子位置估计方法 | |
CN110401393A (zh) | 宽容分层序列模型预测控制方法、设备和存储介质 | |
CN103414422A (zh) | 永磁同步发电机的直接交轴磁链控制装置及其方法 | |
CN106100486A (zh) | 基于双定向电压矢量的永磁同步电机初始位置确定方法 | |
CN103997262B (zh) | 基于无传感器轮毂式电机的电动自行车正弦波控制方法 | |
CN107834914A (zh) | 基于霍尔的无刷电机转子角度计算方法 | |
WO2022120772A1 (zh) | 永磁同步电机的磁场定向校正方法、装置、设备及介质 | |
CN206077268U (zh) | 用于速度和位置反馈的正余弦信号处理电路 | |
CN103346723B (zh) | 一种无位置传感器控制装置及位置检测方法 | |
Wang et al. | Position estimation of outer rotor PMSM using linear Hall effect sensors and neural networks | |
CN106877768A (zh) | 多相永磁电机转子位置辨识方法 | |
CN110086399B (zh) | 一种永磁同步电机转子位置复合检测及起动运行方法 | |
CN112104290B (zh) | 电机转子磁极初始位置辨识方法及装置 | |
Zhao et al. | Research on the decoding method of resolver | |
CN112564572A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191122 |