CN112100569A - 基于频域分析的电机参数追踪方法、装置、设备和介质 - Google Patents
基于频域分析的电机参数追踪方法、装置、设备和介质 Download PDFInfo
- Publication number
- CN112100569A CN112100569A CN202010856065.1A CN202010856065A CN112100569A CN 112100569 A CN112100569 A CN 112100569A CN 202010856065 A CN202010856065 A CN 202010856065A CN 112100569 A CN112100569 A CN 112100569A
- Authority
- CN
- China
- Prior art keywords
- data
- spectrum information
- current
- resistance
- calculating
- 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
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000004458 analytical method Methods 0.000 title claims abstract description 30
- 238000001228 spectrum Methods 0.000 claims abstract description 128
- 230000003068 static effect Effects 0.000 claims abstract description 77
- 230000006698 induction Effects 0.000 claims description 21
- 230000005520 electrodynamics Effects 0.000 claims description 19
- 230000008878 coupling Effects 0.000 claims description 17
- 238000010168 coupling process Methods 0.000 claims description 17
- 238000005859 coupling reaction Methods 0.000 claims description 17
- 230000001965 increasing effect Effects 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000013017 mechanical damping Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 10
- 230000001939 inductive effect Effects 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 16
- 230000000694 effects Effects 0.000 description 6
- 238000010223 real-time analysis Methods 0.000 description 3
- 230000003321 amplification Effects 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000005674 electromagnetic induction Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000009432 framing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Discrete Mathematics (AREA)
- Operations Research (AREA)
- Tests Of Circuit Breakers, Generators, And Electric Motors (AREA)
- Control Of Electric Motors In General (AREA)
Abstract
本发明公开了一种基于频域分析的电机参数追踪方法，该方法包括：采集电机目标帧数量的电流数据及电压数据，根据电流数据计算感应电动势数据；对每一帧的电流数据、电压数据及感应电动势数据分别进行离散傅里叶变换，以得到对应电流数据的电流频谱信息、对应电压数据的电压频谱信息及对应感应电动势数据的电动势频谱信息；根据电流频谱信息及电压频谱信息计算得到每一帧的总电阻，及根据电流频谱信息及电动势频谱信息计算得到每一帧的运动等效电阻；根据总电阻及运动等效电阻计算得到每一帧的平均静态电阻。本发明可精确的分析电机参数，从而实现对电机更精确的控制，以及更好的进行温升保护。此外，还提出了振动播放装置、计算机设备和存储介质。
Description
技术领域
本发明涉及电机参数分析技术领域，尤其是涉及基于频域分析的电机参数追踪方法、装置、设备和介质。
背景技术
线性电机能够提供丰富的振动形式，在触觉反馈领域获得极大的肯定，在日常电子产品中得到广泛应用。而为提供效果较佳的、多样的触觉效果，需要对线性电机的电机参数进行精确的控制，从而为用户带来完美的用户体验。
线性电机作为一种典型的电磁感应设备，具有较明确的机电耦合形式，因此对电机参数进行准确是实现电机精确控制的关键。但电机在实际工作状态下，受到环境的影响，电机的电机参数可能存在一定的波动，从而增大的电机参数控制的难度。特别是电机内部的线圈静态电阻，当电机持续工作时，内部温度升高，而静态电阻也会出现较大的变化。因此如果能精确的追踪电机内部的静态电阻，那么对于电机的精确控制有非常大的帮助。
发明内容
基于此，有必要针对上述问题，提供可精确分析的基于频域分析的电机参数追踪方法、装置、设备和介质。
一种基于频域分析的电机参数追踪方法，所述方法包括：
采集电机目标帧数量的电流数据及电压数据，根据所述电流数据计算感应电动势数据；
对每一帧的所述电流数据、所述电压数据及所述感应电动势数据分别进行离散傅里叶变换，以得到对应所述电流数据的电流频谱信息、对应所述电压数据的电压频谱信息及对应所述感应电动势数据的电动势频谱信息；
根据所述电流频谱信息及所述电压频谱信息计算得到每一帧的总电阻，及根据所述电流频谱信息及所述电动势频谱信息计算得到每一帧的运动等效电阻；
根据所述总电阻及所述运动等效电阻计算得到每一帧的平均静态电阻。
在其中一个实施例中，所述根据所述电流数据计算感应电动势数据，包括：
将所述电流数据代入电机动子运动方程，根据所述电流数据及所述电机动子运动方程计算得到电机动子的运动速度；
将所述运动速度代入感应电动式公式，根据所述运动速度及所述感应电动式公式计算得到所述感应电动势数据。
在其中一个实施例中，所述根据所述电流数据计算感应电动势数据，还包括：
获取所述电机的机械阻尼参数、弹簧系数参数、机电耦合系数及所述电机动子的动子质量；
所述将所述电流数据代入电机动子运动方程，包括：
将所述动子质量、所述机械阻尼参数、所述弹簧系数参数、所述机电耦合系数及所述电流数据代入所述电机动子运动方程；
所述将所述运动速度代入感应电动式公式，包括：
将所述机电耦合系数及所述将所述运动速度代入感应电动式公式。
在其中一个实施例中，所述根据所述电流频谱信息及所述电压频谱信息计算得到每一帧的总电阻，包括：
在同一帧内，将相同频点的所述电压频谱信息与所述电流频谱信息进行点除，得到每个所述频点的总电阻。
在其中一个实施例中，所述根据所述电流频谱信息及所述电动势频谱信息计算得到每一帧的运动等效电阻，包括：
在同一帧内，将相同频点的所述电动势频谱信息与所述电流频谱信息进行点除，得到每个所述频点的运动等效电阻。
在其中一个实施例中，所述根据所述总电阻及所述运动等效电阻计算得到每一帧的平均静态电阻，包括：
在同一帧内，将相同频点的所述总电阻减去所述运动等效电阻，得到每个频点的静态电阻；
对所有频点的所述静态电阻取平均值，得到每一帧的平均静态电阻。
在其中一个实施例中，在所述根据所述总电阻及所述运动等效电阻计算得到平均静态电阻之后，还包括：
当所述平均静态电阻增大且小于预设电阻阈值时，根据所述平均静态电阻增大下一帧的电压数据；
当所述平均静态电阻大于或等于所述预设电阻阈值时，根据所述平均静态电阻减小下一帧的电压数据。
一种振动播放装置，所述装置包括：
数据采集模块，用于采集电机目标帧数量的电流数据及电压数据，根据所述电流数据计算感应电动势数据；
频谱信息获取模块，用于对每一帧的所述电流数据、所述电压数据及所述感应电动势数据分别进行离散傅里叶变换，以得到对应所述电流数据的电流频谱信息、对应所述电压数据的电压频谱信息及对应所述感应电动势数据的电动势频谱信息；
第一计算模块，用于根据所述电流频谱信息及所述电压频谱信息计算得到每一帧的总电阻，及根据所述电流频谱信息及所述电动势频谱信息计算得到每一帧的运动等效电阻；
第二计算模块，用于根据所述总电阻及所述运动等效电阻计算得到每一帧的平均静态电阻。
一种计算机可读存储介质，存储有计算机程序，所述计算机程序被处理器执行时，使得所述处理器执行如下步骤：
采集电机目标帧数量的电流数据及电压数据，根据所述电流数据计算感应电动势数据；
对每一帧的所述电流数据、所述电压数据及所述感应电动势数据分别进行离散傅里叶变换，以得到对应所述电流数据的电流频谱信息、对应所述电压数据的电压频谱信息及对应所述感应电动势数据的电动势频谱信息；
根据所述电流频谱信息及所述电压频谱信息计算得到每一帧的总电阻，及根据所述电流频谱信息及所述电动势频谱信息计算得到每一帧的运动等效电阻；
根据所述总电阻及所述运动等效电阻计算得到每一帧的平均静态电阻。
一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述计算机程序被所述处理器执行时，使得所述处理器执行如下步骤：
采集电机目标帧数量的电流数据及电压数据，根据所述电流数据计算感应电动势数据；
对每一帧的所述电流数据、所述电压数据及所述感应电动势数据分别进行离散傅里叶变换，以得到对应所述电流数据的电流频谱信息、对应所述电压数据的电压频谱信息及对应所述感应电动势数据的电动势频谱信息；
根据所述电流频谱信息及所述电压频谱信息计算得到每一帧的总电阻，及根据所述电流频谱信息及所述电动势频谱信息计算得到每一帧的运动等效电阻；
根据所述总电阻及所述运动等效电阻计算得到每一帧的平均静态电阻。
本发明提供了基于频域分析的电机参数追踪方法、装置、设备和介质，由于是采用频域分析并取均值的方式，降低了采集的数据中误差信息的干扰，也提高了计算的准确性。同时由于考虑了运动等效电阻的干扰，可以实现对电机更精确的控制，以及更好的进行温升保护。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
其中：
图1为一个实施例中基于频域分析的电机参数追踪方法的流程示意图；
图2为一个实施例中电压波形图的示意图；
图3为一个实施例中电流波形图的示意图；
图4为一个实施例中电压频谱信息的示意图；
图5为一个实施例中电流频谱信息的示意图；
图6为一个实施例中平均静态电阻的示意图；
图7为一个实施例中振动播放装置的结构示意图；
图8为一个实施例中计算机设备的结构框图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，图1为一个实施例中基于频域分析的电机参数追踪方法的流程示意图，本基于频域分析的电机参数追踪方法提供的步骤包括：
步骤102，采集电机目标帧数量的电流数据及电压数据，根据电流数据计算感应电动势数据。
本实施例提供的电机参数追踪方法可以对采集到的电流数据及电压数据进行整体分析，或对采集到的电流数据及电压数据进行实时分析。当进行整体分析时，采集较多帧数量的电流数据及电压数据后进行分帧分析。例如从当前开始采集电机的600帧的电流值及电压值，参见图2及图3，其中图2为电压波形图的示意图，图3为电流波形图的示意图。然后图中对每一帧的电流值及电压值进行单独计算。其中每帧的帧长度均相同，通常为10ms-20ms。而当进行实时分析时，从当前开始每采集1帧的电流数据及电压数据便进行后续的计算。
除此之外，为简化计算，也可对多帧数量的电流数据及电压数据进行间隔抽帧后用于后续的计算，或每次采集1帧的电流数据及电压数据后，每隔预设时长后再进行下一帧的电流数据及电压数据的采集。
在一个具体实施例中，当实时采集当前帧的电流数据及电压数据时，通过对当前帧的电流数据进行计算，以得到当前帧的感应电动势数据。具体的，首先将当前帧的电流数据i代入电机动子运动方程。其中，电机动子运动方程如下所示：
mx″+Cx′+kx＝BLi
其中，m表示电机动子的动子质量，c表示电机动子的机械阻尼参数，k表示电机动子的弹簧系数参数，BL表示机电耦合系数，这些参数为线性电机的固有物理参数，可以预先调取得到。而x表示电机动子的运动速度，通过将当前帧的电流数据i代入电机动子运动方程后计算求解得到。进一步的，将电机动子的运动速度x代入感应电动式公式后，可求解得到感应电动式。其中，感应电动式公式如下所示：
ur＝BLv
其中，v为求得的电机动子的运动速度，ur表示感应电动势。
在另一个具体实施例中，当进行整体分析时，在对电流数据进行分帧处理后，对每一帧的电流数据进行单独的电机动子运动方程计算及感应电动式公式计算后，可得到对应每一帧电流数据的感应电动势数据，在此不做赘述。
步骤104，对每一帧的电流数据、电压数据及感应电动势数据分别进行离散傅里叶变换，以得到对应电流数据的电流频谱信息、对应电压数据的电压频谱信息及对应感应电动势数据的电动势频谱信息。
具体的，根据离散傅里叶公式
对每一帧的电压数据u进行离散傅里叶变换，得到电压频谱信息U(f)。参见图4，图4为电压频谱信息U(f)的示意图。对每一帧的电流数据i进行离散傅里叶变换，得到电流频谱信息I(f)。参见图5，图5为电流频谱信息I(f)的示意图。对每一帧的感应电动势数据ur进行离散傅里叶变换，得到电动势频谱信息Ur(f)。
步骤106，根据电流频谱信息及电压频谱信息计算得到每一帧的总电阻，及根据电流频谱信息及电动势频谱信息计算得到每一帧的运动等效电阻。
具体的，在同一帧内，通过将相同频点的电压频谱信息与电流频谱信息进行点除，从而得到频谱内每个频点的总电阻。示例性的，将图4内频率为50Hz处频点的幅值与图5内频率为50Hz的处频点的幅值对应相除，即可得到频率为50Hz处频点的总电阻。用公式表示为Rr(f)＝Ur(f)./I(f)，其中Rr(f)表示总电阻。同理的，在同一帧内，通过将相同频点的电动势频谱信息与电流频谱信息进行点除，从而得到频谱内每个频点的运动等效电阻。用公式表示为Z(f)＝U(f)./I(f)，其中Z(f)为运动等效电阻。
步骤108，根据总电阻及运动等效电阻计算得到每一帧的平均静态电阻。
具体的，在同一帧内，首先将相同频点的总电阻减去运动等效电阻，得到每个频点的静态电阻。示例性的，如将频率为50Hz处频点的Z(f)将去频率为50Hz处频点的Rr(f)，从而得到频率为50Hz处频点的静态电阻，记为Re(f)。进一步的，对同一帧内的所有频点的静态电阻Re(f)取平均值，即将频谱图内所有频点的静态电阻Re(f)相加后除于频点总数，从而得到该帧内的平均静态电阻，记为avgRe(f)。通过上述计算后可得到电机每一帧的平均静态电阻。参见图6，图6为整合后的平均静态电阻的示意图。
在另一个具体实施例中，当步骤102中采用对多帧数量的电流数据及电压数据进行间隔抽帧，或每次采集1帧的电流数据及电压数据后，每隔预设时长后再进行下一帧的电流数据及电压数据采集的方式，采集到的平均静态电阻为存在时间间隔的平均静态电阻。而通过平滑曲线将间隔的平均静态电阻进行顺序连接，可得到连续的平均静态电阻，同样可用于参数的准确追踪分析。
进一步的，计算得到的平均静态电阻用于对下一帧的电压数据进行调整，以求获得较佳的振动播放效果以及对电机进行保护。在其中一个实施例中，当进行整体分析时，该下一帧指的是采集的目标帧数量后的下一帧。当进行实时分析时，该下一帧指的是相对于当前采集帧的下一帧。示例性的，当进行整体分析时，若整体分析时间段内的平均静态电阻逐渐增大且小于预设电压阈值(例如设定为14Ω)，则通过适当提高功放倍数(或在输出前对电压进行更新)，从而增大下一帧的电压数据，以达到期望振动效果。例如，若电机的原电压为6V，原电阻为9Ω，在电机工作一段时间后，计算得到平均静态电阻正逐渐增大且最后一帧的平均静态电阻为12Ω(小于预设电压阈值14Ω)，因此通过调整功放系数(或在输出前对电压进行更新)，使实际输出电压等比放大为9V，可获得较佳的振动播放效果。
在另一个具体实施例中，当平均静态电阻大于或等于预设电阻阈值时，根据平均静态电阻减小下一帧的电压数据。例如，电阻由常温20℃下初始值9Ω，而电机线圈所允许的最高温度为110℃(线圈温升系数30℃/Ω)，因此确定预设电阻阈值为12Ω。当平均静态电阻大于或等于该预设电阻阈值时，适当减小下一帧的电压数据，或直接控制振动播放装置停止工作，从而达到保护电机的效果。
上述基于频域分析的电机参数追踪方法，由于是采用频域分析并取均值的方式，降低了采集的数据中误差信息的干扰，也提高了计算的准确性。同时由于考虑了运动等效电阻的干扰，可以实现对电机更精确的控制，以及更好的进行温升保护。
在一个实施例中，如图7所示，提出了一种振动播放装置，该装置包括：
数据采集模块702，用于采集电机目标帧数量的电流数据及电压数据，根据电流数据计算感应电动势数据；
频谱信息获取模块704，用于对每一帧的电流数据、电压数据及感应电动势数据分别进行离散傅里叶变换，以得到对应电流数据的电流频谱信息、对应电压数据的电压频谱信息及对应感应电动势数据的电动势频谱信息；
第一计算模块706，用于根据电流频谱信息及电压频谱信息计算得到每一帧的总电阻，及根据电流频谱信息及电动势频谱信息计算得到每一帧的运动等效电阻；
第二计算模块708，用于根据总电阻及运动等效电阻计算得到每一帧的平均静态电阻。
上述振动播放装置，由于是采用频域分析并取均值的方式，降低了采集的数据中误差信息的干扰，也提高了计算的准确性。同时由于考虑了运动等效电阻的干扰，可以实现对电机更精确的控制，以及更好的进行温升保护。
在一个实施例中，数据采集模块702，还具体用于：将电流数据代入电机动子运动方程，根据电流数据及电机动子运动方程计算得到电机动子的运动速度；将运动速度代入感应电动式公式，根据运动速度及感应电动式公式计算得到感应电动势数据。
在一个实施例中，数据采集模块702，还具体用于：获取电机的机械阻尼参数、弹簧系数参数、机电耦合系数及电机动子的动子质量；将动子质量、机械阻尼参数、弹簧系数参数、机电耦合系数及电流数据代入电机动子运动方程；将机电耦合系数及将运动速度代入感应电动式公式。
在一个实施例中，第一计算模块706，还具体用于：在同一帧内，将相同频点的电压频谱信息与电流频谱信息进行点除，得到每个频点的总电阻。
在一个实施例中，第一计算模块706，还具体用于：在同一帧内，将相同频点的电动势频谱信息与电流频谱信息进行点除，得到每个频点的运动等效电阻。
在一个实施例中，第二计算模块708，还具体用于：在同一帧内，将相同频点的总电阻减去运动等效电阻，得到每个频点的静态电阻；对所有频点的静态电阻取平均值，得到每一帧的平均静态电阻。
在一个实施例中，振动播放装置还包括：电压调整模块，用于当平均静态电阻增大且小于预设电阻阈值时，根据平均静态电阻增大下一帧的电压数据；当平均静态电阻大于或等于预设电阻阈值时，根据平均静态电阻减小下一帧的电压数据。
图8示出了一个实施例中计算机设备的内部结构图。如图8所示，该计算机设备包括通过系统总线连接的处理器、存储器和网络接口。其中，存储器包括非易失性存储介质和内存储器。该计算机设备的非易失性存储介质存储有操作系统，还可存储有计算机程序，该计算机程序被处理器执行时，可使得处理器实现基于频域分析的电机参数追踪方法。该内存储器中也可储存有计算机程序，该计算机程序被处理器执行时，可使得处理器执行基于频域分析的电机参数追踪方法。本领域技术人员可以理解，图8中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
一种计算机设备，包括存储器、处理器以及存储在该存储器中并可在该处理器上执行的计算机程序，该处理器执行该计算机程序时实现如下步骤：采集电机目标帧数量的电流数据及电压数据，根据电流数据计算感应电动势数据；对每一帧的电流数据、电压数据及感应电动势数据分别进行离散傅里叶变换，以得到对应电流数据的电流频谱信息、对应电压数据的电压频谱信息及对应感应电动势数据的电动势频谱信息；根据电流频谱信息及电压频谱信息计算得到每一帧的总电阻，及根据电流频谱信息及电动势频谱信息计算得到每一帧的运动等效电阻；根据总电阻及运动等效电阻计算得到每一帧的平均静态电阻。
在一个实施例中，根据电流数据计算感应电动势数据，包括：将电流数据代入电机动子运动方程，根据电流数据及电机动子运动方程计算得到电机动子的运动速度；将运动速度代入感应电动式公式，根据运动速度及感应电动式公式计算得到感应电动势数据。
在一个实施例中，根据电流数据计算感应电动势数据，还包括：获取电机的机械阻尼参数、弹簧系数参数、机电耦合系数及电机动子的动子质量；将电流数据代入电机动子运动方程，包括：将动子质量、机械阻尼参数、弹簧系数参数、机电耦合系数及电流数据代入电机动子运动方程；将运动速度代入感应电动式公式，包括：将机电耦合系数及将运动速度代入感应电动式公式。
在一个实施例中，根据电流频谱信息及电压频谱信息计算得到每一帧的总电阻，包括：在同一帧内，将相同频点的电压频谱信息与电流频谱信息进行点除，得到每个频点的总电阻。
在一个实施例中，根据电流频谱信息及电动势频谱信息计算得到每一帧的运动等效电阻，包括：在同一帧内，将相同频点的电动势频谱信息与电流频谱信息进行点除，得到每个频点的运动等效电阻。
在一个实施例中，根据总电阻及运动等效电阻计算得到每一帧的平均静态电阻，包括：在同一帧内，将相同频点的总电阻减去运动等效电阻，得到每个频点的静态电阻；对所有频点的静态电阻取平均值，得到每一帧的平均静态电阻。
在一个实施例中，在根据总电阻及运动等效电阻计算得到平均静态电阻之后，还包括：当平均静态电阻增大且小于预设电阻阈值时，根据平均静态电阻增大下一帧的电压数据；当平均静态电阻大于或等于预设电阻阈值时，根据平均静态电阻减小下一帧的电压数据。
一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现如下步骤：采集电机目标帧数量的电流数据及电压数据，根据电流数据计算感应电动势数据；对每一帧的电流数据、电压数据及感应电动势数据分别进行离散傅里叶变换，以得到对应电流数据的电流频谱信息、对应电压数据的电压频谱信息及对应感应电动势数据的电动势频谱信息；根据电流频谱信息及电压频谱信息计算得到每一帧的总电阻，及根据电流频谱信息及电动势频谱信息计算得到每一帧的运动等效电阻；根据总电阻及运动等效电阻计算得到每一帧的平均静态电阻。
在一个实施例中，根据电流数据计算感应电动势数据，包括：将电流数据代入电机动子运动方程，根据电流数据及电机动子运动方程计算得到电机动子的运动速度；将运动速度代入感应电动式公式，根据运动速度及感应电动式公式计算得到感应电动势数据。
在一个实施例中，根据电流数据计算感应电动势数据，还包括：获取电机的机械阻尼参数、弹簧系数参数、机电耦合系数及电机动子的动子质量；将电流数据代入电机动子运动方程，包括：将动子质量、机械阻尼参数、弹簧系数参数、机电耦合系数及电流数据代入电机动子运动方程；将运动速度代入感应电动式公式，包括：将机电耦合系数及将运动速度代入感应电动式公式。
在一个实施例中，根据电流频谱信息及电压频谱信息计算得到每一帧的总电阻，包括：在同一帧内，将相同频点的电压频谱信息与电流频谱信息进行点除，得到每个频点的总电阻。
在一个实施例中，根据电流频谱信息及电动势频谱信息计算得到每一帧的运动等效电阻，包括：在同一帧内，将相同频点的电动势频谱信息与电流频谱信息进行点除，得到每个频点的运动等效电阻。
在一个实施例中，根据总电阻及运动等效电阻计算得到每一帧的平均静态电阻，包括：在同一帧内，将相同频点的总电阻减去运动等效电阻，得到每个频点的静态电阻；对所有频点的静态电阻取平均值，得到每一帧的平均静态电阻。
在一个实施例中，在根据总电阻及运动等效电阻计算得到平均静态电阻之后，还包括：当平均静态电阻增大且小于预设电阻阈值时，根据平均静态电阻增大下一帧的电压数据；当平均静态电阻大于或等于预设电阻阈值时，根据平均静态电阻减小下一帧的电压数据。
需要说明的是，上述基于频域分析的电机参数追踪方法、装置、设备及计算机可读存储介质属于一个总的发明构思，基于频域分析的电机参数追踪方法、装置、设备及计算机可读存储介质实施例中的内容可相互适用。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，该程序可存储于一非易失性计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本申请专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种基于频域分析的电机参数追踪方法，其特征在于，所述方法包括：
采集电机目标帧数量的电流数据及电压数据，根据所述电流数据计算感应电动势数据；
对每一帧的所述电流数据、所述电压数据及所述感应电动势数据分别进行离散傅里叶变换，以得到对应所述电流数据的电流频谱信息、对应所述电压数据的电压频谱信息及对应所述感应电动势数据的电动势频谱信息；
根据所述电流频谱信息及所述电压频谱信息计算得到每一帧的总电阻，及根据所述电流频谱信息及所述电动势频谱信息计算得到每一帧的运动等效电阻；
根据所述总电阻及所述运动等效电阻计算得到每一帧的平均静态电阻。
2.根据权利要求1所述的基于频域分析的电机参数追踪方法，其特征在于，所述根据所述电流数据计算感应电动势数据，包括：
将所述电流数据代入电机动子运动方程，根据所述电流数据及所述电机动子运动方程计算得到电机动子的运动速度；
将所述运动速度代入感应电动式公式，根据所述运动速度及所述感应电动式公式计算得到所述感应电动势数据。
3.根据权利要求2所述的基于频域分析的电机参数追踪方法，其特征在于，所述根据所述电流数据计算感应电动势数据，还包括：
获取所述电机的机械阻尼参数、弹簧系数参数、机电耦合系数及所述电机动子的动子质量；
所述将所述电流数据代入电机动子运动方程，包括：
将所述动子质量、所述机械阻尼参数、所述弹簧系数参数、所述机电耦合系数及所述电流数据代入所述电机动子运动方程；
所述将所述运动速度代入感应电动式公式，包括：
将所述机电耦合系数及所述将所述运动速度代入感应电动式公式。
4.根据权利要求1所述的基于频域分析的电机参数追踪方法，其特征在于，所述根据所述电流频谱信息及所述电压频谱信息计算得到每一帧的总电阻，包括：
在同一帧内，将相同频点的所述电压频谱信息与所述电流频谱信息进行点除，得到每个所述频点的总电阻。
5.根据权利要求1所述的基于频域分析的电机参数追踪方法，其特征在于，所述根据所述电流频谱信息及所述电动势频谱信息计算得到每一帧的运动等效电阻，包括：
在同一帧内，将相同频点的所述电动势频谱信息与所述电流频谱信息进行点除，得到每个所述频点的运动等效电阻。
6.根据权利要求1所述的基于频域分析的电机参数追踪方法，其特征在于，所述根据所述总电阻及所述运动等效电阻计算得到每一帧的平均静态电阻，包括：
在同一帧内，将相同频点的所述总电阻减去所述运动等效电阻，得到每个频点的静态电阻；
对所有频点的所述静态电阻取平均值，得到每一帧的平均静态电阻。
7.根据权利要求1所述的基于频域分析的电机参数追踪方法，其特征在于，在所述根据所述总电阻及所述运动等效电阻计算得到平均静态电阻之后，还包括：
当所述平均静态电阻增大且小于预设电阻阈值时，根据所述平均静态电阻增大下一帧的电压数据；
当所述平均静态电阻大于或等于所述预设电阻阈值时，根据所述平均静态电阻减小下一帧的电压数据。
8.一种振动播放装置，其特征在于，所述装置包括：
数据采集模块，用于采集电机目标帧数量的电流数据及电压数据，根据所述电流数据计算感应电动势数据；
频谱信息获取模块，用于对每一帧的所述电流数据、所述电压数据及所述感应电动势数据分别进行离散傅里叶变换，以得到对应所述电流数据的电流频谱信息、对应所述电压数据的电压频谱信息及对应所述感应电动势数据的电动势频谱信息；
第一计算模块，用于根据所述电流频谱信息及所述电压频谱信息计算得到每一帧的总电阻，及根据所述电流频谱信息及所述电动势频谱信息计算得到每一帧的运动等效电阻；
第二计算模块，用于根据所述总电阻及所述运动等效电阻计算得到每一帧的平均静态电阻。
9.一种计算机可读存储介质，存储有计算机程序，所述计算机程序被处理器执行时，使得所述处理器执行如权利要求1至7中任一项所述方法的步骤。
10.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述计算机程序被所述处理器执行时，使得所述处理器执行如权利要求1至7中任一项所述方法的步骤。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010856065.1A CN112100569B (zh) | 2020-08-24 | 2020-08-24 | 基于频域分析的电机参数追踪方法、装置、设备和介质 |
| PCT/CN2020/113522 WO2022041309A1 (zh) | 2020-08-24 | 2020-09-04 | 基于频域分析的电机参数追踪方法、装置、设备和介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010856065.1A CN112100569B (zh) | 2020-08-24 | 2020-08-24 | 基于频域分析的电机参数追踪方法、装置、设备和介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112100569A true CN112100569A (zh) | 2020-12-18 |
| CN112100569B CN112100569B (zh) | 2024-04-02 |
Family
ID=73754268
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010856065.1A Active CN112100569B (zh) | 2020-08-24 | 2020-08-24 | 基于频域分析的电机参数追踪方法、装置、设备和介质 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN112100569B (zh) |
| WO (1) | WO2022041309A1 (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113702829A (zh) * | 2021-08-31 | 2021-11-26 | 歌尔股份有限公司 | 振动电机的扫频特性曲线生成方法、装置及存储介质 |
| CN113704966A (zh) * | 2021-07-22 | 2021-11-26 | 西安空间无线电技术研究所 | 一种识别调控囚禁离子振动模式及耦合强度的方法 |
Citations (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090051364A1 (en) * | 2005-03-30 | 2009-02-26 | Pcn Corporatioan | Onboard Battery Management Device |
| TW201018084A (en) * | 2008-10-29 | 2010-05-01 | Ultimax Tw Co Ltd | Electric damper |
| CN102291080A (zh) * | 2011-07-08 | 2011-12-21 | 哈尔滨工业大学 | 基于自适应补偿的异步电机参数辨识方法 |
| US20120116696A1 (en) * | 2009-03-24 | 2012-05-10 | Infinirel Corporation | Systems, devices and methods for predicting power electronics failure |
| US20140275886A1 (en) * | 2013-03-14 | 2014-09-18 | Streamline Automation, Llc | Sensor fusion and probabilistic parameter estimation method and apparatus |
| CN104836486A (zh) * | 2015-04-09 | 2015-08-12 | 哈尔滨工程大学 | 一种基于fpga的无刷直流电机节能控制系统及控制方法 |
| CN106788061A (zh) * | 2017-01-24 | 2017-05-31 | 东南大学 | 一种基于降阶电流环的永磁同步电机转动惯量识别方法 |
| US9954495B1 (en) * | 2016-10-24 | 2018-04-24 | Synaptics Incorporated | PVT invariant peaking stage for continuous time linear equalizer |
| CN108183647A (zh) * | 2017-12-28 | 2018-06-19 | 深圳市英威腾电气股份有限公司 | 一种异步电机离线静态参数辨识方法及其装置 |
| CN108270382A (zh) * | 2016-12-30 | 2018-07-10 | 杭州三花研究院有限公司 | 一种控制方法及装置 |
| CN109581237A (zh) * | 2018-11-15 | 2019-04-05 | 山东金人电气有限公司 | 一种动力电池动态内阻检测方法，设备及计算机可读存储介质 |
| CN109787235A (zh) * | 2019-01-29 | 2019-05-21 | 国网湖南省电力有限公司 | 新能源发电机组全功率范围谐波特性模型获取方法及系统 |
| US20190198303A1 (en) * | 2017-12-22 | 2019-06-27 | Thermo Fisher Scientific (Bremen) Gmbh | Method and Device for Crosstalk Compensation |
| CN110474593A (zh) * | 2018-12-25 | 2019-11-19 | 中国科学院长春光学精密机械与物理研究所 | 一种直流无刷电机电枢等效电阻电感测量方法 |
| CN110736926A (zh) * | 2019-10-14 | 2020-01-31 | 湖南银河电气有限公司 | 电机运行状态的特征参数提取方法和装置 |
| CN110995079A (zh) * | 2019-12-16 | 2020-04-10 | 瑞声科技(新加坡)有限公司 | 电机振动信号的生成方法、装置、终端及存储介质 |
| CN111030546A (zh) * | 2019-12-19 | 2020-04-17 | 深圳市显控科技股份有限公司 | 永磁电机离线参数识别方法和装置 |
| CN111478630A (zh) * | 2019-12-30 | 2020-07-31 | 瑞声科技(新加坡)有限公司 | 一种马达稳态单频失真补偿方法及装置 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8676524B2 (en) * | 2010-05-06 | 2014-03-18 | I-Shou University | System for computing machine parameters of an induction machine |
| CN103376362B (zh) * | 2012-04-25 | 2016-12-14 | 珠海格力电器股份有限公司 | 磁阻电机绕组阻值和换相位置的检测方法及装置 |
| CN109150049B (zh) * | 2018-08-08 | 2020-04-28 | 同济大学 | 一种静止坐标系电机分布式参数模型建立方法 |
| CN110289747A (zh) * | 2019-07-25 | 2019-09-27 | 青岛远洋船员职业学院 | 一种直线振荡电机的等效磁路建模方法 |
-
2020
- 2020-08-24 CN CN202010856065.1A patent/CN112100569B/zh active Active
- 2020-09-04 WO PCT/CN2020/113522 patent/WO2022041309A1/zh active Application Filing
Patent Citations (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090051364A1 (en) * | 2005-03-30 | 2009-02-26 | Pcn Corporatioan | Onboard Battery Management Device |
| TW201018084A (en) * | 2008-10-29 | 2010-05-01 | Ultimax Tw Co Ltd | Electric damper |
| US20120116696A1 (en) * | 2009-03-24 | 2012-05-10 | Infinirel Corporation | Systems, devices and methods for predicting power electronics failure |
| CN102291080A (zh) * | 2011-07-08 | 2011-12-21 | 哈尔滨工业大学 | 基于自适应补偿的异步电机参数辨识方法 |
| US20140275886A1 (en) * | 2013-03-14 | 2014-09-18 | Streamline Automation, Llc | Sensor fusion and probabilistic parameter estimation method and apparatus |
| CN104836486A (zh) * | 2015-04-09 | 2015-08-12 | 哈尔滨工程大学 | 一种基于fpga的无刷直流电机节能控制系统及控制方法 |
| US9954495B1 (en) * | 2016-10-24 | 2018-04-24 | Synaptics Incorporated | PVT invariant peaking stage for continuous time linear equalizer |
| CN108270382A (zh) * | 2016-12-30 | 2018-07-10 | 杭州三花研究院有限公司 | 一种控制方法及装置 |
| CN106788061A (zh) * | 2017-01-24 | 2017-05-31 | 东南大学 | 一种基于降阶电流环的永磁同步电机转动惯量识别方法 |
| US20190198303A1 (en) * | 2017-12-22 | 2019-06-27 | Thermo Fisher Scientific (Bremen) Gmbh | Method and Device for Crosstalk Compensation |
| CN108183647A (zh) * | 2017-12-28 | 2018-06-19 | 深圳市英威腾电气股份有限公司 | 一种异步电机离线静态参数辨识方法及其装置 |
| CN109581237A (zh) * | 2018-11-15 | 2019-04-05 | 山东金人电气有限公司 | 一种动力电池动态内阻检测方法，设备及计算机可读存储介质 |
| CN110474593A (zh) * | 2018-12-25 | 2019-11-19 | 中国科学院长春光学精密机械与物理研究所 | 一种直流无刷电机电枢等效电阻电感测量方法 |
| CN109787235A (zh) * | 2019-01-29 | 2019-05-21 | 国网湖南省电力有限公司 | 新能源发电机组全功率范围谐波特性模型获取方法及系统 |
| CN110736926A (zh) * | 2019-10-14 | 2020-01-31 | 湖南银河电气有限公司 | 电机运行状态的特征参数提取方法和装置 |
| CN110995079A (zh) * | 2019-12-16 | 2020-04-10 | 瑞声科技(新加坡)有限公司 | 电机振动信号的生成方法、装置、终端及存储介质 |
| CN111030546A (zh) * | 2019-12-19 | 2020-04-17 | 深圳市显控科技股份有限公司 | 永磁电机离线参数识别方法和装置 |
| CN111478630A (zh) * | 2019-12-30 | 2020-07-31 | 瑞声科技(新加坡)有限公司 | 一种马达稳态单频失真补偿方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 程啟华: "无刷直流电动机匝间短路故障检测及识别方法", 《微特电机》, vol. 42, no. 02, pages 14 - 17 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113704966A (zh) * | 2021-07-22 | 2021-11-26 | 西安空间无线电技术研究所 | 一种识别调控囚禁离子振动模式及耦合强度的方法 |
| CN113704966B (zh) * | 2021-07-22 | 2024-04-09 | 西安空间无线电技术研究所 | 一种识别调控囚禁离子振动模式及耦合强度的方法 |
| CN113702829A (zh) * | 2021-08-31 | 2021-11-26 | 歌尔股份有限公司 | 振动电机的扫频特性曲线生成方法、装置及存储介质 |
| CN113702829B (zh) * | 2021-08-31 | 2024-03-08 | 歌尔股份有限公司 | 振动电机的扫频特性曲线生成方法、装置及存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2022041309A1 (zh) | 2022-03-03 |
| CN112100569B (zh) | 2024-04-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112100569A (zh) | 基于频域分析的电机参数追踪方法、装置、设备和介质 | |
| Crow et al. | The matrix pencil for power system modal extraction | |
| CN111552371B (zh) | 激励电压生成方法、装置、设备及介质、测试方法及系统 | |
| Wang et al. | Online SVM regression algorithm-based adaptive inverse control | |
| CN110995079B (zh) | 电机振动信号的生成方法、装置、终端及存储介质 | |
| Zhu et al. | H∞ stabilization criterion with less complexity for nonuniform sampling fuzzy systems | |
| CN112650388B (zh) | 马达振动信号生成方法、装置、计算机设备及存储介质 | |
| CN112100884B (zh) | 线性马达的触控装置的振动量确定方法及相关设备 | |
| CN110348055A (zh) | Chaboche粘塑性本构模型材料参数获取与优化方法 | |
| CN111598303A (zh) | 一种基于气象成分分解的夏季短期负荷预测方法 | |
| CN111460639B (zh) | 动态推算电缆导体温度的方法、装置、电子设备及存储介质 | |
| Awoda et al. | Parameter estimation of a permanent magnetic DC motor | |
| CN109102110B (zh) | 一种径流式小水电短期出力预测方法和装置 | |
| Fan | Least squares estimation and Kalman filter based dynamic state and parameter estimation | |
| Takemura et al. | Proposal of novel rolling friction compensation with data-based friction model for ball screw driven stage | |
| Khatibi et al. | Comparative analysis of power system model reduction | |
| CN112198795B (zh) | 机电伺服控制方法、系统、终端设备及存储介质 | |
| Kopchak | Approximation accuracy of electromechanical systems high order objects using different types of fractional order transfer functions | |
| JP5710513B2 (ja) | 信号抽出装置、方法、及びプログラム | |
| KR102226677B1 (ko) | 전력계통의 하모닉 임피던스 측정장치 및 측정방법 | |
| Nguyen et al. | Dynamic state estimation and prediction for real-time control and operation | |
| Larsson et al. | Continuous-time AR parameter estimation by using properties of sampled systems | |
| CN113794198A (zh) | 抑制宽频振荡的方法、装置、终端及存储介质 | |
| CN117970781B (zh) | 一种逆变器智能控制方法及系统 | |
| CN111711388B (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 |