CN108206659B - 基于旋转高频注入算法的永磁同步电机转子位置估计方法 - Google Patents

基于旋转高频注入算法的永磁同步电机转子位置估计方法 Download PDF

Info

Publication number
CN108206659B
CN108206659B CN201810120333.6A CN201810120333A CN108206659B CN 108206659 B CN108206659 B CN 108206659B CN 201810120333 A CN201810120333 A CN 201810120333A CN 108206659 B CN108206659 B CN 108206659B
Authority
CN
China
Prior art keywords
frequency
phase
current
coordinate system
axis
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.)
Active
Application number
CN201810120333.6A
Other languages
English (en)
Other versions
CN108206659A (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.)
Lishui Boyuan Technology Co ltd
Original Assignee
Lishui Boyuan Technology 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 Lishui Boyuan Technology Co ltd filed Critical Lishui Boyuan Technology Co ltd
Priority to CN201810120333.6A priority Critical patent/CN108206659B/zh
Publication of CN108206659A publication Critical patent/CN108206659A/zh
Application granted granted Critical
Publication of CN108206659B publication Critical patent/CN108206659B/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/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • H02P21/18—Estimation of position or speed
    • 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
    • H02P2203/00—Indexing scheme relating to controlling arrangements characterised by the means for detecting the position of the rotor
    • H02P2203/03—Determination of the rotor position, e.g. initial rotor position, during standstill or low speed operation
    • 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
    • H02P2203/00—Indexing scheme relating to controlling arrangements characterised by the means for detecting the position of the rotor
    • H02P2203/11—Determination or estimation of the rotor position or other motor parameters based on the analysis of high frequency signals

Abstract

本发明提出一种基于旋转高频电压注入算法的永磁同步电机转子位置估计方法，主要解决现有技术估计转子位置信号信噪比低，导致估算转子位置精度不高的问题。其实现方案为：首先在静止坐标系注入两个不同的旋转高频电压信号，获取高频电流响应，分别用纯延时滤波获取角频率为ωh高频电流信号响应；用高通滤波器获取角频率为ωh1的高频电流信号响应；根据得到两个的高频信号响应，获得两个估计值，求其平均，得到最终转子位置估计值
Figure DDA0001571856880000011
本发明由于在静止坐标系注入两个不同的旋转高频电压信号，使得估计值更可靠，同时由于利用纯延时滤波器，提高了信噪比，减小了相位滞后的现象，最终提高了估算转子的位置精度。可用于电机控制系统。

Description

基于旋转高频注入算法的永磁同步电机转子位置估计方法
技术领域
本发明属于无位置传感器控制技术领域，特别涉及一种永磁同步电机转子位置的估计方法，可用于电机控制系统。
背景技术
为了实现最佳电压矢量控制或最大的转矩控制，需要精确地获得永磁同步电机转子的位置信息。目前获取永磁同步电机转子位置的方法，主要有利用位置传感器直接检测的方法和无位置传感器估计检测的方法。其中：
位置传感器直接检测的方法，主要是在永磁同步电机中放置机械式传感器，通过传感器直接检测转子位置。但是这种方法的可靠性和控制精确较低，且成本较高。
无位置传感器估计检测的方法，根据不同的原理对永磁同步电机转子位置估算方法可以分为以下几种：
1)利用定子端电压和电流直接计算转子速度和位置。
这种方法时通过检测电动机的端电压、电流、利用磁链计算出转子的位置角和速度。磁链是通过反电动势积分求得，由于积分器有零漂问题，得到的磁链值会有积分误差。当转速较低，无法建立起足够大的反电动势，无法精确计算出转子位置角和速度。
2)模型参考自适应估计法。
这种方法主要是将含有待估计参数的方程作为可调模型，将不含未知参数的方程作为参考模型，两个模型具有相同物理意义的输出量。两个模型同时工作，并利用其输出量的差值，根据自适应算法来调节可调模型的参数，以达到控制对象的输出跟踪参考模型的目的。当差值等于零，表示估计的转子的位置就是转子真实位置，但是电机的参数影响着这种方法的估算精度。
3)观测器基础上的估算方法。
观测器的实质是进行状态重构，即重新构造一个系统，利用原系统中可以测量的变量作为构造系统的输入信号，重构系统的输出信号在一定条件下与原系统的状态等价。将这个用以实现状态重构的系统为观测器。
所述的这种观测器目前主要包括：全阶状态观测器、扩展卡尔曼滤波器、滑模观测器等。其中：
全阶状态观测器，其需要在电动机高速时和低速时采用不同的增益矩阵，而且由于状态观测器受电动机参数变化影响很大，还需要另外一个状态观测器来估计电动机的参数，这样使无传感器传动系统的估计算法更加复杂。同时系统还存在对负载变化比较敏感的不足。
扩展卡尔曼滤波器，其是从随机噪声信号中得到转速的最优观测，不需要电机精确的数学模型，克服了电动机参数的时变对系统的不良影响，可精确估计出转子速度和位置，但是这种方法的计算很大。
滑膜观测器，是利用滑膜变结构控制系统对参数扰动鲁棒性强的特点，把一般状态观测器中的控制回路改成滑膜变结构的形式。滑膜变结构控制的本质是滑膜运动，通过结构变化开关，以很高的频率来回切换，使状态的运动点以很小的幅度在相平面上运动，最终运动到稳定点。由于滑膜变结构控制在本质上是不连续的开关控制，会引起比较大的转矩脉动。
4)人工智能理论基础上的估算方法。
该方法主要是将专家系统、模糊控制、人工神经元网络算法应用到电机调速控制系统。人工神经元网络是对神经元网络进行辨识，先规定网络结构，通过学习系统的输入和输出，使性能指标满足要求，从而归纳出隐含在系统输入输出之间的关系。由于其理论方面还不够成熟，硬件方面也有很大的难度，不易走向实用化。
5)高频注入法。
该方法采用在电动机出线端注入高频电压或者电流信号，根据电机内部固有凸性，检测高频电流或电压响应来获取转子位置和速度信息。高频注入法从注入信号的性质来分，又可分为旋转高频信号和脉振高频信号注入法；从信号前馈的位置来分，又可分为高频电压信号注入法和高频电流信号注入法。
高频脉振注入法，是将高频脉动信号注入到同步坐标系的直轴，由此产生的内置式永磁同步电机的高频响应电流，通过分析内置式永磁同步电机的高频阻抗响应，提取凸极位置信息，最终得到转子的位置。当用于表面式永磁同步电机需要注入的信号需要有比较高的频率，且幅度不能过大。高频脉振注入法可以运行在低速和零速阶段，且同时适用于内嵌式和表贴式的永磁同步电机。但是需要严格保证由水平轴注入，实际中往往存在一定的误差，因此其准确性较差，并且会存在转矩脉动的问题。
高频旋转注入法，是通过将电流信号或者电压信号注入到永磁同步电机的定子中，通过对电机凸极的调制可以获得高频响应电流，经过对该响应电流进行滤波与解调，从而将其中所包含的电机转子位置信息进行提取。该方法通过对其注入的信号不同分为高频旋转电流注入与高频旋转电压注入两种方法。高频旋转电流注入法是将高频旋转电流注入到电机的静止两相旋转坐标系中。该方法的优点是由于高频阻抗，对信号提取方便。所以为了减少对电机造成的损耗需要注入的高频电流的幅值进行限制。缺点是电流调节器需要很大的带宽，需要极高的性能来对其输入的电流信号的高频率和低幅值进行保证，如果电流频率过低会与电机中的基波信号产生干扰。所以对电流调节器的设计难度很大。而高频旋转电压注入法是将高频旋转电压直接从逆变器注入，无需大量运算，与高频电流注入法相比更为简单实用，且高频旋转电压注入法的响应为定子电流，提取方便。但是该高频旋转注入法在估计转子位置信号处理过程中，信噪比较低，并且存在相位滞后现象，导致估算转子位置精度不高。
发明内容
本发明的目的在于针对上述高频旋转电压注入法的不足，提出一种基于旋转高频注入算法的永磁同步电机转子位置估计方法，以减小相位滞后，提高估算转子位置精度。
本发明的技术方案是：首先在静止坐标系注入两个不同的旋转高频电压信号，根据内置式永磁同步电机的不对称性，然后获取高频电流响应，分别用纯延时滤波获取角频率为ωh高频电流信号响应，利用高通滤波器角频率为ωh1的高频电流信号响应，根据两个高频信号的响应得到两个转子位置估计值
Figure BDA0001571856860000031
和
Figure BDA0001571856860000032
其实现步骤包括如下：
(1)定义的坐标系：包括两相静止坐标系、估计两相同步旋转坐标系、和真实两相同步旋转坐标系这三个坐标系，其中：α-β为两相静止坐标系，d-q为实际转子同步旋转坐标系，
Figure BDA0001571856860000041
为估计转子同步旋转坐标系；
(2)在两相静止的坐标系中注入不同信号uαhi和uβhi，即在α轴注入两个高频余弦信号uαhi，在β轴注入两个高频正弦信号uβhi
uαhi＝umhcos(ωht)+umh1cos(ωh1t)，
uβhi＝umhsin(ωht)+umh1sin(ωh1t)，
其中，ωh是第一个高频信号的角频率，ωh1为第二个高频信号角频率，umh是信号角频率为ωh的幅值，umh1是信号角频率为ωh1的幅值，高频电压的角频率远远大于基频信号的角频率，且满足ωh1＝2ωh；
(3)检测出永磁同步电机输出的三相电流ia ib和ic；
(4)将三相电流进行Clarke变换，得到两相静止坐标系α-β下的两相定子电流iα和iβ，该两相电流不仅含有基频电流信号，且含有两个注入高频电压产生的高频电流信号；
(5)对两相静止坐标系α-β下的两相定子电流iα和iβ，进行延时为
Figure BDA0001571856860000042
的纯滤波，获得角频率为ωh的两个高频电流信号iαh和iβh，其中iαh表示为α相高频电流iα经滤波之后的输出信号，iβh表示为β相高频电流iα经滤波之后的输出信号；
(6)将得到的角频率为ωh的高频电流信号iαh和iβh通过同步轴高通滤波器进行滤波，获得角频率为ωh的负序高频电流分量iαhn和iβhn，其中iαhn表示为iαh经过同步轴高通滤波之后的输出信号，iβhn表示为iβh经过同步轴高通滤波之后的输出信号；
(7)定义转子估计位置误差为
Figure BDA0001571856860000043
其中
Figure BDA0001571856860000044
为转子位置估计值，θ为实际转子位置值，对负相序高频电流分量iαhn和iβhn进行锁相环处理，得到估计位置偏差信号εh(Δθ)，其中εh(Δθ)是注入角频率为ωh得到的位置偏差信号；
(8)将估计得到位置偏差信号εh(Δθ)作为PI调节器的输入，PI调节器的输出为估计转子角速度
Figure BDA0001571856860000051
对估计转子角速度
Figure BDA0001571856860000052
进行积分调节得到转子位置估计值
Figure BDA0001571856860000053
其中
Figure BDA0001571856860000054
是注入角频率为ωh得到的估计转子角速度
Figure BDA0001571856860000055
Figure BDA0001571856860000056
是注入角频率为ωh得到的转子位置估计值；
(9)对两相静止坐标系α-β下的两相定子电流iα和iβ,，进行高通滤波，获得角频率为ωh1的高频电流信号iαh1和iβh1，其中iαh1表示为α相高频电流iα经高通滤波之后的输出信号，iβh1表示为β相高频电流iβ经高通滤波之后的输出信号；
(10)对角频率为ωh1的高频电流信号iαh1和iβh1进行同步轴高通滤波器，获得角频率为ωh1的负序高频电流分量iαhn1和iβhn1；
(11)将负相序高频电流分量iαhn1和iβhn1进行锁相环处理，估计位置偏差信号εh1(Δθ)：
(12)将估计得到位置偏差信号εh1(Δθ)作为PI调节器的输入，PI调节器的输出为估计转子角速度
Figure BDA0001571856860000057
对估计转子角速度
Figure BDA0001571856860000058
进行积分调节，得到转子位置估计值
Figure BDA0001571856860000059
(13)将注入角频率为ωh的高频信号得到的转子角速度估计值
Figure BDA00015718568600000520
和注入角频率为ωh1高频信号得到的转子角速度估计值
Figure BDA00015718568600000510
进行相加求平均值，得到最终转子角速度估计值
Figure BDA00015718568600000511
(14)将注入角频率为ωh高频信号得到的转子位置估计值
Figure BDA00015718568600000512
和注入角频率为ωh1的高频信号得到的转子位置估计值
Figure BDA00015718568600000513
进行相加求平均值得到最终转子位置估计值
Figure BDA00015718568600000514
(15)对两相定子电流iα和iβ进行Park变换，通过低通滤波器，输出为两相估计同步旋转坐标系
Figure BDA00015718568600000515
下的两相电流
Figure BDA00015718568600000516
和
Figure BDA00015718568600000517
(16)将最终转子角速度估计值
Figure BDA00015718568600000518
和最终转子位置估计值
Figure BDA00015718568600000519
作为反馈信号，控制电机实际转子位置值与转子位置估计值相差收敛到零，使电机稳定运行。
本发明具有如下优点：
1.本发明通过在两相静止坐标系α-β同时注入角频率为ωh和ωh1的高频信号，且满足ωh1＝2ωh，并通过注入不同的角频率得到不同转子位置估计值，再将两者相加求平均得到最终转子位置估计值
Figure BDA0001571856860000061
提高了转子的估计精度和有效性；
2.本发明通过纯延时滤波器，在提取角频率为ωh的高频电流响应时，对角频率为ωh1的高频电流响应进行全部抑制，使得高频电流响应幅值提高了16倍，且信噪比也得到了明显提高，从而提高转子位置估计值的精度。
附图说明
图1为本发明的原理框图；
图2为本发明中两相静止坐标系、估计两相同步旋转坐标系、真实两相同步旋转坐标系三者之间的关系图；
图3为本发明中纯延时滤波电路原理图。
具体实施方式
下面结合附图对本发明作进一步的描述。
参照附图1，本发明的实现如下：
步骤1，定义坐标系关系图。
参照图2，本步骤定义的坐标系包括两相静止坐标系、估计两相同步旋转坐标系、和真实两相同步旋转坐标系这三个坐标系，其中：α-β为两相静止坐标系，d-q为实际转子同步旋转坐标系，
Figure BDA0001571856860000062
为估计转子同步旋转坐标系；
将转子估计位置误差定义为
Figure BDA0001571856860000063
其中
Figure BDA0001571856860000064
为转子位置估计值，θ为实际转子位置值。
步骤2，在两相静止的坐标系中注入不同信号uαhi和uβhi。
在两相静止的坐标系的α轴注入两个高频余弦信号，即uαhi：
uαhi＝umhcos(ωht)+umh1cos(ωh1t)；
在两相静止的坐标系的β轴注入两个高频正弦信号，即uβhi：
uβhi＝umhsin(ωht)+umh1sin(ωh1t)，
其中，ωh是第一个高频信号的角频率，ωh1为第二个高频信号角频率，umh是信号角频率为ωh的幅值，umh1是信号角频率为ωh1的幅值，高频电压的角频率远远大于基频信号的角频率，且满足ωh1＝2ωh；
步骤3，检测出永磁同步电机输出的三相电流参数。
对永磁同步电机三相进行电流采样，获得永磁同步电机的三相电流值，即A相中的电流为ix，B相中的电流为iy，C相中的电流为iz。
步骤4，对三相电流进行Clarke变换，得到两相静止坐标系α-β下的两相定子电流iα和iβ，该两相电流不仅含有基频电流信号，且含有两个注入高频电压产生的高频电流信号。
4a)通过如下公式对永磁同步电机三相电流ix，iy和iz进行Clarke变换,得到两相静止坐标系α-β下的两相定子电流iα和iβ：
Figure BDA0001571856860000071
其中，α相定子电流为iα，β相定子电流为iβ，分别表示为：
Figure BDA0001571856860000072
Figure BDA0001571856860000073
其中t是时间变量，iαb为基频余弦电流响应，iαi是角频率为ωh余弦电流响应，iαi1是角频率为ωh1余弦电流响应，iβb为基频正弦电流响应，iβi是角频率为ωh正弦电流响应，iβi1是角频率为ωh1正弦电流响应，IA表示为基频电流幅值，ωr为转子的速度，
Figure BDA0001571856860000081
为基频电流初始相位，I11表示为角频率为ωh正相幅值，I12表示为角频率为ωh负相幅值，I21表示为角频率为ωh1正相幅值，I22表示为角频率为ωh1负相幅值，
Figure BDA0001571856860000082
Figure BDA0001571856860000083
Ld为两相旋转坐标系d轴直轴电感，Lq为两相旋转坐标系q轴的交轴电感。
步骤5，对两相静止坐标系α-β下的两相定子电流iα和iβ，进行延时纯滤波，获得角频率为ωh的两个高频电流信号iαh和iβh。
参照附图3，本步骤具体实施，参考如下：
5a)对两相静止坐标系α-β下的两相定子电流iα和iβ进行一次纯延时，α轴电流iα延时
Figure BDA0001571856860000084
后得到的信号为iαu，β轴电流iβ延时
Figure BDA0001571856860000085
后得到的信号为iβu,其可表示如下：
Figure BDA0001571856860000086
Figure BDA0001571856860000087
5b)对两相静止坐标系α-β下的两相定子电流iα和iβ进行滤波，得到α轴电流iα经过一次纯延时滤波后的电流iαoi，β轴电流iβ经过一次纯延时滤波的电流iβoi，其表示如下：
Figure BDA0001571856860000091
Figure BDA0001571856860000092
由于高频电压的角频率ωi远远大于基频信号的角频率ωr，即得到：
Figure BDA0001571856860000093
故可将上式简化为：
Figure BDA0001571856860000094
Figure BDA0001571856860000095
5c)对5b)的结果iαoi,iβoi再进行延时为
Figure BDA0001571856860000096
的二次纯延时，得到α轴一次延迟电流iαoi的二次延迟信号iαui，β轴一次延迟电流iβoi的二次延迟信号iβui，其分别表示如下：
Figure BDA0001571856860000097
Figure BDA0001571856860000098
5d)对5b)的结果iαoi和iβoi进行纯延时滤波，得到α轴的二次纯延时滤波后的电流iαh，β轴的二次纯延时滤波后的电流iβh，其分别表示如下：
Figure BDA0001571856860000099
Figure BDA00015718568600000910
步骤6，对5d)的结果进行同步轴高通滤波器滤波。
6a)对角频率为ωh的α轴的二次纯延时滤波后的电流iαh进行同步轴高通滤波，得到角频率为ωh的α相高频电流负序分量iαhn：
Figure BDA0001571856860000101
6b)对角频率为ωh的β轴的二次纯延时滤波后的电流iβh进行同步轴高通滤波，得到角频率为ωh的β相高频电流负序分量iβhn：
Figure BDA0001571856860000102
步骤7，根据6a)和6b)的结果，得到估计位置偏差信号εh(Δθ)。
7a)设永磁同步电机转子位置估计值
Figure BDA0001571856860000103
的初始值设为0，即
Figure BDA0001571856860000104
7b)设角频率为ωh的余弦信号为io，设角频率为ωh的正弦信号为iq，并将转子位置估计值
Figure BDA0001571856860000105
作为这两个信号的相位角，分别表示为：
Figure BDA0001571856860000106
Figure BDA0001571856860000107
7c)将6a)中得到的结果iαhn与io相乘，得到α轴解调后的信号ix：
Figure BDA0001571856860000108
7d)将6b)中得到结果iβh与iq相乘，得到β轴解调后的信号iy：
Figure BDA0001571856860000109
7e)将7c)和7d)得到的结果ix，iy两者相加，得到估计位置偏差信号εh(Δθ)。
Figure BDA00015718568600001010
步骤8，根据估计位置偏差信号εh(Δθ)，得到转子角速度估计值
Figure BDA00015718568600001011
根据估计位置偏差信号εh(Δθ)计算转子角速度估计值
Figure BDA0001571856860000111
时，需要对估计位置偏差信号εh(Δθ)进行线性调节，本实例通过引入PI调节器进行调节，由于PI调节器中含有比例放大和积分两个部分，故要分别通过这两部分进行调节，其实现如下：
8a)将位置偏差信号εh(Δθ)作为PI调节器中的放大部分的输入，其输出为εh(Δθ)的放大值
Figure BDA0001571856860000112
8b)将位置偏差信号εh(Δθ)作为PI调节器中的积分部分的输入，其输出为εh(Δθ)的积分值
Figure BDA0001571856860000113
8c)将εh(Δθ)的放大值
Figure BDA0001571856860000114
和εh(Δθ)的积分值
Figure BDA0001571856860000115
相加，得到转子角速度估计值
Figure BDA0001571856860000116
即：
Figure BDA0001571856860000117
步骤9，根据转子角速度估计值
Figure BDA0001571856860000118
得到转子位置估计值
Figure BDA0001571856860000119
通过对转子角速度估计值
Figure BDA00015718568600001110
进行积分，获得转子位置估计值
Figure BDA00015718568600001111
步骤10，对步骤4得到的电流iα和iβ，进行高通滤波。
10a)对α轴电流iα进行高通滤波，得到角频率为ωh1的高频电流信号iαh1：
Figure BDA00015718568600001112
10b)对β轴电流iβ进行高通滤波，得到角频率为ωh1的高频电流信号iβh1：
Figure BDA00015718568600001113
步骤11，对10a)和10b)得到的结果iαh1和iβh1进行同步轴高通滤波。
11a)对α轴角频率为ωh1的高频电流信号iαh1进行同步轴高通滤波，得到角频率为ωh1的负序高频电流分量iαhn1：
Figure BDA00015718568600001114
11b)对β轴角频率为ωh1的高频电流信号iβh1进行同步轴高通滤波，得到角频率为ωh1的负序高频电流分量iβhn1：
Figure BDA0001571856860000121
步骤12，根据11a)和11b)的结果，得到估计位置偏差信号εh1(Δθ)。
12a)设永磁同步电机转子位置估计值
Figure BDA0001571856860000122
的初始值设为0，即
Figure BDA0001571856860000123
12b)设角频率为ωh1的余弦信号为im，设角频率为ωh1的正弦信号为ig，并将转子位置估计值
Figure BDA0001571856860000124
作为这两个信号的相位角，分别表示为：
Figure BDA0001571856860000125
Figure BDA0001571856860000126
12c)将11a)中得到的结果iαhn1与im相乘，得到α轴解调后的信号ir：
Figure BDA0001571856860000127
12d)将11b)中得到结果iβh1与ig相乘，得到β轴解调后的信号is：
Figure BDA0001571856860000128
12e)将12c)和12d)得到的结果ir，is两者相加，得到估计位置偏差信号εh1(Δθ)。
Figure BDA0001571856860000129
步骤13，根据估计位置偏差信号εh1(Δθ)，得到转子角速度估计值
Figure BDA00015718568600001210
根据估计位置偏差信号εh1(Δθ)计算转子角速度估计值
Figure BDA00015718568600001211
时，需要对估计位置偏差信号εh1(Δθ)进行线性调节，本实例通过引入PI调节器进行调节，由于PI调节器中含有比例放大和积分两个部分，故要分别通过这两部分进行调节，其实现如下：
13a)将位置偏差信号εh1(Δθ)作为PI调节器中的放大部分的输入，其输出为εh1(Δθ)的放大值
Figure BDA00015718568600001212
13b)将位置偏差信号εh1(Δθ)作为PI调节器中的积分部分的输入，其输出为εh1(Δθ)的积分值
Figure BDA0001571856860000131
13c)将εh1(Δθ)的放大值
Figure BDA0001571856860000132
和εh1(Δθ)的积分值
Figure BDA0001571856860000133
相加，得到转子角速度估计值
Figure BDA0001571856860000134
即：
Figure BDA0001571856860000135
步骤14，根据转子角速度估计值
Figure BDA0001571856860000136
得到转子位置估计值
Figure BDA0001571856860000137
通过对转子角速度估计值
Figure BDA0001571856860000138
进行积分，获得转子位置估计值
Figure BDA0001571856860000139
步骤15，根据步骤8和步骤13的结果，求得最终角速度估计值
Figure BDA00015718568600001310
将步骤8因注入角频率为ωh信号而得到的角速度估计值
Figure BDA00015718568600001311
与步骤13因注入角频率为ωh1信号得到的角速度估计值
Figure BDA00015718568600001312
求平均，得到最终角速度估计值
Figure BDA00015718568600001313
其可表示为：
Figure BDA00015718568600001314
步骤16，将步骤9和步骤14的结果，求得最终转子位置估计值
Figure BDA00015718568600001315
将步骤9因注入角频率为ωh信号得到的角速度估计值
Figure BDA00015718568600001316
与步骤14因注入角频率为ωh1信号得到的角速度估计值
Figure BDA00015718568600001317
求平均，得到最终转子位置估计值
Figure BDA00015718568600001318
其可表示为：
Figure BDA00015718568600001319
步骤17，根据步骤16的结果，对两相定子电流iα和iβ进行Park变换，获得两相估计同步旋转坐标系
Figure BDA00015718568600001320
下的电流
Figure BDA00015718568600001321
和
Figure BDA00015718568600001322
17a)对两相静止坐标系α-β下的两相定子电流iα和iβ进行Park变换，获得两相估计同步旋转坐标系
Figure BDA00015718568600001323
下的
Figure BDA00015718568600001324
轴电流信号
Figure BDA00015718568600001325
和
Figure BDA00015718568600001326
轴电流信号
Figure BDA00015718568600001327
17b)对同步旋转坐标系
Figure BDA00015718568600001328
下的
Figure BDA00015718568600001329
轴电流信号
Figure BDA00015718568600001330
和
Figure BDA00015718568600001331
轴电流信号
Figure BDA00015718568600001332
进行低通滤波，分别得到
Figure BDA00015718568600001333
轴的反馈电流
Figure BDA00015718568600001334
和
Figure BDA00015718568600001335
轴的反馈电流
Figure BDA00015718568600001336
步骤18，根据步骤15和步骤17的结果，获得
Figure BDA00015718568600001337
轴的参考值
Figure BDA00015718568600001338
电压
Figure BDA00015718568600001339
和
Figure BDA00015718568600001340
轴的电压
Figure BDA0001571856860000141
18a)将最终转子角速度估计值
Figure BDA0001571856860000142
作为反馈信号，并将目标角速度
Figure BDA0001571856860000143
与该转子角速度估计值
Figure BDA0001571856860000144
做差，再将差值输入到PI速度调节器，获得
Figure BDA0001571856860000145
轴电流的参考值
Figure BDA0001571856860000146
18b)对
Figure BDA0001571856860000147
轴电流的参考值
Figure BDA0001571856860000148
与
Figure BDA0001571856860000149
轴电流反馈电流
Figure BDA00015718568600001410
做差，并将该差值再经过PI调节器调节，获得
Figure BDA00015718568600001411
轴的电压
Figure BDA00015718568600001412
18c)对
Figure BDA00015718568600001413
轴的电流与
Figure BDA00015718568600001414
轴电流反馈电流
Figure BDA00015718568600001415
做差，并将差值作为PI调节器的输入，获得
Figure BDA00015718568600001416
轴的电压
Figure BDA00015718568600001417
其中
Figure BDA00015718568600001419
轴的电流为0。
步骤19，根据转子位置估计值
Figure BDA00015718568600001420
Figure BDA00015718568600001421
轴的电压
Figure BDA00015718568600001422
和
Figure BDA00015718568600001423
轴的电压
Figure BDA00015718568600001424
获得两相静止坐标系α-β中α轴的电压值
Figure BDA00015718568600001425
和β轴的电压值为
Figure BDA00015718568600001426
对两相估计同步旋转坐标系
Figure BDA00015718568600001427
下的
Figure BDA00015718568600001428
轴电压值
Figure BDA00015718568600001429
和
Figure BDA00015718568600001430
轴电压值
Figure BDA00015718568600001431
进行逆Park变换，获得两相静止坐标系α-β中的α轴的电压值
Figure BDA00015718568600001432
和β轴的电压值为
Figure BDA00015718568600001433
当实际转子位置值θ，与转子位置估计值
Figure BDA00015718568600001434
相差收敛到零时，电机稳定运行。
以上描述仅是本发明的一个具体实例，并未构成对本发明的任何限制，显然对于本领域的专业人员来说，在了解了本发明内容和原理后，都可能在不背离本发明原理、结构的情况下，进行形式和细节上的各种修改和改变，但是这些基于本发明思想的修正和改变仍在本发明的权利要求保护范围之内。

Claims (10)

1.一种基于旋转高频注入算法的永磁同步电机转子位置估计方法，其特征在于，包括：
(1)定义的坐标系：包括两相静止坐标系、估计两相同步旋转坐标系、和真实两相同步旋转坐标系这三个坐标系，其中：α-β为两相静止坐标系，d-q为实际转子同步旋转坐标系，
Figure FDA0002247915220000011
为估计转子同步旋转坐标系；
(2)在两相静止的坐标系中注入不同信号uαhi和uβhi，即在α轴注入两个高频余弦信号uαhi，在β轴注入两个高频正弦信号uβhi
uαhi＝umhcos(ωht)+umh1cos(ωh1t)，
uβhi＝umhsin(ωht)+umh1sin(ωh1t)，
其中，ωh是第一个高频信号的角频率，ωh1为第二个高频信号角频率，umh是信号角频率为ωh的幅值，umh1是信号角频率为ωh1的幅值，高频电压的角频率远远大于基频信号的角频率，且满足ωh1＝2ωh；
(3)检测出永磁同步电机输出的三相电流ia ib和ic；
(4)将三相电流进行Clarke变换，得到两相静止坐标系α-β下的两相定子电流iα和iβ，该两相电流不仅含有基频电流信号，且含有两个注入高频电压产生的高频电流信号；
(5)对两相静止坐标系α-β下的两相定子电流iα和iβ，进行延时为
Figure FDA0002247915220000012
的纯滤波，获得角频率为ωh的两个高频电流信号iαh和iβh，其中iαh表示为α相高频电流iα经滤波之后的输出信号，iβh表示为β相高频电流iα经滤波之后的输出信号；
(6)将得到的角频率为ωh的高频电流信号iαh和iβh通过同步轴高通滤波器进行滤波，获得角频率为ωh的负序高频电流分量iαhn和iβhn，其中iαhn表示为iαh经过同步轴高通滤波之后的输出信号，iβhn表示为iβh经过同步轴高通滤波之后的输出信号；
(7)定义转子估计位置误差为
Figure FDA0002247915220000021
其中
Figure FDA0002247915220000022
为转子位置估计值，θ为实际转子位置值，对负相序高频电流分量iαhn和iβhn进行锁相环处理，得到估计位置偏差信号εh(Δθ)，其中εh(Δθ)是注入角频率为ωh得到的位置偏差信号；
(8)将估计得到位置偏差信号εh(Δθ)作为PI调节器的输入，PI调节器的输出为估计转子角速度
Figure FDA0002247915220000023
对估计转子角速度
Figure FDA0002247915220000024
进行积分调节得到转子位置估计值
Figure FDA0002247915220000025
其中
Figure FDA0002247915220000026
是注入角频率为ωh得到的估计转子角速度
Figure FDA0002247915220000027
Figure FDA0002247915220000028
是注入角频率为ωh得到的转子位置估计值；
(9)对两相静止坐标系α-β下的两相定子电流iα和iβ，进行高通滤波，获得角频率为ωh1的高频电流信号iαh1和iβh1，其中iαh1表示为α相高频电流iα经高通滤波之后的输出信号，iβh1表示为β相高频电流iβ经高通滤波之后的输出信号；
(10)对角频率为ωh1的高频电流信号iαh1和iβh1进行同步轴高通滤波器，获得角频率为ωh1的负序高频电流分量iαhn1和iβhn1；
(11)将负相序高频电流分量iαhn1和iβhn1进行锁相环处理，估计位置偏差信号εh1(Δθ)：
(12)将估计得到位置偏差信号εh1(Δθ)作为PI调节器的输入，PI调节器的输出为估计转子角速度
Figure FDA0002247915220000029
对估计转子角速度
Figure FDA00022479152200000210
进行积分调节，得到转子位置估计值
Figure FDA00022479152200000211
(13)将注入角频率为ωh的高频信号得到的转子角速度估计值
Figure FDA00022479152200000212
和注入角频率为ωh1高频信号得到的转子角速度估计值
Figure FDA00022479152200000213
进行相加求平均值，得到最终转子角速度估计值
Figure FDA00022479152200000214
(14)将注入角频率为ωh高频信号得到的转子位置估计值
Figure FDA00022479152200000215
和注入角频率为ωh1的高频信号得到的转子位置估计值
Figure FDA00022479152200000216
进行相加求平均值得到最终转子位置估计值
Figure FDA00022479152200000217
(15)对两相定子电流iα和iβ进行Park变换，通过低通滤波器，输出为两相估计同步旋转坐标系
Figure FDA0002247915220000031
下的两相电流
Figure FDA0002247915220000032
和
Figure FDA0002247915220000033
(16)将最终转子角速度估计值
Figure FDA0002247915220000034
和最终转子位置估计值
Figure FDA0002247915220000035
作为反馈信号，控制电机实际转子位置值与转子位置估计值相差收敛到零，使电机稳定运动。
2.根据权利要求1所述的方法，其特征在于，步骤(4)中两相静止坐标系α-β下的两相定子电流iα和iβ，分别表示如下：
Figure FDA0002247915220000036
Figure FDA0002247915220000037
其中t是时间变量，iαb为基频余弦电流响应，iαi是角频率为ωh余弦电流响应，iαi1是角频率为ωh1余弦电流响应，iβb为基频正弦电流响应，iβi是角频率为ωh正弦电流响应，iβi1是角频率为ωh1正弦电流响应，IA表示为基频电流幅值，ωr为转子的速度，
Figure FDA0002247915220000038
为基频电流初始相位，I11表示为角频率为ωh正相幅值，I12表示为角频率为ωh负相幅值，I21表示为角频率为ωh1正相幅值，I22表示为角频率为ωh1负相幅值，
Figure FDA0002247915220000039
Figure FDA00022479152200000310
Ld为两相旋转坐标系d轴直轴电感，Lq为两相旋转坐标系q轴的交轴电感。
3.根据权利要求1所述的方法，其特征在于，步骤(5)中α相高频电流iα经纯延时滤波之后的输出信号iαh和β相高频电流iβ经纯延时滤波之后的输出信号iβh，分别表示如下：
Figure FDA0002247915220000041
Figure FDA0002247915220000042
其中t是时间变量，ωr为转子的频率，I11表示为角频率为ωh正相幅值，I12表示为角频率为ωh负相幅值，
Figure FDA0002247915220000043
Figure FDA0002247915220000044
Ld为两相旋转坐标系d轴直轴电感，Lq为两相旋转坐标系q轴的交轴电感。
4.根据权利要求1所述的方法，其特征在于，步骤(6)中α相高频电流iαh经同步轴高通滤波后的输出信号iαhn和β相高频电流iβh经同步轴高通滤波后的输出信号iβhn，分别表示如下：
Figure FDA0002247915220000045
Figure FDA0002247915220000046
其中t是时间变量，ωr为转子的频率，I12表示为角频率为ωh负相幅值，
Figure FDA0002247915220000047
Figure FDA0002247915220000048
Ld为两相旋转坐标系d轴直轴电感，Lq为两相旋转坐标系q轴的交轴电感。
5.根据权利要求1所述的方法，其特征在于，步骤(7)中的估计位置偏差信号εh(Δθ)，表示如下:
Figure FDA0002247915220000051
其中
Figure FDA0002247915220000052
为转子位置的估计值，θr表示真实转子位置的值，I12表示为角频率为ωh负相幅值，Δθ为真实转子位置与转子位置的估计值的差值，
Figure FDA0002247915220000053
Figure FDA0002247915220000054
Ld为两相旋转坐标系d轴直轴电感，Lq为两相旋转坐标系q轴的交轴电感。
6.根据权利要求1所述的方法，其特征在于，步骤(9)中，α相高频电流iα进行高通滤波之后的输出信号iαh1和β相高频电流进行高通滤波之后的输出信号iβh1，分别表示如下：
Figure FDA0002247915220000055
Figure FDA0002247915220000056
其中t是时间变量，ωr为转子的频率，I21表示为角频率为ωh1正相幅值，I22表示为角频率为ωh1负相幅值，
Figure FDA0002247915220000057
Ld为两相旋转坐标系d轴直轴电感，Lq为两相旋转坐标系q轴的交轴电感。
7.根据权利要求1所述的方法，其特征在于，步骤(10)中，α相高频电流iαh1经同步轴高通滤波后的输出信号iαhn1和β相高频电流iβh1经同步轴高通滤波后的输出信号iβhn1，分别表示如下:
Figure FDA0002247915220000061
Figure FDA0002247915220000062
其中t是时间变量，ωr为转子的频率，I22表示为角频率为ωh1负相幅值，
Figure FDA0002247915220000063
Ld为两相旋转坐标系d轴直轴电感，Lq为两相旋转坐标系q轴的交轴电感。
8.根据权利要求1所述的方法，其特征在于，步骤(11)中的估计位置偏差信号εh1(Δθ)，表示如下：
Figure FDA0002247915220000064
其中
Figure FDA0002247915220000065
为注入角频率为ωh1的高频信号得到的转子位置估计值
Figure FDA0002247915220000066
θr表示真实转子位置的值，I12表示为角频率为ωh负相幅值，Δθ为真实转子位置值与转子位置的估计值的差值，
Figure FDA0002247915220000067
Ld为两相旋转坐标系d轴直轴电感，Lq为两相旋转坐标系q轴的交轴电感。
9.根据权利要求1所述的方法，其特征在于：
步骤(13)中的最终转子角速度估计值
Figure FDA0002247915220000068
表示为：
Figure FDA0002247915220000069
步骤(14)中的到最终转子位置估计值
Figure FDA00022479152200000610
表示为：
Figure FDA00022479152200000611
10.根据权利要求1所述的方法，其特征在于，所述步骤(16)的具体实现如下：
16a)利用最终转子位置估计值
Figure FDA00022479152200000612
对两相静止坐标系α-β下的两相定子电流iα和iβ进行Park变换，再分别进行低通滤波，获得两相同步旋转坐标系
Figure FDA00022479152200000613
下的
Figure FDA00022479152200000614
轴电流反馈信号
Figure FDA00022479152200000615
和
Figure FDA00022479152200000616
轴电流反馈信号
Figure FDA00022479152200000617
16b)将最终转子角速度估计值
Figure FDA0002247915220000071
作为反馈信号，并将目标角速度
Figure FDA0002247915220000072
与该转子角速度估计值
Figure FDA0002247915220000073
做差，再将差值输入到PI速度调节器，获得q轴电流的参考值
Figure FDA0002247915220000074
16c)将该参考值
Figure FDA0002247915220000075
与
Figure FDA0002247915220000076
轴电流反馈信号
Figure FDA0002247915220000077
做差，将该差值再经过PI调节器，获得q轴电压
Figure FDA0002247915220000078
16d)将d轴的电流0与
Figure FDA0002247915220000079
轴电流反馈信号
Figure FDA00022479152200000710
做差，将差值作为PI调节器的输入，获得d轴的电压
Figure FDA00022479152200000711
16e)利用转子位置估计值
Figure FDA00022479152200000712
两相同步旋转坐标系d-q下的d轴电压
Figure FDA00022479152200000713
和q轴电压
Figure FDA00022479152200000714
进行逆Park变换，获得两相静止坐标系α-β，其中α轴的电压值为
Figure FDA00022479152200000715
和β轴的电压值为
Figure FDA00022479152200000716
CN201810120333.6A 2018-02-07 2018-02-07 基于旋转高频注入算法的永磁同步电机转子位置估计方法 Active CN108206659B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810120333.6A CN108206659B (zh) 2018-02-07 2018-02-07 基于旋转高频注入算法的永磁同步电机转子位置估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810120333.6A CN108206659B (zh) 2018-02-07 2018-02-07 基于旋转高频注入算法的永磁同步电机转子位置估计方法

Publications (2)

Publication Number Publication Date
CN108206659A CN108206659A (zh) 2018-06-26
CN108206659B true CN108206659B (zh) 2020-04-07

Family

ID=62605682

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810120333.6A Active CN108206659B (zh) 2018-02-07 2018-02-07 基于旋转高频注入算法的永磁同步电机转子位置估计方法

Country Status (1)

Country Link
CN (1) CN108206659B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109150052B (zh) * 2018-08-13 2020-04-24 哈尔滨工业大学 一种考虑数字控制延迟的随机正弦注入永磁电机转子位置观测器的高频噪音抑制方法
CN109358228B (zh) * 2018-11-09 2020-12-15 哈工大（张家口）工业技术研究院 基于双增强型锁相环的电网电压正负序分量实时估计方法
CN109921712A (zh) * 2019-02-26 2019-06-21 浙江大学 基于注入高频脉振电压的永磁同步电机双闭环i/f控制方法
CN109889117B (zh) * 2019-04-04 2020-09-01 合肥工业大学 基于旋转高频注入法的ipmsm位置观测方法、系统及驱动系统
CN110460270B (zh) * 2019-07-24 2021-03-16 东南大学盐城新能源汽车研究院 一种考虑交叉饱和电感的无位置传感器控制方法
CN112564572B (zh) * 2020-11-20 2022-05-03 江苏科技大学 基于电流偏置的永磁同步电机转子位置检测系统及方法
CN113517837B (zh) * 2021-04-08 2023-07-18 华南理工大学 基于高频脉振注入法的pmsm转速提取方法
CN113467418B (zh) * 2021-06-25 2022-06-28 河北工业大学 一种控制环路性能指标测量方法
CN114123901B (zh) * 2021-11-24 2023-12-01 南京航空航天大学 一种抑制转矩脉动的逆二倍频同步坐标脉振高频注入的六相永磁电机无位置传感器控制方法
CN114050755B (zh) * 2022-01-12 2022-03-29 希望森兰科技股份有限公司 基于高频旋转电压注入的永磁同步电机位置观测改进算法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014200148A (ja) * 2013-03-29 2014-10-23 株式会社デンソー 回転機の制御装置
JP2016025714A (ja) * 2014-07-18 2016-02-08 富士電機株式会社 永久磁石形同期電動機の制御装置
CN105450127A (zh) * 2015-11-23 2016-03-30 南京航空航天大学 基于高频信号注入的永磁同步电机转子位置检测方法
CN105450125A (zh) * 2015-11-20 2016-03-30 南京航空航天大学 一种旋转高频电压信号注入法永磁同步电机转子位置估算误差的补偿方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014200148A (ja) * 2013-03-29 2014-10-23 株式会社デンソー 回転機の制御装置
JP2016025714A (ja) * 2014-07-18 2016-02-08 富士電機株式会社 永久磁石形同期電動機の制御装置
CN105450125A (zh) * 2015-11-20 2016-03-30 南京航空航天大学 一种旋转高频电压信号注入法永磁同步电机转子位置估算误差的补偿方法
CN105450127A (zh) * 2015-11-23 2016-03-30 南京航空航天大学 基于高频信号注入的永磁同步电机转子位置检测方法

