一种基于参数回归的平面电机控制方法及终端设备与流程

本发明涉及平面电机技术领域，特别涉及一种基于参数回归的平面电机控制方法及终端设备。

背景技术：

平面电机在集成电路、超精密加工等领域极具应用前景。根据电磁力产生原理，平面电机主要分为平面步进电机、平面感应电机、平面永磁同步电机、平面直流电机和平面开关磁阻电机5类。其中，平面开关磁阻电机因其结构简单、制造容易、精度高、成本低、热耗低、可靠性高等优点，在平面电机领域已成为重要发展趋势。

目前平面开关磁阻电机的位置控制方法有：1.比例微分积分控制：实现简单，点对点定位精度为5μm，动态跟踪精度为37.4μm；2.鲁棒控制：增强了系统的鲁棒性，但未给出电机的位置控制精度；3.自抗扰控制：抗干扰能力强，但未给出电机的位置控制精度；4.模型参考控制：自适应能力强，定位精度小于0.5μm，但较难实现高精度动态跟踪控制；5.滑模变结构控制：增强了系统的鲁棒性，点对点定位精度为15.5μm，动态跟踪精度为79.2μm；6.极点配置自校正控制：自调整能力强，动态跟踪精度为58.2μm。虽然上述方法可以保证平面电机点对点定位精度较高，但是，无法保证平面电机动态跟踪精度。

技术实现要素：

鉴于现有技术的不足，本发明旨在提供一种基于参数回归的平面电机控制方法及终端设备。

本发明所采用的技术方案如下：

一种基于参数回归的平面电机控制方法，其包括：

根据平面电机的输入输出信号建立脉冲传递函数，并根据所述脉冲传递函数确定平面电机的参数回归模型；

根据所述参数回归模型确定用于预测所述平面电机的第一预设数量时域的预测位置的预测模型；

基于所述预测模型确定的预测位置建立代价函数，并根据所述代价函数计算平面电机的推力指令信号；

将所述推力指令信号作用于平面电机的驱动装置，以控制所述平面电机运动。

所述参数回归的平面电机控制方法，其中，所述根据平面电机的输入输出信号建立脉冲传递函数，并根据所述脉冲传递函数确定平面电机的参数回归模型具体包括：

以推力信号为输入信号、位置信号为输出信号建立平面电机j轴的脉冲传递函数，其中，j轴为x轴或y轴；

获取所述脉冲传递函数的差分方程，并根据所述差分方程确定所述平面电机的参数回归模型。

所述参数回归的平面电机控制方法，其中，所述以推力信号为输入信号、位置信号为输出信号建立平面电机j轴的脉冲传递函数具体包括：

以推力信号为输入信号、位置信号为输出信号建立平面电机j轴的脉冲传递函数；

采用具有遗忘因子的递归最小二乘法辨识所述脉冲传递函数的系数，以得到所述脉冲传递函数。

所述的参数回归的平面电机控制方法，其中，所述脉冲传递函数的表达式为

其中，fj(z)为输入信号，pj(z)为输出信号，aj0，aj1，bj0，bj1和bj2均为系数。

所述的参数回归的平面电机控制方法，其中，所述参数回归模型的表达式为：

其中，cj＝[10]，xj_1(k)＝pj(k)-bj2fj(k)，xj_2(k)＝-aj0pj(k-1)+bj0fj(k-1)。

所述参数回归的平面电机控制方法，其中，所述基于所述预测模型确定的预测位置建立代价函数，并根据所述代价函数计算平面电机的推力指令信号具体包括：

基于所述预测模型确定的预测位置建立代价函数，其中，所述代价函数包括预测位置序列和控制序列；

根据所述代价函数计算所述平面电机的控制序列，并根据所述控制序列确定下一时刻的水平推力，以根据所述水平推力确定推力指令信号。

所述参数回归的平面电机控制方法，其中，所述根据所述代价函数计算所述平面电机的控制序列，并根据所述控制序列确定下一时刻的水平推力，以根据所述水平推力确定推力指令信号具体包括：

以所述代价函数最小为依据计算最优控制序列；

选取所述最优控制序列的第二分量，并将所述第二分量作为下一时刻的水平推力，以根据所述水平推力确定推力指令信号

所述参数回归的平面电机控制方法，其中，所述基于所述预测模型确定的预测位置建立代价函数具体包括：

获取前一时刻的控制序列，并根据所述前一时刻的控制序列确定未来预设数量时域的预测位置序列；

根据所述预测位置序列以及当前时刻的控制序列建立代价函数。

一种计算机可读存储介质，所述计算机可读存储介质存储有一个或者多个程序，所述一个或者多个程序可被一个或者多个处理器执行，以实现如上任一所述的参数回归的平面电机控制方法中的步骤。

一种终端设备，其包括：处理器、存储器及通信总线；所述存储器上存储有可被所述处理器执行的计算机可读程序；

所述通信总线实现处理器和存储器之间的连接通信；

所述处理器执行所述计算机可读程序时实现如上所述的参数回归的平面电机控制方法中的步骤。

有益效果：与现有技术相比，本发明提供了一种基于参数回归的平面电机控制方法及终端设备，所述方法包括：根据平面电机的输入输出信号建立脉冲传递函数，并根据所述脉冲传递函数确定平面电机的参数回归模型；根据所述参数回归模型确定用于预测所述平面电机的第一预设数量时域的预测位置的预测模型；基于所述预测模型确定的预测位置建立代价函数，并根据所述代价函数计算平面电机的推力指令信号；将所述推力指令信号作用于平面电机的驱动装置，以控制所述平面电机运动。本发明基于参数回归的方法确定以水平推力为输入信号、位置为输出信号的预测模型，根据所述预测模型预测的位置确定代价函数，并根据代价函数计算平面电机的推力指令信号，以将所述推力指令信号作用于电机，提高了平面电机的动态跟踪精度。

附图说明

图1为本发明提供的基于参数回归的平面电机控制方法的系统结构框图。

图2为本发明提供的基于参数回归的平面电机控制方法的流程图。

图3为本发明提供的基于参数回归的平面电机控制方法的一个实施例的流程图。

图4为本发明提供的基于参数回归的平面电机控制方法中步骤s10的流程图。

图5为本发明提供的基于参数回归的平面电机控制方法中步骤s30的流程图。

图6为本发明提供的一种终端设备的结构原理图。

具体实施方式

本发明提供一种基于参数回归的平面电机控制方法及终端设备，为使本发明的目的、技术方案及效果更加清楚、明确，以下参照附图并举实施例对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。

本技术领域技术人员可以理解，除非特意声明，这里使用的单数形式“一”、“一个”、“所述”和“该”也可包括复数形式。应该进一步理解的是，本发明的说明书中使用的措辞“包括”是指存在所述特征、整数、步骤、操作、元件和/或组件，但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件、组件和/或它们的组。应该理解，当我们称元件被“连接”或“耦接”到另一元件时，它可以直接连接或耦接到其他元件，或者也可以存在中间元件。此外，这里使用的“连接”或“耦接”可以包括无线连接或无线耦接。这里使用的措辞“和/或”包括一个或更多个相关联的列出项的全部或任一单元和全部组合。

本技术领域技术人员可以理解，除非另外定义，这里使用的所有术语(包括技术术语和科学术语)，具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是，诸如通用字典中定义的那些术语，应该被理解为具有与现有技术的上下文中的意义一致的意义，并且除非像这里一样被特定定义，否则不会用理想化或过于正式的含义来解释。

本发明提供了一种基于参数回归的平面电机控制方法，所述方法应用于如图1所述控制系统，如图2所示，所述方法具体是预先根据平面电机的输入输出信号建立脉冲传输函数，并确定相应的参数回归模型，再基于所述参数回归模型确定位置预测模型，根据位置预测模型计算位置预测序列，再根据位置预测序列以及参考轨迹确定代价函数，再通过对代价函数求极小值得到最优控制量，之后根据最优控制量确定当前时刻的控制量，并将所述控制量作用于平面电机，同时，在获取到当前时刻的控制量后，判断当前时刻是否为最后时刻，当不为最后时刻时，重复上述过程直到最后时刻。

下面为了对发明内容作进一步说明，结合一个具体实施例对本发明做详细说明。

如图3所示，本实施例提供了一种基于参数回归的平面电机控制方法，所述方法包括：

s10、根据平面电机的输入输出信号建立脉冲传递函数，并根据所述脉冲传递函数确定平面电机的参数回归模型；

s20、根据所述参数回归模型确定用于预测所述平面电机的第一预设数量时域的预测位置的预测模型；

s30、基于所述预测模型确定的预测位置建立代价函数，并根据所述代价函数计算平面电机的推力指令信号；

s40、将所述推力指令信号作用于平面电机的驱动装置，以控制所述平面电机运动。

本实施例根据平面电机的输入输出信号建立平面电机的参数回归模型，并根据参数回归模型确定用于预测位置的预测模型，根据所述预测模型得到的预测位置确定代价函数，并根据代价函数计算平面电机的水平推力信号，再将所述推力指令信号作用于电机，提高了平面电机的动态跟踪精度。

