CN103364579B - 一种用于乒乓球机器人的乒乓球自旋角速度的测量方法和测量装置 - Google Patents
一种用于乒乓球机器人的乒乓球自旋角速度的测量方法和测量装置 Download PDFInfo
- Publication number
- CN103364579B CN103364579B CN201310272639.0A CN201310272639A CN103364579B CN 103364579 B CN103364579 B CN 103364579B CN 201310272639 A CN201310272639 A CN 201310272639A CN 103364579 B CN103364579 B CN 103364579B
- Authority
- CN
- China
- Prior art keywords
- omega
- prime
- rightarrow
- partiald
- ping
- 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.)
- Active
Links
Landscapes
- Image Analysis (AREA)
Abstract
本发明是一种用于乒乓球机器人的乒乓球自旋角速度的预测方法,建立了乒乓球飞行轨迹的动力学模型,利用该模型分析由视觉捕捉系统得到的乒乓球飞行轨迹,并由此来计算乒乓球的自旋角速度,完善现有的乒乓球飞行轨迹的预测方法,使乒乓球机器人可以在进行乒乓球对打时提高成功率。
Description
技术领域
本发明涉及一种用于乒乓球机器人的乒乓球飞行自旋角速度的计算,属于机器人技术领域。
背景技术
机器人进行快速击球作业时,需要针对事先难以预料的来球轨迹与环境因素进行瞬时交互。为了使机器人能完成有效的快速回球,使球越过球网装置,再触及对方台区,需要尽可能地快速、准确获取环境信息。在乒乓球作业过程中,中等飞行速度的球,从对手击球至球飞抵机器人能击球空间的时间约为0.5秒左右,而机器人手臂驱动球拍从静止加速到一定的速度,与球轨迹相交接触需要一定的时间,因此机器人必须尽可能早地检测到运动球的轨迹,尽早地对球将来的空中飞行、碰撞和反弹的轨迹进行预测,而不能等到球飞行至手臂的作业空间内再开始动作,这些是实现机器人有效回球的必要前提和基础。
对于乒乓球轨迹的预测,乒乓球的旋转测量是一个比较重要并且比较难的问题。在现有的测量方法中,精确度高的方法是通过在乒乓球上画网格,然后用高速摄像机采集乒乓球在飞行中的网格速度的变化,通过对所采集的数据进行处理得到乒乓球的旋转速度。这种方法虽然精度较高,但是不具有实时性,不能满足乒乓球机器人对于乒乓球轨迹快速识别的需求。因而,现有的乒乓球机器人在进行打乒乓球作业时,无法处理旋转中的乒乓球。不考虑接发弧旋球。如果对方打弧旋球,通常采用光木板球拍,对旋转球不敏感。
现有专利201110136109.4提供了一种获取乒乓球机器人球拍击球姿态和击球速度的方法,该方法能够通过预测乒乓球落点位置及乒乓球落点速度来计算机器人击球姿态,使机器人实现击打乒乓球。并不考虑到乒乓球的旋转问题。
现有专利200910095291.6公开了一种精确测量和预测乒乓球轨迹系统与系统运行方法,该方法可以测量并能够实时的在显示屏上回放乒乓球的轨迹,该方法只是实时的测量乒乓球的运行轨迹,也没有给出乒乓球飞行中旋转速度的相关参数。
本专利提出的飞行中的乒乓球自旋角速的检测方法,可以通过机器人的视觉捕捉系统来获得乒乓球飞行的轨迹,建立乒乓球飞行的动力学模型,并利用该模型分析乒乓球飞行的轨迹,进而机算出乒乓球在飞行中的自旋角速度,使预测出乒乓球的轨迹更加精确,从而使得乒乓球机器人快速击球的成功率大大的提高。
发明内容
本发明提出了一种用于测量飞行中的对象的自旋角速度的方法,为乒乓球机器人进行打乒乓球作业提供更加精确的乒乓球轨迹预测,使乒乓球机器人提高快速击球的成功率。所述的技术方案如下:
一方面,本发明提供了一种自旋角速度的测量方法,包括:
(1)获取一系列的被测量对象在绝对坐标系中的四维坐标,包括3个方向的位置信息和一个时间信息;
(2)利用所获取的数据得到所述对象飞行的速度,并用最小二乘法拟合所述对象沿3个方向飞行轨迹的3次多项式;
(3)从上一步骤得到的3次多项式和速度信息,得到所述对象的加速度分量;
(4)对加速度方程式两边求导,并联合速度方程的行列式得到所述对象的自旋角速度ωx,ωy,ωz。
另一方面,本发明提供了一种自旋角速度的测量装置,包括:
视觉捕捉系统,能够获取一系列的被测量对象在绝对坐标系中的四维坐标,包括3个方向的位置信息和一个时间信息;
视觉处理系统,能够利用所获取的数据得到所述对象飞行的速度,并用最小二乘法拟合所述对象沿3个方向飞行轨迹的3次多项式;
所述视觉处理系统还能够利用得到的3次多项式和速度信息,得到所述对象的加速度分量;以及
对加速度方程式两边求导,并联合速度方程的行列式得到所述对象的自旋角速度ωx,ωy,ωz。
本发明提供的技术方案的有益效果是:
通过视觉捕捉系统得到的乒乓球飞行的轨迹,对该轨迹的简单分析,即可以有效地预测出乒乓球飞行的自旋角速度,可以实时的为机器人提供乒乓球轨迹数据,是机器人更加顺利的完成打乒乓球的作业。
附图说明
图1是本专利所提出的计算运行中乒乓球轨迹的方法流程图。
图2是本专利所用到的乒乓球台的详细参数以及模型的坐标系的建立位置的示意图。
图3是本发明所用到的乒乓球机器人进行乒乓球作业时的系统示意图。
图4是本发明所用到的乒乓球机器人的视觉识别系统的示意图。
图5是本发明施例中机器人完成击打乒乓球作业的系统流程图。
具体实施方式:
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
本发明所用的乒乓球、球台、球拍、球网采用与《乒乓球竞赛规则》中规定的相同。其详细参数如下:
(1)球台外形为与水平面平行的长方形,长2.74m,宽1.525m,台面离地面高0.76m。台面具有一致的弹性,当标准球从离台面0.30m高处落至台面时,弹起高度约为0.23m,参见图2。
(2)球网悬挂在一根绳子上,整个球网的顶端距离比赛台面0.1525m。
(3)乒乓球直径为40mm,球的质量为2.7g。
(4)球拍为椭圆形,短轴长约为0.15m,长轴长约为17cm。选用普通球拍。
建立机器人作业系统的世界坐标系,坐标原点设在球台的一角,沿着乒乓球台“长”的方向为+X(朝着机器人方向为正方向,离球桌中线300mm),沿着乒乓球台“宽”的方向为+Y,由右手法则,确定+Z(正方向与重力加速度方向相反)。
机器人的进行打乒乓球作业的整体系统组成如下:
天花板上有高速立体视觉,处理速度100幅/秒,对5米/秒运动的球,分辨率为5厘米一个球图,高速图像处理系统可以实时计算出球的位置与速度,并根据动力学方程计算未来轨迹、落点与反弹轨迹,将计算结果通过无线通讯传给机器人,机器人预测球未来的轨迹。机器人根据自己位置、姿态,计算拦截点,进行运动规划,控制全身上下肢运动击球,参见图3。
机器人的视觉捕捉系统主要组成如下所示:
系统包括两个同型号的黑白智能相机,分别固定于标准乒乓球台(2.74m×1.525m)左上方和右上方对称位置的天花板上,它的视频输出线直接连接到显示器,可以实时看到相机采集的图像。智能相机集成了高性能数字处理芯片DSP与FPGA,具有较强的图像处理能力,最 高采集处理能力能够达到250帧/秒。而且具有串口、以太网接口等,并且支持TCP/IP协议。
两台智能摄像机A、B与个人电脑通过以太网接口构成局域网。其中智能相机A为服务器,智能相机B与个人电脑为客户端。当对面球台的球飞来时,两台智能相机同步采集乒乓球的轨迹,并将处理结果发送给个人电脑。个人电脑利用离线标定好的信息,完成乒乓球的三维轨迹计算。参见图4。
该方法得到乒乓球轨迹的主要过程为:
机器人的视觉捕捉系统从球飞过预测线时起,获取一系列的乒乓球在绝对坐标系中的四维坐标(xi,yi,zi,ti),i为序列号。当获取了6组以上乒乓球的位置坐标后,分别对坐标的X、Y、Z分量用最小二乘法拟合乒乓球沿X、Y、Z轴飞行轨迹的3次多项式(即取m=3,n≥6):
根据空气动力学,在乒乓球飞行阶段,乒乓球受主要受重力、空气阻力、马格努斯力(Magus force)的作用。空气阻力与球速的平方关系,马格努斯力与自旋速度和球速的平方成正比。一般地,只要球拍对乒乓球球的作用力不通过球的质心,球就会产生自旋,产生马格努斯效应。记球的速度矢量为自旋矢量重力加速度为则在飞行阶段乒乓球加速度的公式为:
式中,Cd是空气阻力系数,Cm是马格努斯系数。在室内条件情况下,对于乒乓球而言,Cd的取值大约为0.12,Cm约为4×10-5。式(1.2)右端第二项,矢量与是叉乘关系。式(1.1)是一个非线性耦合的微分方程,不易获得解析解。
对于给定速度分量vx,vy,vz及自旋向量分量ωx,ωy,ωz,依据式(1.2)可得作用于乒乓 球的加速度分量ax,ay,az。其表达式如下式:
如果直接将拟合得到的初始速度及加速度代入式(1.3),不能求解出初始自旋角速度ωx0、ωy0、ωz0。方程式(1.3)是奇异的,因为下面的行列式等于零。
实际上,可以利用式(1.3)所得的加速度率来求解角速度ωx0、ωy0、ωz0。
对等式ax=-Cdvvx+Cmv(ωyvz-vyωz)两边时间t求导,整理后得:
其中,ax'=jx0,
联合式(1.4)和(1.5),得:
将ωz代入av=-Cdvvy+Cmv(ωzvx-vzωx),求解得:
乒乓球在飞行期间,作用在球上的力矩可以忽略,因此,认为乒乓球在整个飞行期间自旋角速度不变。至此我们利用乒乓球飞行的轨迹得到了乒乓球飞行中的三个方向上的自旋角速度。乒乓球在飞行中的基本参数都已经通过我们的专利提出的方法得到了。
将得到的乒乓球飞行的预测轨迹以及自旋角速度发送给乒乓球机器人,乒乓球机器人将其和自己的身体位置与手臂位置进行对比,计算球拍最快可以达到的地点,将其作为预定拦截点。同时根据在拦截点的球速、球向,我方击球后预定的落点,计算自己球拍以什么速度与角度到达拦截点。并在运动过程中不停的根据自己的实际位置与球的实际位置,调节拦截点,最后精确到达拦截点,将球击回。机器人完成击打乒乓球作业的系统流程图参见图5。
Claims (2)
1.一种自旋角速度的测量方法,包括:
(1)获取一系列的被测量对象在绝对坐标系中的四维坐标(xi,yi,zi,ti),其中i为序列号,包括3个方向的位置信息和一个时间信息;
(2)当获取了6组以上所述对象的位置坐标后,利用所获取的数据得到所述对象飞行的速度v,并用最小二乘法拟合所述对象沿3个方向飞行轨迹的3次多项式:
所述被测对象的加速度公式为:
其中为所述对象的速度矢量,为自旋矢量,为重力加速度,Cd是空气阻力系数,Cm是马格努斯系数;
(3)对于给定速度分量vx,vy,vz及自旋角速度ωx,ωy,ωz,得到作用于所述对象的加速度分量ax,ay,az,其表达式如下式:
(4)对加速度方程式两边求导,并联合速度方程的行列式得到所述对象的自旋角速度ωx,ωy,ωz,过程如下:
对等式ax=-Cdvvx+Cmv(ωyvz-vyωz)两边时间t求导,整理后得:
其中,ax'=jx0,
进一步求解得:
即
将ωz代入av=-Cdvvy+Cmv(ωzvx-vzωx),求解得:
2.一种自旋角速度的测量装置,包括:
视觉捕捉系统,能够获取一系列的被测量对象在绝对坐标系中的四维坐标(xi,yi,zi,ti),其中i为序列号,包括3个方向的位置信息和一个时间信息;
视觉处理系统,能够利用所获取的数据得到所述对象飞行的速度v,并用最小二乘法拟合所述对象沿3个方向飞行轨迹的3次多项式:
记所述对象的速度矢量为自旋矢量重力加速度为Cd是空气阻力系数,Cm是马格努斯系数,则所述对象加速度的公式为:
所述视觉处理系统还能够利用得到的3次多项式和速度信息,对于给定速度分量vx,vy,vz及自旋角速度ωx,ωy,ωz,得到作用于所述对象的加速度分量ax,ay,az:
以及
对加速度方程式两边求导,并联合速度方程的行列式得到所述对象的自旋角速度ωx,ωy,ωz,过程如下:
对等式ax=-Cdvvx+Cmv(ωyvz-vyωz)两边时间t求导,整理后得:
其中,ax'=jx0,
进一步求解得:
即
将ωz代入av=-Cdvvy+Cmv(ωzvx-vzωx),求解得:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310272639.0A CN103364579B (zh) | 2013-07-02 | 2013-07-02 | 一种用于乒乓球机器人的乒乓球自旋角速度的测量方法和测量装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310272639.0A CN103364579B (zh) | 2013-07-02 | 2013-07-02 | 一种用于乒乓球机器人的乒乓球自旋角速度的测量方法和测量装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103364579A CN103364579A (zh) | 2013-10-23 |
CN103364579B true CN103364579B (zh) | 2015-04-15 |
Family
ID=49366388
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310272639.0A Active CN103364579B (zh) | 2013-07-02 | 2013-07-02 | 一种用于乒乓球机器人的乒乓球自旋角速度的测量方法和测量装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103364579B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105373128B (zh) * | 2015-06-17 | 2018-06-12 | 电子科技大学 | 羽毛球机器人及其全场定位方法 |
CN106325302B (zh) * | 2015-06-17 | 2020-03-13 | 电子科技大学 | 一种移动物体位置跟踪平台、跟踪装置及跟踪方法 |
CN106390409B (zh) * | 2016-11-17 | 2018-10-02 | 中科新松有限公司 | 一种乒乓球机器人的击球方法及装置 |
CN109807888A (zh) * | 2019-01-24 | 2019-05-28 | 哈尔滨理工大学 | 一种基于并联机械臂的类冰球对决方法 |
CN109901581A (zh) * | 2019-03-15 | 2019-06-18 | 智久(厦门)机器人科技有限公司上海分公司 | 一种agv车自旋角的标定方法及自旋运动控制方法 |
CN110941795B (zh) * | 2019-12-16 | 2023-05-12 | 上海创屹科技有限公司 | 乒乓球旋转角度获取方法、获取装置及存储介质 |
CN111721962B (zh) * | 2020-07-17 | 2021-05-14 | 厦门大学 | 一种基于马格努斯效应的流速测量方法 |
CN113362366B (zh) * | 2021-05-21 | 2023-07-04 | 上海奥视达智能科技有限公司 | 一种球体转速的确定方法、装置、终端和存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101458434A (zh) * | 2009-01-08 | 2009-06-17 | 浙江大学 | 精确测量和预测乒乓球轨迹系统与系统运行方法 |
CN102200760A (zh) * | 2011-05-25 | 2011-09-28 | 中国科学院自动化研究所 | 一种获取乒乓球机器人球拍击球姿态和击球速度的方法 |
CN103019024A (zh) * | 2012-11-29 | 2013-04-03 | 浙江大学 | 实时精确观测和分析乒乓球旋转系统与系统运行方法 |
-
2013
- 2013-07-02 CN CN201310272639.0A patent/CN103364579B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101458434A (zh) * | 2009-01-08 | 2009-06-17 | 浙江大学 | 精确测量和预测乒乓球轨迹系统与系统运行方法 |
CN102200760A (zh) * | 2011-05-25 | 2011-09-28 | 中国科学院自动化研究所 | 一种获取乒乓球机器人球拍击球姿态和击球速度的方法 |
CN103019024A (zh) * | 2012-11-29 | 2013-04-03 | 浙江大学 | 实时精确观测和分析乒乓球旋转系统与系统运行方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103364579A (zh) | 2013-10-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103364579B (zh) | 一种用于乒乓球机器人的乒乓球自旋角速度的测量方法和测量装置 | |
CN104504694B (zh) | 一种获取运动球体三维信息的方法 | |
CN103389738B (zh) | 一种用于乒乓球机器人预测乒乓球轨迹的方法和装置 | |
CN106390409B (zh) | 一种乒乓球机器人的击球方法及装置 | |
CN101964047B (zh) | 一种基于多跟踪点的人体动作识别方法 | |
CN110610547B (zh) | 基于虚拟现实的座舱实训方法、系统及存储介质 | |
CN104548555A (zh) | 运动分析方法、运动分析装置 | |
CN106582012A (zh) | 一种vr场景下的攀爬操作处理方法和装置 | |
CN105536233B (zh) | 一种用于田径场上跑步速度的采集测量系统 | |
Huang et al. | Trajectory prediction of spinning ball for ping-pong player robot | |
CN108888918B (zh) | 一种用于复杂路径下多目标运动速度测量系统及方法 | |
Goff et al. | Creating drag and lift curves from soccer trajectories | |
CN105069829A (zh) | 一种基于多目视频的人体动画生成方法 | |
CN113386126B (zh) | 旋转推断装置及方法、轨迹预测装置及方法、回球机器人 | |
Liu et al. | Racket control and its experiments for robot playing table tennis | |
Liu et al. | Accurate real-time ball trajectory estimation with onboard stereo camera system for humanoid ping-pong robot | |
KR101462079B1 (ko) | 골프 시뮬레이터 | |
CN105879349A (zh) | 显示屏上显示高尔夫球落在果岭上的位置的方法及系统 | |
TWI764768B (zh) | 桌球軌跡預測方法及其系統 | |
CN113041587B (zh) | 一种基于感应台面和球拍面采集兵乓球击打轨迹的方法、装置、计算机设备及存储介质 | |
CN107754286A (zh) | 一种基于惯性导航与触感校正定位的乒乓球运动系统 | |
Sun et al. | A novel ball-racket rebound model for table tennis robot based on continuous contact force | |
Lluna et al. | Measurement of aerodynamic coefficients of spherical objects using an electro-optic device | |
CN107823864B (zh) | 一种基于运动轨迹的乒乓球对局发球识别系统及方法 | |
Zhang et al. | The design and physical implementation of tennis training robots |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |