CN109100532B - 基于自适应交互双模算法的测速发电机的滤波测速方法 - Google Patents
基于自适应交互双模算法的测速发电机的滤波测速方法 Download PDFInfo
- Publication number
- CN109100532B CN109100532B CN201810640792.7A CN201810640792A CN109100532B CN 109100532 B CN109100532 B CN 109100532B CN 201810640792 A CN201810640792 A CN 201810640792A CN 109100532 B CN109100532 B CN 109100532B
- Authority
- CN
- China
- Prior art keywords
- model
- algorithm
- filtering
- speed measurement
- adaptive
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P3/00—Measuring linear or angular speed; Measuring differences of linear or angular speeds
- G01P3/42—Devices characterised by the use of electric or magnetic means
- G01P3/44—Devices characterised by the use of electric or magnetic means for measuring angular speed
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Control Of Electric Motors In General (AREA)
Abstract
发明提供了一种基于自适应交互双模算法的测速发电机的滤波测速方法。所述基于自适应交互双模算法的测速发电机的滤波测速方法包括如下步骤:步骤1:针对测速发电机不对称的特性,建立类爪极永磁式交流测速电机测速模型;步骤2:根据步骤1建立的测速发电机数学模型,提出滤波测速算法,并应用于该测速发电机数学模型中,从而精确计算得出测速发电机的转速。本发明的有益效果是:所述基于自适应交互双模算法的测速发电机的滤波测速方法提出的宽范围滤波测速算法对于结构不对称扰动具有良好的鲁棒性,能够适用于转速高机动调节,实时性和精度高,调速范围宽。
Description
技术领域
本发明属于测速发电机技术领域,具体地涉及一种基于自适应交互双模算法的测速发电机的滤波测速方法。
背景技术
近年来,各类伺服控制系统广泛应用于智能制造、工业机器人、精密数控装备等行业领域。在高性能伺服系统中,电机速度反馈的实时性和精度是影响伺服系统控制性能的决定性因素。在当前的工程应用中,受限于测速传感器成本和分辨率以及无传感器技术依赖于复杂、高阶的电机模型和负载特性,环境扰动、复杂工况等因素的影响,宽范围、高精度电机转速测量问题仍是学术和工程界具有挑战的研究方向,受到了广泛的关注。
目前,对于宽范围、高精度电机转速测量问题的解决方案有两种:基于对各类测速传感器的测速性能进行优化、及基于各型电机模型参数的无传感器测速技术。
尽管基于各型电机模型参数的无传感器测速技术虽然取得了长足的进步,但由于都需要依赖于各型电机精确的模型参数,受限于系统环境,与基于传感器的测速方法相比精度较低,抗干扰能力较差,主要应用于对空间、密封性、环境工况等有苛刻要求的特定场景。
当前,综合考虑性价比、可靠性和测速性能,基于传感器的测速方法仍是普通环境下测速应用中的首选方法。由于测速传感器的精度(线数)与价格成正比,高精度、高性能测速传感器高昂的价格让普通用户望而却步。因此针对常规精度测速传感器,设计在保证高、中转速区测速精度,并能有效改善低速区测速精度的宽范围、高可靠、实时滤波测速算法对提高普通测速传感器的应用范围具有重要的工程意义。
发明内容
本发明的目的在于针对现有技术的缺陷或问题,提供一种基于自适应交互双模算法的测速发电机的滤波测速方法。
本发明的技术方案如下:一种基于自适应交互双模算法的测速发电机的滤波测速方法包括如下步骤:步骤1:针对测速发电机不对称的特性,建立类爪极永磁式交流测速电机测速模型;步骤2:根据步骤1建立的测速发电机数学模型,提出滤波测速算法,并应用于该测速发电机数学模型中,从而精确计算得出测速发电机的转速。
优选地,在步骤一中,具体包括如下步骤:
步骤1.1:由于结构不对称参数Δθ的存在使得测速发电机的导磁爪距呈不均匀分布特性,即出现机械弧度为的长导磁爪距和机械弧度为的短导磁爪距;假定永磁转子以电角速度ω相对定子绕组顺时针方向旋转,在k时刻,永磁转子的N极或S极旋转至和段的电角速度和经历的电角度分别为:
步骤1.2:假定永磁转子以磁爪D的中心作为永磁转子旋转的起始参考点,则永磁转子定子绕组基波电动势的离散测速方程为:
优选地,在步骤二中,考虑到测速发电机永磁转子的转动状态一般为恒速或变速运动的组合,则将测速发电机永磁转子(被测电机)的运动状态对应为两个运动模型:恒速模型M1和变速模型M2,设Mj,j=1,2的离散状态方程为:
优选地,在步骤二中,所述滤波测速算法为无迹卡尔曼滤波算法,并包括如下步骤
步骤2.1:根据UT变换,基于τ阶模型Mj计算由2τ+1个sigma点组成的采样点集χk|k:
式中:λ=α2(τ+κ)-τ,是一个比例因子,α和κ都是正常数,α决定x估值m周围sigma点的分布状态,通常选择0≤α≤1;适当调节α、κ可以提高估计均值的精度,调节β可以提高方差的精度。表示矩阵的第i列,i∈[1,2τ];
步骤2.2:将2τ+1个sigma点[χk|k]i(χk|k的第i列),i∈[1,2τ]依次带入Mj,j=1,2的离散状态方程计算采样点集χk|k的预测值χk+1|k=f(χk|k);
式中:
分别为均值和协方差加权值;
步骤2.4:根据UT变换步骤2.3中的预测公式,则预测估计采样点集为:
步骤2.5:将2τ+1个sigma点预测向量[χk+1|k]i(χk+1|k的第i列),i∈[1,2τ]依次带入步骤1.2中的离散测速方程预测观测值ek+1|k,进而计算测量均值协方差状态协方差卡尔曼增益矩阵Kk+1如下:
步骤2.6:最后,更新模型Mj的状态和协方差:
优选地,在步骤二中,基于所述恒速模型M1和所述变速模型M2的自适应交互双模算法包括如下步骤:
一、模型交互
二、滤波
三、更新模型概率
根据贝叶斯(Bayes)公式,模型Mj在k时刻的概率更新为:
其中ρ为总采样次数;
四、联合估计
根据模型概率,加权合并计算全部模型Mj的滤波值,可得联合估计和协方差:
五、自适应估计转移概率矩阵
在线自适应更新转移概率矩阵;
六、自适应降阶
定义滤波误差的均值和标准差为:
式中:M为采样次数;
优选地,在步骤六中,自适应降阶运行策略包括如下步骤:
步骤6.1:设定自适应交互双模算法中各变量的初值以及ε,进入步骤6.2。
步骤6.3:并行运行降阶和全阶自适应交互双模算法,在采用降阶估计值继续保障系统实时性的同时,在计算耗时较长的全阶算法运行完成一次后修正一次降阶估计值,提高测速精度;若返回步骤6.2,若重复步骤6.3。
本发明提供的技术方案具有如下有益效果:
所述基于自适应交互双模算法的测速发电机的滤波测速方法在结构不对称条件下建立了一类爪极永磁式交流测速电机测速模型,提出了一种面向工程应用、计算量小的交互双模自适应降阶无迹卡尔曼滤波算法实时估计永磁转子转速;而且,所述基于自适应交互双模算法的测速发电机的滤波测速方法提出的宽范围滤波测速算法对于结构不对称扰动具有良好的鲁棒性,能够适用于转速高机动调节,实时性和精度高,调速范围宽。
附图说明
图1是爪极永磁式交流测速电机的工作原理示意图;
图2是爪极永磁式交流测速电机结构不对称条件下转速模型原理图;
图3是爪极永磁式交流测速电机结构不对称条件下匀转速电动势波形图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
除非上下文另有特定清楚的描述,本发明中的元件和组件,数量既可以单个的形式存在,也可以多个的形式存在,本发明并不对此进行限定。本发明中的步骤虽然用标号进行了排列,但并不用于限定步骤的先后次序,除非明确说明了步骤的次序或者某步骤的执行需要其他步骤作为基础,否则步骤的相对次序是可以调整的。可以理解,本文中所使用的术语“和/或”涉及且涵盖相关联的所列项目中的一者或一者以上的任何和所有可能的组合。
接下来以爪极永磁式交流测速电机为例对本发明进行阐述。当然,应当理解,本发明实施例提供的基于自适应交互双模算法的测速发电机的滤波测速方法不仅适用于爪极永磁式交流测速电机的测速,也可以适用于其他合适的测速电机,本发明对此不做限定。
而且,爪极永磁式交流测速电机的工作原理如图1所示。永磁转子轴向充磁,磁极的极性用N和S表示。当前瞬时,永磁转子N极、导磁爪A、上爪极(白色)、下爪极(黑色)、导磁爪D、永磁转子S极构成磁回路如图1(a)所示,磁场方向垂直纸面,生成的感应电动势按照正弦规律变化如图1(b)所示。
用电动机驱动永磁转子以恒定转速n(单位:r/min)相对于定子顺时针方向旋转,与定子绕组之间产生相对运动。根据电磁感应定律,定子绕组中会产生感应电动势,其瞬时值:
e=eN-eS=Em sinωt (1)
其中:eN和eS分别为N极和S极的电动势,它们大小相等,方向相反,即在时间相位上彼此相差π弧度电角度;ω为空间电角速度(单位:rad/s),Em为e的最大值(单位:V),可以表示为:
若已知ω和导磁爪对数p,则永磁转子顺时针方向旋转的转速n可以表示为:
如图2所示,考虑到工程中应用的爪极永磁式交流测速电机常在结构不对称条件下运行,感应电动势e将发生畸变,不再是标准的正弦特性,转速n将与有效值E也不再成简单的正比关系,需要重新建立新的测速模型。
例如,当上爪极在安装时与下极爪偏差电角度θ,则导磁爪D与永磁转子S极对齐时,导磁爪A与对称位置A1偏差电机械角度Δθ,会导致N极和S极分别产生的电动势eN和eS不对称;若导磁爪A的长度与标准长度相差δ,则气隙磁通量将随之变化。结构不对称因子Δθ、δ以及转速波动及干扰噪声等随机因素的存在是使得定子绕组的感应电动势e产生不对称畸变问题的重要因素,如图3所示。例如A、B、C半波的周期和幅值均不一致,在C半波的D处还产生了毛刺(导磁爪D结构瑕疵引起)。建立结构不对称条件下爪极永磁式交流测速电机数学模型,并在其基础上采用滤波方法估计对称电角速度ω,将有助于提高永磁转子(原动机)的转速n的测量精度和鲁棒性。
如此,则本发明实施例提供的基于自适应交互双模算法的测速发电机的滤波测速方法包括如下步骤:
步骤1:针对测速发电机不对称的特性,建立类爪极永磁式交流测速电机测速模型;
步骤2:根据步骤1建立的测速发电机数学模型,提出滤波测速算法,并应用于该测速发电机数学模型中,从而精确计算得出测速发电机的转速。
具体地,在步骤1中,包括如下步骤:
步骤1.2:假定永磁转子以磁爪D的中心作为永磁转子旋转的起始参考点,则永磁转子定子绕组基波电动势的离散测速方程为:
其中:υk为零均值白噪声序列,协方差矩阵为R;而且,在利用所述离散测速方程估计目标转速时,为了避免电角度数据累加溢出,依据或的周期性,将Θk的值域限定在一个周期内,即在算法运行过程中,当Θk从零递增至2πp就将其清零,使得0≤Θk≤2πp。
在步骤二中,考虑到测速发电机永磁转子的转动状态一般为恒速或变速运动的组合,则将测速发电机永磁转子(被测电机)的运动状态对应为两个运动模型:恒速模型M1和变速模型M2,设Mj,j=1,2的离散状态方程为:
实际上,步骤1.2中永磁转子定子绕组基波电动势的离散测速方程为非线性的,在本实施例中,在步骤二中的滤波测速算法为交互双模自适应降阶无迹卡尔曼滤波算法,从而能够实时、宽范围估计被测转速,并提高测速的精度和鲁棒性。
具体地,所述滤波测速算法为无迹卡尔曼滤波算法,并包括如下步骤:
步骤2.1:根据UT变换,基于τ阶模型Mj计算由2τ+1个sigma点组成的采样点集χk|k:
式中:λ=α2(τ+κ)-τ,是一个比例因子,α和κ都是正常数,α决定x估值m周围sigma点的分布状态,通常选择0≤α≤1;适当调节α、κ可以提高估计均值的精度,调节β可以提高方差的精度。表示矩阵的第i列,i∈[1,2τ];
步骤2.2:将2τ+1个sigma点[χk|k]i(χk|k的第i列),i∈[1,2τ]依次带入Mj,j=1,2的离散状态方程计算采样点集χk|k的预测值χk+1|k=f(χk|k);
式中:
分别为均值和协方差加权值;
步骤2.4:根据UT变换步骤2.3中的预测公式,则预测估计采样点集为:
步骤2.5:将2τ+1个sigma点预测向量[χk+1|k]i(χk+1|k的第i列),i∈[1,2τ]依次带入步骤1.2中的离散测速方程预测观测值ek+1|k,进而计算测量均值协方差状态协方差卡尔曼增益矩阵Kk+1如下:
步骤2.6:最后,更新模型Mj的状态和协方差:
而且,在步骤二中,基于所述恒速模型M1和所述变速模型M2的自适应交互双模算法包括如下步骤:
一、模型交互
二、滤波
三、更新模型概率
根据贝叶斯(Bayes)公式,模型Mj在k时刻的概率更新为:
其中ρ为总采样次数;
四、联合估计
根据模型概率,加权合并计算全部模型Mj的滤波值,可得联合估计和协方差:
五、自适应估计转移概率矩阵
在线自适应更新转移概率矩阵;
六、自适应降阶
定义滤波误差的均值和标准差为:
式中:M为采样次数;
此外,在步骤六中,自适应降阶运行策略包括如下步骤:
步骤6.1:设定自适应交互双模算法中各变量的初值以及ε,进入步骤6.2。
步骤6.3:并行运行降阶和全阶自适应交互双模算法,在采用降阶估计值继续保障系统实时性的同时,在计算耗时较长的全阶算法运行完成一次后修正一次降阶估计值,提高测速精度;若返回步骤6.2,若重复步骤6.3。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
Claims (2)
1.一种基于自适应交互双模算法的测速发电机的滤波测速方法,其特征在于:包括如下步骤:
步骤1:针对测速发电机不对称的特性,建立类爪极永磁式交流测速电机测速模型;
步骤2:根据步骤1建立的测速发电机数学模型,提出滤波测速算法,并应用于该测速发电机数学模型中,从而精确计算得出测速发电机的转速;其中:
在步骤一中,具体包括如下步骤:
步骤1.1:由于结构不对称参数Δθ的存在使得测速发电机的导磁爪距呈不均匀分布特性,即出现机械弧度为的长导磁爪距和机械弧度为的短导磁爪距;假定永磁转子以电角速度ω相对定子绕组顺时针方向旋转,在k时刻,永磁转子的N极或S极旋转至和段的电角速度和经历的电角度分别为:
步骤1.2:假定永磁转子以磁爪D的中心作为永磁转子旋转的起始参考点,则永磁转子定子绕组基波电动势的离散测速方程为:
在步骤二中,考虑到测速发电机永磁转子的转动状态一般为恒速或变速运动的组合,则将测速发电机永磁转子的运动状态对应为两个运动模型:恒速模型M1和变速模型M2,设Mj,j=1,2的离散状态方程为:
在步骤二中,所述滤波测速算法为无迹卡尔曼滤波算法,并包括如下步骤
步骤2.1:根据UT变换,基于τ阶模型Mj计算由2τ+1个sigma点组成的采样点集χk|k:
式中:λ=α2(τ+κ)-τ,是一个比例因子,α和κ都是正常数,α决定x估值m周围sigma点的分布状态,通常选择0≤α≤1;适当调节α、κ可以提高估计均值的精度,调节β可以提高方差的精度;表示矩阵的第i列,i∈[1,2τ];
步骤2.2:将2τ+1个sigma点[χk|k]i(χk|k的第i列),i∈[1,2τ]依次带入Mj,j=1,2的离散状态方程计算采样点集χk|k的预测值χk+1|k=f(χk|k);
式中:
分别为均值和协方差加权值;
步骤2.4:根据UT变换步骤2.3中的预测公式,则预测估计采样点集为:
步骤2.5:将2τ+1个sigma点预测向量[χk+1|k]i(χk+1|k的第i列),i∈[1,2τ]依次带入步骤1.2中的离散测速方程预测观测值ek+1|k,进而计算测量均值协方差状态协方差卡尔曼增益矩阵Kk+1如下:
步骤2.6:最后,更新模型Mj的状态和协方差:
在步骤二中,基于所述恒速模型M1和所述变速模型M2的自适应交互双模算法包括如下步骤:
一、模型交互
二、滤波
三、更新模型概率
根据贝叶斯(Bayes)公式,模型Mj在k时刻的概率更新为:
其中ρ为总采样次数;
四、联合估计
根据模型概率,加权合并计算全部模型Mj的滤波值,可得联合估计和协方差:
五、自适应估计转移概率矩阵
在线自适应更新转移概率矩阵;
六、自适应降阶
定义滤波误差的均值和标准差为:
式中:M为采样次数;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810640792.7A CN109100532B (zh) | 2018-06-21 | 2018-06-21 | 基于自适应交互双模算法的测速发电机的滤波测速方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810640792.7A CN109100532B (zh) | 2018-06-21 | 2018-06-21 | 基于自适应交互双模算法的测速发电机的滤波测速方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109100532A CN109100532A (zh) | 2018-12-28 |
CN109100532B true CN109100532B (zh) | 2020-08-11 |
Family
ID=64844958
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810640792.7A Active CN109100532B (zh) | 2018-06-21 | 2018-06-21 | 基于自适应交互双模算法的测速发电机的滤波测速方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109100532B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109962647B (zh) * | 2019-03-21 | 2020-11-10 | 上海交通大学 | 带有双模耦合结构观测器的电机控制系统及状态估计方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE2314952A1 (de) * | 1973-03-26 | 1974-10-31 | Abakus Gmbh | Drehzahlregelung fuer gleichspannungsscheibenlaeufer-motoren |
WO1985002591A1 (en) * | 1983-12-16 | 1985-06-20 | Robert Bosch Gmbh | Process for transforming the output signal of a rotational speed measuring sensor into an interference-poor signal |
CN101005263A (zh) * | 2007-01-25 | 2007-07-25 | 上海交通大学 | 交流电机伺服系统速度控制方法 |
CN101741302A (zh) * | 2008-11-21 | 2010-06-16 | 上海电机学院 | 一种用于永磁无刷直流电机的转子位置检测方法 |
CN102176653A (zh) * | 2011-01-19 | 2011-09-07 | 哈尔滨工业大学 | 带指数渐消因子的卡尔曼滤波器的感应电机转速观测方法 |
CN104483502B (zh) * | 2014-12-23 | 2017-09-29 | 三明恒顺印刷机械有限公司 | 一种基于单片机的电机实时转速宽范围精确测速方法 |
CN104914269B (zh) * | 2015-05-20 | 2017-12-08 | 三明学院 | 测速发电机不对称运行下的宽范围实时测速方法 |
CN107037235B (zh) * | 2016-11-28 | 2023-08-01 | 东南大学 | 一种制动下滑量软测量方法及测量装置 |
-
2018
- 2018-06-21 CN CN201810640792.7A patent/CN109100532B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109100532A (zh) | 2018-12-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107070341B (zh) | 基于鲁棒迭代学习控制的永磁同步电机转矩脉动抑制方法 | |
Xu et al. | Very-low speed control of PMSM based on EKF estimation with closed loop optimized parameters | |
KR101087581B1 (ko) | 영구자석형 동기 전동기의 센서리스 제어방법 | |
CN110289795A (zh) | 一种电动汽车用永磁同步电机控制系统及控制方法 | |
CN107065551A (zh) | 一种基于模型参数精确辨识的仿真转台自校正控制方法 | |
CN108258967A (zh) | 一种基于新型磁链观测器的永磁电机无位置直接转矩控制方法 | |
CN110739893A (zh) | 一种改进自适应无轨迹卡尔曼滤波转动惯量辨识方法 | |
CN111371362B (zh) | 高频注入法永磁直线电机动子位置估计的补偿方法 | |
CN111504360A (zh) | 基于时间坐标磁电编码器角度精分方法及装置 | |
CN109510541A (zh) | 基于分段式永磁同步电机滑模伺服控制的方法 | |
CN107395083A (zh) | Pmlsm伺服控制系统位置与速度估计方法及装置 | |
CN110995102A (zh) | 一种永磁同步电机直接转矩控制方法及系统 | |
EP2493067B1 (en) | Method and apparatus for estimating rotor angle of synchronous reluctance motor | |
Zhao et al. | Back EMF-based dynamic position estimation in the whole speed range for precision sensorless control of PMLSM | |
CN109194208A (zh) | 无速度传感器控制方法 | |
CN111193448A (zh) | 基于扩展卡尔曼滤波器的表贴式永磁同步电机负载转矩观测方法 | |
CN111817633A (zh) | 一种永磁同步电机的机械参数辨识方法 | |
CN107681937B (zh) | 一种基于神经网络的超高速永磁同步电机速度观测方法 | |
CN109100532B (zh) | 基于自适应交互双模算法的测速发电机的滤波测速方法 | |
CN106602952A (zh) | 一种pmsm永磁体磁链满秩辨识方法 | |
CN110768600A (zh) | 一种pmsm无速度传感器转子检测方法 | |
CN112149274A (zh) | 一种带有死区输入非线性的多轴雕刻机系统在线建模方法 | |
CN105871277A (zh) | 一种基于最小方差的永磁伺服系统非线性模型预测控制器设计方法 | |
CN110995098B (zh) | 一种永磁同步电机的反演控制方法 | |
Lesan et al. | Comparative study of speed estimation techniques for sensorless vector control of induction machine |
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 | ||
CB03 | Change of inventor or designer information |
Inventor after: Wen Xia Inventor after: Ren Wen Inventor after: Lai Sencai Inventor before: Ren Wen Inventor before: Lai Sencai |
|
CB03 | Change of inventor or designer information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |