CN103956956A - 一种基于扩展卡尔曼滤波器的无刷直流电机状态估计方法 - Google Patents
一种基于扩展卡尔曼滤波器的无刷直流电机状态估计方法 Download PDFInfo
- Publication number
- CN103956956A CN103956956A CN201410200184.6A CN201410200184A CN103956956A CN 103956956 A CN103956956 A CN 103956956A CN 201410200184 A CN201410200184 A CN 201410200184A CN 103956956 A CN103956956 A CN 103956956A
- Authority
- CN
- China
- Prior art keywords
- interval
- kalman filter
- phase
- apex coordinate
- extended kalman
- 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
Landscapes
- Control Of Motors That Do Not Use Commutators (AREA)
Abstract
本发明公开了一种基于扩展卡尔曼滤波器的无位置传感器无刷直流电机状态估计方法。使用本发明能够在***噪声和测量噪声为非零均值白噪声以及***模型存在误差的情况降低调试难度,提高观测精度。本发明通过模糊规则调整测量误差协方差矩阵R,从而大幅度降低了***调试难度,并提高了精度;同时在状态更新方程中加入一个衰减因子,调整了新旧观测量对估计值的修正权重,提高了跟踪性能;此外,通过将p(0),Q的值做适当分组调整的方法,抑制严重失真情况并提高扩展卡尔曼滤波器的估计准确性。
Description
技术领域
本发明涉及无刷直流电机的无位置传感器控制技术领域,具体涉及一种基于扩展卡尔曼滤波器的无位置传感器无刷直流电机状态估计方法。
背景技术
永磁无刷直流电机(Brushless DC Motor,BLDCM)融合了直流电动机控制方式简单,转矩特性好,调速性能好,以及交流电动机制造简单,无励磁损耗,功率密度高等优点,所以在各个领域都有应用,并引起国内外学者的广泛关注与研究。但是无刷直流电机的运行需要转子位置信号,传统的方法是采用位置传感器。随着对额定功率更小的电机的研究与需求,位置传感器的体积占整个电机***体积的百分比越来越大,这就使得整个***难以实现小型化,限制了其在精密微小场合的应用。此外,由于位置传感器的存在,使得***复杂程度加大、成本增加、可靠性降低、增加了生产和维护难度并且对电机的制造工艺也带来不利的影响。
针对位置传感器对无刷直流电机控制***带来的上述不利影响,比较好的解决办法是采用无位置传感器控制技术。所谓的无位置传感器控制,实质上是无机械位置传感器的控制,即在电机的运转过程中,由新的位置信号检测方法(一般通过软件与硬件相互结合的技术)代替机械的位置传感器来提供逆变桥功率器件换相导通的时序信号。无位置传感器控制不仅可以克服机械位置传感器带来的问题,还可以充分利用软件技术改善并提高电机的动、静态性能。
所以在无刷直流电机的无位置传感器控制技术中,状态估计是最关键的核心问题,也就是如何不通过机械的位置传感器,而只通过三相电压与电流,准确、快速、可靠的得到转子位置与速度信号。
Gamazo-Real JC等人在文献(Position and Speed Control of Brushless DCMotors Using Sensorless Techniques and Application Trends[J].SENSORS,2010,10(7):6901-6947.)中总结,现有的无位置传感器无刷直流电机的状态估计方法主要有传统方法和状态观测器方法两种。其中传统方法主要有:反电势过零法、三次谐波电压积分法、续流二级管法、反电势积分法等;状态观测器法主要有:扩展卡尔曼滤波法、参考模型自适应法、神经网络法和滑模观测器法等。
传统方法主要存在的问题是在某一特定的速度范围内可以较好的估计出转子的速度、位置信息,但是在电机低速或高速运行时,观测效果大幅下降甚至无法使用。
近年来的主要研究和应用方向都集中在利用状态观测器对无刷直流电机***进行状态估计,无刷直流电机***是一个典型的非线性的***,结合扩展卡尔曼滤波器,可以对其进行状态估计。
卡尔曼滤波器有一系列递推数学公式描述。他们提供了一种高效可计算的方法来估计过程的状态,并使估计均方误差最小。卡尔曼滤波器应用广泛且功能强大;它可以估计信号的过去和当前状态,甚至能估计将来的状态,即使并不知道模型的确切性质。但是卡尔曼滤波器用于无刷直流电机***的估计存在一些问题。首先,由于卡尔曼滤波器的参数是在某一特定转速(如额定转速)下设置的,当速度变化时,会导致观测精度下降。其次,在扩展卡尔曼滤波器的设计过程中,假设了***噪声和测量噪声为精确已知的零均值白噪声,而当***噪声与测量噪声不满足此条件时,可能会导致误差协方差矩阵变大,降低了估计精度,如果实时人工调整滤波器的参数,调试难度很大。再次,扩展卡尔曼滤波的速度跟踪性能不是很好,主要体现在估计的转速比实际转速存在一个明显的相位延迟,造成这种现象的原因主要是新旧观测量对估计值的修正权重不一样。
发明内容
有鉴于此,本发明提供了一种基于扩展卡尔曼滤波器的无位置传感器无刷直流电机状态估计方法,能够在***噪声和测量噪声为非零均值白噪声以及***模型存在误差的情况降低调试难度,提高观测精度。
本发明的基于扩展卡尔曼滤波器的无刷直流电机状态估计方法包括如下步骤:
步骤一,对传统的扩展卡尔曼滤波器进行改进,采用模糊算法调整测量误差协方差矩阵Rk,Rk=R+ΔR,其中,R为***噪声wk的协方差矩阵,ΔR为根据模糊规则获得的调整量;
其中,模糊规则为:
以测量变量与预测值之差的理论值与实际值之差ek作为模糊输入,测量噪声方差阵Rk的调整量ΔR作为模糊输出;
当ek≈0时,ΔR为0;
当ek>0时,ΔR为负值;
当ek<0时,ΔR为正值;
步骤二,采用步骤一改进后的的扩展卡尔曼滤波器对无刷直流电机的状态进行估计。
其中,所述模糊算法中,ek和ΔR的隶属度函数均分为负大NB、负中NM、负小NS、零ZO、正小PS、正中PM和正大PB7个区间,每个区间的隶属度函数均为三角形,其中,模糊输入ek隶属度函数中,NB区间三角形三个顶点坐标为(-∞,1)、(-4.8,1)、(-2.5,0);NM区间三角形三个顶点坐标为(-4.8,0)、(-3,1)、(-2,0);NS区间三角形三个顶点坐标为(-3.2,0)、(-2,1)、(0,0);ZO区间三角形三个顶点坐标为(-2.6,0)、(0,1)、(2.6,0);PS区间三角形三个顶点坐标为(0,0)、(2,1)、(3.2,0);PM区间三角形三个顶点坐标为(2,0)、(3,1)、(4.8,0);PB区间三角形三个顶点坐标为(2.5,0)、(4.8,1)、(+∞,1);模糊输出ΔR隶属度函数中,NB区间三角形三个顶点坐标为(-∞,1)、(-2.8,1)、(-1.8,0);NM区间三角形三个顶点坐标为(-2.8,0)、(-2,1)、(-1.3,0);NS区间三角形三个顶点坐标为(-.2,0)、(-1,1)、(0,0);ZO区间三角形三个顶点坐标为(-1.3,0)、(0,1)、(1.3,0);PS区间三角形三个顶点坐标为(0,0)、(1.3,1)、(2,0);PM区间三角形三个顶点坐标为(1.3,0)、(2,1)、(2.8,0);PB区间三角形三个顶点坐标为(1.8,0)、(2.8,1)、(+∞,1);
模糊规则为:
(1)如果ek∈ZO,则令ΔR∈ZO;
(2)如果ek∈NS,则令ΔR∈PS;
(3)如果ek∈NM,则令ΔR∈PM;
(4)如果ek∈NB,则令ΔR∈PB;
(5)如果ek∈PS,则令ΔR∈NS;
(6)如果ek∈PM,则令ΔR∈NM;
(7)如果ek∈PB,则令ΔR∈NB。
进一步地,可以在扩展卡尔曼滤波器的状态更新方程中加入一个衰减因子γ,即扩展卡尔曼滤波器的状态更新方程为
其中,为先验估计误差的协方差;为后验估计误差的协方差;为已知测量变量zk时k时刻的后验状态估计;Qk为过程噪声wk的协方差矩阵;
ia、ib、ic为定子三相的绕组电流,ea、eb、ec为定子三相绕组的反电动势,ω为转子速度,θ为转子位置,uab为a相与b相的相电压,ubc为b相与c相的相电压,uca为c相与a相的相电压,L为电机每相绕组的自感,M为电机每两相绕组间的互感,R为电机每相定子绕组的电阻,T为采样周期,衰减因子γ在1.05~1.15之间。
进一步地,可以在状态更新方程中根据电机的转子速度对***误差协方差矩阵Q、滤波器初值p(0)进行分组,其中,分别将转子速度小于额定转速30%、转子速度为额定转速30%-60%和转子速度为额定转速60%-120%分为三个区间,根据仿真实验结果分别选取p(0),Q的值。
有益效果:
(1)本发明针对***噪声和测量噪声非零均值白噪声,以及***模型存在误差造成的估计精度下降,调试难度大的问题,通过模糊规则调整测量误差协方差矩阵R,从而大幅度降低了***调试难度,并提高了精度。
(2)本发明针对速度跟踪性能不好的问题,在状态更新方程中加入一个衰减因子,调整了新旧观测量对估计值的修正权重,提高了跟踪性能。
(3)本发明针对速度变化影响卡尔曼滤波器估计精度的问题,通过将p(0),Q的值做适当分组调整的方法,可以抑制严重失真情况并提高扩展卡尔曼滤波器的估计准确性。
(4)本发明可以对以上三种改进进行有机结合,应用到无刷直流电机***中取得了较好的观测效果。
附图说明
图1为本发明流程图。
图2为模糊输入ek隶属度函数。
图3为模糊输出ΔR隶属度函数。
图4为扩展卡尔曼滤波器无刷直流电机状态观测的仿真结构图。
图5为本发明中的方法与传统的扩展卡尔曼滤波器方法的调速对比图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提供了一种基于扩展卡尔曼滤波器的无位置传感器无刷直流电机状态估计方法,对传统的扩展卡尔曼滤波器进行了改进,采用以模糊规则调整测量误差协方差矩阵R,以速度为准对滤波器初值p(0)、***误差协方差矩阵Q分组并加入衰减因子等方法,在几乎不增加卡尔曼滤波器结构复杂度的前提下,提高了扩展卡尔曼滤波器在无刷直流电机存在建模误差与非理想噪声时的估计效果。具体实现步骤如下:
步骤一,控制对象建模。
本发明控制对象为无刷直流电机。以两相导通星形三相六状态的无刷直流电机为例,无刷直流电机的电压平衡方程如下:
其中,ia、ib、ic为定子三相的绕组电流(单位为A);ua、ub、uc为定子三相的绕组电压(单位为V);ea、eb、ec为定子三相绕组的反电动势(单位为V);L为电机每相绕组的自感(单位为H);M为电机每两相绕组间的互感(单位为H),R为电机每相定子绕组的电阻(单位为Ω),表示微分。
无刷直流电机的随机离散模型为:
式中 ω为转子速度,θ为转子位置,uab为a相与b相的相电压,ubc为b相与c相的相电压,uca为c相与a相的相电压,wk、vk分别为***噪声和测量噪声,这里假设两者均为零均值的高斯白噪声,协方差矩阵为R,Q。
ea=KeT(a11+a12(θ-θkom))
eb=KeT(a21+a22(θ-θkom))
ec=KeT(a31+a32(θ-θkom))
式中Ke为反电势系数,T为采样周期,aij、θkom的取值如表1所示。
表1反电势表达式系数表
步骤二,构建扩展卡尔曼滤波器。
传统的扩展卡尔曼滤波器的方程为:
其中,
其中,(“^”代表估计)为在已知k时刻状态的情况下对k+1时刻的先验状态估计;为已知测量变量zk时k时刻的后验状态估计,zk∈Rn是k时刻的测量值;uk为k时刻***的控制量。
为后验估计误差的协方差;为先验估计误差的协方差:
其中,为先验估计误差, 为后验估计误差,
Q为过程噪声wk(也称为***噪声)的协方差矩阵;R为观测噪声vk(也称为测量噪声或量测噪声)的协方差矩阵;I为单位矩阵;Kk为k时刻测量变量与预测之差(即定义为测量过程的残余,或称为革新)的增益(或称为混合因数),使得后验估计误差最小。
传统的扩展卡尔曼滤波器在无刷直流电机***应用中存在三个问题:第一,在扩展卡尔曼滤波器的设计过程中,假设了***噪声和测量噪声为精确已知的零均值白噪声,而在现实***中,由于受温度、使用环境、负载变化等因素的影响,都会引起电机的参数产生变化,在扩展卡尔曼滤波器的设计过程中所采用的电机模型,也不可避免的由于参数变化或未建模态的存在,这些都会导致***噪声与测量噪声不再为零均值白噪声,会降低扩展卡尔曼滤波器的估计精度,如果实时调整***参数,难度很大;第二,由于卡尔曼滤波器的参数是在某一特定转速(如额定转速)下设置的,当速度变化较大时(如速度变化超过30%额定转速的范围),会导致观测精度大幅下降;第三,传统扩展卡尔曼滤波的速度跟踪性能不是很好,主要体现在估计的转速比实际转速存在一个明显的相位延迟,而且估计转速峰值小于实际转速峰值,如果闭环则相当于加入一个较大的惯性环节,造成这种现象的原因主要是新旧观测量对估计值的修正权重不一样。因此,本发明针对以上三个问题,结合模糊控制理论,强跟踪算法等方法,提出了扩展卡尔曼滤波器。
针对实际电机***中测量噪声不是零均值高斯白噪声的情况,结合模糊控制理论对卡尔曼滤波的数学模型进行处理,通过模糊控制中蕴含的专家***弥补噪声的不准确性,不仅提高了估计精度,还可以通过模糊***实时调整***参数,大幅度降低了调试难度。专家***主要由模糊化、知识库、模糊推理和反模糊化4个模块构成的。
本发明中,隶属度函数选择为三角形。模糊规则的形式采用if A is a then B isb,其中A与B是语言变量(linguistic variable),a和b是由隶属函数映射到的语言值(linguistic values),控制规则根据控制专家的经验或大量实验总结得出。
根据残余的定义(即测量变量与预测值之差),为这一时刻的实际量减去上一时刻对这一时刻的预测量。在理想状况下,残余应该是一个零均值白噪声。如果残余不满足这个条件,则说明***噪声或者测量噪声发生了变化,这是需要对测量噪声方差阵Rk和残余增益Kk进行调整,而由于Rk包含在Kk的求解过程中,所以只需要对Rk进行调整即可。(在传统卡尔曼滤波器中,测量噪声方差阵R为定值,即为常数阵,而在本发明中,将R调整为根据模糊规则变化的在不同时间数值不同的矩阵,故表示为Rk)。
以残余的理论值与实际值的差作为模糊输入,测量噪声方差阵Rk的调整量ΔR作为模糊输出。设新息的理论值为Lk,残余的实际值为Sk,残余的理论与实际值的差为ek=Lk-Sk。
1)当ek≈0时,说明Lk与Sk几乎相等,Rk不用调整,ΔR为0。
2)当ek>0时,说明Lk大于Sk,需要减小Rk,ΔR为负值。
3)当ek<0时,说明Lk小于Sk,需要增大Rk,ΔR为正值。
经过大量实验,调整隶属度函数规则的数量,以及函数的顶点,起点与形状,最终根据扩展卡尔曼滤波器的观测效果,选取最理想的数据,建立如下的隶属度函数,模糊输入ek隶属度函数如图2所示,模糊输出ΔR隶属度函数如图3所示。图中,NB表示负大,NM表示负中,NS表示负小,ZO表示零,PS表示正小,PM表示正中,PB表示正大。模糊输入ek隶属度函数的各个区间均为三角形,NB区间三角形三个顶点坐标为(-∞,1)、(-4.8,1)、(-2.5,0);NM区间三角形三个顶点坐标为(-4.8,0)、(-3,1)、(-2,0);NS区间三角形三个顶点坐标为(-3.2,0)、(-2,1)、(0,0);ZO区间三角形三个顶点坐标为(-2.6,0)、(0,1)、(2.6,0);PS区间三角形三个顶点坐标为(0,0)、(2,1)、(3.2,0);PM区间三角形三个顶点坐标为(2,0)、(3,1)、(4.8,0);PB区间三角形三个顶点坐标为(2.5,0)、(4.8,1)、(+∞,1)。模糊输出ΔR隶属度函数的各个区间均为三角形,NB区间三角形三个顶点坐标为(-∞,1)、(-2.8,1)、(-1.8,0);NM区间三角形三个顶点坐标为(-2.8,0)、(-2,1)、(-1.3,0);NS区间三角形三个顶点坐标为(-.2,0)、(-1,1)、(0,0);ZO区间三角形三个顶点坐标为(-1.3,0)、(0,1)、(1.3,0);PS区间三角形三个顶点坐标为(0,0)、(1.3,1)、(2,0);PM区间三角形三个顶点坐标为(1.3,0)、(2,1)、(2.8,0);PB区间三角形三个顶点坐标为(1.8,0)、(2.8,1)、(+∞,1)。
本发明建立模糊规则的基本思想为(1)当误差为负大时,若误差变化为负,这时误差有增大的趋势,为尽快消除已有的负大误差并抑制误差变大,所以控制量的取负大。当误差为负大而误差变化为正时,***本身已有减小误差的趋势,为尽快消除误差且又不引起超调,应取较小的控制量。(2)当误差为负中时,控制量应使误差尽快消除,取值与误差为负大时相同。(3)当误差为负小时,***接近稳态,若误差变化为负,选取控制量为负中,以抑制误差往负方向变化,若误差变化为正时,***本身已有趋势消除负小的偏差,选取控制量为零或正小即可。(4)当误差为正时,控制思想与此基本相同,仅符号相反。
根据上述规则结合本文中的无刷直流电机***,最终的模糊规则表示为:
(1)IF ek∈ZO,THENΔR∈ZO
(2)IF ek∈NS,THENΔR∈PS
(3)IF ek∈NM,THENΔR∈PM
(4)IF ek∈NB,THENΔR∈PB
(5)IF ek∈PS,THENΔR∈NS
(6)IF ek∈PM,THENΔR∈NM
(7)IF ek∈PB,THENΔR∈NB
同时,以速度为界限将滤波器初值p(0),***噪声Q进行分组。通过大量实验发现,滤波器初值p(0),***噪声Q在不同速度下的观测结果差异较大。在某一特定速度下,调整好p(0)与Q的值,可以得到比较准确的结果。但是当给定速度发生较大变化时(例如超过30%额定转速的速度变化),速度估计的误差会增大,有时甚至出现严重失真。电机的实际应用中,往往要求电机有较广的调速范围,而固定的p(0)与Q的值,在不同速度下的估计结果差异较大,显然不能满足需求。考虑到不再增加扩展卡尔曼滤波器的复杂度,以及p(0)与Q的值对较小速度变化影响较小的情况,采用了将p(0)、Q的值做适当分组调整的方法,可以抑制严重失真情况并提高扩展卡尔曼滤波器的估计准确性。
本发明将速度分为小于额定转速30%、额定转速30%-60%和额定转速60%-120%三个区间,分别选取适合的p(0),Q的值,可以使估计结果在各个速度时都保证一定的精确性。将p(0),Q按照额定转速0-30%,30%-60%,60%-120%分为3组,p(0)1、p(0)2、p(0)3与Qk1、Qk2和Qk3,其值根据仿真实验获取。
经试验总结出规律,在高速运行时,p(0),Q值选取过小,则无法得出正确的估计结果,在较低速运行时,p(0),Q选取过大,可以得出估计结果,但是稳态误差较大,故应该随着速度的增大,适当增大p(0),Q的值。
传统扩展卡尔曼滤波的速度跟踪性能不是很好,主要体现在估计的转速比实际转速存在一个明显的相位延迟,而且估计转速峰值小于实际转速峰值,如果闭环则相当于加入一个较大的惯性环节,造成这种现象的原因主要是新旧观测量对估计值的修正权重不一样。借鉴于强跟踪算法,可以在状态更新方程中加入一个衰减因子,使新测量值对估计值的修正作用比旧测量值大。但是加入衰减因子,也同时加重了***的不稳定性,而且衰减因子过大会加大速度估计的稳态误差。根据经验,选取衰减因子在1.05~1.15之间,可以在减少相位延迟的情况下,保证***稳定性。
综上所述,本发明结合扩展卡尔曼滤波器原理,经过以模糊规则调整测量误差协方差矩阵R,以速度为界限对滤波器初值p(0)、***误差协方差矩阵Q分组并加入衰减因子等方法予以改进,本发明的扩展卡尔曼滤波器为:
式中γ为衰减因子。带“^”为其观测值。
根据式(8)-(12)和递推初值p(0)i,就可以进行扩展卡尔曼迭代运算了。与传统方法的区别在于:第一,相比于式(4),在式(9)中加入衰减因子,并且将常数矩阵Q根据速度分组改变为矩阵Qki;第二,相比于式(5),在式(10)中将常数矩阵R根据模糊规则调整修改为矩阵Rk;第三,递推初值由常值p(0)修改为根据速度分组的p(0)i。
步骤三,采用步骤二中改进的扩展卡尔曼滤波器对无位置传感器无刷直流电机的状态进行估计。
下面进行仿真实验,电机参数:额定电压U=160V,极对数p=3,定子电阻R=0.086Ω,定子电感L=0.0012mh,定子自感M=0.00037mH,转子转动惯量J=0.00021405kgm2。反电势系数Ke=0.1805V/(r·min-1)。
扩展卡尔曼滤波器无刷直流电机无位置传感器控制***如图4所示,主要分为以下三个模块:电机***模块(图中的电机***BLDCM),速度控制模块(图中的PID控制器),扩展卡尔曼滤波器模块(图中的EKF)。其中扩展卡尔曼滤波器模块通过S函数实现。电机本体模块即为现实中无刷直流电机的仿真,扩展卡尔曼滤波器模块通过三相电压与电流,估计出电机的转子位置,速度,作为速度控制器的反馈,以及参考电流模块的输入。EKF输入为电机***BLDCM的三相电压与电流,输出为位置与速度估计值,速度估计值与参考速度值做差,作为PID控制器的输入,电机***BLDCM的输入为PID控制器的输入,EKF的估计位置,以及负载转矩。
图5为调速实验图,0-0.4s给定速度为w=50rad/s,0.4s-0.8s时给定速度为w=100rad/s,0.8到1.2s时给定速度为w=200rad/s,图中实线为扩展卡尔曼滤波器的估计结果,虚线为传统扩展卡尔曼滤波器的估计结果。可见,扩展卡尔曼滤波器在调速的过程中,也有比传统卡尔曼滤波器更优越的估计性能,稳态误差缩小三分之二以上,达到3%以内,上升时间缩短0.02秒以上。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种基于扩展卡尔曼滤波器的无刷直流电机状态估计方法,其特征在于包括如下步骤:
步骤一,对传统的扩展卡尔曼滤波器进行改进,采用模糊算法调整测量误差协方差矩阵Rk,Rk=R+ΔR,其中,R为***噪声wk的协方差矩阵,ΔR为根据模糊规则获得的调整量;
其中,模糊规则为:
以测量变量与预测值之差的理论值与实际值之差ek作为模糊输入,测量噪声方差阵Rk的调整量ΔR作为模糊输出;
当ek≈0时,ΔR为0;
当ek>0时,ΔR为负值;
当ek<0时,ΔR为正值;
步骤二,采用步骤一改进后的的扩展卡尔曼滤波器对无刷直流电机的状态进行估计。
2.如权利要求1所述的无位置传感器无刷直流电机的扩展卡尔曼滤波器,其特征在于,所述模糊算法中,ek和ΔR的隶属度函数均分为负大NB、负中NM、负小NS、零ZO、正小PS、正中PM和正大PB7个区间,每个区间的隶属度函数均为三角形,其中,模糊输入ek隶属度函数中,NB区间三角形三个顶点坐标为(-∞,1)、(-4.8,1)、(-2.5,0);NM区间三角形三个顶点坐标为(-4.8,0)、(-3,1)、(-2,0);NS区间三角形三个顶点坐标为(-3.2,0)、(-2,1)、(0,0);ZO区间三角形三个顶点坐标为(-2.6,0)、(0,1)、(2.6,0);PS区间三角形三个顶点坐标为(0,0)、(2,1)、(3.2,0);PM区间三角形三个顶点坐标为(2,0)、(3,1)、(4.8,0);PB区间三角形三个顶点坐标为(2.5,0)、(4.8,1)、(+∞,1);模糊输出ΔR隶属度函数中,NB区间三角形三个顶点坐标为(-∞,1)、(-2.8,1)、(-1.8,0);NM区间三角形三个顶点坐标为(-2.8,0)、(-2,1)、(-1.3,0);NS区间三角形三个顶点坐标为(-.2,0)、(-1,1)、(0,0);ZO区间三角形三个顶点坐标为(-1.3,0)、(0,1)、(1.3,0);PS区间三角形三个顶点坐标为(0,0)、(1.3,1)、(2,0);PM区间三角形三个顶点坐标为(1.3,0)、(2,1)、(2.8,0);PB区间三角形三个顶点坐标为(1.8,0)、(2.8,1)、(+∞,1);
模糊规则为:
(1)如果ek∈ZO,则令ΔR∈ZO;
(2)如果ek∈NS,则令ΔR∈PS;
(3)如果ek∈NM,则令ΔR∈PM;
(4)如果ek∈NB,则令ΔR∈PB;
(5)如果ek∈PS,则令ΔR∈NS;
(6)如果ek∈PM,则令ΔR∈NM;
(7)如果ek∈PB,则令ΔR∈NB。
3.如权利要求1或2所述的无位置传感器无刷直流电机的扩展卡尔曼滤波器,其特征在于,在扩展卡尔曼滤波器的状态更新方程中加入一个衰减因子γ,即扩展卡尔曼滤波器的状态更新方程为
其中,为先验估计误差的协方差;为后验估计误差的协方差;为已知测量变量zk时k时刻的后验状态估计;Qk为过程噪声wk的协方差矩阵;
ia、ib、ic为定子三相的绕组电流,ea、eb、ec为定子三相绕组的反电动势,ω为转子速度,θ为转子位置,uab为a相与b相的相电压,ubc为b相与c相的相电压,uca为c相与a相的相电压,L为电机每相绕组的自感,M为电机每两相绕组间的互感,R为电机每相定子绕组的电阻,T为采样周期,衰减因子γ在1.05~1.15之间。
4.如权利要求1或2所述的无位置传感器无刷直流电机的扩展卡尔曼滤波器,其特征在于,在状态更新方程中根据电机的转子速度对***误差协方差矩阵Q、滤波器初值p(0)进行分组,其中,分别将转子速度小于额定转速30%、转子速度为额定转速30%-60%和转子速度为额定转速60%-120%分为三个区间,根据仿真实验结果分别选取p(0),Q的值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410200184.6A CN103956956B (zh) | 2014-05-13 | 2014-05-13 | 一种基于扩展卡尔曼滤波器的无刷直流电机状态估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410200184.6A CN103956956B (zh) | 2014-05-13 | 2014-05-13 | 一种基于扩展卡尔曼滤波器的无刷直流电机状态估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103956956A true CN103956956A (zh) | 2014-07-30 |
CN103956956B CN103956956B (zh) | 2017-01-04 |
Family
ID=51334199
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410200184.6A Active CN103956956B (zh) | 2014-05-13 | 2014-05-13 | 一种基于扩展卡尔曼滤波器的无刷直流电机状态估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103956956B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104993765A (zh) * | 2015-08-04 | 2015-10-21 | 重庆大学 | 一种无刷直流电机的转速估计方法 |
CN106019045A (zh) * | 2016-05-16 | 2016-10-12 | 安徽大学 | 一种永磁同步电机匝间短路故障诊断的方法 |
CN106533297A (zh) * | 2016-12-22 | 2017-03-22 | 东南大学 | 一种基于卡尔曼滤波的定子永磁型记忆电机调磁方法 |
CN108490472A (zh) * | 2018-01-29 | 2018-09-04 | 哈尔滨工程大学 | 一种基于模糊自适应滤波的无人艇组合导航方法 |
CN108649848A (zh) * | 2018-05-18 | 2018-10-12 | 湖南大学 | 基于扩展卡尔曼滤波方法的磁齿轮电机失步预测控制方法 |
CN110556927A (zh) * | 2019-06-19 | 2019-12-10 | 国家电网有限公司 | 一种感应式无线电能传输***的状态观测的方法及*** |
CN111308522A (zh) * | 2020-03-27 | 2020-06-19 | 上海天好信息技术股份有限公司 | 一种基于新型Kalman滤波的多维时空数据估计方法 |
CN112858929A (zh) * | 2021-03-16 | 2021-05-28 | 上海理工大学 | 一种基于模糊逻辑与扩展卡尔曼滤波的电池soc估计方法 |
CN113381665A (zh) * | 2021-06-09 | 2021-09-10 | 上海晟矽微电子股份有限公司 | 电机控制方法及装置、电机和电动工具 |
CN114257150A (zh) * | 2022-01-26 | 2022-03-29 | 合肥倍豪海洋装备技术有限公司 | 一种基于改进卡尔曼滤波器的永磁同步电机死区补偿方法 |
CN114744938A (zh) * | 2022-04-11 | 2022-07-12 | 华中科技大学 | 一种基于卡尔曼滤波的全参数观测器及全参数辨识方法 |
CN116137505A (zh) * | 2023-04-18 | 2023-05-19 | 深圳市浮思特科技有限公司 | 一种无刷直流电机的卡尔曼滤波矢量控制***及方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5940814A (en) * | 1997-05-01 | 1999-08-17 | Nec Research Institute, Inc. | Defuzzification strategy and parameter learning process |
US20100121573A1 (en) * | 2008-11-11 | 2010-05-13 | Seiko Epson Corporation | Position calculating method and position calculating device |
CN103034249A (zh) * | 2012-12-27 | 2013-04-10 | 中国电子科技集团公司第二十六研究所 | 一种基于模糊决策的双轴光电转台复合控制***和控制方法 |
CN103684183A (zh) * | 2013-11-14 | 2014-03-26 | 西安理工大学 | 异步电机转速估计方法 |
-
2014
- 2014-05-13 CN CN201410200184.6A patent/CN103956956B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5940814A (en) * | 1997-05-01 | 1999-08-17 | Nec Research Institute, Inc. | Defuzzification strategy and parameter learning process |
US20100121573A1 (en) * | 2008-11-11 | 2010-05-13 | Seiko Epson Corporation | Position calculating method and position calculating device |
CN103034249A (zh) * | 2012-12-27 | 2013-04-10 | 中国电子科技集团公司第二十六研究所 | 一种基于模糊决策的双轴光电转台复合控制***和控制方法 |
CN103684183A (zh) * | 2013-11-14 | 2014-03-26 | 西安理工大学 | 异步电机转速估计方法 |
Non-Patent Citations (3)
Title |
---|
D.LOEBIS,R.SUTTON,J.CHUDLEY,W.NAEEM: "《Adaptive tuning of a Kalman filter via fuzzy logic for an intelligent AUV navigation system》", 《CONTROL ENGINEERING PRACTICE 12(2004)》 * |
MOULAY RACHID DOUIRI,MOHAMED CHERKAOUI,SIDI MOHANMED DOURI: "《Rotor resistance and speed identification using extended Kalman filter and fuzzy logic controller for induction machine drive》", 《INTERNATIONAL CONFERENCE ON MULTIMEDIA COMPUTING & SYSTEM》 * |
周健,王辉: "《基于扩展卡尔曼滤波器的无刷直流电动机控制》", 《微特电机》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104993765A (zh) * | 2015-08-04 | 2015-10-21 | 重庆大学 | 一种无刷直流电机的转速估计方法 |
CN106019045A (zh) * | 2016-05-16 | 2016-10-12 | 安徽大学 | 一种永磁同步电机匝间短路故障诊断的方法 |
CN106533297A (zh) * | 2016-12-22 | 2017-03-22 | 东南大学 | 一种基于卡尔曼滤波的定子永磁型记忆电机调磁方法 |
CN108490472B (zh) * | 2018-01-29 | 2021-12-03 | 哈尔滨工程大学 | 一种基于模糊自适应滤波的无人艇组合导航方法 |
CN108490472A (zh) * | 2018-01-29 | 2018-09-04 | 哈尔滨工程大学 | 一种基于模糊自适应滤波的无人艇组合导航方法 |
CN108649848A (zh) * | 2018-05-18 | 2018-10-12 | 湖南大学 | 基于扩展卡尔曼滤波方法的磁齿轮电机失步预测控制方法 |
CN108649848B (zh) * | 2018-05-18 | 2021-03-30 | 湖南大学 | 基于扩展卡尔曼滤波方法的磁齿轮电机失步预测控制方法 |
CN110556927A (zh) * | 2019-06-19 | 2019-12-10 | 国家电网有限公司 | 一种感应式无线电能传输***的状态观测的方法及*** |
CN111308522A (zh) * | 2020-03-27 | 2020-06-19 | 上海天好信息技术股份有限公司 | 一种基于新型Kalman滤波的多维时空数据估计方法 |
CN112858929A (zh) * | 2021-03-16 | 2021-05-28 | 上海理工大学 | 一种基于模糊逻辑与扩展卡尔曼滤波的电池soc估计方法 |
CN113381665A (zh) * | 2021-06-09 | 2021-09-10 | 上海晟矽微电子股份有限公司 | 电机控制方法及装置、电机和电动工具 |
CN114257150A (zh) * | 2022-01-26 | 2022-03-29 | 合肥倍豪海洋装备技术有限公司 | 一种基于改进卡尔曼滤波器的永磁同步电机死区补偿方法 |
CN114744938A (zh) * | 2022-04-11 | 2022-07-12 | 华中科技大学 | 一种基于卡尔曼滤波的全参数观测器及全参数辨识方法 |
CN114744938B (zh) * | 2022-04-11 | 2024-06-04 | 华中科技大学 | 一种基于卡尔曼滤波的全参数观测器及全参数辨识方法 |
CN116137505A (zh) * | 2023-04-18 | 2023-05-19 | 深圳市浮思特科技有限公司 | 一种无刷直流电机的卡尔曼滤波矢量控制***及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103956956B (zh) | 2017-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103956956A (zh) | 一种基于扩展卡尔曼滤波器的无刷直流电机状态估计方法 | |
CN103997272B (zh) | 永磁同步电机的负载扰动补偿装置及方法 | |
CN109194219B (zh) | 一种基于无模型非奇异终端滑模控制永磁同步电机的方法及*** | |
CN104682805B (zh) | 基于扩张状态观测器的永磁同步电机全阶滑模变结构位置伺服控制方法 | |
CN103956953B (zh) | 基于滑模观测器的无刷直流电机状态估计方法 | |
DE102014106668A1 (de) | Optimierte steuerung für synchronmotoren | |
EP2544357A2 (en) | Control of electrical machines | |
CN106533295A (zh) | 永磁同步电机无位置传感器控制方法和装置 | |
CN105577058A (zh) | 基于新型模糊自抗扰控制器的五相容错永磁电机速度控制方法 | |
CN102969968A (zh) | 一种永磁同步电机控制方法 | |
CN106911282B (zh) | 一种改进模糊控制的永磁电机无速度测速*** | |
CN105337546B (zh) | 基于变阶次分数阶滑模的永磁同步电机控制装置及方法 | |
CN103051274A (zh) | 基于变阻尼的二自由度永磁同步电机的无源性控制方法 | |
CN106911281B (zh) | 一种基于模糊控制和mras的永磁同步电机无速度传感器测速方法 | |
CN106059423A (zh) | 一种基于fc和smo的无速度传感器控制*** | |
CN102647144B (zh) | 用于估计同步磁阻电动机的转子角度的方法和设备 | |
CN106788051A (zh) | 一种无轴承异步电机转速估计方法 | |
CN111193448A (zh) | 基于扩展卡尔曼滤波器的表贴式永磁同步电机负载转矩观测方法 | |
Krzeminski | Observer of induction motor speed based on exact disturbance model | |
CN109194224B (zh) | 基于扩张状态观测器的永磁同步电机无传感器控制方法 | |
CN107395080B (zh) | 基于级联非奇异终端滑模观测器的无速度传感器转矩控制***及方法 | |
CN108448983B (zh) | 一种抑制伺服***极低速时非线性干扰的方法 | |
CN108574440A (zh) | 一种基于滑模参考自适应的永磁同步电机状态估计方法 | |
CN109687792A (zh) | 面向矢量控制***的牵引电机转子参数在线辨识优化方法 | |
CN112436774A (zh) | 一种无速度传感器驱动的异步电机控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |