CN116442240B - 一种基于高通滤波解耦的机器人零力控制方法及装置 - Google Patents

Abstract

本发明公开了一种基于高通滤波解耦的机器人零力控制方法及装置，方法包括：获取机械臂端部所受的合力；合力为机械臂抓握物体时机械臂端部所受的第一作用力与操作者对机械臂端部的第二作用力的合力；根据合力、机器人的质量特性系数、机器人的阻尼特性系数、机器人的弹簧特性系数以及机械臂的电机关节位置信息，计算得到待处理的输入信号；根据预先构建的双输入单输出状态空间表达式确定对待处理的输入信号的高通滤波截止频率；根据高通滤波截止频率对待处理的输入信号进行高通滤波，以消除待处理的输入信号中的第一作用力，得到第二作用力。本发明将机械臂端部所受合力解耦并得到操作者对机械臂端部的作用力，可广泛应用于机器人控制领域。

Description

一种基于高通滤波解耦的机器人零力控制方法及装置

技术领域

本发明涉及机器人控制领域，尤其是一种基于高通滤波解耦的机器人零力控制方法及装置。

背景技术

机器人学作为一门迅速发展的综合性的前沿科学，机器人的理论和技术在不断的更新与迭代。传统的工业中，为了能让机器人更好的帮助人类完成任务，出现了机器人的示教需求，主要分成在线和离线两类，前者通过操作者使用示教工具来控制机器人的末端位姿到达期望位姿，后者则是通过编程实现工件的轨迹规划。然而两者都要求操作者具有特定的机器人学知识，导致示教效率较低，但协作机器人的诞生改变了这一现状。研究人员在研究人机交互的过程中，出现了一种新兴的示教模式，即直接示教技术。人与机器人的直接接触来牵引机器人的末端来达到期望位姿，进而实现对机器人的运动控制。这样一来便可大大提高示教的效率。

在拖动示教中，由于打破了传统的人与机器单独分离的空间间隔，因此不能像传统机器人的控制方式那样，仅仅依靠位置模式就实现机器人的运动与轨迹规划等。由于协作机器人处于复杂的无约束的环境中，相较于工业机器人而言缺乏一定的感知能力，所以需要解决的关键问题，即零力控制。

目前，零力控制包括将机器人与人交互表现中的运动参数实时的传递给控制层，然后通过补偿等方法控制机器人与外界力的感知。但机械臂端部抓握的物体重力与人施加于机械臂端部的作用力将耦合成一个合力作为信号输入到零力控制的算法当中，导致现有零力控制方法无法灵活运用。

发明内容

有鉴于此，本发明实施例提供一种基于高通滤波解耦的机器人零力控制方法及装置，用于将机械臂端部所受合力解耦并从解耦的合力中得到操作者对机械臂端部的作用力。

本发明实施例的一方面提供了一种基于高通滤波解耦的机器人零力控制方法，包括：

获取机械臂端部所受的合力；其中，所述合力为机械臂抓握物体时所述机械臂端部所受的第一作用力与操作者对所述机械臂端部的第二作用力的合力；

根据所述合力、机器人的质量特性系数、机器人的阻尼特性系数、机器人的弹簧特性系数以及所述机械臂的电机关节位置信息，计算得到待处理的输入信号；

根据预先构建的双输入单输出状态空间表达式确定对所述待处理的输入信号的高通滤波截止频率；

根据所述高通滤波截止频率对所述待处理的输入信号进行高通滤波，以消除所述待处理的输入信号中的所述第一作用力，得到所述第二作用力。

可选地，所述获取机械臂端部所受的合力，包括：

通过设置于所述机械臂端部的六轴传感器获取所述机械臂端部所受的合力。

可选地，所述根据所述合力、机器人的质量特性系数、机器人的阻尼特性系数、机器人的弹簧特性系数以及所述机械臂的电机关节位置信息，计算得到待处理的输入信号，包括：

根据第一表达式计算得到待处理的输入信号，所述第一表达式为：

其中，F_N表示所述待处理的输入信号，F_c表示所述合力，M表示所述质量特性系数，B表示所述阻尼特性系数，K表示所述弹簧特性系数，x表示所述电机关节位置信息，与/>分别表示所述电机关节位置信息的一阶微分与二阶微分。

可选地，所述根据预先构建的双输入单输出状态空间表达式确定对所述待处理的输入信号的高通滤波截止频率，包括：

根据所述双输入单输出状态空间表达式确定所述总输入信号中位置输入系统的系统矩阵和所述总输入信号中力矩输入系统的系统矩阵；

所述双输入单输出状态空间表达式为：

其中，x_R和x_N分别表示所述位置输入系统和所述力矩输入系统的状态向量，A_R和A_N分别表示所述位置输入系统和所述力矩输入系统的系统矩阵，b_R和b_N分别表示所述位置输入系统和所述力矩输入系统的控制矩阵，c_R和c_N表示所述位置输入系统和所述力矩输入系统的输出矩阵，y为所述双输入单输出状态空间表达式的输出；

根据由所述系统矩阵的特征值确定的最小根，确定所述高通滤波截止频率。

可选地，所述系统矩阵、所述控制矩阵以及所述输出矩阵分别表示为：

其中，各个所述矩阵的参数如下：

A_i＝(m+M)L,i＝{R N}

B_i＝(f+mω_c+B)L+(m+M)R,i＝{R N}

C_i＝(fω_c+K)L+(mω_c+f+B)R+C_kK_t,i＝{R N}

D_i＝(fω_c+K)R+C_kK_tω_c+K_tK_p,i＝{R N}

E_i＝K_tK_pω_c,i＝{R N}

F_R＝K_tK_p

F_N＝-L

G_R＝K_tK_pω_c

G_N＝-R

其中，m表示转子质量，M表示所述质量特性系数，B表示所述阻尼特性系数，K表示所述弹簧特性系数，f表示摩擦系数，K_t表示推力常数，K_p表示PID控制器系数，C_k表示反电势常数，R表示电阻，L表示电感。

可选地，所述根据由所述系统矩阵的特征值确定的最小根，确定所述高通滤波截止频率，包括：

将所述最小根表示为x₁，确定所述高通滤波截止频率的取值范围为0<ω_c<x₁，其中，ω_c表示所述高通滤波截止频率。

可选地，所述根据所述高通滤波截止频率对所述待处理的输入信号进行高通滤波，以消除所述待处理的输入信号中的所述第一作用力，得到所述第二作用力，包括：

根据高通滤波器的第二表达式对所述待处理的输入信号进行高通滤波；

所述第二表达式为：

其中，s表示拉氏变量，ω_c表示所述高通滤波截止频率。

本发明实施例的另一方面还提供了一种基于高通滤波解耦的机器人零力控制装置，包括：

合力获取单元，用于获取机械臂端部所受的合力；其中，所述合力为机械臂抓握物体时所述机械臂端部所受的第一作用力与操作者对所述机械臂端部的第二作用力的合力；

信号计算单元，用于根据所述合力、机器人的质量特性系数、机器人的阻尼特性系数、机器人的弹簧特性系数以及所述机械臂的电机关节位置信息，计算得到待处理的输入信号；

频率确定单元，用于根据预先构建的双输入单输出状态空间表达式确定对所述待处理的输入信号的高通滤波截止频率；

合力解耦单元，用于根据所述高通滤波截止频率对所述待处理的输入信号进行高通滤波，以消除所述待处理的输入信号中的所述第一作用力，得到所述第二作用力。

本发明实施例的另一方面还提供了一种电子设备，包括处理器以及存储器；

所述存储器用于存储程序；

所述处理器执行所述程序实现所述的一种基于高通滤波解耦的机器人零力控制方法。

本发明实施例的另一方面还提供了一种计算机可读存储介质，所述存储介质存储有程序，所述程序被处理器执行实现所述的一种基于高通滤波解耦的机器人零力控制方法。

本发明实施例还公开了一种计算机程序产品或计算机程序，该计算机程序产品或计算机程序包括计算机指令，该计算机指令存储在计算机可读存储介质中。计算机设备的处理器可以从计算机可读存储介质读取该计算机指令，处理器执行该计算机指令，使得该计算机设备执行前面的方法。

本发明先根据机械臂端部所受的合力、机器人的质量特性系数、机器人的阻尼特性系数、机器人的弹簧特性系数以及机械臂的电机关节位置信息，计算得到待处理的输入信号，再根据预先构建的双输入单输出状态空间表达式确定对待处理的输入信号的高通滤波截止频率，然后根据高通滤波截止频率对待处理的输入信号进行高通滤波，将机械臂端部抓握物体的重力，即第一作用力作为低频信号进行消除，从而将合力解耦分开，实现了在协作机器人零力控制过程中，只将操作者对机械臂端部的作用力作为零力控制系统的唯一输入。

附图说明

为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例提供的一种基于高通滤波解耦的机器人零力控制方法的流程示意图；

图2为本发明实施例提供的一种对输入信号进行高通滤波的示例图；

图3为本发明实施例提供的一种基于高通滤波解耦的机器人零力控制原理框图；

图4为本发明实施例提供的一种基于高通滤波解耦的机器人零力控制装置的结构框图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

参照图1，本发明实施例提供了一种基于高通滤波解耦的机器人零力控制方法，具体包括以下步骤：

S100：获取机械臂端部所受的合力；其中，所述合力为机械臂抓握物体时所述机械臂端部所受的第一作用力与操作者对所述机械臂端部的第二作用力的合力。

具体的，本发明实施例可以通过设置于所述机械臂端部的六轴传感器获取所述机械臂端部所受的合力。

S110：根据所述合力、机器人的质量特性系数、机器人的阻尼特性系数、机器人的弹簧特性系数以及所述机械臂的电机关节位置信息，计算得到待处理的输入信号。

具体的，根据第一表达式计算得到待处理的输入信号，所述第一表达式为：

其中，F_N表示所述待处理的输入信号，F_c表示所述合力，M表示所述质量特性系数，B表示所述阻尼特性系数，K表示所述弹簧特性系数，x表示所述电机关节位置信息，与/>分别表示所述电机关节位置信息的一阶微分与二阶微分。

S120：根据预先构建的双输入单输出状态空间表达式确定对所述待处理的输入信号的高通滤波截止频率。

具体的，确定高通滤波截止频率的过程包括：

S1、根据所述双输入单输出状态空间表达式确定所述总输入信号中位置输入系统的系统矩阵和所述总输入信号中力矩输入系统的系统矩阵；

所述双输入单输出状态空间表达式为：

其中，x_R和x_N分别表示所述位置输入系统和所述力矩输入系统的状态向量，A_R和A_N分别表示所述位置输入系统和所述力矩输入系统的系统矩阵，b_R和b_N分别表示所述位置输入系统和所述力矩输入系统的控制矩阵，c_R和c_N表示所述位置输入系统和所述力矩输入系统的输出矩阵，y为所述双输入单输出状态空间表达式的输出。

具体的，所述系统矩阵、所述控制矩阵以及所述输出矩阵分别表示为：

其中，各个所述矩阵的参数如下：

A_i＝(m+M)L,i＝{R N}

B_i＝(f+mω_c+B)L+(m+M)R,i＝{R N}

C_i＝(fω_c+K)L+(mω_c+f+B)R+C_kK_t,i＝{R N}

D_i＝(fω_c+K)R+C_kK_tω_c+K_tK_p,i＝{R N}

E_i＝K_tK_pω_c,i＝{R N}

F_R＝K_tK_p

F_N＝-L

G_R＝K_tK_pω_c

G_N＝-R

其中，m表示转子质量，M表示所述质量特性系数，B表示所述阻尼特性系数，K表示所述弹簧特性系数，f表示摩擦系数，K_t表示推力常数，K_p表示PID控制器系数，C_k表示反电势常数，R表示电阻，L表示电感。

S2、根据由所述系统矩阵的特征值确定的最小根，确定所述高通滤波截止频率。

具体的，本发明实施例计算该最小根时，可以根据李雅普诺夫稳定判据，确定系统矩阵A的所有特征值小于0，进而设定一个tolerance，当sum(A_k-A_k-1(:))≤tolerance时结束迭代，由此得到的上三角矩阵的对角元素即为系统矩阵A的特征值λ_i，令λ_i＝0，得到最小与最大的两个根，分别表示为x₁和x₂。

然后，确定所述高通滤波截止频率的取值范围，本发明实施例的取值范围可以设置为0<ω_c<x₁，其中，ω_c表示所述高通滤波截止频率。

一种可选的实施方式下，取ω_c＝0.1x₁。

S130：根据所述高通滤波截止频率对所述待处理的输入信号进行高通滤波，以消除所述待处理的输入信号中的所述第一作用力，得到所述第二作用力。

具体的，根据高通滤波器的第二表达式对所述待处理的输入信号进行高通滤波；

所述第二表达式为：

其中，s表示拉氏变量，ω_c表示所述高通滤波截止频率。

接下来将以具体实例说明本发明的应用过程。

具体的，本发明实施例可以包括以下步骤：

步骤1、测量协作机械臂抓握物体重力G与人对机械臂端部作用力F_h的合力F_c。

具体的，所述步骤1通过在协作机械臂端部安装的六轴力传感器实现对两者合力F_c的测量。

步骤2、将合力F_c与系统“质量-阻尼-弹簧”动态特性结合，作为需要处理的总输入信号F_N。

具体的，所述步骤2的具体计算方法为：

其中，F_N表示待处理的输入信号，F_c表示合力，M、B、K分别表示系统的“质量、阻尼和弹簧”特性系数，x表示电机关节的位置信息，与/>表示位置信息的一阶微分与二阶微分。

步骤3、引入高通滤波器对待处理的输入信号F_N进行滤波处理。

具体的，所述步骤3将待处理的输入信号F_N，经过如图2所示的高通滤波器，得到处理后的输入信号F′_N。

步骤4、建立系统的一个双输入单输出状态空间表达式y。

具体的，所述步骤4建立关节电机引入高通滤波器之后包含控制器的整个系统的一个双输入单输出的状态空间表达式，表达式如下：

其中，x_R和x_N分别表示位置输入系统和力矩输入系统的状态向量，A_R和A_N分别表示位置输入系统和力矩输入系统的系统矩阵，b_R和b_N分别表示位置输入系统和力矩输入系统的控制矩阵，c_R和c_N则表示位置输入系统和力矩输入系统的输出矩阵，y为系统输出。其中，各个矩阵表示如下：

上式各个矩阵的参数如下：

A_i＝(m+M)L,i＝{R N}

B_i＝(f+mω_c+B)L+(m+M)R,i＝{R N}

C_i＝(fω_c+K)L+(mω_c+f+B)R+C_kK_t,i＝{R N}

D_i＝(fω_c+K)R+C_kK_tω_c+K_tK_p,i＝{R N}

E_i＝K_tK_pω_c,i＝{R N}

F_R＝K_tK_p

F_N＝-L

G_R＝K_tK_pω_c

G_N＝-R

其中，m表示转子质量，M，B，K分别表示关节电机系统的质量，阻尼，弹簧特性系数，f表示摩擦系数，K_t表示推力常数，K_p表示PID控制器系数，C_k表示反电势常数，R表示电阻，L表示电感。具体的机器人零力控制原理框图请参照图3。

步骤5、分析高通滤波器的可行与最佳截止频率ω_c。

具体的，所述步骤5由上述构建的双输入单输出状态空间表达式，可以得到机器人的系统矩阵A，控制矩阵B以及输出矩阵C。若需要机器人的系统稳定，则可以根据李雅普诺夫稳定判据，确定系统矩阵A的所有特征值小于0。

本发明实施例的系统矩阵A可以是8阶矩阵，因此，可以使用QR迭代法计算矩阵A的特征值：

设A∈C^n×n,记A₁＝A；

对A_t进行QR分解得：A₁＝Q₁×R₁,记A₂＝Q₁×R₁；

对A₂进行QR分解得：A₂＝Q₂×R₂,记A₃＝Q₁×R₁；

对A_k进行QR分解得：A_k＝Q_k×R_k,记A_k+1＝Q_k×R_k；

本发明实施例可以设定一个tolerance，当sum(A_k-A_k-1(:))≤tolerance时结束迭代，由此得到的上三角矩阵的对角元素即为系统矩阵A的特征值λ_i，令λ_i＝0，得到最小与最大的两个根，分别表示为x₁和x₂。

当高通滤波器的截至频率取值取ω_c>x₂时，在理论上机器人的系统是稳定的，但是高通滤波截至频率过大，会使关节电机的“质量-阻尼-弹簧”动态特性被消除，从而导致机器人的系统出现刚度过高问题。

当高通滤波器的截至频率取值取0<ω_c<x₁时，系统稳定，可选地，本发明实施例取ω_c＝0.1x₁。该取值可以使得关节电机系统收敛速度快，系统稳定。需要注意的是，在0<ω_c<x₁取值范围内，当ω_c取值过低，高通滤波几乎不起作用，收敛速度慢，当ω_c取值过高，系统震荡严重。

本发明实施例可以将六轴力传感器测量得到的合力解耦分离，实现零力控制中机器人仅受操作者作用于机械臂端部的作用力，从而驱动并产生运动。同时经过仿真数据分析，证明本发明实施例具有较好效果的高通滤波截至频率，降低了零力控制系统设计的难度。

参照图4，本发明实施例提供了一种基于高通滤波解耦的机器人零力控制装置，包括：

合力获取单元，用于获取机械臂端部所受的合力；其中，所述合力为机械臂抓握物体时所述机械臂端部所受的第一作用力与操作者对所述机械臂端部的第二作用力的合力；

信号计算单元，用于根据所述合力、机器人的质量特性系数、机器人的阻尼特性系数、机器人的弹簧特性系数以及所述机械臂的电机关节位置信息，计算得到待处理的输入信号；

频率确定单元，用于根据预先构建的双输入单输出状态空间表达式确定对所述待处理的输入信号的高通滤波截止频率；

合力解耦单元，用于根据所述高通滤波截止频率对所述待处理的输入信号进行高通滤波，以消除所述待处理的输入信号中的所述第一作用力，得到所述第二作用力。

本发明实施例还公开了一种计算机程序产品或计算机程序，该计算机程序产品或计算机程序包括计算机指令，该计算机指令存储在计算机可读存储介质中。计算机设备的处理器可以从计算机可读存储介质读取该计算机指令，处理器执行该计算机指令，使得该计算机设备执行图1所示的方法。

在一些可选择的实施例中，在方框图中提到的功能/操作可以不按照操作示图提到的顺序发生。例如，取决于所涉及的功能/操作，连续示出的两个方框实际上可以被大体上同时地执行或所述方框有时能以相反顺序被执行。此外，在本发明的流程图中所呈现和描述的实施例以示例的方式被提供，目的在于提供对技术更全面的理解。所公开的方法不限于本文所呈现的操作和逻辑流程。可选择的实施例是可预期的，其中各种操作的顺序被改变以及其中被描述为较大操作的一部分的子操作被独立地执行。

此外，虽然在功能性模块的背景下描述了本发明，但应当理解的是，除非另有相反说明，所述的功能和/或特征中的一个或多个可以被集成在单个物理装置和/或软件模块中，或者一个或多个功能和/或特征可以在单独的物理装置或软件模块中被实现。还可以理解的是，有关每个模块的实际实现的详细讨论对于理解本发明是不必要的。更确切地说，考虑到在本文中公开的装置中各种功能模块的属性、功能和内部关系的情况下，在工程师的常规技术内将会了解该模块的实际实现。因此，本领域技术人员运用普通技术就能够在无需过度试验的情况下实现在权利要求书中所阐明的本发明。还可以理解的是，所公开的特定概念仅仅是说明性的，并不意在限制本发明的范围，本发明的范围由所附权利要求书及其等同方案的全部范围来决定。

所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

在流程图中表示或在此以其他方式描述的逻辑和/或步骤，例如，可以被认为是用于实现逻辑功能的可执行指令的定序列表，可以具体实现在任何计算机可读介质中，以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用，或结合这些指令执行系统、装置或设备而使用。就本说明书而言，“计算机可读介质”可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。

计算机可读介质的更具体的示例(非穷尽性列表)包括以下：具有一个或多个布线的电连接部(电子装置)、便携式计算机盘盒(磁装置)、随机存取存储器(RAM)、只读存储器(ROM)、可擦除可编辑只读存储器(EPROM或闪速存储器)、光纤装置以及便携式光盘只读存储器(CDROM)。另外，计算机可读介质甚至可以是可在其上打印所述程序的纸或其他合适的介质，因为可以例如通过对纸或其他介质进行光学扫描，接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得所述程序，然后将其存储在计算机存储器中。

应当理解，本发明的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中，多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。例如，如果用硬件来实现，和在另一实施方式中一样，可用本领域公知的下列技术中的任一项或他们的组合来实现：具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路，具有合适的组合逻辑门电路的专用集成电路，可编程门阵列(PGA)，现场可编程门阵列(FPGA)等。

在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。

尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。

以上是对本发明的较佳实施进行了具体说明，但本发明并不限于所述实施例，熟悉本领域的技术人员在不违背本发明精神的前提下还可做出种种的等同变形或替换，这些等同的变形或替换均包含在本发明权利要求所限定的范围内。

Claims (5)

1.一种基于高通滤波解耦的机器人零力控制方法，其特征在于，包括：

获取机械臂端部所受的合力；其中，所述合力为机械臂抓握物体时所述机械臂端部所受的第一作用力与操作者对所述机械臂端部的第二作用力的合力；

根据所述合力、机器人的质量特性系数、机器人的阻尼特性系数、机器人的弹簧特性系数以及所述机械臂的电机关节位置信息，计算得到待处理的输入信号；

根据预先构建的双输入单输出状态空间表达式确定对所述待处理的输入信号的高通滤波截止频率；

根据所述高通滤波截止频率对所述待处理的输入信号进行高通滤波，以消除所述待处理的输入信号中的所述第一作用力，得到所述第二作用力；

所述根据所述合力、机器人的质量特性系数、机器人的阻尼特性系数、机器人的弹簧特性系数以及所述机械臂的电机关节位置信息，计算得到待处理的输入信号，包括：

根据第一表达式计算得到待处理的输入信号，所述第一表达式为：

其中，F_N表示所述待处理的输入信号，F_c表示所述合力，M表示所述质量特性系数，B表示所述阻尼特性系数，K表示所述弹簧特性系数，x表示所述电机关节位置信息，与/>分别表示所述电机关节位置信息的一阶微分与二阶微分；

所述根据预先构建的双输入单输出状态空间表达式确定对所述待处理的输入信号的高通滤波截止频率，包括：

