阅读说明：本技术 一种基于自适应扩展卡尔曼滤波的电机参数估计方法 (Motor parameter estimation method based on adaptive extended Kalman filtering ) 是由 丁洁 曹正鑫 林金星 于 2020-06-23 设计创作，主要内容包括：本发明公开了一种自适应扩展卡尔曼滤波的电机参数估计方法,包括：将在αβ静止参考坐标系下的电流方程作为电机的连续状态空间表达式描述电机的内部动态特性；离散化状态空间方程同时将所需辨识参数增广到系统状态中；仿真电机采用双闭环控制方式,通过从电流检测单元和电压检测单元获得相电流、相电压；用当前时刻状态估计值代替系统真实值,并考虑满足白噪声时刻无关的性质,估计出当前时刻的系统噪声协方差矩阵；对噪声协方差估计值进行加权,得到当前时刻估计的噪声协方差矩阵；通过获得的相电流、相电压,利用已计算得到噪声协方差矩阵的扩展卡尔曼滤波算法估计反电动势同时辨识电机参数。该方法可提升电机参数估计的精度。(The invention discloses a motor parameter estimation method of adaptive extended Kalman filtering, which comprises the following steps: taking a current equation under an alpha beta static reference coordinate system as a continuous state space expression of the motor to describe the internal dynamic characteristics of the motor; the discretization state space equation simultaneously expands the needed identification parameters into the system state; the simulation motor adopts a double closed-loop control mode, and phase current and phase voltage are obtained from a current detection unit and a voltage detection unit; replacing a system true value with a current time state estimation value, and estimating a system noise covariance matrix at the current time by considering the property of meeting the white noise moment independence; weighting the noise covariance estimated value to obtain a noise covariance matrix estimated at the current moment; and estimating the back electromotive force and identifying the motor parameters by utilizing an extended Kalman filtering algorithm of the noise covariance matrix obtained by calculation through the obtained phase current and phase voltage. The method can improve the precision of motor parameter estimation.)

一种基于自适应扩展卡尔曼滤波的电机参数估计方法

技术领域

本发明涉及电机参数估计方法，特别涉及一种基于自适应扩展卡尔曼滤波的电机参数估计方法。

背景技术

无刷直流电机以其使用寿命长，控制简单，运行可靠等优点被广泛应用于各种工业设备中，如：工业机器人。目前用于控制无刷直流电机的控制器有很多，如各种型号的单片机、DSP、FPGA等，且大多采用PWM和PID的控制方式。而电感、电阻参数不只是实现电机精确控制的前提条件，还是变频器控制算法实现、电机性能分析和优化设计的重要依据，所以为了保证控制精度，提高电机性能，需要在工程实际中对电机参数进行辨识。

无刷直流电机(BLDCM)以其寿命长、控制简单、运行可靠等优点，广泛应用于工业机器人、数控机床等工业设备中。目前，对无刷直流电机的热点主要是其设计与控制，包括无刷电机本体设计、转矩脉动、无传感器控制、转矩控制等。大多实现无传感器控制和转矩控制的关键是获得电机的准确、实时的反电动势，一般认为电机的反电势是理想的梯形波，但会造成控制精度低；或是通过查表计算控制方案中的反电势值，然而，这增加了操作步骤。

当将电机反电动势当作状态变量时，可以用状态观测器或滤波器进行估计。问题是状态估计方法依赖于精确的电机参数。对于无刷直流电动机，换流过程中电流变化较大，电感和电阻参数的偏差对反电势估计有很大影响。因此，为了保证控制精度，提高电机性能，在工程实践中对电机参数进行辨识是十分必要的。参数识别方法主要包括最小二乘法、滑模辨识、模型参考自适应法、人工神经网络和卡尔曼滤波器等。卡尔曼滤波对于多维状态和非平稳过程是一种有效的滤波方法，即使在有噪声干扰的情况下也具有很高的精度。由于非线性，它是电机参数辨识中最常用的方法之一。扩展卡尔曼滤波算法计算简单，易于跟踪。但是，由于只保留了Taylor展开的一次项，所以该算法得到的结果可能存在误差，且系统的非线性程度越强，预测结果的误差就越大。未知噪声或估计精度差可能影响电机的精度控制。

发明内容

本发明的目的是提供一种精确度高、滤波不会产生较大的偏差、发散的基于自适应扩展卡尔曼滤波的电机参数估计方法。

本发明的基于自适应扩展卡尔曼滤波的电机参数估计方法，包括如下步骤：

(1)将在αβ静止参考坐标系下的电流方程作为电机的连续状态空间表达式描述电机的内部动态特性；

(2)离散化状态空间方程同时将所需辨识参数增广到系统状态中；

(3)仿真电机采用双闭环控制方式，通过从电流检测单元和电压检测单元获得相电流、相电压；

(4)用当前时刻状态估计值代替系统真实值，并考虑满足白噪声时刻无关的性质，估计出当前时刻的系统噪声协方差矩阵；对噪声协方差估计值进行加权，得到当前时刻估计的噪声协方差矩阵；

(5)通过获得的相电流、相电压，利用己计算得到噪声协方差矩阵的扩展卡尔曼滤波算法估计反电动势同时辨识电机参数。

进一步地，所述步骤(1)中在αβ静止参考坐标系下的电流方程如下：

其中i_α，i_β是定子电流，e_α，e_β是绕组反电势，u_α，u_β是αβ坐标系中的定子电压；

R为定子电阻，L为定子电感，等于绕组自感与互感的差，

状态空间表达式表示为：

其中x₀＝[i_α，i_α，e_α，e_β]^T为系统状态向量，z＝[i_α，i_β]^T为系统输出，u＝[u_α，u_β]^T为系统输入，以及

进一步地，将L和R增广到系统状态向量中，得到以下六阶模型：

其中x＝[x₀ ^T，L，R]^T和

采用一阶欧拉离散化方法得到一个离散时间模型，以T_s为采样周期，

A_k＝(I+T_sA)，B_k＝T_sB，H_k＝T_sH

由于L和R是未知的，所以A_kx_k+B_ku_k是一个非线性函数，用f_k(x_k，u_k)表示非线性函数，

x_k＝f_k-1(x_k-1，u_k-1)+w_k-1

Z_k＝H_kx_k+v_k

其中w_k和v_k是系统噪声和测量噪声。

进一步地，所述步骤(3)的仿真电机采用两两120°导通模式控制，在双闭环控制模式下，电机转速的给定值与霍尔信号计算的实际转速值之差作为输入速度控制器输入，速度控制器的输出和电流检测单元采集到的当前电流反馈量的差为电流控制器的输入，PWM控制信号产生器，根据霍尔传感器判读当前的电机转子位置，然后将PWM控制信号发生器连接到需要打开的功率管上，以完成电机的调速，同时，通过从电流检测单元和电压检测单元获得相电流和相电压。

进一步地，所述步骤(4)的具体过程为：

为了估计噪声协方差矩阵Q_k-1，将x_k用代替，可得：

进而可得：

同时，希望过程噪声估计不同时刻无关，即满足：

最终得噪声协方差的估计式：

为了避免某些系统状态在H中对应的系数全为0，从而使得对应的的值为0的情况，

其中，λ为(0，1)的常数。

进一步地，估计过程如下：

K_k＝P_k/k-1H_k ^T(H_kP_k/k-1H_k ^T+R_k)^-1 (3)

将得到的

替换原来的

回代到(2)-(4)式

P_k＝(I-K_kH_k)P_k/k-1。

有益效果：本发明充分利用了电机模型的输入输出数据，提升了参数估计的精度，简单易于实现，滤波偏差较小，不易于发散。

附图说明

图1为本发明实施例中所述基于自适应扩展卡尔曼滤波的无刷直流电机的参数估计方法的流程框图示意；

图2为本发明仿真中无刷直流电机控制框架图；

图3为本发明仿真中EKF和AEKF对电机状态e_α估计值的比较图；

图4为本发明仿真中EKF和AEKF对电机参数L估计值和误差图；

图5为本发明仿真中EKF和AEKF对电机参数R估计值和误差图。

具体实施方式

本实施例的基于自适应扩展卡尔曼滤波的无刷直流电机的参数估计方法。

参阅图1，所述方法具体包括如下步骤：

1、利用在αβ静止参考坐标系下电流方程作为电机的连续状态空间表达式描述电机的内部动态特性。

2、离散化状态空间方程同时将所需辨识参数增广到系统状态中；首先，将电机参数L，R当作电机状态增广到原电机状态中，得到六阶模型：其中x＝[x₀ ^T，L，R]^T和

H＝[H₀，0_2×2]；其次利用一阶欧拉离散化得到一个离散时间模型，它以T_s为采样周期：A_k＝(I+T_sA)，B_k＝T_sB，H_k＝T_sH。由于L和R是未知的，所以A_kx_k+B_ku_k是一个非线性函数。用f_k(x_k，u_k)表示非线性函数：x_k＝f_k-1(x_k-1，u_k-1)+w_k-1，z_k＝H_kx_k+v_k

其中w_k和v_k是系统噪声和测量噪声。

3、仿真系统中永磁同步直流电机实验电机为表贴式永磁无刷直流电机，采用两两120°导通模式控制，使用双闭环控制方式，如图2所示，电机转速的给定值与霍尔信号计算的实际转速值之差作为输入速度控制器输入。速度控制器的输出和电流检测单元采集到的当前电流反馈量的差为电流控制器的输入。PWM控制信号产生器，根据霍尔传感器可以判读当前的电机转子位置，然后将PWM控制信号发生器连接到需要打开的功率管上，以完成电机的调速。同时，通过从电流检测单元和电压检测单元获得相电流和相电压。

4、由传统的扩展卡尔曼滤波算法可知：

为了估计Q_k-1，将x_k用

代替，可得：

进而可得：

同时，希望过程噪声估计

不同时刻无关，即满足：

等价于：(I-K_kH_k)F_k-1P_k-1F_k-1 ^T(H_k ^TK_k ^T-I)+(I-K_kH_k)Q_k-1H_k ^TK_k ^T-K_kR_kK_k ^T＝0可得：使用代替Q_k-1，则可以将修正后的噪声协方差修改为

为了避免某些系统状态在H中对应的系数全为0，从而使得

对应的的值为0的情况：其中，λ为(0，1)的常数

5、设采样周期T_s＝2×10^-6，λ＝0.7，各初值设为x₀＝[0，0，0，0，0.01，0.5]^T，P₀＝diag[1，1，1，1，1，1]^T，Q₀＝diag(0.01，0.01，10，10，0，0)R＝diag(1000，1000)。通过获得的相电流和相电压，分别使用扩展卡尔曼滤波算法和自适应扩展卡尔曼滤波算法，估计电机反电动势并同时辨识电机参数。通过图3能够看出，两种算法估计的反电动势的误差逐渐减小，但扩展卡尔曼滤波算法估计初始误差相对较大，自适应扩展卡尔曼滤波算法估计的误差较前者略小；由图4可知，两种算法估计的电机参数L误差均较小，且收敛速度快，但自适应扩展卡尔曼滤波算法估计的误差更小；由图5可以看出，使用自适应扩展卡尔曼滤波算法估计R精度较传统的扩展卡尔曼滤波方法有了很大的提升，并且收敛速度更快。对比可知，本发明方法估计无刷直流电机参数具有更好的精确度，从而起到提高电机反电动势估计的准确性，提高电机的控制精度。