CN105424277A - 一种获得电机转动惯量的方法和装置 - Google Patents

Abstract

本申请实施例提供了一种获得电机转动惯量的方法和装置，用于获得电机的转动惯量。所述方法包括：获得电机在第n-1周期的第n-1实际转速；n为正整数；获得与所述电机对应的可调模型输出的第n-1推定转速；所述第n-1推定转速为所述可调模型基于第n-1轮计算结果输出的转速；根据所述第n-1实际转速和所述第n-1推定转速，获得第n轮计算结果；所述第n轮计算结果中包括第n转动惯量相关量；判断所述第n转动惯量相关量和所述第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系；当所述第n转动惯量相关量和所述第n-1转动惯量相关量满足所述预设关系时，确定所述第n转动惯量相关量对应的第n转动惯量为所述电机的实际转动惯量。

Description

一种获得电机转动惯量的方法和装置

技术领域

本发明涉及电子技术领域，尤其涉及一种获得电机转动惯量的方法和装置。

背景技术

转动惯量是电机的一个重要参数。目前，在对电机进行控制时，通常直接采用电机生产商提供的规格参数当中的转动惯量。然而，转动惯量并不是一个恒定的值，它与电机实际运行状况有关系。例如电机带风机叶片和不带风机叶片运行，两种情况下的转动惯量差别非常大。所以，测试规格参数的环境与电机实际应用环境的不同，将导致规格参数中给定的转动惯量和实际转动惯量之间存在或大或小的误差。而给定转动惯量与实际转动惯量之间的误差，将会导致对电机的控制不准确的技术问题。

发明内容

本申请实施例提供了一种获得电机转动惯量的方法和装置，用于获得电机的转动惯量。

第一方面，本申请提供了一种获得电机转动惯量的方法，包括：

获得电机在第n-1周期的第n-1实际转速；n为正整数；

获得与所述电机对应的可调模型输出的第n-1推定转速；所述第n-1推定转速为所述可调模型基于第n-1轮计算结果输出的转速；

根据所述第n-1实际转速和所述第n-1推定转速，获得第n轮计算结果；所述第n轮计算结果中包括第n转动惯量相关量；

判断所述第n转动惯量相关量和所述第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系；

当所述第n转动惯量相关量和所述第n-1转动惯量相关量满足所述预设关系时，确定所述第n转动惯量相关量对应的第n转动惯量为所述电机的实际转动惯量。

可选的，在判断所述第n转动惯量相关量和所述第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系之后，还包括：

当所述第n转动惯量相关量和所述第n-1转动惯量相关量不满足所述预设关系时，将所述第n轮计算结果输入所述可调模型，以使所述可调模型基于所述第n轮计算结果在第n+1周期输出与所述第n-1推定转速不同的第n推定转速。

可选的，根据所述第n-1实际转速和所述第n-1推定转速，获得第n轮计算结果，包括：

根据下列公式获得所述第n轮计算结果：