Also Published As

Publication number Publication date
CN108206659A (zh) 2018-06-26

Similar Documents

Publication Publication Date Title
CN108206659B (zh) 基于旋转高频注入算法的永磁同步电机转子位置估计方法
CN107482977B (zh) 一种永磁同步电机转子位置和转速检测方法
CN110198150B (zh) 一种永磁同步电机多参数在线辨识方法
CN104393811B (zh) 一种提高永磁同步电机转子位置检测精度的方法
CN112737450B (zh) 一种用于spmsm转子位置估计的高频注入补偿方法
EP1852967B1 (en) Apparatus for controlling high speed operation of motor and method thereof
CN106788054B (zh) 一种基于旋转高频注入法的无速度传感器控制方法
CN106533300B (zh) 一种基于速度环模糊控制和高频注入法的无传感器控制系统
CN112713834B (zh) 一种永磁同步电机无位置传感器控制方法及系统
CN111464105B (zh) 一种永磁同步电机初始位置检测方法
CN114598206B (zh) 永磁同步电机宽速域转子位置观测器设计方法
CN110620533A (zh) 一种表贴式永磁同步电机无传感器控制方法
CN115276487A (zh) 一种永磁同步电机控制方法
CN108512475A (zh) 一种基于脉振高频电压注入法的转子位置检测方法
CN115765563A (zh) 一种表贴式永磁同步电机转子位置及转速信息检测方法
Kang et al. Improved sliding mode observer based sensorless control for PMSM
CN110768590A (zh) 一种采用改进滑模观测器的永磁同步电机转子位置和速度估算系统及方法
CN105932924B (zh) 异步电机转子电角速度观测方法
CN107769655B (zh) 永磁同步电机转速估算方法、装置、计算设备及存储介质
CN114337416A (zh) 电机控制方法、装置、压缩机、存储介质及空调器
CN111224600B (zh) 永磁同步电机速度控制方法、装置、计算机设备及介质
CN113965129A (zh) 一种永磁同步电机控制系统电流测量偏移误差的补偿方法
CN110995092A (zh) 一种应用于pmsm无位置传感器控制的磁极判断方法
CN111900908B (zh) 一种基于无差拍反电动势预测器的永磁同步电机转子位置和速度估算方法
CN116979838A (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
PE01 Entry into force of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: Rotor position estimation of permanent magnet synchronous motor based on rotating high frequency injection algorithm

Effective date of registration: 20221104

Granted publication date: 20200407

Pledgee: CITIC Bank Limited by Share Ltd. Lishui branch

Pledgor: LISHUI BOYUAN TECHNOLOGY Co.,Ltd.

Registration number: Y2022980020826