根据所述双输入单输出状态空间表达式确定所述待处理的输入信号中位置输入系统的系统矩阵和所述待处理的输入信号中力矩输入系统的系统矩阵；

所述双输入单输出状态空间表达式为：

其中，x_R和x_N分别表示所述位置输入系统和所述力矩输入系统的状态向量，A_R和A_N分别表示所述位置输入系统和所述力矩输入系统的系统矩阵，b_R和b_N分别表示所述位置输入系统和所述力矩输入系统的控制矩阵，c_R和c_N表示所述位置输入系统和所述力矩输入系统的输出矩阵，y为所述双输入单输出状态空间表达式的输出；

根据由所述系统矩阵的特征值确定的最小根，确定所述高通滤波截止频率；

所述系统矩阵、所述控制矩阵以及所述输出矩阵分别表示为：

其中，各个所述矩阵的参数如下：

A_i＝(m+M)L,i＝{R N}

B_i＝(f+mω_c+B)L+(m+M)R,i＝{R N}

C_i＝(fω_c+K)L+(mω_c+f+B)R+C_kK_t,i＝{R N}

D_i＝(fω_c+K)R+C_kK_tω_c+K_tK_p,i＝{R N}

E_i＝K_tK_pω_c,i＝{R N}

F_R＝K_tK_p

F_N＝-L

G_R＝K_tK_pω_c

G_N＝-R

其中，m表示转子质量，M表示所述质量特性系数，B表示所述阻尼特性系数，K表示所述弹簧特性系数，f表示摩擦系数，K_t表示推力常数，K_p表示PID控制器系数，C_k表示反电势常数，R表示电阻，L表示电感；

所述根据由所述系统矩阵的特征值确定的最小根，确定所述高通滤波截止频率，包括：

将所述最小根表示为x₁，确定所述高通滤波截止频率的取值范围为0<ω_c<x₁，其中，ω_c表示所述高通滤波截止频率；

所述根据所述高通滤波截止频率对所述待处理的输入信号进行高通滤波，以消除所述待处理的输入信号中的所述第一作用力，得到所述第二作用力，包括：

根据高通滤波器的第二表达式对所述待处理的输入信号进行高通滤波；

所述第二表达式为：

其中，s表示拉氏变量，ω_c表示所述高通滤波截止频率。

2.根据权利要求1所述的一种基于高通滤波解耦的机器人零力控制方法，其特征在于，所述获取机械臂端部所受的合力，包括：

通过设置于所述机械臂端部的六轴传感器获取所述机械臂端部所受的合力。

3.一种基于高通滤波解耦的机器人零力控制装置，其特征在于，应用于如权利要求1或2所述的一种基于高通滤波解耦的机器人零力控制方法，所述装置包括：

合力获取单元，用于获取机械臂端部所受的合力；其中，所述合力为机械臂抓握物体时所述机械臂端部所受的第一作用力与操作者对所述机械臂端部的第二作用力的合力；

信号计算单元，用于根据所述合力、机器人的质量特性系数、机器人的阻尼特性系数、机器人的弹簧特性系数以及所述机械臂的电机关节位置信息，计算得到待处理的输入信号；

频率确定单元，用于根据预先构建的双输入单输出状态空间表达式确定对所述待处理的输入信号的高通滤波截止频率；

合力解耦单元，用于根据所述高通滤波截止频率对所述待处理的输入信号进行高通滤波，以消除所述待处理的输入信号中的所述第一作用力，得到所述第二作用力。

4.一种电子设备，其特征在于，包括处理器以及存储器；

所述存储器用于存储程序；

所述处理器执行所述程序实现如权利要求1或2所述的一种基于高通滤波解耦的机器人零力控制方法。

5.一种计算机可读存储介质，其特征在于，所述存储介质存储有程序，所述程序被处理器执行实现如权利要求1或2所述的一种基于高通滤波解耦的机器人零力控制方法。