具体地，在所述步骤s10中，所述平面电机的输入输出信号指的是平面电机的输入信号，以及平面电机的输出信号。其中，所述输入信号为水平推力信号fj，所述输出信号为位置信号pj，所述j表示j轴，j轴为x轴或者y轴。在本实施例中，所述平面电机为平面开关磁阻电机(psrm)，所述平面开关磁阻电机的脉冲传递函数gj(z)的表达式可以为：

其中，fj(z)为输入信号，pj(z)为输出信号，aj0，aj1，bj0，bj1和bj2均为系数。

相应的，如图4所示，所述根据平面电机的输入输出信号建立脉冲传递函数，并根据所述脉冲传递函数确定平面电机的参数回归模型具体包括：

s11、以推力信号为输入信号、位置信号为输出信号建立平面电机j轴的脉冲传递函数，其中，j轴为x轴或者y轴；

s12、获取所述脉冲传递函数的差分方程，并根据所述差分方程确定所述平面电机的参数回归模型。

具体地，所述差分方程是对所述脉冲传递函数进行差分化得到，根据所述脉冲传递函数得到的差分方程为：

pj(k+2)+aj1pj(k+1)+aj0pj(k)＝bj2fj(k+2)+bj1fj(k+1)+bj0fj(k)

其中，k表示k时刻，k+1表示k时刻的下一时刻，k+2表示k+1的下一时刻，所述aj0，aj1，bj0，bj1和bj2为脉冲传递函数的系数。此外，在所述脉冲传递函数gj(z)中，所述系数aj0，aj1，bj0，bj1和bj2均为通过具有遗忘因子的递归最小二乘法辨识得到的系数。相应的，所述以推力信号为输入信号、位置信号为输出信号建立平面电机j轴的脉冲传递函数具体包括：以推力信号为输入信号、位置信号为输出信号建立平面电机j轴的脉冲传递函数，其中，j轴为x轴或者y轴；采用具有遗忘因子的递归最小二乘法辨识所述脉冲传递函数的系数，以得到所述脉冲传递函数。

进一步，在获取到差分方程后，根据所述差分方程可以获取到参数回归模型。在获取所述参数回归模型时，可以定义离散状态变量xj_1(k)和xj_2(k)，所述xj_1(k)和xj_2(k)的表达式分别为：

xj_1(k)＝pj(k)-bj2fj(k)

xj_2(k)＝-aj0pj(k-1)+bj0fj(k-1)

那么，根据差分方程和离散状态变量可以得到k+1时刻的离散状态变量xj_1(k+1)和xj_2(k+1)，所述xj_1(k+1)和xj_2(k+1)的表达式分别为：

xj_1(k+1)＝-aj1xj_1(k)+xj_2(k)+(bj1-aj1bj2)fj(k)

xj_2(k+1)＝-aj0xj_1(k)+(bj0-aj0bj2)fj(k)

进一步，根据所述xj_1(k+1)和xj_2(k+1)可以确定所述参数回归模型的表达式为：

其中，cj＝[10]，yj(k)为k时刻的预测位置。

进一步，在所述步骤s20中，所述第一预设数量为预先设定，这里用p表示第一预设数量。也就是说，所述预测模型可以预测未来p个时刻的电机位置并输出电机位置对应的水平推力。在确定参数回归模型后，取预测时域等于控制时域为nt，t为采样周期，根据参数回归模型确定k+1时刻的预测模型为

yj_n(k+1|k-1)＝gjxj(k-1)+hjuj(k)

其中，n为整数且n＞2，yj_n(k+1|k-1)，xj(k-1)，uj(k)，gj和hj分别为

此外，根据所述预测模型可以得出预测位置序列yj_n(k+1|k-1)是控制序列uj(k)的函数，(k+1)t，(k+2)t到(k+n)t时刻的预测位置可由序列yj_n(k+1|k-1)得出。

进一步，在实施例的一个变形实施例中，本发明中所述脉冲传递函数gj(z)的分子分母均为z的二阶多项式，若bj2＝0，则离散状态变量可变形为x'j_1(k)和x'j_2(k)，分别为

x′j_1(k)＝pj(k)

x′j_2(k)＝-aj0pj(k-1)+bj0fj(k-1)

所述平面开关磁阻电机j轴运动的状态空间方程可变形为

式中，

那么，所述预测模型为

y′j_n(k+1|k)＝gjx′j(k)+h′ju′j(k)

其中，yj_n(k+1|k)，x'j(k)，u'j(k)和h'j分别为

