CN106452247A - 一种永磁同步电机转动惯量辨识的方法及装置 - Google Patents
一种永磁同步电机转动惯量辨识的方法及装置 Download PDFInfo
- Publication number
- CN106452247A CN106452247A CN201611139450.4A CN201611139450A CN106452247A CN 106452247 A CN106452247 A CN 106452247A CN 201611139450 A CN201611139450 A CN 201611139450A CN 106452247 A CN106452247 A CN 106452247A
- Authority
- CN
- China
- Prior art keywords
- formula
- synchronous motor
- permagnetic synchronous
- value
- moment
- 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
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/0003—Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
- H02P21/0017—Model reference adaptation, e.g. MRAS or MRAC, useful for control or parameter estimation
-
- 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/143—Inertia or moment of inertia estimation
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Ac Motors In General (AREA)
- Control Of Electric Motors In General (AREA)
Abstract
本发明实施例公开了一种永磁同步电机转动惯量辨识的方法及装置,用于解决目前转动惯量离线辨识是在系统模型已知的情况下,获得所有的数据,然后匹配相应的算法,对数据进行估计处理,得到辨识目标参数的估计值,现有的离线辨识方法,如加减速法、转矩限幅加速法、积分辨识法等等,它们原理简单,算法简单,操作起来方便易行,但是其需要进行大量的数据存储和运算,辨识收敛时间长,离线辨识方法具有对计算机要求高,辨识精度不高,不具有实时性的技术问题。
Description
技术领域
本发明涉及永磁同步电机领域,尤其涉及一种永磁同步电机转动惯量辨识的方法及装置。
背景技术
随着现代工业自动化技术的飞速发展,永磁同步电机因其体积小、效率高、过载能力大并且可靠性高和响应速度快方面的优点,以及高启动力矩、高功率密度、转子转动惯量小等特性,在伺服领域中占据重要位置。同时,永磁同步电机也广泛应用于家电设备,交通运输中。伴随着永磁同步电机的伺服系统范围的不断扩大,必然也会使人们对其性能方面的要求不断提高,尤其是跟踪精度和鲁棒性方面。针对这种情况下,国内外大量学者先后提出了许多先进的控制策略,但这其中大多数控制策略都是基于永磁同步电机精确参数的基础之上,这就奠定了精确电机参数在整个控制系统设计中的重要基础位置。
一般来说,永磁同步电机的参数是可以直接从电机的名牌上获得,但是这些参数是电机在特定运行条件下(多为固定情况)测试得到的,而永磁同步电机是一个强耦合、非线性系统,电机的转子永磁体会因为电机老化而发生磁性退化的现象,当电机的相电流增大到额定电流以上的时候,电机定子磁芯会饱和,因此,要获得实时的准确的电机参数变得异常困难,通过准确辨识电机的参数,可以避免参数误差对系统性能的影响,实现调节器参数的重新整定,提高系统的鲁棒性。
转动惯量的变化对伺服系统的机械性能影响较大,根据电机的运动方程,我们可以知道转动惯量随着负载变化,惯量过大,使系统响应变慢,导致速度不稳,精度下降;惯量过小,会加快系统的响应速度,容易使速度产生超调和震荡,直接进行测量永磁同步电机的转动惯量并不容易。
发明内容
本发明实施例提供了一种永磁同步电机转动惯量辨识的方法及装置,解决了目前转动惯量离线辨识是在系统模型已知的情况下,获得所有的数据,然后匹配相应的算法,对数据进行估计处理,得到辨识目标参数的估计值,现有的离线辨识方法,如加减速法、转矩限幅加速法、积分辨识法等等,它们原理简单,算法简单,操作起来方便易行,但是其需要进行大量的数据存储和运算,辨识收敛时间长,离线辨识方法具有对计算机要求高,辨识精度不高,不具有实时性的技术问题。
本发明实施例提供的一种永磁同步电机转动惯量辨识的方法,其特征在于,包括:
S1:通过递推最小二乘法计算辨识参数最小二乘的递推公式;
S2:通过预置的性能指标,引入遗传因子对所述辨识参数最小二乘的递推公式进行数据饱和的改进,计算遗忘因子最小二乘参数估计公式;
S3:通过电机转矩平衡方程,计算出最小二乘形式,并通过所述遗忘因子最小二乘参数估计公式,计算出k时刻的永磁同步电机的转动惯量J的值J(k);
S4:接收预置的满足条件,更新k时刻的永磁同步电机的转动惯量J的值J(k),直至k时刻的永磁同步电机的转动惯量J的值J(k)满足预置的满足条件,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k)。
优选地,
所述步骤S4具体包括:
接收预置的遗传因子β的初始值0和预置的当前时刻转动惯量的辨识结果和当前时刻的前一时刻的辨识结果差的绝对值E=|J(k)-J(k-1)小于等于0.000001的满足条件,更新k时刻的永磁同步电机的转动惯量J的值J(k),若E≤0.000001,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k),若E≥0.000001,遗传因子β自动加0.01,继续更新k时刻的永磁同步电机的转动惯量J的值J(k),直至E≤0.000001,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k)。
优选地,
所述步骤S1具体包括:
接收预置的非奇异方阵A,(A+BC),(I+CA-1B),(A+BC)的逆A-1+X,计算出X=-A-1B(I+CA-1B)-1CA-1,(A+BC)-1=A-1-A-1B(I+CA-1B)-1CA-1,根据 通过递推最小二乘法,计算出
接收预置的A=P-1(k-1),计算出辨识参数最小二乘的递推公式:
优选地,
所述步骤S2具体包括:
通过预置的性能指标式中β为遗忘因子,0≤β≤1,对所述辨识参数最小二乘的递推公式进行数据饱和的改进,计算遗忘因子最小二乘参数估计公式为:
优选地,
所述步骤S3具体包括:
通过永磁同步电机转矩平衡方程:
其中Tl是永磁同步电机的负载转矩,J为转动惯量,B为运动阻尼系数,ωm为转子机械速度,计算出
再进行离散化处理:
计算出最小二乘形式:
接收预置的
y(k)=ωm(k)-2ωm(k-1)+ωm(k-2)
即最小二乘形式为通过遗忘因子最小二乘参数估计式,计算出k时刻的永磁同步电机的转动惯量J的值J(k)。
本发明实施例中提供的一种永磁同步电机转动惯量辨识的装置,包括:
第一计算单元,用于通过递推最小二乘法计算辨识参数最小二乘的递推公式;
第二计算单元,用于通过预置的性能指标,引入遗传因子对所述辨识参数最小二乘的递推公式进行数据饱和的改进,计算遗忘因子最小二乘参数估计公式;
第三计算单元,用于通过电机转矩平衡方程,计算出最小二乘形式,并通过所述遗忘因子最小二乘参数估计公式,计算出k时刻的永磁同步电机的转动惯量J的值J(k);
更新单元,用于接收预置的满足条件,更新k时刻的永磁同步电机的转动惯量J的值J(k),直至k时刻的永磁同步电机的转动惯量J的值J(k)满足预置的满足条件,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k)。
优选地,
所述更新单元具体包括:
接收子单元,用于接收预置的遗传因子β的初始值0和预置的当前时刻转动惯量的辨识结果和当前时刻的前一时刻的辨识结果差的绝对值E=|J(k)-J(k-1)小于等于0.000001的满足条件;
更新子单元,用于更新k时刻的永磁同步电机的转动惯量J的值J(k);
输出子单元,用于若E≤0.000001,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k);
循环子单元,用于若E≥0.000001,遗传因子β自动加0.01;继续更新k时刻的永磁同步电机的转动惯量J的值J(k),直至E≤0.000001,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k)。
优选地,
所述第一计算单元具体包括:
第一计算子单元,用于接收预置的非奇异方阵A,(A+BC),(I+CA-1B),(A+BC)的逆A-1+X,计算出X=-A-1B(I+CA-1B)-1CA-1,(A+BC)-1=A-1-A-1B(I+CA-1B)-1CA-1;
第二计算子单元,用于根据
通过递推最小二乘法,计算出
第三计算子单元,用于接收预置的A=P-1(k-1), 计算出辨识参数最小二乘的递推公式:
优选地,
所述第二计算单元具体包括:
改进子单元,用于通过预置的性能指标式中β为遗忘因子,0≤β≤1,对所述辨识参数最小二乘的递推公式进行数据饱和的改进;
第四计算子单元,用于计算遗忘因子最小二乘参数估计公式为:
优选地,
所述第三计算单元具体包括:
第五计算子单元,用于通过永磁同步电机转矩平衡方程:
其中Tl是永磁同步电机的负载转矩,J为转动惯量,B为运动阻尼系数,ωm为转子机械速度,计算出
处理子单元,用于对
进行离散化处理:
计算出最小二乘形式:
第六计算子单元,用于接收预置的
y(k)=ωm(k)-2ωm(k-1)+ωm(k-2)
即最小二乘形式为通过遗忘因子最小二乘参数估计式,计算出k时刻的永磁同步电机的转动惯量J的值J(k)。
从以上技术方案可以看出,本发明实施例具有以下优点:
本发明实施例中提供的一种永磁同步电机转动惯量辨识的方法及装置,其中一种永磁同步电机转动惯量辨识的方法,包括:S1:通过递推最小二乘法计算辨识参数最小二乘的递推公式;S2:通过预置的性能指标,引入遗传因子对所述辨识参数最小二乘的递推公式进行数据饱和的改进,计算遗忘因子最小二乘参数估计公式;S3:通过电机转矩平衡方程,计算出最小二乘形式,并通过所述遗忘因子最小二乘参数估计公式,计算出k时刻的永磁同步电机的转动惯量J的值J(k);S4:接收预置的满足条件,更新k时刻的永磁同步电机的转动惯量J的值J(k),直至k时刻的永磁同步电机的转动惯量J的值J(k)满足预置的满足条件,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k)。本实施例中,通过在递推最小二乘法的基础上加入遗忘因子,并对遗忘因子的变化,设计了循环环节,减少输出震荡,保证了输出的平稳性,具有实时辨识,辨识收敛时间短的优点,解决了目前转动惯量离线辨识是在系统模型已知的情况下,获得所有的数据,然后匹配相应的算法,对数据进行估计处理,得到辨识目标参数的估计值,现有的离线辨识方法,如加减速法、转矩限幅加速法、积分辨识法等等,它们原理简单,算法简单,操作起来方便易行,但是其需要进行大量的数据存储和运算,辨识收敛时间长,离线辨识方法具有对计算机要求高,辨识精度不高,不具有实时性的技术问题。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其它的附图。
图1本发明实施例提供的一种磁同步电机转动惯量辨识的方法的一个实施例的流程示意图;
图2本发明实施例提供的一种循环环节的流程示意图;
图3本发明实施例提供的一种磁同步电机转动惯量辨识的装置的一个实施例的结构示意图;
图4本发明实施例提供的一种磁同步电机转动惯量辨识的装置的另一个实施例的结构示意图。
具体实施方式
本发明实施例提供了一种永磁同步电机转动惯量辨识的方法及装置,用于解决目前转动惯量离线辨识是在系统模型已知的情况下,获得所有的数据,然后匹配相应的算法,对数据进行估计处理,得到辨识目标参数的估计值,现有的离线辨识方法,如加减速法、转矩限幅加速法、积分辨识法等等,它们原理简单,算法简单,操作起来方便易行,但是其需要进行大量的数据存储和运算,辨识收敛时间长,离线辨识方法具有对计算机要求高,辨识精度不高,不具有实时性的技术问题。
为使得本发明的发明目的、特征、优点能够更加的明显和易懂,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,下面所描述的实施例仅仅是本发明一部分实施例,而非全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
请参阅图1,本发明实施例提供的一种磁同步电机转动惯量辨识的方法的一个实施例包括:
101、通过递推最小二乘法计算辨识参数最小二乘的递推公式;
102、通过预置的性能指标,引入遗传因子对所述辨识参数最小二乘的递推公式进行数据饱和的改进,计算遗忘因子最小二乘参数估计公式;
103、通过电机转矩平衡方程,计算出最小二乘形式,并通过所述遗忘因子最小二乘参数估计公式,计算出k时刻的永磁同步电机的转动惯量J的值J(k);
104、接收预置的满足条件,更新k时刻的永磁同步电机的转动惯量J的值J(k),直至k时刻的永磁同步电机的转动惯量J的值J(k)满足预置的满足条件,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k)。
上面是对一种永磁同步电机转动惯量辨识的方法进行详细的描述,下面将对一种永磁同步电机转动惯量辨识的方法的具体过程进行更详细的描述,本发明实施例提供的一种永磁同步电机转动惯量辨识的方法的另一个实施例包括:
201、接收预置的非奇异方阵A,(A+BC),(I+CA-1B),(A+BC)的逆A-1+X,计算出X=-A-1B(I+CA-1B)-1CA-1,(A+BC)-1=A-1-A-1B(I+CA-1B)-1CA-1,根据 通过递推最小二乘法,计算出
接收预置的A=P-1(k-1),计算出辨识参数最小二乘的递推公式:
取非奇异方阵A,(A+BC),(I+CA-1B),根据矩阵的求逆引理,可逆矩阵A+BCD的逆可以表示成A-1+X,其中X表示未知矩阵,得:
(A+BC)-1=A-1+X (1)
X=-A-1B(I+CA-1B)-1CA-1 (2)
则:(A+BC)-1=A-1-A-1B(I+CA-1B)-1CA-1 (3)
由批处理最小二乘得到下面两式:
最后得到k时刻的递推最小二次估计:
令A=P-1(k-1),代入式(1)得:
将式(6)代入中可以得到:
将(7)式代入(6)式得到协方差矩阵:
综上所述,得到参数最小二乘估计的递推公式为:
202、通过预置的性能指标式中β为遗忘因子,0≤β≤1,对所述辨识参数最小二乘的递推公式进行数据饱和的改进,计算遗忘因子最小二乘参数估计公式为:
为了解决数据饱和的问题,引入遗忘因子对递推公式进行数据饱和的改进,设置性能指标为:
式(13)中β为遗忘因子,0≤β≤1,可以得到遗忘因子最小二乘参数估计公式为:
203、通过永磁同步电机转矩平衡方程:
其中Tl是永磁同步电机的负载转矩,J为转动惯量,B为运动阻尼系数,ωm为转子机械速度,计算出
再进行离散化处理:
计算出最小二乘形式:
接收预置的
y(k)=ωm(k)-2ωm(k-1)+ωm(k-2)
即最小二乘形式为通过遗忘因子最小二乘参数估计式,计算出k时刻的永磁同步电机的转动惯量J的值J(k);
永磁同步电机转矩平衡方程:
式中Tl是永磁同步电机的负载转矩,J为转动惯量,B为运动阻尼系数,ωm为转子机械速度,在控制系统中,B可以忽略,可以得到:
对上式进行离散化处理得:
在采样周期Ts极短的情况下,可以将转矩Tl忽略变化,视为常数,即Tl(k-1)=Tl(k),得最小二乘形式:
令:
y(k)=ωm(k)-2ωm(k-1)+ωm(k-2) (21)
最小二乘形式:将其带入遗忘因子最小二乘参数估计式中,得J(k)。
204、接收预置的遗传因子β的初始值0和预置的当前时刻转动惯量的辨识结果和当前时刻的前一时刻的辨识结果差的绝对值E=|J(k)-J(k-1)小于等于0.000001的满足条件,通过采样,更新k时刻的永磁同步电机的转动惯量J的值J(k),若E≤0.000001,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k),若E≥0.000001,遗传因子β自动加0.01,继续更新k时刻的永磁同步电机的转动惯量J的值J(k),直至E≤0.000001,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k)。
设计一个循环环节,使得输出平稳,针对遗忘因子的变化,定义差值E是辨识算法当前时刻转动惯量的辨识结果和前一时刻的辨识结果差的绝对值,E=|J(k)-J(k-1),为了得到平稳的输出,设E≥0.000001,β加0.01(β的初始值为0,0≤β≤1),算法进行更新J(k),直到E≤0.000001,输出J(k)。具体流程图如图2所示。
本发明对永磁同步电机转动惯量进行实时辨识,以便能够根据因外界负载的变化而带来的转动惯量的变化进行调节器参数的整定,提高伺服系统的精确性和鲁棒性。本发明在最小二乘法中加入了一个遗忘因子,以抵抗数据饱和,并针对遗忘因子,加入了一个遗忘因子变化的约束方法,使遗忘因子的取值适应系统的变化,并很好的辨识出永磁同步电机的转动惯量。
下面对涉及的专用名词的解释
1、永磁同步电机(PMSM)
永磁同步电机是采用永磁体取代其电励磁的同步电机。相对于其他的电机,永磁同步电机有着以下显著的优势:A.永磁同步电机在基速以下不需要励磁电流,在稳定运行时没有转子电阻损耗,可以显著提高功率因数;B.永磁同步电机没有电刷和滑环,结构简单,使用方便,可靠性高;C.永磁同步电机转子结构多样,灵活,可根据不同需要选择不同转子结构。
2、参数辨识
参数辨识是根据实验数据和建立的模型来确定一个或一组参数,使得由模型计算得到的数值能最好的拟合测试数据,进而为生产过程进行预测,提供一定的理论指导。参数辨识主要解决当系统模型已知的条件下,确定模型中的一些未知参数的问题。辨识有3个要素,即数据、模型类和准则,其中数据是指系统的输入输出数据,模型类则定义了模型的基本结构类型,准则即为评价模型与输入输出数据拟合程度的量度标准。
3、最小二乘法
最小二乘法的基本原理如下:用给定的参数计算模型的观测向量和状态变量,将实际系统输出值减去观测值,将差值进行平方和计算,采用某种算法持续改变参数数值,当差值的平方和为最小时,则已经辨识出了系统的参数矩阵,且辨识出的参数值与实际值非常接近。最小二乘法种类很多,有遗忘因子递推算法、加权最小二乘法、递推最小二乘法、多级二乘法等等,算法的核心思想是将电机非线性模型转化为线性模型,得到线性化后的电机模型后与电机参数有着直接联系。在参数辨识中,最小二乘法作为一种在线辨识方法被应用来对所需辨识的参数作实时跟踪估计。
本发明是在递推最小二乘法辨识算法的基础上进行改进的,递推最小二乘法适用于处理数据量大,实时控制中不断有数据输入输出的情况。递推最小二乘算法首先根据输入输出序列列出估计的观测矩阵,给需要辨识的参数和协方差矩阵附初值,计算增益矩阵,最后算出辨识参数。递推最小二乘辨识算法可以比较准确的辨识出转动惯量,但是,它存在的不可缺点,就算数据更新能力弱,出现数据饱和现象,辨识精度不高。转动惯量的辨识方法基本可以归结为两大类:离线辨识和在线辨识。离线辨识是在系统模型已知的情况下,获得所有的数据,然后匹配相应的算法,对数据进行估计处理,得到辨识目标参数的估计值。现有的离线辨识方法,如加减速法、转矩限幅加速法、积分辨识法等等,它们原理简单,算法简单,操作起来方便易行,但是其需要进行大量的数据存储和运算,辨识收敛时间长。离线辨识方法具有对计算机要求高,辨识精度不高,不具有实时性等缺点。针对现有方法的缺点,本发明在递推小二乘法的基础上加入遗忘因子,并对遗忘因子的变化提出相应的约束条件,对电机的转动惯量实现实时辨识,缩短辨识收敛时间,提高辨识精度,解决了现有的离线辨识需要进行大量数据的存储和运算,对计算机要求高,不具有实时性的技术问题。
请参阅图3,本发明实施例提供的一种磁同步电机转动惯量辨识的装置的一个实施例包括:
第一计算单元301,用于通过递推最小二乘法计算辨识参数最小二乘的递推公式;
第二计算单元302,用于通过预置的性能指标,引入遗传因子对所述辨识参数最小二乘的递推公式进行数据饱和的改进,计算遗忘因子最小二乘参数估计公式;
第三计算单元303,用于通过电机转矩平衡方程,计算出最小二乘形式,并通过所述遗忘因子最小二乘参数估计公式,计算出k时刻的永磁同步电机的转动惯量J的值J(k);
更新单元304,用于接收预置的满足条件,更新k时刻的永磁同步电机的转动惯量J的值J(k),直至k时刻的永磁同步电机的转动惯量J的值J(k)满足预置的满足条件,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k)。
上面是对一种永磁同步电机转动惯量辨识的装置各单元进行详细的描述,下面将对一种永磁同步电机转动惯量辨识的装置各附加单元进行详细的描述,请参阅图4,本发明实施例提供的一种磁同步电机转动惯量辨识的装置的另一个实施例包括:
第一计算单元401,用于通过递推最小二乘法计算辨识参数最小二乘的递推公式;
第一计算单元401具体包括:
第一计算子单元4011,用于接收预置的非奇异方阵A,(A+BC),(I+CA-1B),(A+BC)的逆A-1+X,计算出X=-A-1B(I+CA-1B)-1CA-1,(A+BC)-1=A-1-A-1B(I+CA-1B)-1CA-1;
第二计算子单元4012,用于根据
通过递推最小二乘法,计算出
第三计算子单元4013,用于接收预置的A=P-1(k-1), 计算出辨识参数最小二乘的递推公式:
第二计算单元402,用于通过预置的性能指标,引入遗传因子对所述辨识参数最小二乘的递推公式进行数据饱和的改进,计算遗忘因子最小二乘参数估计公式;
第二计算单元402具体包括:
改进子单元4021,用于通过预置的性能指标式中β为遗忘因子,0≤β≤1,对所述辨识参数最小二乘的递推公式进行数据饱和的改进;
第四计算子单元4022,用于计算遗忘因子最小二乘参数估计公式为:
第三计算单元403,用于通过电机转矩平衡方程,计算出最小二乘形式,并通过所述遗忘因子最小二乘参数估计公式,计算出k时刻的永磁同步电机的转动惯量J的值J(k);
第三计算单元403具体包括:
第五计算子单元4031,用于通过永磁同步电机转矩平衡方程:
其中Tl是永磁同步电机的负载转矩,J为转动惯量,B为运动阻尼系数,ωm为转子机械速度,计算出
处理子单元4032,用于对
进行离散化处理:
计算出最小二乘形式:
第六计算子单元4033,用于接收预置的
y(k)=ωm(k)-2ωm(k-1)+ωm(k-2)
即最小二乘形式为通过遗忘因子最小二乘参数估计式,计算出k时刻的永磁同步电机的转动惯量J的值J(k)。
更新单元404,用于接收预置的满足条件,更新k时刻的永磁同步电机的转动惯量J的值J(k),直至k时刻的永磁同步电机的转动惯量J的值J(k)满足预置的满足条件,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k)。
更新单元404具体包括:
接收子单元4041,用于接收预置的遗传因子β的初始值0和预置的当前时刻转动惯量的辨识结果和当前时刻的前一时刻的辨识结果差的绝对值E=|J(k)-J(k-1)小于等于0.000001的满足条件;
更新子单元4042,用于通过采样,更新k时刻的永磁同步电机的转动惯量J的值J(k);
输出子单元4043,用于若E≤0.000001,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k);
循环子单元4044,用于若E≥0.000001,遗传因子β自动加0.01;继续更新k时刻的永磁同步电机的转动惯量J的值J(k),直至E≤0.000001,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k)。
本发明实施例中提供的一种永磁同步电机转动惯量辨识的方法及装置,其中一种永磁同步电机转动惯量辨识的方法,包括:S1:通过递推最小二乘法计算辨识参数最小二乘的递推公式;S2:通过预置的性能指标,引入遗传因子对所述辨识参数最小二乘的递推公式进行数据饱和的改进,计算遗忘因子最小二乘参数估计公式;S3:通过电机转矩平衡方程,计算出最小二乘形式,并通过所述遗忘因子最小二乘参数估计公式,计算出k时刻的永磁同步电机的转动惯量J的值J(k);
S4:接收预置的满足条件,更新k时刻的永磁同步电机的转动惯量J的值J(k),直至k时刻的永磁同步电机的转动惯量J的值J(k)满足预置的满足条件,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k)。本实施例中,通过在递推最小二乘法的基础上加入遗忘因子,并对遗忘因子的变化,设计了循环环节,减少输出震荡,保证了输出的平稳性,具有实时辨识,辨识收敛时间短的优点,解决了目前转动惯量离线辨识是在系统模型已知的情况下,获得所有的数据,然后匹配相应的算法,对数据进行估计处理,得到辨识目标参数的估计值,现有的离线辨识方法,如加减速法、转矩限幅加速法、积分辨识法等等,它们原理简单,算法简单,操作起来方便易行,但是其需要进行大量的数据存储和运算,辨识收敛时间长,离线辨识方法具有对计算机要求高,辨识精度不高,不具有实时性的技术问题。
本发明目的在于实现永磁同步电机的转动惯量的在线辨识,其关键点在于抵抗数据饱和,加入了遗忘因子,并且针对遗忘因子,设计了一个循环环节,保证了输出的稳定性。
现有的离线辨识需要进行大量数据的存储和运算,对计算机要求高,不具有实时性。本发明是在递推最小二乘法的基础上加入遗忘因子,并对遗忘因子的变化,设计了循环环节,减少输出震荡,保证了输出的平稳性。本专利提出的转动惯量辨识的方法属于在线辨识,具有实时辨识,辨识收敛时间短的优点。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的系统,装置和单元的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
在本申请所提供的几个实施例中,应该理解到,所揭露的系统,装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-OnlyMemory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种永磁同步电机转动惯量辨识的方法,其特征在于,包括:
S1:通过递推最小二乘法计算辨识参数最小二乘的递推公式;
S2:通过预置的性能指标,引入遗传因子对所述辨识参数最小二乘的递推公式进行数据饱和的改进,计算遗忘因子最小二乘参数估计公式;
S3:通过电机转矩平衡方程,计算出最小二乘形式,并通过所述遗忘因子最小二乘参数估计公式,计算出k时刻的永磁同步电机的转动惯量J的值J(k);
S4:接收预置的满足条件,更新k时刻的永磁同步电机的转动惯量J的值J(k),直至k时刻的永磁同步电机的转动惯量J的值J(k)满足预置的满足条件,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k)。
2.根据权利要求1所述的永磁同步电机转动惯量辨识的方法,其特征在于,所述步骤S4具体包括:
接收预置的遗传因子β的初始值0和预置的当前时刻转动惯量的辨识结果和当前时刻的前一时刻的辨识结果差的绝对值E=|J(k)-J(k-1)|小于等于0.000001的满足条件,通过采样,更新k时刻的永磁同步电机的转动惯量J的值J(k),若E≤0.000001,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k),若E≥0.000001,遗传因子β自动加0.01,继续更新k时刻的永磁同步电机的转动惯量J的值J(k),直至E≤0.000001,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k)。
3.根据权利要求2所述的永磁同步电机转动惯量辨识的方法,其特征在于,所述步骤S1具体包括:
接收预置的非奇异方阵A,(A+BC),(I+CA-1B),(A+BC)的逆A-1+X,计算出X=-A-1B(I+CA- 1B)-1CA-1,(A+BC)-1=A-1-A-1B(I+CA-1B)-1CA-1,根据 通过递推最小二乘法,计算出第一公式第二公式第三公式
接收预置的计算出辨识参数最小二乘的递推公式:第四公式第五公式第六公式
4.根据权利要求3所述的永磁同步电机转动惯量辨识的方法,其特征在于,所述步骤S2具体包括:
通过预置的性能指标式中β为遗忘因子,0≤β≤1,对所述辨识参数最小二乘的递推公式进行数据饱和的改进,计算遗忘因子最小二乘参数估计公式的第七公式、第八公式、第九公式分别为:
5.根据权利要求4所述的永磁同步电机转动惯量辨识的方法,其特征在于,所述步骤S3具体包括:
通过永磁同步电机转矩平衡方程:
其中Tl是永磁同步电机的负载转矩,J为转动惯量,B为运动阻尼系数,ωm为转子机械速度,计算出第十公式:
再进行离散化处理计算出第十一公式:
计算出最小二乘形式的第十二公式
接收预置的第十三公式、第十四公式、第十五公式分别为:
y(k)=ωm(k)-2ωm(k-1)+ωm(k-2)
即最小二乘形式为通过遗忘因子最小二乘参数估计式,计算出k时刻的永磁同步电机的转动惯量J的值J(k)。
6.一种永磁同步电机转动惯量辨识的装置,其特征在于,包括:
第一计算单元,用于通过递推最小二乘法计算辨识参数最小二乘的递推公式;
第二计算单元,用于通过预置的性能指标,引入遗传因子对所述辨识参数最小二乘的递推公式进行数据饱和的改进,计算遗忘因子最小二乘参数估计公式;
第三计算单元,用于通过电机转矩平衡方程,计算出最小二乘形式,并通过所述遗忘因子最小二乘参数估计公式,计算出k时刻的永磁同步电机的转动惯量J的值J(k);
更新单元,用于接收预置的满足条件,更新k时刻的永磁同步电机的转动惯量J的值J(k),直至k时刻的永磁同步电机的转动惯量J的值J(k)满足预置的满足条件,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k)。
7.根据权利要求6所述的永磁同步电机转动惯量辨识的方法,其特征在于,所述更新单元具体包括:
接收子单元,用于接收预置的遗传因子β的初始值0和预置的当前时刻转动惯量的辨识结果和当前时刻的前一时刻的辨识结果差的绝对值E=|J(k)-J(k-1)|小于等于0.000001的满足条件;
更新子单元,用于通过采样,更新k时刻的永磁同步电机的转动惯量J的值J(k);
输出子单元,用于若E≤0.000001,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k);
循环子单元,用于若E≥0.000001,遗传因子β自动加0.01;继续更新k时刻的永磁同步电机的转动惯量J的值J(k),直至E≤0.000001,输出平稳的k时刻的永磁同步电机的转动惯量J的值J(k)。
8.根据权利要求7所述的永磁同步电机转动惯量辨识的方法,其特征在于,所述第一计算单元具体包括:
第一计算子单元,用于接收预置的非奇异方阵A,(A+BC),(I+CA-1B),(A+BC)的逆A-1+X,计算出X=-A-1B(I+CA-1B)-1CA-1,(A+BC)-1=A-1-A-1B(I+CA-1B)-1CA-1;
第二计算子单元,用于根据
通过递推最小二乘法,计算出第一公式、第二公式、第三公式分别为:
第三计算子单元,用于接收预置的 计算出辨识参数最小二乘的递推公式:第四公式、第五公式、第六公式分别为:
9.根据权利要求8所述的永磁同步电机转动惯量辨识的方法,其特征在于,所述第二计算单元具体包括:
改进子单元,用于通过预置的性能指标式中β为遗忘因子,0≤β≤1,对所述辨识参数最小二乘的递推公式进行数据饱和的改进;
第四计算子单元,用于计算遗忘因子最小二乘参数估计公式为:第七公式、第八公式、第九公式分别为:
10.根据权利要求9所述的永磁同步电机转动惯量辨识的方法,其特征在于,所述第三计算单元具体包括:
第五计算子单元,用于通过永磁同步电机转矩平衡方程:
其中Tl是永磁同步电机的负载转矩,J为转动惯量,B为运动阻尼系数,ωm为转子机械速度,计算出第十公式:
处理子单元,用于对第十公式:
进行离散化处理计算出第十一公式:
计算出最小二乘形式的第十二公式:
第六计算子单元,用于接收预置的第十三公式、第十四公式、第十五公式分别为:
y(k)=ωm(k)-2ωm(k-1)+ωm(k-2)
即最小二乘形式为通过遗忘因子最小二乘参数估计式,计算出k时刻的永磁同步电机的转动惯量J的值J(k)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611139450.4A CN106452247B (zh) | 2016-12-12 | 2016-12-12 | 一种永磁同步电机转动惯量辨识的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611139450.4A CN106452247B (zh) | 2016-12-12 | 2016-12-12 | 一种永磁同步电机转动惯量辨识的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106452247A true CN106452247A (zh) | 2017-02-22 |
CN106452247B CN106452247B (zh) | 2020-06-09 |
Family
ID=58217778
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611139450.4A Expired - Fee Related CN106452247B (zh) | 2016-12-12 | 2016-12-12 | 一种永磁同步电机转动惯量辨识的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106452247B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106787992A (zh) * | 2017-03-24 | 2017-05-31 | 哈尔滨理工大学 | 永磁无刷电机霍尔传感器容错控制方法 |
CN107124132A (zh) * | 2017-05-22 | 2017-09-01 | 珠海格力节能环保制冷技术研究中心有限公司 | 一种转动惯量辨识方法及辨识器 |
CN111464097A (zh) * | 2020-04-03 | 2020-07-28 | 华南理工大学 | 一种永磁同步电机转动惯量离线辨识方法 |
CN112781789A (zh) * | 2020-12-29 | 2021-05-11 | 深圳市雷赛软件技术有限公司 | 转动惯量的识别方法及装置、转动惯量识别系统 |
CN112960501A (zh) * | 2021-02-23 | 2021-06-15 | 杭州优迈科技有限公司 | 一种电梯运行控制方法、装置及电子设备 |
CN114337427A (zh) * | 2021-12-17 | 2022-04-12 | 南京理工大学 | 一种带遗忘因子的递推最小二乘法的转动惯量辨识方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103178758A (zh) * | 2011-12-20 | 2013-06-26 | 上海新世纪机器人有限公司 | 一种永磁同步电动机负载参数辨识方法 |
CN103338003A (zh) * | 2013-06-28 | 2013-10-02 | 西安交通大学 | 一种电机负载转矩及惯量同时在线辨识的方法 |
CN104639003A (zh) * | 2013-11-06 | 2015-05-20 | 沈阳高精数控技术有限公司 | 一种交流伺服系统的转动惯量辨识方法 |
-
2016
- 2016-12-12 CN CN201611139450.4A patent/CN106452247B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103178758A (zh) * | 2011-12-20 | 2013-06-26 | 上海新世纪机器人有限公司 | 一种永磁同步电动机负载参数辨识方法 |
CN103338003A (zh) * | 2013-06-28 | 2013-10-02 | 西安交通大学 | 一种电机负载转矩及惯量同时在线辨识的方法 |
CN104639003A (zh) * | 2013-11-06 | 2015-05-20 | 沈阳高精数控技术有限公司 | 一种交流伺服系统的转动惯量辨识方法 |
Non-Patent Citations (1)
Title |
---|
桑妤: "永磁同步电机转动惯量辨识算法研究", 《中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑)》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106787992A (zh) * | 2017-03-24 | 2017-05-31 | 哈尔滨理工大学 | 永磁无刷电机霍尔传感器容错控制方法 |
CN106787992B (zh) * | 2017-03-24 | 2019-08-02 | 哈尔滨理工大学 | 永磁无刷电机霍尔传感器容错控制方法 |
CN107124132A (zh) * | 2017-05-22 | 2017-09-01 | 珠海格力节能环保制冷技术研究中心有限公司 | 一种转动惯量辨识方法及辨识器 |
CN107124132B (zh) * | 2017-05-22 | 2019-06-14 | 珠海格力电器股份有限公司 | 一种转动惯量辨识方法及辨识器 |
CN111464097A (zh) * | 2020-04-03 | 2020-07-28 | 华南理工大学 | 一种永磁同步电机转动惯量离线辨识方法 |
CN112781789A (zh) * | 2020-12-29 | 2021-05-11 | 深圳市雷赛软件技术有限公司 | 转动惯量的识别方法及装置、转动惯量识别系统 |
CN112781789B (zh) * | 2020-12-29 | 2023-04-18 | 深圳市雷赛软件技术有限公司 | 转动惯量的识别方法及装置、转动惯量识别系统 |
CN112960501A (zh) * | 2021-02-23 | 2021-06-15 | 杭州优迈科技有限公司 | 一种电梯运行控制方法、装置及电子设备 |
CN114337427A (zh) * | 2021-12-17 | 2022-04-12 | 南京理工大学 | 一种带遗忘因子的递推最小二乘法的转动惯量辨识方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106452247B (zh) | 2020-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106452247A (zh) | 一种永磁同步电机转动惯量辨识的方法及装置 | |
Sakthivel et al. | Multi-objective parameter estimation of induction motor using particle swarm optimization | |
Zerdali et al. | The comparisons of optimized extended Kalman filters for speed-sensorless control of induction motors | |
CN107645257B (zh) | 机械学习器、电动机控制装置以及机械学习方法 | |
Bai et al. | A novel parameter identification approach via hybrid learning for aggregate load modeling | |
CN108896924B (zh) | 基于深度置信网络和卡尔曼滤波融合的锂离子电池荷电状态估计方法 | |
CN108763831B (zh) | 一种水泵水轮机线性模型参数辨识方法 | |
CN104639003B (zh) | 一种交流伺服系统的转动惯量辨识方法 | |
CN103219939B (zh) | 交流永磁同步电机伺服系统的转动惯量在线辨识方法 | |
CN103944481B (zh) | 一种交流异步电机矢量控制系统模型参数在线修正方法 | |
CN110932274A (zh) | 一种电力系统测量与负荷参数分析辨识方法 | |
ÇANAKOĞLU et al. | Induction motor parameter estimation using metaheuristic methods | |
CN110988674A (zh) | 永磁同步电机的健康状态监测方法、系统及移动终端 | |
CN109066651B (zh) | 风电-负荷场景的极限传输功率的计算方法 | |
CN104009696A (zh) | 一种基于滑模控制的交互式模型参考自适应速度与定子电阻的辨识方法 | |
Rutczyńska-Wdowiak | Replacement strategies of genetic algorithm in parametric identification of induction motor | |
CN106059412B (zh) | 基于信度规则库推理的他励直流电动机转速控制方法 | |
Nikranajbar et al. | Parameter identification of a cage induction motor using particle swarm optimization | |
CN106849801A (zh) | 一种无轴承异步电机转速估计方法 | |
Bila | Power system dynamic state estimation and load modeling | |
CN112564557B (zh) | 一种永磁同步电机的控制方法、装置、设备及存储介质 | |
CN103986400B (zh) | 基于二维模糊控制的模型参考自适应系统参数自整定方法 | |
CN110927597B (zh) | 一种确定电池放电曲线的方法 | |
CN112711542B (zh) | 一种用于手机客户端进行测试环境切换的方法 | |
CN112149896A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200609 Termination date: 20201212 |