$\left\{\begin{array}{c}{a^{\′}}_n=a0-\underset{i=0}{\overset{n-1}{\Σ}}\left(k_1{{\mathrm{\Ω}}^{\′}}_i\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\\ {b^{\′}}_n=b0+\underset{i=0}{\overset{n-1}{\Σ}}\left(k_2{T}_{ei}\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\\ {c^{\′}}_n=c0-\underset{i=0}{\overset{n-1}{\Σ}}\left(k_3\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\end{array}\right.$

其中，B_mn'为第n粘滞摩擦系数，J_n'为第n转动惯量，T_ln'为第n负载转矩，B_m0为所述可调模型的初始粘滞摩擦系数，J₀为所述可调模型的初始转动惯量，T_l0为所述可调模型的初始负载转矩，Ω_i为第i实际转速，Ω'_i为第i推定转速，k₁k₂k₃≠0，T_ei为第i转矩，△T为周期。

可选的，当所述第n轮计算结果至少包括J_n'、B_mn'和T_ln'，所述第n转动惯量相关量为J_n'时，判断所述第n转动惯量相关量和所述第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系，包括：

判断第n-1转动惯量J_n-1'与J_n'的比值是否在预设范围内；当J_n-1'与J_n'的所述比值在所述预设范围内时，表示所述第n转动惯量相关量和所述第n-1转动惯量相关量满足所述预设关系；其中，J_n-1'为所述第n-1转动惯量相关量。

可选的，当所述第n轮计算结果至少包括a_n'、b_n'和c_n'，所述第n转动惯量相关量为b_n'时，判断所述第n转动惯量相关量和所述第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系，包括：

判断b_n-1'与b_n'之差的绝对值是否小于阈值；当所述绝对值小于所述阈值时，表示所述第n转动惯量相关量和所述第n-1转动惯量相关量满足所述预设关系；其中，b_n-1'为所述第n-1转动惯量相关量。

可选的，所述可调模型按照下列公式，基于所述第n轮计算结果在第n+1周期输出与所述第n-1推定转速不同的第n推定转速：

${J_n}^{\′}\frac{{{\mathrm{\Ω}}^{\′}}_n-{{\mathrm{\Ω}}^{\′}}_{n-1}}{\mathrm{\Δ}T}=T_{en}-{T_{ln}}^{\′}-{B_{mn}}^{\′}{{\mathrm{\Ω}}^{\′}}_n$

其中，Ω_n'为所述第n推定转速，Ω_n-1'为所述第n-1推定转速，T_en为第n转矩。

另一方面，本申请提供了一种获得电机转动惯量的装置，包括：

实际转速获得单元，用于获得电机在第n-1周期的第n-1实际转速；n为正整数；

推定转速获得单元，用于获得与所述电机对应的可调模型输出的第n-1推定转速；所述第n-1推定转速为所述可调模型基于第n-1轮计算结果输出的转速；

计算单元，用于根据所述第n-1实际转速和所述第n-1推定转速，获得第n轮计算结果；所述第n轮计算结果中包括第n转动惯量相关量；

判断单元，用于判断所述第n转动惯量相关量和所述第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系；

确定单元，用于当所述第n转动惯量相关量和所述第n-1转动惯量相关量满足所述预设关系时，确定所述第n转动惯量相关量对应的第n转动惯量为所述电机的实际转动惯量。

可选的，所述装置还包括：

输入单元，用于在判断所述第n转动惯量相关量和所述第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系之后，当所述第n转动惯量相关量和所述第n-1转动惯量相关量不满足所述预设关系时，将所述第n轮计算结果输入所述可调模型，以使所述可调模型基于所述第n轮计算结果在第n+1周期输出与所述第n-1推定转速不同的第n推定转速。

可选的，所述计算单元用于根据下列公式获得所述第n轮计算结果：

$\left\{\begin{array}{c}{a^{\′}}_n=a0-\underset{i=0}{\overset{n-1}{\Σ}}\left(k_1{{\mathrm{\Ω}}^{\′}}_i\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\\ {b^{\′}}_n=b0+\underset{i=0}{\overset{n-1}{\Σ}}\left(k_2{T}_{ei}\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\\ {c^{\′}}_n=c0-\underset{i=0}{\overset{n-1}{\Σ}}\left(k_3\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\end{array}\right.$

其中，B_mn'为第n粘滞摩擦系数，J_n'为第n转动惯量，T_ln'为第n负载转矩，B_m0为所述可调模型的初始粘滞摩擦系数，J₀为所述可调模型的初始转动惯量，T_l0为所述可调模型的初始负载转矩，Ω_i为第i实际转速，Ω'_i为第i推定转速，k₁k₂k₃≠0，T_ei为第i转矩，△T为周期。

可选的，当所述第n轮计算结果至少包括J_n'、B_mn'和T_ln'，所述第n转动惯量相关量为J_n'时，所述判断单元用于判断第n-1转动惯量J_n-1'与J_n'的比值是否在预设范围内；当J_n-1'与J_n'的所述比值在所述预设范围内时，表示所述第n转动惯量相关量和所述第n-1转动惯量相关量满足所述预设关系；其中，J_n-1'为所述第n-1转动惯量相关量。

可选的，当所述第n轮计算结果至少包括a_n'、b_n'和c_n'，所述第n转动惯量相关量为b_n'时，所述判断单元用于判断b_n-1'与b_n'之差的绝对值是否小于阈值；当所述绝对值小于所述阈值时，表示所述第n转动惯量相关量和所述第n-1转动惯量相关量满足所述预设关系；其中，b_n-1'为所述第n-1转动惯量相关量。

可选的，所述可调模型按照下列公式，基于所述第n轮计算结果在第n+1周期输出与所述第n-1推定转速不同的第n推定转速：

${J_n}^{\′}\frac{{{\mathrm{\Ω}}^{\′}}_n-{{\mathrm{\Ω}}^{\′}}_{n-1}}{\mathrm{\Δ}T}=T_{en}-{T_{ln}}^{\′}-{B_{mn}}^{\′}{{\mathrm{\Ω}}^{\′}}_n$

其中，Ω_n'为所述第n推定转速，Ω_n-1'为所述第n-1推定转速，T_en为第n转矩。

本申请实施例中的上述一个或多个技术方案，至少具有如下一种或多种技术效果：

在本申请实施例的技术方案中，首先获得电机在第n-1周期的第n-1实际转速，以及与电机对应的可调模型输出的第n-1推定转速。其中，第n-1推定转速为可调模型基于第n-1轮计算结果输出的转速；n为正整数。然后，根据第n-1实际转速和第n-1推定转速，获得第n轮计算结果，第n轮计算结果中包括第n转动惯量相关量。接下来，判断第n转动惯量相关量和第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系，当满足预设关系时，确定第n转动惯量相关量对应的第n转动惯量为电机的实际转动惯量。所以，通过本申请实施例中的技术方案就获得了电机的实际转动惯量，那么，根据实际转动惯量对电机进行控制就会更加准确。

附图说明

图1为本申请实施例中获得电机转动惯量的方法流程图；

图2为本申请实施例中获得电机转动惯量的模型框架的示意图；

图3为本申请实施例中获得电机转动惯量的装置结构示意图。

具体实施方式

本申请实施例提供了一种获得电机转动惯量的方法和装置，用于获得电机的转动惯量。

本申请提供的技术方案总体思路如下：

在本申请实施例的技术方案中，首先获得电机在第n-1周期的第n-1实际转速，以及与电机对应的可调模型输出的第n-1推定转速。其中，第n-1推定转速为可调模型基于第n-1轮计算结果输出的转速；n为正整数。然后，根据第n-1实际转速和第n-1推定转速，获得第n轮计算结果，第n轮计算结果中包括第n转动惯量相关量。接下来，判断第n转动惯量相关量和第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系，当满足预设关系时，确定第n转动惯量相关量对应的第n转动惯量为电机的实际转动惯量。所以，通过本申请实施例中的技术方案就获得了电机的实际转动惯量，那么，根据实际转动惯量对电机进行控制就会更加准确。

下面通过附图以及具体实施例对本发明技术方案做详细的说明，应当理解本申请实施例以及实施例中的具体特征是对本申请技术方案的详细的说明，而不是对本申请技术方案的限定，在不冲突的情况下，本申请实施例以及实施例中的技术特征可以相互组合。

本文中术语“和/或”，仅仅是一种描述关联对象的关联关系，表示可以存在三种关系，例如，A和/或B，可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。另外，本文中字符“/”，一般表示前后关联对象是一种“或”的关系。

本申请第一方面提供了一种获得电机转动惯量的方法，请参考图1，包括如下步骤：

S101：获得电机在第n-1周期的第n-1实际转速。

S102：获得与所述电机对应的可调模型输出的第n-1推定转速。

S103：根据所述第n-1实际转速和所述第n-1推定转速，获得第n轮计算结果。

S104：判断所述第n转动惯量相关量和所述第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系。

S105：当所述第n转动惯量相关量和所述第n-1转动惯量相关量满足所述预设关系时，确定所述第n转动惯量相关量对应的第n转动惯量为所述电机的实际转动惯量。

在S101中，获得电机在第n-1周期的第n-1实际转速，具体为通过转速传感器检测或者控制软件获得电机在第n-1周期的实际转速。在本申请实施例中，周期是指每轮计算的周期，周期长度△T可以与脉冲调制信号的周期相一致，例如1.4×10^-4s到2.5×10^-4s中的任意长度。

在S102中，获得与电机对应的可调模型输出的第n-1推定转速。具体来讲，电机在运行中，可以得到电机的实际模型，即电机传递函数。在本申请实施例中，可调模型的传递函数为与电机的实际模型形式一致，而参数可调的函数。所以，可调模型是与电机对应的。推定转速是根据可调模型的函数和输入而得到的转速。举例来说，例如电机的实际模型为其中，J为电机的实际转动惯量，Ω为输出的实际转速，T_e为输入的转矩命令值，T_l为实际负载转矩，B_m为实际的粘滞摩擦系数。那么可调模型为与电机实际模型形式一致的其中，J'为可调模型的转动惯量，Ω'为可调模型输出的转速，T_e为输入可调模型的转矩，与输入电机的转矩命令值相同，T_l'为可调模型的负载转矩，B_m'为可调模型的粘滞摩擦系数。

第n-1推定转速为可调模型基于上一轮计算结果，即第n-1轮计算结果输出的转速。具体来讲，可调模型基于第n-1轮计算结果调整参数，例如可调模型的转动惯量，可调模型的负载转矩和/或可调模型的粘滞摩擦系数，然后根据输入而输出第n-1推定转速。

在具体实现过程中，为获得电机的实际转动惯量而进行第n轮计算时，S101可以先于S102执行，S102也可以先于S101执行，还可以同时执行S101和S102，本申请不做具体限制。

获得第n-1实际转速和第n-1推定转速后，S103中，基于第n-1实际转速和第n-1推定转速进行第n轮计算，进而获得第n轮计算结果。如图2所示，在本申请实施例中，每一轮计算，从电机获得的实际转速Ω和可调模型输出的推定转速Ω'均输入到计算模块中进行计算。计算模块按照预设的算法处理，进而输出本轮计算结果。

具体来讲，本申请实施例中的第n轮计算结果中包括第n转动惯量相关量。在本申请实施例中，第n转动惯量相关量可以为第n转动惯量本身，也可以为第n转动惯量的应变量，例如第n转动惯量的倒数值或第n转动惯量的3倍值，本申请不做具体限制。

获得第n转动惯量相关量后，在S104中，判断第n转动惯量相关量和第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系。在本申请实施例中，预设关系表示可调模型的参数与电机实际模型的参数一致或接近。进一步，当第n转动惯量相关量和第n-1转动惯量相关量满足预设关系时，S105中确定第n转动惯量相关量对应的第n转动惯量为电机的实际转动惯量。

具体来讲，由于可调模型的函数形式与电机实际模型的形式是一致的，而电机实际模型又表示电机的实际状态，所以，当第n转动惯量相关量和第n-1转动惯量相关量满足预设关系时，表示可调模型的参数已调整至与电机实际模型的参数一致或接近，进而表示此时的可调模型与电机实际模型是相同或接近相同。因此，此时可调模型的第n转动惯量相关量对应的第n转动惯量就与电机的实际转动惯量相同或接近，故可以作为电机的实际转动惯量。

所以，由上述描述可知，当第n转动惯量相关量和第n-1转动惯量相关量满足预设关系时，可调模型相当于电机实际模型，可调模型的第n转动惯量就相当于电机的实际转动惯量，所以，将第n转动惯量作为实际转动惯量，就获得了电机的实际转动惯量。

在本申请另一实施例中，在S104之后，还可以包括：

当所述第n转动惯量相关量和所述第n-1转动惯量相关量不满足所述预设关系时，将所述第n轮计算结果输入所述可调模型，以使所述可调模型基于所述第n轮计算结果在第n+1周期输出与所述第n-1推定转速不同的第n推定转速。

具体来讲，当S104中判断出第n转动惯量相关量和第n-1转动惯量相关量不满足预设关系时，表示此时可调模型的参数还没有调整到与电机一致或者接近一致，所以，可调模型仍然需要继续调整参数。所以，在本申请实施例中，当第n转动惯量相关量和第n-1转动惯量相关量不满足预设关系时，进一步将计算模块输出的第n轮计算结果输入可调模型，进而使可调模型调整参数，从而在第n+1周期输出与第n-1推定转速不同的第n推定转速。

在本申请实施例中，如果第n转动惯量相关量和第n-1转动惯量相关量不满足预设关系时，将第n轮计算结果输入可调模型，然后再次执行S101至S104，直到第n转动惯量相关量和第n-1转动惯量相关量满足预设关系，确定出电机的实际转动惯量才结束。

接下来，对如何获得第n轮计算结果进行详细介绍。

在本申请实施例中，根据下列公式(1)获得第n轮计算结果：

$\left\{\begin{array}{c}{a^{\′}}_n=a0-\underset{i=0}{\overset{n-1}{\Σ}}\left(k_1{{\mathrm{\Ω}}^{\′}}_i\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\\ {b^{\′}}_n=b0+\underset{i=0}{\overset{n-1}{\Σ}}\left(k_2{T}_{ei}\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\\ {c^{\′}}_n=c0-\underset{i=0}{\overset{n-1}{\Σ}}\left(k_3\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\end{array}\right.;$ 公式(1)

具体来讲，B_mn'为第n粘滞摩擦系数，J_n'为第n转动惯量，T_ln'为第n负载转矩。B_m0为可调模型的初始粘滞摩擦系数，J₀为可调模型的初始转动惯量，T_l0为可调模型的初始负载转矩。B_m0、J₀和T_l0由本领域技术人员任意设置。其中，由于J₀将作为分母，所以在设置时J₀具体可以设置为任意非零实数，本申请不做具体限制。其中Ω_i为电机在第i周期输出的第i实际转速，Ω'_i为可调模型在第i周期输出的第i推定转速。

在初始时刻，将a0,b0,c0作为初始值，或者将B_m0、J₀和T_l0作为初始值。当n＝1，即进行第1轮计算时，将B_m0'＝B_m0、J₀'＝J₀和T_l0'＝T_l0代入或者将a'₀＝a0,b'₀＝b0和c'₀＝c0代入获得可调模型的第0推定转速Ω'₀。检测电机输出的第0实际转速Ω₀，进而，将Ω'₀和Ω₀，代入公式(1)获得

${a^{\′}}_1=a0-\underset{i=0}{\overset{0}{\Σ}}\left(k_1\right({\mathrm{\Ω}}_0-{{\mathrm{\Ω}}^{\′}}_0\left){{\mathrm{\Ω}}^{\′}}_0\mathrm{\Δ}T\right)=a0-k_1({\mathrm{\Ω}}_0-{{\mathrm{\Ω}}^{\′}}_0){{\mathrm{\Ω}}^{\′}}_0\mathrm{\Δ}T.$

所以本申请实施例中的n从1开始取值。获得b'₁和c'₁的过程类似。

当n＝2，即进行第2轮计算时:将a'₁、b'₁和c'₁代入获得第1推定转速Ω'₁。并检测第1实际转速Ω₁，代入公式(1)获得

${a^{\′}}_2=a0-\underset{i=0}{\overset{1}{\Σ}}\left(k_1{{\mathrm{\Ω}}^{\′}}_0\right({\mathrm{\Ω}}_0-{{\mathrm{\Ω}}^{\′}}_0\left)\mathrm{\Δ}T\right)=a0-k_1{{\mathrm{\Ω}}^{\′}}_0({\mathrm{\Ω}}_0-{{\mathrm{\Ω}}^{\′}}_0)\mathrm{\Δ}T-k_1{{\mathrm{\Ω}}^{\′}}_1({\mathrm{\Ω}}_1-{{\mathrm{\Ω}}^{\′}}_1)\mathrm{\Δ}T.$

后续过程以此类推，这里就不一一赘述了。

在本申请实施例中，k₁具体为使a_n'收敛的任意值，k₂具体为使b_n'收敛的任意值，k₃具体为使c_n'收敛的任意值。同时，k₁k₂k₃≠0。在具体实现过程中，k₁、k₂和k₃可以相同也可以不完全相同。Ω_n-1为第n-1实际转速，Ω'_n-1为第n-1推定转速。T_ei为第i转矩，即第i周期输入可调模型和电机的转矩命令值，T_en-1为第n-1转矩，即第n-1周期输入电机和可调模型的转矩命令值。△T为周期。

在本申请实施例中，第n轮计算结果至少包括J_n'、B_mn'和T_ln'，或者至少包括a_n'、b_n'和c_n'。所以，当需要执行第n轮计算时，将Ω_n-1和Ω'_n-1输入计算模块，进而计算模块基于公式(1)执行计算，就可以获得J_n'、B_mn'和T_ln'，或者a_n'、b_n'和c_n'。

在本申请实施例中，电机实际模型为：进而可调模型为其中，J为电机的实际转动惯量，Ω为输出的实际转速，T_e为输入的转矩，B_m为实际粘滞摩擦系数，T_l为实际的负载转矩。J'为可调模型的转动惯量，Ω'为可调模型输出的转速，T_e为输入的转矩，B_m'为可调模型的粘滞摩擦系数，T_ln'为可调模型的负载转矩。

在具体实现过程中，基于第n-1实际转速和第n-1推定转速获得第n转动惯量的方法有多种，例如基于李雅普诺夫函数、波波夫不等式或最小二乘法等，本申请不做具体限制。在本申请实施例中，将以李雅普诺夫函数为例详细获得公式(1)的推导过程。

首先，令电机实际模型可以变形为：

$\stackrel{\·}{\mathrm{\Ω}}=-a\mathrm{\Ω}+{\mathrm{bT}}_e-c.$ 公式(2)

其中，公式(2)中的“·”表示为微分算符。类似地，可以将可调模型变形为：

${\stackrel{\·}{\mathrm{\Ω}}}^{\′}=-a^{\′}{\mathrm{\Ω}}^{\′}+b^{\′}{T}_e-c^{\′},$ 公式(3)

B_m'为待确定的可调模型的粘滞摩擦系数，J'为待确定的可调模型的转动惯量，T_l'为待确定的可调模型的负载转矩。

接下来，公式(2)-公式(3)，得

$\stackrel{\·}{e}=-ae-(a-a^{\′}){\mathrm{\Ω}}^{\′}+(b-b^{\′})T_e-(c-c^{\′}),$ 公式(4)

其中，e＝Ω-Ω'。

进一步，令a1＝a-a',a2＝b-b',a3＝c-c'，那么公式(4)可以简化为：

$\stackrel{\·}{e}=-ae-a1{\mathrm{\Ω}}^{\′}+a2T_e-a3.$ 公式(5)

再进一步，令则，A＝[-a]，φ^T＝[a1,a2,a3],s＝[-Ω',T_e,-1]^T。

由于李雅普诺夫函数形式为其中任意p和Γ为对称正定矩阵。所以，在本申请实施例中，

$V(\stackrel{\·}{e},t)=\frac{1}{2}(\stackrel{\·}{e^T}pe+e^{T}p\stackrel{\·}{e}+\stackrel{\·}{{\mathrm{\φ}}^T}\mathrm{\Γ}\mathrm{\φ}+{\mathrm{\φ}}^T\mathrm{\Γ}\stackrel{\·}{\mathrm{\φ}}).$ 公式(6)

由于 $\stackrel{\·}{e}=Ae+{\mathrm{\φ}}^{T}s,$ 所以 $\stackrel{\·}{e^T}pe={(Ae+{\mathrm{\φ}}^{T}s)}^{T}pe=e^T{A}^{T}pe+s^T\mathrm{\φ}e,$ 公式(7)

以及 $e^{T}p\stackrel{\·}{e}=e^{T}p(Ae+{\mathrm{\φ}}^{T}s)=e^{T}pAe+e^T{\mathrm{p\φ}}^{T}s.$ 公式(8)

那么结合公式(7)和公式(8)可得

$\stackrel{\·}{e^T}pe+e^{T}p\stackrel{\·}{e}=(e^T{A}^{T}pe+s^T\mathrm{\φ}pe)+(e^{T}pAe+e^T{\mathrm{p\φ}}^{T}s)=e^T(A^{T}p+pA)e+(s^T\mathrm{\φ}pe+e^T{\mathrm{p\φ}}^{T}s).$

公式(9)

而对于公式(9)，由于φ^T＝[a1,a2,a3],s＝[-Ω',T_e,-1]^T，所以

s^Tφpe+e^Tpφ^Ts＝[-Ω',T_e,-1][a1,a2,a3]^Tpe+e^Tp[a1,a2,a3][-Ω',T_e,-1]^T。

公式(10)

＝(-a1Ω'+a2T_e-a3)pe+e^Tp(-a1Ω'+a2T_e-a3)

在本申请实施例中，选择p为单位矩阵，即p＝[1]，那么

s^Tφpe+e^Tpφ^Ts＝2(-a1Ω'+a2T_e-a3)e。公式(11)

另外，令本申请实施例中的 $\mathrm{\Γ}=\left[\begin{array}{c}{g}_1,0,0\\ 0,g_2,0\\ 0,0,g_3\end{array}\right],$ g₁、g₂和g₃为正数，所以

$\stackrel{\·}{{\mathrm{\φ}}^T}\mathrm{\Γ}\mathrm{\φ}=\[\stackrel{\·}{a1},\stackrel{\·}{a2},\stackrel{\·}{a3}\]\left[\begin{array}{c}{g}_1,0,0\\ 0,g_2,0\\ 0,0,g_3\end{array}\right]{\[a1,a2,a3\]}^T=g_{1}a1\stackrel{\·}{a1}+g_{2}a2\stackrel{\·}{a2}+g_{3}a3\stackrel{\·}{a3}.$ 公式(12)

同理可得 ${\mathrm{\φ}}^T\mathrm{\Γ}\stackrel{\·}{\mathrm{\φ}}=g_{1}a1\stackrel{\·}{a1}+g_{2}a2\stackrel{\·}{a2}+g_{3}a3\stackrel{\·}{a3}.$ 公式(13)

所以，最终公式(6)中的李雅普诺夫函数为

$V(\stackrel{\·}{e},t)=\frac{1}{2}{e}^T(A^{T}p+pA)e+(-a1{\mathrm{\Ω}}^{\′}+a2T_e-a3)e+(g_1\stackrel{\·}{a1}a1+g_2\stackrel{\·}{a2}a2+g_3\stackrel{\·}{a3}a3).$ 公式(14)

而根据李雅普诺夫稳定条件可知，要使李雅普诺夫函数稳定，需负定。而A^Tp+pA＝-2a，所以A^Tp+pA负定。

再次整理公式(14)得到

$V(\stackrel{\·}{e},t)=\frac{1}{2}{e}^T(A^{T}p+pA)e+a1\[g_1\stackrel{\·}{a1}-{\mathrm{\Ω}}^{\′}e\]+a2\[g_2\stackrel{\·}{a2}+T_{e}e\]+a3\[g_3\stackrel{\·}{a3}-e\].$ 公式(15)

由此可见，当

$\left\{\begin{array}{c}{g}_1\stackrel{\·}{a1}-{\mathrm{\Ω}}^{\′}e=0\\ g_2\stackrel{\·}{a2}+T_{e}e=0\\ g_3\stackrel{\·}{a3}-e=0\end{array}\right.$ 公式(16)

成立时，负定。

而a1＝a-a',a2＝b-b',a3＝c-c'中，J，B_m和T_l都是实际参数，是定值，所以a，b和c的导数均为0。进而公式(16)进一步可化简为

$\left\{\begin{array}{c}-g_1\stackrel{\·}{a^{\′}}-{\mathrm{\Ω}}^{\′}e=0\\ -g_2\stackrel{\·}{b^{\′}+T_{e}e=0}\\ -g\stackrel{\·}{c^{\′}}-e=0\end{array}\right..$ 公式(17)

故可得到，

$\left\{\begin{array}{c}{a}^{\′}=a0-{\∫}_0^t{k}_1{\mathrm{\Ω}}^{\′}(\mathrm{\Ω}-{\mathrm{\Ω}}^{\′})dt\\ b^{\′}=b0+{\∫}_0^t{k}_2{T}_e(\mathrm{\Ω}-{\mathrm{\Ω}}^{\′})dt\\ c^{\′}=c0-{\∫}_0^t{k}_3(\mathrm{\Ω}-{\mathrm{\Ω}}^{\′})dt\end{array}\right..$ 公式(19)

其中，t表示进行计算时刻与初始时刻直接的时间差，具体为t＝n△T。

最后，为便于计算机处理，将上述公式(19)离散化，得到

$\left\{\begin{array}{c}{a^{\′}}_n=a0-\underset{i=0}{\overset{n-1}{\Σ}}\left(k_1{{\mathrm{\Ω}}^{\′}}_i\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\\ {b^{\′}}_n=b0+\underset{i=0}{\overset{n-1}{\Σ}}\left(k_2{T}_{ei}\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\\ {c^{\′}}_n=c0-\underset{i=0}{\overset{n-1}{\Σ}}\left(k_3\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\end{array}\right.;$ 公式(1)

在本申请实施例中，第n轮计算结果有两种可能。第一种，至少将J_n'、B_mn'和T_ln'作为第n轮计算结果。在具体实现过程中，第n轮计算结果进一步还可以包括a_n'、b_n'和/或c_n'等，本申请不做具体限制。第二种，至少将a_n'、b_n'和c_n'作为第n轮计算结果。在具体实现过程中，第n轮计算结果进一步还可以包括J_n'、B_mn'和/或T_ln'。

由于第n轮计算结果有两种可能，所以S104中判断第n转动惯量相关量与第n-1转动惯量相关量是否满足预设关系，也有两种具体实现方式。

第一种：

当第n轮计算结果至少包括J_n'、B_mn'和T_ln'，且第n转动惯量相关量为J_n'时，S104具体包括：

判断第n-1转动惯量J_n-1'与J_n'的比值是否在预设范围内。

当第n转动惯量相关量为第n转动惯量J_n'本身时，获得第n-1轮计算结果中的第n-1转动惯量J_n-1'。在第一种实现方式中，J_n-1'就是第n-1转动惯量相关量。然后进一步获得

在第一种实现方式中，判断第n转动惯量相关量与第n-1转动惯量相关量是否满足预设关系，具体为判断是否在预设范围内。在本申请实施例中，预设范围为1附近的小范围，例如[0.995，1.005]，或[0.998，1.002]等。在具体实现过程中，预设范围越小，最终作为实际转动惯量的第n转动惯量越接近实际转动惯量。本申请所属领域的普通技术人员可以根据实际进行设置，本申请不做具体限制。

具体来讲，当在预设范围内时，表示等于1或者接近于1。而等于1或者接近于1，表明Ω'_n-1等于或接近于Ω_n-1，进而表明此时可调模型的参数等于或接近于电机实际模型，所以可以将J_n'视为电机实际转动惯量J。所以，当在预设范围内时，表示第n转动惯量相关量与第n-1转动惯量相关量满足预设关系。

进一步，在第一种实现方式中，由于第n转动惯量相关量就是第n转动惯量本身，所以在S105中直接确定J_n'为实际转动惯量。

第二种：

当第n轮计算结果至少包括a_n'、b_n'和c_n'，且第n转动惯量相关量为b_n'时，S104具体包括：

判断b_n-1'与b_n'之差的绝对值是否小于阈值。

当第n转动惯量相关量为b_n'时，获得第n-1轮计算结果中的第n-1转动惯量b_n-1'。在第二种实现方式中，b_n-1'就是第n-1转动惯量相关量。然后进一步获得|b'_n-b'_n-1|或|b'_n-1-b'_n|。

在第二种实现方式中，判断第n转动惯量相关量与第n-1转动惯量相关量是否满足预设关系，具体为判断|b'_n-b'_n-1|或|b'_n-1-b'_n|是否小于阈值。在本申请实施例中，阈值为0或接近于0的数，例如0，0.05，或0.1等。在具体实现过程中，阈值越小，越接近0，最终作为实际转动惯量的第n转动惯量越接近实际转动惯量。本申请所属领域的普通技术人员可以根据实际进行设置，本申请不做具体限制。

具体来讲，|b'_n-b'_n-1|或|b'_n-1-b'_n|小于阈值时，表示|b'_n-b'_n-1|或|b'_n-1-b'_n|等于0或者接近于0。而|b'_n-b'_n-1|或|b'_n-1-b'_n|等于0或者接近于0，表明Ω'_n-1等于或接近于Ω_n-1，进而表明此时可调模型的参数等于或接近于电机实际模型，所以可以将J_n'视为电机实际转动惯量J。所以，当|b'_n-b'_n-1|或|b'_n-1-b'_n|小于阈值时，表示第n转动惯量相关量与第n-1转动惯量相关量满足预设关系。

进一步，由于所以在第二种实现方式中，S105中确定的第n转动惯量相关量对应的第n转动惯量就是

进一步，可调模型为所以按照公式(20)输出第n-1推定转速。

${J_{n-1}}^{\′}\frac{{{\mathrm{d\Ω}}^{\′}}_{n-1}}{dt}=T_{en-1}-{T_{ln-1}}^{\′}-{B_{mn-1}}^{\′}{{\mathrm{\Ω}}^{\′}}_{n-1}$ 公式(20)

其中，J_n-1'为第n-1轮计算出的第n-1转动惯量，T_ln-1'为第n-1负载转矩，B_mn-1'为第n-1粘滞摩擦系数，T_en-1为第n-1转矩。

以及，可调模型按照公式(21)输出第n推定转速。

${J_n}^{\′}\frac{{{\mathrm{d\Ω}}^{\′}}_n}{dt}=T_{en}-{T_{ln}}^{\′}-{B_{mn}}^{\′}{{\mathrm{\Ω}}^{\′}}_n$ 公式(21)

其中，Ω_n'为第n推定转速，T_en为第n转矩，即第n周期输入电机和可调模型的转矩值。

进一步，在具体实现过程中，若利用计算机获得推定转速，则计算机具体为利用将离散化后的公式进行计算。例如利用公式(21)离散化后的公式获得Ω_n'，具体为

${J_n}^{\′}\frac{{{\mathrm{\Ω}}^{\′}}_n-{{\mathrm{\Ω}}^{\′}}_{n-1}}{\mathrm{\Δ}T}=T_{en}-{T_{ln}}^{\′}-{B_{mn}}^{\′}{{\mathrm{\Ω}}^{\′}}_n.$ 公式(22)

基于与前述实施例中获得电机转动惯量同样的发明构思，本申请第二方面还提供一种获得电机转动惯量的装置，如图3所示，包括：

实际转速获得单元301，用于获得电机在第n-1周期的第n-1实际转速；n为正整数；

推定转速获得单元302，用于获得与电机对应的可调模型输出的第n-1推定转速；第n-1推定转速为可调模型基于第n-1轮计算结果输出的转速；

计算单元303，用于根据第n-1实际转速和第n-1推定转速，获得第n轮计算结果；第n轮计算结果中包括第n转动惯量相关量；

判断单元304，用于判断第n转动惯量相关量和第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系；

确定单元305，用于当第n转动惯量相关量和第n-1转动惯量相关量满足预设关系时，确定第n转动惯量相关量对应的第n转动惯量为电机的实际转动惯量。

进一步，本申请实施例中的装置还包括：

输入单元，用于在判断第n转动惯量相关量和第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系之后，当第n转动惯量相关量和第n-1转动惯量相关量不满足预设关系时，将第n轮计算结果输入可调模型，以使可调模型基于第n轮计算结果在第n+1周期输出与第n-1推定转速不同的第n推定转速。

具体来讲，计算单元303用于根据下列公式获得第n轮计算结果：

$\left\{\begin{array}{c}{a^{\′}}_n=a0-\underset{i=0}{\overset{n-1}{\Σ}}\left(k_1{{\mathrm{\Ω}}^{\′}}_i\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\\ {b^{\′}}_n=b0+\underset{i=0}{\overset{n-1}{\Σ}}\left(k_2{T}_{ei}\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\\ {c^{\′}}_n=c0-\underset{i=0}{\overset{n-1}{\Σ}}\left(k_3\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\end{array}\right.$

其中，B_mn'为第n粘滞摩擦系数，J_n'为第n转动惯量，T_ln'为第n负载转矩，B_m0为可调模型的初始粘滞摩擦系数，J₀为可调模型的初始转动惯量，T_l0为可调模型的初始负载转矩，Ω_i为第i实际转速，Ω'_i为第i推定转速，k₁k₂k₃≠0，T_ei为第i转矩，△T为周期。

而当第n轮计算结果至少包括a_n'、b_n'和c_n'，第n转动惯量相关量为b_n'时，判断单元304用于判断b_n-1'与b_n'之差的绝对值是否小于阈值；当绝对值小于阈值时，表示第n转动惯量相关量和第n-1转动惯量相关量满足预设关系；其中，b_n-1'为第n-1转动惯量相关量。

具体来讲，可调模型按照下列公式，基于第n轮计算结果在第n+1周期输出与第n-1推定转速不同的第n推定转速：

${J_n}^{\′}\frac{{{\mathrm{d\Ω}}_n}^{\′}}{dt}=T_{en}-{T_{ln}}^{\′}-{B_{mn}}^{\′}{{\mathrm{\Ω}}_n}^{\′}$

其中，Ω_n'为所述第n推定转速，Ω_n-1'为所述第n-1推定转速，T_en为第n转矩。

前述图1-图2实施例中的获得电机转动惯量的方法各种变化方式和具体实例同样适用于本实施例的获得电机转动惯量的装置，通过前述对获得电机转动惯量的方法的详细描述，本领域技术人员可以清楚的知道本实施例中获得电机转动惯量的装置的实施方法，所以为了说明书的简洁，在此不再详述。

本申请实施例中的上述一个或多个技术方案，至少具有如下一种或多种技术效果：

在本申请实施例的技术方案中，首先获得电机在第n-1周期的第n-1实际转速，以及与电机对应的可调模型输出的第n-1推定转速。其中，第n-1推定转速为可调模型基于第n-1轮计算结果输出的转速；n为正整数。然后，根据第n-1实际转速和第n-1推定转速，获得第n轮计算结果，第n轮计算结果中包括第n转动惯量相关量。接下来，判断第n转动惯量相关量和第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系，当满足预设关系时，确定第n转动惯量相关量对应的第n转动惯量为电机的实际转动惯量。所以，通过本申请实施例中的技术方案就获得了电机的实际转动惯量，那么，根据实际转动惯量对电机进行控制就会更加准确。

显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。

Claims (12)

1.一种获得电机转动惯量的方法，其特征在于，包括：

获得电机在第n-1周期的第n-1实际转速；n为正整数；

获得与所述电机对应的可调模型输出的第n-1推定转速；所述第n-1推定转速为所述可调模型基于第n-1轮计算结果输出的转速；

根据所述第n-1实际转速和所述第n-1推定转速，获得第n轮计算结果；所述第n轮计算结果中包括第n转动惯量相关量；

判断所述第n转动惯量相关量和所述第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系；

当所述第n转动惯量相关量和所述第n-1转动惯量相关量满足所述预设关系时，确定所述第n转动惯量相关量对应的第n转动惯量为所述电机的实际转动惯量。

2.如权利要求1所述的方法，其特征在于，在判断所述第n转动惯量相关量和所述第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系之后，还包括：

当所述第n转动惯量相关量和所述第n-1转动惯量相关量不满足所述预设关系时，将所述第n轮计算结果输入所述可调模型，以使所述可调模型基于所述第n轮计算结果在第n+1周期输出与所述第n-1推定转速不同的第n推定转速。

3.如权利要求2所述的方法，其特征在于，根据所述第n-1实际转速和所述第n-1推定转速，获得第n轮计算结果，包括：

根据下列公式获得所述第n轮计算结果：

$\left\{\begin{array}{c}{a^{\′}}_n=a0-\underset{i=0}{\overset{n-1}{\mathrm{\Σ}}}\left(k_1{{\mathrm{\Ω}}^{\′}}_i\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\\ {b^{\′}}_n=b0+\underset{i=0}{\overset{n-1}{\mathrm{\Σ}}}\left(k_2{T}_{ei}\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\\ {c^{\′}}_n=c0-\underset{i=0}{\overset{n-1}{\mathrm{\Σ}}}\left(k_3\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\end{array}\right.$

其中，B_mn'为第n粘滞摩擦系数，J_n'为第n转动惯量，T_ln'为第n负载转矩，B_m0为所述可调模型的初始粘滞摩擦系数，J₀为所述可调模型的初始转动惯量，T_l0为所述可调模型的初始负载转矩，Ω_i为第i实际转速，Ω'_i为第i推定转速，k₁k₂k₃≠0，T_ei为第i转矩，ΔT为周期。

4.如权利要求3所述的方法，其特征在于，当所述第n轮计算结果至少包括J_n'、B_mn'和T_ln'，所述第n转动惯量相关量为J_n'时，判断所述第n转动惯量相关量和所述第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系，包括：

判断第n-1转动惯量J_n-1'与J_n'的比值是否在预设范围内；当J_n-1'与J_n'的所述比值在所述预设范围内时，表示所述第n转动惯量相关量和所述第n-1转动惯量相关量满足所述预设关系；其中，J_n-1'为所述第n-1转动惯量相关量。

5.如权利要求3所述的方法，其特征在于，当所述第n轮计算结果至少包括a_n'、b_n'和c_n'，所述第n转动惯量相关量为b_n'时，判断所述第n转动惯量相关量和所述第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系，包括：

判断b_n-1'与b_n'之差的绝对值是否小于阈值；当所述绝对值小于所述阈值时，表示所述第n转动惯量相关量和所述第n-1转动惯量相关量满足所述预设关系；其中，b_n-1'为所述第n-1转动惯量相关量。

6.如权利要求3所述的方法，其特征在于，所述可调模型按照下列公式，基于所述第n轮计算结果在第n+1周期输出与所述第n-1推定转速不同的第n推定转速：

${J_n}^{\′}\frac{{{\mathrm{\Ω}}^{\′}}_n-{{\mathrm{\Ω}}^{\′}}_{n-1}}{\mathrm{\Δ}T}=T_{en}-{T_{ln}}^{\′}-{B_{mn}}^{\′}{{\mathrm{\Ω}}^{\′}}_n$

其中，Ω_n'为所述第n推定转速，Ω_n-1'为所述第n-1推定转速，T_en为第n转矩。

7.一种获得电机转动惯量的装置，其特征在于，包括：

实际转速获得单元，用于获得电机在第n-1周期的第n-1实际转速；n为正整数；

推定转速获得单元，用于获得与所述电机对应的可调模型输出的第n-1推定转速；所述第n-1推定转速为所述可调模型基于第n-1轮计算结果输出的转速；

计算单元，用于根据所述第n-1实际转速和所述第n-1推定转速，获得第n轮计算结果；所述第n轮计算结果中包括第n转动惯量相关量；

判断单元，用于判断所述第n转动惯量相关量和所述第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系；

确定单元，用于当所述第n转动惯量相关量和所述第n-1转动惯量相关量满足所述预设关系时，确定所述第n转动惯量相关量对应的第n转动惯量为所述电机的实际转动惯量。

8.如权利要求7所述的装置，其特征在于，所述装置还包括：

输入单元，用于在判断所述第n转动惯量相关量和所述第n-1轮计算结果中的第n-1转动惯量相关量是否满足预设关系之后，当所述第n转动惯量相关量和所述第n-1转动惯量相关量不满足所述预设关系时，将所述第n轮计算结果输入所述可调模型，以使所述可调模型基于所述第n轮计算结果在第n+1周期输出与所述第n-1推定转速不同的第n推定转速。

9.如权利要求8所述的装置，其特征在于，所述计算单元用于根据下列公式获得所述第n轮计算结果：

$\left\{\begin{array}{c}{a^{\′}}_n=a0-\underset{i=0}{\overset{n-1}{\mathrm{\Σ}}}\left(k_1{{\mathrm{\Ω}}^{\′}}_i\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\\ {b^{\′}}_n=b0+\underset{i=0}{\overset{n-1}{\mathrm{\Σ}}}\left(k_2{T}_{ei}\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\\ {c^{\′}}_n=c0-\underset{i=0}{\overset{n-1}{\mathrm{\Σ}}}\left(k_3\right({\mathrm{\Ω}}_i-{{\mathrm{\Ω}}^{\′}}_i\left)\mathrm{\Δ}T\right)\end{array}\right.$

其中，B_mn'为第n粘滞摩擦系数，J_n'为第n转动惯量，T_ln'为第n负载转矩，B_m0为所述可调模型的初始粘滞摩擦系数，J₀为所述可调模型的初始转动惯量，T_l0为所述可调模型的初始负载转矩，Ω_i为第i实际转速，Ω'_i为第i推定转速，k₁k₂k₃≠0，T_ei为第i转矩，ΔT为周期。

10.如权利要求9所述的装置，其特征在于，当所述第n轮计算结果至少包括J_n'、B_mn'和T_ln'，所述第n转动惯量相关量为J_n'时，所述判断单元用于判断第n-1转动惯量J_n-1'与J_n'的比值是否在预设范围内；当J_n-1'与J_n'的所述比值在所述预设范围内时，表示所述第n转动惯量相关量和所述第n-1转动惯量相关量满足所述预设关系；其中，J_n-1'为所述第n-1转动惯量相关量。

11.如权利要求9所述的装置，其特征在于，当所述第n轮计算结果至少包括a_n'、b_n'和c_n'，所述第n转动惯量相关量为b_n'时，所述判断单元用于判断b_n-1'与b_n'之差的绝对值是否小于阈值；当所述绝对值小于所述阈值时，表示所述第n转动惯量相关量和所述第n-1转动惯量相关量满足所述预设关系；其中，b_n-1'为所述第n-1转动惯量相关量。

12.如权利要求9所述的装置，其特征在于，所述可调模型按照下列公式，基于所述第n轮计算结果在第n+1周期输出与所述第n-1推定转速不同的第n推定转速：

${J_n}^{\′}\frac{{{\mathrm{\Ω}}^{\′}}_n-{{\mathrm{\Ω}}^{\′}}_{n-1}}{\mathrm{\Δ}T}=T_{en}-{T_{ln}}^{\′}-{B_{mn}}^{\′}{{\mathrm{\Ω}}^{\′}}_n$

其中，Ω_n'为所述第n推定转速，Ω_n-1'为所述第n-1推定转速，T_en为第n转矩。