进一步，在所述步骤s30中，所述代价函数是以参考位置与预测位置的误差和输入控制量为目标，并通过所述代价函数确定最优输入控制序列。所述预测模型确定的预测位置建立代价函数是根据当前时刻的前一时刻的控制序列确定预测位置序列，即所述基于所述预测模型确定的预测位置建立代价函数具体为获取前一时刻的控制序列，并根据所述前一时刻的控制序列确定未来预设数量时域的预测位置序列；根据所述预测位置序列以及当前时刻的控制序列建立代价函数。在本实施例中，所述代价函数的表达式可以为：

其中，r(k)为k时刻的参考位置，q1和q2是正对角的加权矩阵，分别为：

r(k)＝[r(k)lr(k+n-1)]^t1×n

此外，由于预测位置序列yj_n(k+1|k-1)是控制序列uj(k)的函数，那么将代价函数中的yj_n(k+1|k-1)采用uj(k)替换，并计算替换后的代价函数的极小值可以得到最优控制序列。所述最优控制序列uj(k)的表达式可以为：

同时在本实施例中，如图5所示，所述基于所述预测模型确定的预测位置建立代价函数，并根据所述代价函数计算平面电机的推力指令信号具体包括：

s31、基于所述预测模型确定的预测位置建立代价函数，其中，所述代价函数包括预测位置序列和控制序列；

s32、根据所述代价函数计算所述平面电机的控制序列，并根据所述控制序列确定下一时刻的水平推力，以根据所述水平推力确定推力指令信号。

具体地，所述预测位置序列是根据预测模型得到，即根据所述预测模型获取预测位置序列，再根据预测位置序列以及控制序列确定代价函数。所述代价函数包括两个目标优化量，其中，一个目标约束条件为参考位置与预测位置的差值，另一个是控制序列。此外，根据代价函数计算得到最优控制序列，再根据所述最优控制序列得到下一时刻的水平推力，其中，当前时刻为k时刻，下一时刻为k+1时刻。

示例性地，所述根据所述代价函数计算所述平面电机的控制量序列，并根据所述控制量序列确定下一时刻的水平推力具体包括：

以所述代价函数最小为依据计算最优控制序列；

选取所述最优控制序列的第二分量，并将所述第二分量作为下一时刻的水平推力，以根据所述水平推力确定推力指令信号。

具体地，所述下一时刻的水平推力为输入控制量序列的第二分量，即在计算得到输入控制量序列时，选取所述输入控制量序列的第二分量，并将所述第二分量作为下一刻的水平推力；基于所述水平推力确定推力指令信号，并将所述确定的水平推力指令作用于平面电机的驱动装置，以控制所述平面电机运动，这样可以提高平面电机运动的精确性。

进一步，在步骤s40中，在获取到k时刻的水平推力后，在k+1时刻重复上述过程可以获取k+2时刻的输入控制量，依次类推直至k+n时刻。此外在实际应用中，在k+1时刻获取k+2时刻的输入控制量之前，可以先判断所述k+1时刻是否为最后时刻，当k+1不为最后时刻时，执行上步骤获取k+2时刻的输入控制量，当k+1为最后时刻时，则完成输入控制量获取操作。

基于上述参数回归的平面电机控制方法，本发明还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有一个或者多个程序，所述一个或者多个程序可被一个或者多个处理器执行，以实现如上任一所述的参数回归的平面电机控制方法中的步骤。

基于上述参数回归的平面电机控制方法，本发明还提供了一种移动终端，如图6所示，其包括至少一个处理器(processor)20；显示屏21；以及存储器(memory)22，还可以包括通信接口(communicationsinterface)23和总线24。其中，处理器20、显示屏21、存储器22和通信接口23可以通过总线24完成相互间的通信。显示屏21设置为显示初始设置模式中预设的用户引导界面。通信接口23可以传输信息。处理器20可以调用存储器22中的逻辑指令，以执行上述实施例中的方法。

此外，上述的存储器22中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。

存储器22作为一种计算机可读存储介质，可设置为存储软件程序、计算机可执行程序，如本公开实施例中的方法对应的程序指令或模块。处理器20通过运行存储在存储器22中的软件程序、指令或模块，从而执行功能应用以及数据处理，即实现上述实施例中的方法。

存储器22可包括程序存储区和数据存储区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序；存储数据区可存储根据终端设备的使用所创建的数据等。此外，存储器22可以包括高速随机存取存储器，还可以包括非易失性存储器。例如，u盘、移动硬盘、只读存储器(read-onlymemory，rom)、随机存取存储器(randomaccessmemory，ram)、磁碟或者光盘等多种可以存储程序代码的介质，也可以是暂态存储介质。

此外，上述存储介质以及移动终端中的多条指令处理器加载并执行的具体过程在上述方法中已经详细说明，在这里就不再一一陈述。

最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。