CN109167546B - 一种基于数据生成模型的异步电机参数在线辨识方法 - Google Patents
一种基于数据生成模型的异步电机参数在线辨识方法 Download PDFInfo
- Publication number
- CN109167546B CN109167546B CN201811119061.4A CN201811119061A CN109167546B CN 109167546 B CN109167546 B CN 109167546B CN 201811119061 A CN201811119061 A CN 201811119061A CN 109167546 B CN109167546 B CN 109167546B
- Authority
- CN
- China
- Prior art keywords
- model
- identification
- motor
- stage
- parameter
- 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
- 238000000034 method Methods 0.000 title claims abstract description 75
- 238000012549 training Methods 0.000 claims abstract description 25
- 230000008569 process Effects 0.000 claims abstract description 14
- 238000009826 distribution Methods 0.000 claims description 25
- 230000008859 change Effects 0.000 claims description 8
- 238000013499 data model Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 2
- 230000005284 excitation Effects 0.000 claims description 2
- 230000006870 function Effects 0.000 claims description 2
- 101000744515 Homo sapiens Ras-related protein M-Ras Proteins 0.000 description 6
- 102100039789 Ras-related protein M-Ras Human genes 0.000 description 6
- 238000010586 diagram Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005389 magnetism Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
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
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Electric Motors In General (AREA)
- Control Of Ac Motors In General (AREA)
Abstract
本发明公开了一种基于数据生成模型的异步电机参数在线辨识方法。该方法融合了基于模型的辨识方法和基于数据的辨识方法的优点,整体思路为:①,首先通过离线数据建立输出转矩与状态间的概率模型,②,在线辨识过程中,将基于模型的辨识方法作为“数据产生器”,再采用基于数据的方法,结合①中的转矩概率模型,进行电机参数的在线辨识,最终得到电机参数辨识值的概率模型。本发明提出的方法的优点在于:辨识过程中,训练数据集很容易获得;辨识的参数不受模型误差的影响;辨识的参数使电机在运行中保持最优转矩电流比;辨识的参数考虑了电机运行过程中的不确定性,具有很高的鲁棒性。
Description
技术领域
本发明涉及电机控制技术领域,特别涉及一种基于数据生成模型的异步电机参数在线辨识方法。
背景技术
异步电机在运行过程中,输出特性会随着内部参数的变化而改变。因此,有必要在电机运行的过程中,在线的对其参数进行辨识。传统的异步电机参数方法大多是基于电机的电压-电流动力学模型或等效电路模型,称为基于模型的方法,如模型参考自适应法(MRAS)、滑模观测器法(SMO)以及线性最小二乘法(LSM)等等。然而,基于模型的方法具有如下缺点:
1)参数的辨识精度和辨识收敛性极大的依赖于电机的等效模型,容易受到噪声和模型稳定性的影响。因此,在某些极端工况下(例如极低速、高速弱磁等),存在着辨识精度低,鲁棒性差等缺点。
2)只能辨识出参数的物理值,然而,在一些情况下,这些参数的物理值并不能使电机运行在最优状态。例如,在使用间接矢量控制方法对电机进行控制时,为了保证转矩电流比最优,有时会人为的过估计或欠估计电机的转子时间常数,而并非使用转子时间常数的物理值。
3)电机在运行的过程中,参数的变化会经常受到各种不确定因素的影响,例如温度,湿度,振动,老化等等。然而,基于模型的辨识方法大多使用简化的等效模型,并没有考虑上述不确定性因素,因此,随着电机运行过程中不确定性因素的增多,辨识值与实际值之间的偏差会越来越大。
针对基于模型方法的上述缺点,学者和工程人员提出了一系列基于实测数据的电机参数辨识方法。如基于支持向量机的辨识方法(SVM)、基于粒子群优化的辨识方法(PSO)以及基于深度Q学习的辨识方法(DQL)等。此类方法不使用电机的等效模型,而是完全依赖于实测数据,因此不会受到模型误差的制约,具有更高的精确性和鲁棒性。然而,这类方法同样也存在一些缺陷:
1)基于数据的方法需要大量的带有标签的数据进行训练,此类数据的获取较为困难。
2)训练的过程为离线运行,和基于模型的方法一样,无法消除在线运行时各种不确定性因素的影响。因此,仅仅适用于电机参数的离线辨识。
发明内容
针对上述两类方法的缺陷,本发明提出了一种基于数据生成模型的异步电机参数在线辨识方法,该方法的思路为:在电机的运行过程中使用基于模型的方法实时产生训练数据,再使用该训练数据实时地生成一个概率模型。最后,通过生成的概率模型进行电机参数的辨识。
本发明采用的技术方案为:一种基于数据生成模型的异步电机参数在线辨识方法,本发明所述的方法分为两个阶段:
阶段1,离线训练阶段。此阶段的目标为训练出一个转矩预测概率模型,使得电机可以在任意给定状态下预测出输出转矩的大小。将训练好的模型用于在线辨识阶段,从而保证了在线辨识阶段可以辨识出使得电机保持转矩电流比最大的参数。在此阶段下电机需在装有转矩传感器的对拖台架上带载运行。运行时,记录当前时刻d-q轴电流、d-q轴电压,记为{id *,iq *,ud *,uq *},作为训练集的输入数据x*;记录当前时刻输出转矩,作为训练集的标签数据在不同时刻t记录x*和记为考虑不确定性,假设服从以x*为自变量的高斯分布,因此可以高斯过程回归,训练出在未知输入x下输出转矩Te的概率分布模型,记为:
阶段2,在线辨识阶段,此阶段的目的是使电机在运行过程中可以在线的调整参数,使得电机永远保持最优转矩电流比。在此阶段下电机运行于实际工况中,首先采用传统基于模型的方法进行参数辨识,记录参数变化前的d-q轴电流、d-q轴电压,记为xt={idt,iqt,udt,uqt},当参数辨识值变化后,记录变化后的参数辨识值θt,以及参数变化后的d-q轴电流、d-q轴电压,记为xt+1={idt+1,iqt+1,udt+1,uqt+1}。然后建立在θt的动作下,从xt到xt+1的状态转移模型,即在xt以及θt条件下,对下一时刻状态xt+1的预测模型。考虑不确定性,假设xt与xt+1均服从高斯分布,使用高斯过程回归,生成在xt以及θt条件下xt+1出现的概率模型p(xt+1|xt,θt)。同样p(xt+1|xt,θt)也服从高斯分布,即
最终目标问题是是在t时刻,观测为xt的条件下,如何选择θ,使得在t+1时刻,转矩Te最大,记为argmaxθ(Te|xt,θt)。则首先要求出p(Te|xt,θt)的分布。由统计学知识得:
p(Te|xt,θt)=∫p(Te|xt+1)p(xt+1|xt,θt)dxt+1 (2)
其中p(Te|xt+1)由阶段1算得,p(xt+1|xt,θt)由阶段2算得。
计算出p(Te|xt,θ)后,取p(Te|xt,θ)均值μ(xt,θt),并令:
即可计算出参数辨识值。
经过若干次计算后,将传统基于模型方法切换为数据生成模型,即可完成在线辨识。
本发明与现有技术相比的优势有:
相对于基于数据的方法,其优势在于:
1)本发明使用传统的基于模型的方法生成训练数据,解决了训练数据难以获取的问题。
2)数据的生成和训练都是实时的,适用于电机参数的在线辨识问题。
相对于基于模型的方法,其优势在于:
3)该方法并辨识的参数具有使电机可以保持转矩电流比最大的特性。而传统的基于模型的方法只能辨识电机参数的物理值,不具有使得电机保持转矩电流比最大的特性。
4)本发明考虑了电机运行过程中的各种不确定性问题,使用贝叶斯推断生成一个概率模型,相比传统的基于模型的方法本发明的模型更具鲁棒性。
附图说明
图1为离线训练阶段所需环境;
图2为在线辨识框图;
图3为将数据模型替换为MRAS模型框图。
具体实施方式
下面结合附图以及具体实施方式进一步说明本发明。
本发明公开了一种基于数据生成模型的异步电机参数在线辨识方法,具体如下:
1、离线训练阶段
离线训练阶段所需环境如图1所示,被测电机安装在电机对拖台架上,使用电机控制器使得被测电机运行于转矩模式,测功机电机运行于转速模式。使用转矩传感器测量被测电机的输出转矩Te。
①电机运行过程中,每隔5秒将被测电机的d-q轴电流、d-q轴电压数据x={id,iq,ud,uq}以及被测电机的转矩数据Te传输给数据采集器,然后数据采集器通过CAN-BUS总线将数据传输给计算机作为训练集数据,同时每隔30秒随机调节一次被测电机的转矩和测功机电机的转速。整个过程若干分钟。最终,采集n组x作为训练集的输入数据,记为 同时采集n组Te作为训练集的标签数据,记为则最终训练数据集为:
②考虑到数据采集时的不确定性,设x和Te均服从高斯分布,可使用高斯过程回归生成Te在x的条件下的预测模型,记为:
Te=f(x)+ε (4)
式(4)的具体求解步骤为:
给定输入x’={x’1,…,x’m}时,f(x’)与训练集的联合分布为:
利用高斯分布的性质,得到:
其中:
考虑噪声,则最终可得转矩的估计值Te’的概率分布:
式(4)求解完成。
2、在线辨识阶段
在线辨识框图如图2所示,电机运行在实际工况下,控制算法为有速度传感器的间接矢量控制,其中xt为电机参数变化前状态{idt,iqt,udt,uqt},θt为变化后的辨识参数值,辨识的参数可以是定子电阻,转子电阻,定转子漏感以及激磁电感。xt+1为已知xt和θt的条件下,电机参数变化后的状态{idt+1,iqt+1,udt+1,uqt+1}。传统的基于模型的辨识方法包括但不限于使用MRAS,SMO,LSM等方法。为表述方便,图2中以模型参考自适应(MRAS)为例,作为传统的基于模型的辨识方法进行说明。
先使用MRAS方法进行参数辨识,每隔1分钟,采集xt,θt作为训练数据。并把当前的xt数据设为上一分钟的xt+1,考虑到实际运行中的不确定因素,设xt和xt+1之间并非一一对应的映射关系,而是xt+1服从以xt为条件的高斯分布,则可使用贝叶斯推断生成电机的状态转移模型,即:在已知观测为xt,选择辨识参数的值为θt的条件下,状态xt+1发生的概率,记为:
式(10)的求解步骤如下:
使用全概率公式以及贝叶斯推断可得:
p(xt+1|xt,θt)=∫p(xt+1|xt,θt,ht)p(ht|xt,xt-1,θt-1)dht (11)
其中ht为状态转移概率模型的参数,满足:
p(ht|xt,xt-1,θt-1)=p(xt|xt-1,θt-1,ht)p(ht-1|xt-1,xt-2,θt-2) (12)
由式(12)可知当前概率模型的参数ht与上一时刻的ht-1,xt-1,θt-1有关,设初始时刻x0、θ0为随机值,则初始时刻模型参数将式(12)代入(11),通过链式法则可以求得p(xt+1|xt,θt)。
最终可得,在xt,θt的条件下,下一时刻转矩Te的概率模型为:
p(Te|xt,θt)=∫p(Te|xt+1)p(xt+1|xt,θt)dxt+1 (13)
根据高斯分布性质可知p(Te|xt,θt)同样服从高斯分布,记为:
其中μ(xt,θt)表示分布的均值为xt,θt的函数,同理,Σ(xt,θt)表示分布的协方差为xt,θt的函数。
式(10)求解完成。
令:
计算得到θtmax,即为当前状态xt时,能使得转矩Te为最大的参数辨识值。易知θtmax能使电机保持最优转矩电流比。
经过若干次计算后,可认为数据模型已经收敛,将数据模型替换为MRAS模型,如图3所示。此时,只需记录当前状态xt,便可直接预测下一时刻转矩为最大的参数辨识值,从而实现在线辨识。
Claims (3)
1.一种基于数据生成模型的异步电机参数在线辨识方法,其特征在于:所述的方法分为两个阶段:
阶段1,离线训练阶段,此阶段的目标为训练出一个转矩预测模型,使得电机可以在任意给定状态下预测出输出转矩的大小,从而保证了在线辨识阶段,辨识出的参数值可以使得电机保持最优转矩电流比;在此阶段下电机需在装有转矩传感器的对拖台架上带载运行;运行时,记录当前时刻d-q轴电流id,iq、d-q轴电压ud,uq,记为{id *,iq *,ud *,uq *},作为训练集的输入数据x*;记录当前时刻输出转矩Te,作为训练集的标签数据在不同时刻t记录x*和记为考虑不确定性,假设服从以x*为自变量的高斯分布,因此可以高斯过程回归,训练出在未知输入x下输出转矩Te的概率分布模型,记为:
阶段2,在线辨识阶段,此阶段的目的是使电机在运行过程中可以在线的调整参数,使得电机保持最优转矩电流比;在此阶段下电机运行于实际工况中,首先采用基于模型的方法进行参数辨识,记录t时刻参数变化前的d-q轴电流、d-q轴电压,记为xt={idt,iqt,udt,uqt},当参数辨识值变化后,记录变化后的参数辨识值θt,辨识的参数可以是定子电阻,转子电阻,定转子漏感以及激磁电感,以及参数变化后t+1时刻的d-q轴电流、d-q轴电压,记为xt+1={idt+1,iqt+1,udt+1,uqt+1};假设xt与xt+1均服从高斯分布,使用高斯过程回归,生成在xt以及θt条件下xt+1出现的概率模型p(xt+1|xt,θt),同样p(xt+1|xt,θt)也服从高斯分布,即(μ(xt,θt),Σ(xt,θt));
最终目标问题是在t时刻,观测为xt的条件下,如何选择参数辨识值θ,使得在t+1时刻,转矩Te最大,记为argmaxθ(Te|xt,θt),则首先要求出p(Te|xt,θt)的分布;由统计学知识得:
p(Te|xt,θt)=∫p(Te|xt+1)p(xt+1|xt,θt)dxt+1 (2)
其中p(Te|xt+1)由阶段1算得,p(xt+1|xt,θt)由阶段2算得;
计算出p(Te|xt,θ)后,取p(Te|xt,θ)均值μ(xt,θt),并令:
即可计算出参数辨识值;
经过若干次计算后,将基于模型方法切换为数据生成模型,即可完成在线辨识。
2.根据权利要求1所述的一种基于数据生成模型的异步电机参数在线辨识方法,其特征在于:①阶段1中,只需采集{id,iq,ud,uq}作为状态量x,无需再采集其他信号;②阶段1中,考虑到运行中的不确定性因素,假定状态变量x和输出转矩Te并非一一映射关系,而是Te服从以x为自变量的高斯分布,可以用高斯过程回归对Te和x的映射进行概率建模。
3.根据权利要求1所述的一种基于数据生成模型的异步电机参数在线辨识方法,其特征在于:①,阶段2中,初始状态下,使用基于模型的方法进行在线参数辨识,同时,使用基于模型的方法产生的数据来生成数据模型;②,阶段2中,先使用基于模型的方法进行参数动作,再采集参数变化前的状态xt,用来计算生成数据模型的先验概率,同时采集参数变化后的状态xt+1,用来计算参数生成数据模型的似然概率,最终可以根据贝叶斯推断生成在线参数辨识的数据模型;③,阶段2中,最终生成的模型仍然是一个服从高斯分布的概率模型,并且该模型的均值和方差均为辨识参数θ的函数,只需取该模型的均值,并对其均值求极值,便可得到辨识参数θ的解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811119061.4A CN109167546B (zh) | 2018-09-25 | 2018-09-25 | 一种基于数据生成模型的异步电机参数在线辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811119061.4A CN109167546B (zh) | 2018-09-25 | 2018-09-25 | 一种基于数据生成模型的异步电机参数在线辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109167546A CN109167546A (zh) | 2019-01-08 |
CN109167546B true CN109167546B (zh) | 2021-08-03 |
Family
ID=64880425
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811119061.4A Active CN109167546B (zh) | 2018-09-25 | 2018-09-25 | 一种基于数据生成模型的异步电机参数在线辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109167546B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113459824B (zh) * | 2020-03-31 | 2023-04-28 | 安徽威灵汽车部件有限公司 | 电动汽车抖动抑制方法、装置以及电动汽车、存储介质 |
CN113459822B (zh) * | 2020-03-31 | 2023-04-28 | 安徽威灵汽车部件有限公司 | 电动汽车抖动抑制方法、装置以及电动汽车、存储介质 |
CN111814408B (zh) * | 2020-08-25 | 2024-02-09 | 安徽大学 | 一种基于条件变分自编码器的电机数据生成模型建模方法 |
CN112886890B (zh) * | 2021-02-07 | 2022-11-11 | 安徽大学 | 一种数据驱动的异步电机动力学模型建模方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106705999A (zh) * | 2016-12-21 | 2017-05-24 | 南京航空航天大学 | 一种无人机陀螺仪故障的诊断方法 |
CN107147342A (zh) * | 2017-05-31 | 2017-09-08 | 合肥申芯电子技术有限责任公司 | 一种感应电机参数辨识系统及方法 |
CN107861061A (zh) * | 2017-10-27 | 2018-03-30 | 安徽大学 | 一种数据驱动的感应电机参数在线辨识方法 |
-
2018
- 2018-09-25 CN CN201811119061.4A patent/CN109167546B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106705999A (zh) * | 2016-12-21 | 2017-05-24 | 南京航空航天大学 | 一种无人机陀螺仪故障的诊断方法 |
CN107147342A (zh) * | 2017-05-31 | 2017-09-08 | 合肥申芯电子技术有限责任公司 | 一种感应电机参数辨识系统及方法 |
CN107861061A (zh) * | 2017-10-27 | 2018-03-30 | 安徽大学 | 一种数据驱动的感应电机参数在线辨识方法 |
Non-Patent Citations (2)
Title |
---|
Parameter Identification of Linear Induction Motor Model in Extended Range of Operation by Means of Input-Output Data;Francesco Alonge等;《IEEE Transactions on Industry Applications》;20140430;第50卷(第2期);第959-972页 * |
基于模糊模型的交流异步电力测功机转矩和转速辨识;陈鸿蔚等;《中国机械工程》;20090125;第20卷(第2期);第162-167页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109167546A (zh) | 2019-01-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109167546B (zh) | 一种基于数据生成模型的异步电机参数在线辨识方法 | |
CN113343633B (zh) | 动力锂电池热失控故障分类及风险预测方法、系统 | |
Chen et al. | Combining empirical mode decomposition and deep recurrent neural networks for predictive maintenance of lithium-ion battery | |
CN109800537B (zh) | 一种基于深度神经网络和蒙特卡洛法的机床热误差模型可靠度计算方法 | |
Wu et al. | Modeling and analysis of tool wear prediction based on SVD and BiLSTM | |
Berghout et al. | A semi-supervised deep transfer learning approach for rolling-element bearing remaining useful life prediction | |
Chen et al. | Deep attention relation network: A zero-shot learning method for bearing fault diagnosis under unknown domains | |
CN108375474B (zh) | 一种航空发动机过渡态关键性能参数预测方法 | |
CN110763997A (zh) | 一种同步电机定子早期故障预警方法 | |
CN116261690A (zh) | 提供用于高炉热控制的操作指令的计算机系统和方法 | |
US11853047B2 (en) | Sensor-agnostic mechanical machine fault identification | |
CN111241748A (zh) | 基于长短期记忆模型循环神经网络的风力发电机故障诊断 | |
CN103649858A (zh) | 用于预测组件或系统状况的方法和设备及计算机程序产品 | |
CN111680875A (zh) | 基于概率基线模型的无人机状态风险模糊综合评价方法 | |
CN115859077A (zh) | 一种变工况下多特征融合的电机小样本故障诊断方法 | |
Lin et al. | Master general parking skill via deep learning | |
CN111537888A (zh) | 一种数据驱动的梯次电池soc预测方法 | |
CN113485261A (zh) | 一种基于CAEs-ACNN的软测量建模方法 | |
Lu et al. | GAN-LSTM predictor for failure prognostics of rolling element bearings | |
CN111075647B (zh) | 一种基于elm的变速风电机组最大风能捕获方法 | |
CN108959787A (zh) | 考虑实际工况的宏宏双驱动系统的热变形预测方法及系统 | |
Yang et al. | A balanced deep transfer network for bearing fault diagnosis | |
CN114330571A (zh) | 基于迁移学习的电力系统主导失稳模式识别方法及系统 | |
Melendez et al. | Self-supervised Multi-stage Estimation of Remaining Useful Life for Electric Drive Units | |
CN117113870A (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 |