CN113110542B - 一种基于单轴云台的四旋翼无人机目标追踪系统 - Google Patents

一种基于单轴云台的四旋翼无人机目标追踪系统 Download PDF

Info

Publication number
CN113110542B
CN113110542B CN202110411015.7A CN202110411015A CN113110542B CN 113110542 B CN113110542 B CN 113110542B CN 202110411015 A CN202110411015 A CN 202110411015A CN 113110542 B CN113110542 B CN 113110542B
Authority
CN
China
Prior art keywords
coordinate system
aerial vehicle
unmanned aerial
target
control
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.)
Expired - Fee Related
Application number
CN202110411015.7A
Other languages
English (en)
Other versions
CN113110542A (zh
Inventor
王雪飞
马宏军
陈豹
宋大雷
白桦
徐少杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northeastern University China
Original Assignee
Northeastern University China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northeastern University China filed Critical Northeastern University China
Priority to CN202110411015.7A priority Critical patent/CN113110542B/zh
Publication of CN113110542A publication Critical patent/CN113110542A/zh
Application granted granted Critical
Publication of CN113110542B publication Critical patent/CN113110542B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
    • G05D1/08Control of attitude, i.e. control of roll, pitch, or yaw
    • G05D1/0808Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/101Simultaneous control of position or course in three dimensions specially adapted for aircraft

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明公开了一种基于单轴云台的四旋翼无人机目标追踪系统,包括机架,所述机架的上端活动安装有电机,所述电机的上端活动安装有螺旋桨,所述机架的上端固定安装有JetsonTX2板载计算机,所述JetsonTX2板载计算机的上端固定安装有pixhawk飞行控制器,所述机架的下端固定安装有单轴云台,所述单轴云台的下端活动安装有摄像机,该基于单轴云台的四旋翼无人机目标追踪系统,在算法中设置一个虚拟平面估计目标倾角,准确照在目标平面上使目标平面和摄像头成像平面平行,在当目标平面发生倾斜的时候,无人机可以自行判定目标倾角的大小然后调整云台的角度对目标进行追踪,以获得同等高度下目标的最大图像。在追踪侦测领域有巨大作用。也可在航拍领域进行广泛推广。

Description

一种基于单轴云台的四旋翼无人机目标追踪系统
技术领域
本发明涉及无人机技术领域,具体为一种基于单轴云台的四旋翼无人机目标追踪系统。
背景技术
随着嵌入式技术,视觉处理技术和空气动力学领域的不断发展,无人机迎来了一个蓬勃发展的黄金时期,由仅仅用于飞行试验和航拍等任务逐渐扩展到了货物配送,电路设施检测,人员营救,环境探测,军事活动等领域,无人机的优势在于其活动灵活,具有空中优势还可搭配传感器获取更多信息。
无人机搭配视觉传感器和激光传感器等可以充分发挥其本身机动性强,获取信息多的特点。旋翼无人机较固定翼无人机而言,其机动性和灵活性都有了很大的提升,除此之外,它能够实现自由悬停、垂直起降等比较难的任务,旋翼无人机还可以实现低速低空平稳飞行,能够在隐蔽的空间里面执行任务、在很小的范围进行检测,因此旋翼无人机有更多的优势。
无人机最大的优势是搭配视觉传感器获取视觉信息对周围环境进行分析定位。视觉追踪是无人机应用的一个比较盛行的方向,可以跟踪目标进行拍摄或者抓取。
发明内容
本发明的目的在于提供一种基于单轴云台的四旋翼无人机目标追踪系统,以解决上述背景技术中提出针对追踪技术不完善,只能框选出目标对其位置进行追踪无法判别目标平面是否发生变化,对斜面上的目标追踪效果不佳的问题。
为实现上述目的,本发明提供如下技术方案:一种基于单轴云台的四旋翼无人机目标追踪系统,其步骤如下:
步骤S1:在此方法中我们一共使用四个坐标系,分别是世界坐标系(惯性坐标系)I={Oi,Xi,Yi,Zi},机体坐标系B={Ob,Xb,Yb,Zb},相机坐标系C={Oc,Xc,Yc,Zc},虚拟坐标系V={Ov,Xv,Yv,Zv},我们假设相机坐标系和机体坐标系的原点重合。虚拟坐标系下的虚拟图像平面是始终水平的。在机体坐标系下的无人机的运动学方程如下:
Figure GDA0003680379960000021
Figure GDA0003680379960000022
Figure GDA0003680379960000023
Figure GDA0003680379960000024
F=-U1E3+mgRTe3
其中ζ代表的是世界坐标系下的位移,R表示的是从机体坐标系到世界坐标系的转换矩阵。上述的第二个公式表示的是转换矩阵的求导公式。第三个公式表示机体坐标系下无人机的加速度
Figure GDA0003680379960000025
m表示的是无人机的质量,Ω表示无人机的角速度,F是机体坐标系下无人机所受的合力,J是无人机的转动惯量。第四个公式是无人机旋转方向上的运动公式,τ是作用在无人机上的转矩。第五个公式表示的是无人机在机体坐标系下所受到的合力。
步骤S2:在虚拟坐标系计算完偏角并且云台调整完之后在相机坐标系下基于backstepping控制法的无人机控制率如下:
Figure GDA0003680379960000026
一共经过三层的推导,第一层是以速度为控制量,第二层是以力为控制量,第三层是以力的导为控制量。上式是经过backstepping推导计算得到的最终结果。其中
Figure GDA0003680379960000027
表示的是推力的导数,sk(·)表示反对称矩阵e3表示的是三维列向量[0,0,1]T
Figure GDA0003680379960000028
并且q1,q2,q3是backstepping方法中每一步的控制量误差。
Figure GDA0003680379960000029
是构建速度观测器时q1的测量值和估计值之差。z*是高度的期望值。c1,k1,k2,k3,k4,k5是控制参数。偏航控制的控制率为
Figure GDA0003680379960000031
其中Ω2和Ω3都是
步骤S3:所使用的虚拟平面是水平的一个平面,世界坐标系,相机坐标系,虚拟坐标系,机体坐标系的关系如下:
I=RPRψV=RψC=RψRθRφB
其中I是世界坐标系,B是机体坐标系,C是相机坐标系,V是虚拟坐标系。相机坐标系和世界坐标系之间平行只相差一个偏航的误差,RP是虚拟坐标系与相机坐标系之间的变换矩阵只有一个俯仰变化。
步骤S4:判别世界坐系下目标平面的偏角其方法如下:
如图1,追踪的目标至少可以是四个特征点,假设这四个特征点分别是位于正方形的四个角上。当目标处于水平的平面上,虚拟坐标系下四个特征点也听该是两两之间距离相等的。目标点在虚拟平面下的横纵坐标分别是u1,v1,u2,v2,u3,v3,u4,v4。首先求虚拟平面下四个边的长度:
Figure GDA0003680379960000032
Figure GDA0003680379960000033
然后比较相对的一对边与另一对边的长短,即比较A+C和B+D。当A+C大于B+D则表示目标平面是在俯仰方向产生了倾斜,反之则是图像平面在横滚方向产生了倾斜。至于目标倾斜角的正负,由于通过摄像机透视法则,离摄像头近的在图像平面中会比里的远的大。则是需要直接比较相对的两个边的大小。当A>C时表示目标的偏转角为正,此时控制单轴云台向机头方向偏转,相机平面与目标平面平行。角度的计算是
Figure GDA0003680379960000034
步骤S5对于无人机控制信息获取如下:在相机坐标系下的特征点位置信息是:
Figure GDA0003680379960000041
对于位置信息的微分如下:
Figure GDA0003680379960000042
其中
Figure GDA0003680379960000043
是相机坐标系相对于惯性坐标系的线速度。v=[cvx cvy cvz]是在相机坐标系下的线速度。根据相机透视法得到像素坐标系下点的坐标:
Figure GDA0003680379960000044
可以根据上式得到像素坐标的导如下:
Figure GDA0003680379960000045
然后计算特征点坐标的平均值,并且使用图像矩来控制无人机的追踪动作。特征点平均坐标是:
Figure GDA0003680379960000046
a=cμ20+cμ02,设置a*是a的理想值。控制位移和偏航的图像矩为:
Figure GDA0003680379960000047
Figure GDA0003680379960000048
步骤S6在backstepping控制法下进行仿真实验,验证追踪方法的有效性。
一种基于单轴云台的四旋翼无人机,包括机架(4),所述机架(4)的上端活动安装有电机(5),所述电机(5)的上端活动安装有螺旋桨(6),所述机架(4)的上端固定安装有Jetson TX2板载计算机(2),所述Jetson TX2板载计算机(2)的上端固定安装有pixhawk飞行控制器(1),所述机架(4)的下端固定安装有单轴云台(7),所述单轴云台(7)的下端活动安装有摄像机(3)。
优选的,所述无人机机体坐标系下的运动方程其特征在于为保证准确追踪无人机需尽量保证飞行的稳定性。在水平的虚拟坐标系下所受的合力为f=-RφθU1E3+mge3,其中U1是所有螺旋桨产生的无人机的推力,Rφθ是机体和摄像头一起产生的偏转造成的。E3和e3都是三维列向量[0 0 1]T。在相机坐标系下的合力为
Figure GDA0003680379960000051
其中R是机体的偏转矩阵,Rp是测得的目标平面的倾角。
优选的,所述在S2解算出来的推力约束是如下
Figure GDA0003680379960000052
可以对合力
Figure GDA0003680379960000053
求导来获得的
Figure GDA0003680379960000054
可以得到如下:
Figure GDA0003680379960000055
其中Ω1和Ω2是无人机的俯仰、横滚角速度。根据求得的姿态角速度Ω1、Ω2然后根据公式:
Figure GDA0003680379960000056
得到惯性坐标系下角度的导
Figure GDA0003680379960000057
然后使用PD控制求得角度
Figure GDA0003680379960000058
对于位置控制使用如下的无人机动力学模型:
Figure GDA0003680379960000059
Figure GDA00036803799600000510
Figure GDA00036803799600000511
首先通过PD控制求得线速度,然后再通过PD控制求得位置。
优选的,所述对于偏航的控制还是在水平虚拟平面中进行偏航控制。定义偏航的图像矩为如下:
Figure GDA0003680379960000061
对其求导
Figure GDA0003680379960000062
设偏航角的理想值是
Figure GDA0003680379960000063
定义特征误差
Figure GDA0003680379960000064
对其求导有
Figure GDA0003680379960000065
选择偏航角的控制输入
Figure GDA0003680379960000066
使用如下公式:
Figure GDA0003680379960000067
获得最终的控制
Figure GDA0003680379960000068
优选的,所述在出事状态下单轴云台是垂直向下的即角度为0,每一次的调整量是通过当前时刻计算出的角度和前一时刻的角度相减得到控制量。
优选的,所述单轴云台使用的是步进电机,由pixhawk飞行控制器的pwm输出口进行控制。
与现有技术相比,本发明的有益效果是:该基于单轴云台的四旋翼无人机目标追踪系统,该方法在当目标平面发生倾斜的时候无人机可以自行判定目标倾角的大小然后调整云台的角度对目标进行追踪,以获得同等高度下目标的最大图像,在追踪侦测等领域有巨大作用,也可在航拍领域进行广泛的推广,直接控制无人机姿态和位置的是pixhawk飞行控制器,由pixhawk飞行控制器直接连接电调控制四个电机的转速进而控制无人机的姿态和位置,USB摄像头和TX2板载计机直接连接获取图像信息,TX2板载计算机负责特征点提取和图像的处理,计算控制量传达给pixhawk飞行控制器,TX2板载计算机和pixhawk飞行控制器通过mavros进行通信,在算法中设置一个虚拟平面估计目标倾角,将估计到的倾角信息传达到pixhawk飞行控制器控制步进电机调整摄像头角度,准确照在目标平面上使目标平面和摄像头成像平面平行,可以拍摄到更加准确形变更少的图像。
附图说明
图1为本发明流程框图;
图2为本发明整体结构示意图;
图3为本发明在追踪过程中qx,qy,qz,qψ四个控制量的变化示意图;
图4为本发明在追踪过程中无人机自身位置x,y,z的变化示意图;
图5为本发明在追踪过程中识别的目标角度的变化示意图。
图中:1、pixhawk飞行控制器;2、Jetson TX2板载计算机;3、摄像机;4、机架;5、电机;6、螺旋桨;7、单轴云台。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图1-5,本发明提供一种技术方案:一种基于单轴云台的四旋翼无人机目标追踪系统,其特征在于:其步骤如下:
步骤S1:在此方法中我们一共使用四个坐标系,分别是世界坐标系(惯性坐标系)I={Oi,Xi,Yi,Zi},机体坐标系B={Ob,Xb,Yb,Zb},相机坐标系C={Oc,Xc,Yc,Zc},虚拟坐标系V={Ov,Xv,Yv,Zv},我们假设相机坐标系和机体坐标系的原点重合。虚拟坐标系下的虚拟图像平面是始终水平的。在机体坐标系下的无人机的运动学方程如下:
Figure GDA0003680379960000081
Figure GDA0003680379960000082
Figure GDA0003680379960000083
Figure GDA0003680379960000084
F=-U1E3+mgRTe3
其中ζ代表的是世界坐标系下的位移,R表示的是从机体坐标系到世界坐标系的转换矩阵。上述的第二个公式表示的是转换矩阵的求导公式。第三个公式表示机体坐标系下无人机的加速度
Figure GDA0003680379960000085
m表示的是无人机的质量,Ω表示无人机的角速度,F是机体坐标系下无人机所受的合力,J是无人机的转动惯量。第四个公式是无人机旋转方向上的运动公式,τ是作用在无人机上的转矩。第五个公式表示的是无人机在机体坐标系下所受到的合力。
步骤S2:在虚拟坐标系计算完偏角并且云台调整完之后在相机坐标系下基于backstepping控制法的无人机控制率如下:
Figure GDA0003680379960000086
一共经过三层的推导,第一层是以速度为控制量,第二层是以力为控制量,第三层是以力的导为控制量。上式是经过backstepping推导计算得到的最终结果。其中
Figure GDA0003680379960000087
表示的是推力的导数,sk(·)表示反对称矩阵e3表示的是三维列向量[0,0,1]T
Figure GDA0003680379960000088
并且q1,q2,q3是backstepping方法中每一步的控制量误差。
Figure GDA0003680379960000089
是构建速度观测器时q1的测量值和估计值之差。z*是高度的期望值。c1,k1,k2,k3,k4,k5是控制参数。偏航控制的控制率为
Figure GDA00036803799600000810
其中Ω2和Ω3都是
步骤S3:所使用的虚拟平面是水平的一个平面,世界坐标系,相机坐标系,虚拟坐标系,机体坐标系的关系如下:
I=RPRψV=RψC=RψRθRφB
其中I是世界坐标系,B是机体坐标系,C是相机坐标系,V是虚拟坐标系。相机坐标系和世界坐标系之间平行只相差一个偏航的误差,RP是虚拟坐标系与相机坐标系之间的变换矩阵只有一个俯仰变化。
步骤S4:判别世界坐系下目标平面的偏角其方法如下:
如图1,追踪的目标至少可以是四个特征点,假设这四个特征点分别是位于正方形的四个角上。当目标处于水平的平面上,虚拟坐标系下四个特征点也听该是两两之间距离相等的。目标点在虚拟平面下的横纵坐标分别是u1,v1,u2,v2,u3,v3,u4,v4。首先求虚拟平面下四个边的长度:
Figure GDA0003680379960000091
Figure GDA0003680379960000092
然后比较相对的一对边与另一对边的长短,即比较A+C和B+D。当A+C大于B+D则表示目标平面是在俯仰方向产生了倾斜,反之则是图像平面在横滚方向产生了倾斜。至于目标倾斜角的正负,由于通过摄像机透视法则,离摄像头近的在图像平面中会比里的远的大。则是需要直接比较相对的两个边的大小。当A>C时表示目标的偏转角为正,此时控制单轴云台向机头方向偏转,相机平面与目标平面平行。角度的计算是
Figure GDA0003680379960000093
步骤S5对于无人机控制信息获取如下:在相机坐标系下的特征点位置信息是:
Figure GDA0003680379960000094
对于位置信息的微分如下:
Figure GDA0003680379960000095
其中
Figure GDA0003680379960000101
是相机坐标系相对于惯性坐标系的线速度。v=[cvx cvy cvz]是在相机坐标系下的线速度。根据相机透视法得到像素坐标系下点的坐标:
Figure GDA0003680379960000102
可以根据上式得到像素坐标的导如下:
Figure GDA0003680379960000103
然后计算特征点坐标的平均值,并且使用图像矩来控制无人机的追踪动作。特征点平均坐标是:
Figure GDA0003680379960000104
a=cμ20+cμ02,设置a*是a的理想值。控制位移和偏航的图像矩为:
Figure GDA0003680379960000105
Figure GDA0003680379960000106
步骤S6在backstepping控制法下进行仿真实验,验证追踪方法的有效性。
一种基于单轴云台的四旋翼无人机,包括机架4,其特征在于:机架4的上端活动安装有电机5,电机5的上端活动安装有螺旋桨6,机架4的上端固定安装有Jetson TX2板载计算机2,Jetson TX2板载计算机2的上端固定安装有pixhawk飞行控制器,机架4的下端固定安装有单轴云台7,单轴云台的下端活动安装有摄像机3。
进一步的,无人机机体坐标系下的运动方程其特征在于为保证准确追踪无人机需尽量保证飞行的稳定性。在水平的虚拟坐标系下所受的合力为f=-RφθU1E3+mge3,其中U1是所有螺旋桨产生的无人机的推力,Rφθ是机体和摄像头一起产生的偏转造成的。E3和e3都是三维列向量[0 0 1]T。在相机坐标系下的合力为
Figure GDA0003680379960000111
其中R是机体的偏转矩阵,Rp是测得的目标平面的倾角;
进一步的,在S2解算出来的推力约束是如下
Figure GDA0003680379960000112
可以对合力
Figure GDA0003680379960000113
求导来获得的
Figure GDA0003680379960000114
可以得到如下:
Figure GDA0003680379960000115
其中Ω1和Ω2是无人机的俯仰、横滚角速度。根据求得的姿态角速度Ω1、Ω2然后根据公式:
Figure GDA0003680379960000116
得到惯性坐标系下角度的导
Figure GDA0003680379960000117
然后使用PD控制求得角度
Figure GDA0003680379960000118
对于位置控制使用如下的无人机动力学模型:
Figure GDA0003680379960000119
Figure GDA00036803799600001110
Figure GDA00036803799600001111
首先通过PD控制求得线速度,然后再通过PD控制求得位置;
进一步的,对于偏航的控制还是在水平虚拟平面中进行偏航控制。定义偏航的图像矩为如下:
Figure GDA00036803799600001112
对其求导
Figure GDA00036803799600001113
设偏航角的理想值是
Figure GDA00036803799600001114
定义特征误差
Figure GDA00036803799600001115
对其求导有
Figure GDA00036803799600001116
选择偏航角的控制输入
Figure GDA00036803799600001117
使用如下公式:
Figure GDA0003680379960000121
获得最终的控制
Figure GDA0003680379960000122
进一步的,在出事状态下单轴云台是垂直向下的即角度为0,每一次的调整量是通过当前时刻计算出的角度和前一时刻的角度相减得到控制量,即相机坐标系C和机体坐标系B是保持一致的[x y z]T的方向表示的非别是前右下。每一次的调整量是通过当前时刻计算出的调整角度和前一时刻的角度相减得到控制量;
进一步的,单轴云台使用的是步进电机,由pixhawk飞行控制器的pwm输出口进行控制,然后配合编码器可以测得改变后的角度对云台进行闭环控制。本发明提供了一种基于单轴云台的无人机目标追踪系统。首先对无人机进行建模,分析其动力系统,保证追踪过程中无人机的姿态保持稳定。其次设置一个虚拟平面原点与相机平面重合,虚拟平面通过惯导模块(IMU)测定的无人机的姿态和当前相机的姿态调整自身保持水平。在水平的虚拟坐标系下通过目标图像得到变化测定其偏转的角度,并且偏航控制在虚拟坐标系下完成,同时调整单轴云台使相机平面与目标平面平行,获得同等分飞行高度下最大目标图像。在虚拟平面和相机平面下都对图像矩进行计算利用图像矩进行目标的追踪。因为无人机具有非线性欠驱动的特点,本专利使用反步控制法解算出其控制量
Figure GDA0003680379960000123
Figure GDA0003680379960000124
最终使用PD控制实现无人机的位置和姿态调整。
工作原理:首先明确无人机在直线和角度上的动力学模型,然后建立世界坐标系,机体坐标系,相机坐标系,虚拟坐标系之间的关系,得到目标点的像素之后,计算图像矩当作被控量,然后使用反步控制计算出控制量,进而算得角速度控制量和加速度控制量,使用PD控制来控制无人机的姿态赫尔位置进行目标追踪,根据S1中,无人机机体坐标系下的运动方程其特征在于为保证准确追踪无人机需尽量保证飞行的稳定性。在水平的虚拟坐标系下所受的合力为f=-RφθU1E3+mge3,其中U1是所有螺旋桨产生的无人机的推力,Rφθ是机体和摄像头一起产生的偏转造成的。E3和e3都是三维列向量[0 0 1]T。在相机坐标系下的合力为
Figure GDA0003680379960000131
其中R是机体的偏转矩阵,Rp是测得的目标平面的倾角。
根据S2解算出来的推力约束是如下
Figure GDA0003680379960000132
可以对合力
Figure GDA0003680379960000133
求导来获得的
Figure GDA0003680379960000134
可以得到如下:
Figure GDA0003680379960000135
其中Ω1和Ω2是无人机的俯仰、横滚角速度。根据求得的姿态角速度Ω1、Ω2然后根据公式:
Figure GDA0003680379960000136
得到惯性坐标系下角度的导
Figure GDA0003680379960000137
然后使用PD控制求得角度
Figure GDA0003680379960000138
对于位置控制使用如下的无人机动力学模型:
Figure GDA0003680379960000139
Figure GDA00036803799600001310
Figure GDA00036803799600001311
首先通过PD控制求得线速度,然后再通过PD控制求得位置。
根据所述步骤S5中,对于偏航的控制还是在水平虚拟平面中进行偏航控制。定义偏航的图像矩为如下:
Figure GDA0003680379960000141
对其求导
Figure GDA0003680379960000142
设偏航角的理想值是
Figure GDA0003680379960000143
定义特征误差
Figure GDA0003680379960000144
对其求导有
Figure GDA0003680379960000145
选择偏航角的控制输入
Figure GDA0003680379960000146
使用如下公式:
Figure GDA0003680379960000147
获得最终的控制
Figure GDA0003680379960000148
根据所述步骤S3中基于单轴云台的无人机目标追踪系统,在初始状态下单轴云台是垂直向下的即角度为0,即相机坐标系C和机体坐标系B是保持一致的[x y z]T的方向表示的非别是前右下。每一次的调整量是通过当前时刻计算出的调整角度和前一时刻的角度相减得到控制量。
为更好的说明本系统对于无人机目标追踪具有较好的控制效果和对控制理论进行了仿真:设定无人机质量m=2Kg,重力加速度g=9.81m/s2。无人机的初始位置是[x y z]=[4,3,-7]。期望的高度为4m。转动惯量I1=0.0090kg.m2/rad2,I2=0.0090kg.m2/rad2,I2=0.0176kg.m2/rad2。相机的焦距f=0.0032mm,所追踪的四个点在世界坐标系下设为:a1=[0.25 0.5 0]T,a2=[-0.25 0.5 0]T,a3=[-0.25 -0.5 0.5]T,a4=[0.25 -0.5 0.5]T单位是m。所使用到的参数分别是c1=2,k1=4,k2=8,k3=2,k4=2,k5=0.7。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。

Claims (6)

1.一种基于单轴云台的四旋翼无人机目标追踪系统,包括机架(4),其特征在于:所述机架(4)的上端活动安装有电机(5),所述电机(5)的上端活动安装有螺旋桨(6),所述机架(4)的上端固定安装有Jetson TX2板载计算机(2),所述Jetson TX2板载计算机(2)的上端固定安装有pixhawk飞行控制器(1),所述机架(4)的下端固定安装有单轴云台(7),所述单轴云台(7)的下端活动安装有摄像机(3);
步骤S1:一共使用四个坐标系,分别是世界坐标系I={Oi,Xi,Yi,Zi},机体坐标系B={Ob,Xb,Yb,Zb},相机坐标系C={Oc,Xc,Yc,Zc},虚拟坐标系V={Ov,Xv,Yv,Zv},假设相机坐标系和机体坐标系的原点重合,虚拟坐标系下的虚拟图像平面是始终水平的,在机体坐标系下的无人机的运动学方程如下:
Figure FDA0003680379950000011
Figure FDA0003680379950000012
Figure FDA0003680379950000013
Figure FDA0003680379950000014
F=-U1E3+mgRTe3
其中ζ代表的是世界坐标系下的位移,R表示的是从机体坐标系到世界坐标系的转换矩阵,上述的第二个公式表示的是转换矩阵的求导公式,第三个公式表示机体坐标系下无人机的加速度
Figure FDA0003680379950000015
m表示的是无人机的质量,Ω表示无人机的角速度,F是机体坐标系下无人机所受的合力,J是无人机的转动惯量,第四个公式是无人机旋转方向上的运动公式,τ是作用在无人机上的转矩,第五个公式表示的是无人机在机体坐标系下所受到的合力;
步骤S2:在虚拟坐标系计算完偏角并且云台调整完之后在相机坐标系下基于backstepping控制法的无人机控制率如下:
Figure FDA0003680379950000016
一共经过三层的推导,第一层是以速度为控制量,第二层是以力为控制量,第三层是以力的导为控制量,上式是经过backstepping推导计算得到的最终结果,其中
Figure FDA0003680379950000021
表示的是推力的导数,sk(·)表示反对称矩阵,e3表示三维列向量[0,0,1]T
Figure FDA0003680379950000022
并且q1,q2,q3是backstepping方法中每一步的控制量误差,
Figure FDA0003680379950000023
是构建速度观测器时q1的测量值和估计值之差,z*是高度的期望值,c1,k1,k2,k3,k4,k5是控制参数,偏航控制的控制率为
Figure FDA0003680379950000024
其中Ω1和Ω2是无人机的俯仰、横滚角速度;
步骤S3:所使用的虚拟平面是水平的一个平面,世界坐标系,相机坐标系,虚拟坐标系,机体坐标系的关系如下:
I=RPRψV=RψC=RψRθRφB
其中I是世界坐标系,B是机体坐标系,C是相机坐标系,V是虚拟坐标系,相机坐标系和世界坐标系之间平行只相差一个偏航的误差,RP是虚拟坐标系与相机坐标系之间的变换矩阵只有一个俯仰变化;
步骤S4:判别世界坐系下目标平面的偏角其方法如下:
追踪的目标至少可以是四个特征点,假设这四个特征点分别是位于正方形的四个角上,当目标处于水平的平面上,虚拟坐标系下四个特征点也应该是两两之间距离相等的,目标点在虚拟平面下的横纵坐标分别是u1,v1,u2,v2,u3,v3,u4,v4,首先求虚拟平面下四个边的长度:
Figure FDA0003680379950000025
Figure FDA0003680379950000026
然后比较相对的一对边与另一对边的长短,即比较A+C和B+D,当A+C大于B+D则表示目标平面是在俯仰方向产生了倾斜,反之则是图像平面在横滚方向产生了倾斜,至于目标倾斜角的正负,由于通过摄像机透视法则,离摄像头近的在图像平面中会比里的远的大,则是需要直接比较相对的两个边的大小,当A>C时表示目标的偏转角为正,此时控制单轴云台向机头方向偏转,相机平面与目标平面平行,角度的计算是
Figure FDA0003680379950000031
步骤S5:对于无人机控制信息获取如下:在相机坐标系下的特征点位置信息是:
Figure FDA0003680379950000032
对于位置信息的微分如下:
Figure FDA0003680379950000033
其中
Figure FDA0003680379950000034
是相机坐标系相对于惯性坐标系的线速度,v=[cvx cvy cvz]是在相机坐标系下的线速度,根据相机透视法得到像素坐标系下点的坐标:
Figure FDA0003680379950000035
可以根据上式得到像素坐标的导如下:
Figure FDA0003680379950000036
然后计算特征点坐标的平均值,并且使用图像矩来控制无人机的追踪动作,特征点平均坐标是:
Figure FDA0003680379950000037
a=cμ20+cμ02
设置a*是a的理想值,控制位移和偏航的图像矩为:
Figure FDA0003680379950000038
Figure FDA0003680379950000041
步骤S6:在backstepping控制法下进行仿真实验,验证追踪方法的有效性。
2.根据权利要求1所述的一种基于单轴云台的四旋翼无人机目标追踪系统,其特征在于:所述无人机机体坐标系下的运动方程其特征在于为保证准确追踪无人机需尽量保证飞行的稳定性,在水平的虚拟坐标系下所受的合力为f=-RφθU1E3+mge3,其中U1是所有螺旋桨产生的无人机的推力,Rφθ是机体和摄像头一起产生的偏转造成的,E3和e3都是三维列向量[0 0 1]T,在相机坐标系下的合力为
Figure FDA0003680379950000042
其中R是机体的偏转矩阵,Rp是测得的目标平面的倾角。
3.根据权利要求2所述的一种基于单轴云台的四旋翼无人机目标追踪系统,其特征在于:在S2解算出来的推力约束是如下
Figure FDA0003680379950000043
可以对合力
Figure FDA0003680379950000044
求导来获得的
Figure FDA0003680379950000045
可以得到如下:
Figure FDA0003680379950000046
根据求得的姿态角速度Ω1、Ω2然后根据公式:
Figure FDA0003680379950000047
得到惯性坐标系下角度的导
Figure FDA0003680379950000048
然后使用PD控制求得角度
Figure FDA0003680379950000049
对于位置控制使用如下的无人机动力学模型:
Figure FDA0003680379950000051
Figure FDA0003680379950000052
Figure FDA0003680379950000053
首先通过PD控制求得线速度,然后再通过PD控制求得位置。
4.根据权利要求3所述的一种基于单轴云台的四旋翼无人机目标追踪系统,其特征在于:所述对于偏航的控制还是在水平虚拟平面中进行偏航控制,定义偏航的图像矩为如下:
Figure FDA0003680379950000054
对其求导
Figure FDA0003680379950000055
设偏航角的理想值是
Figure FDA0003680379950000056
定义特征误差
Figure FDA0003680379950000057
对其求导有
Figure FDA0003680379950000058
选择偏航角的控制输入
Figure FDA0003680379950000059
使用如下公式:
Figure FDA00036803799500000510
获得最终的控制
Figure FDA00036803799500000511
5.根据权利要求1所述的一种基于单轴云台的四旋翼无人机目标追踪系统,其特征在于:在初始状态下单轴云台是垂直向下的即角度为0,每一次的调整量是通过当前时刻计算出的角度和前一时刻的角度相减得到控制量。
6.根据权利要求1所述的一种基于单轴云台的四旋翼无人机目标追踪系统,其特征在于:所述单轴云台使用的是步进电机,由pixhawk飞行控制器的pwm输出口进行控制。
CN202110411015.7A 2021-04-16 2021-04-16 一种基于单轴云台的四旋翼无人机目标追踪系统 Expired - Fee Related CN113110542B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110411015.7A CN113110542B (zh) 2021-04-16 2021-04-16 一种基于单轴云台的四旋翼无人机目标追踪系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110411015.7A CN113110542B (zh) 2021-04-16 2021-04-16 一种基于单轴云台的四旋翼无人机目标追踪系统

Publications (2)

Publication Number Publication Date
CN113110542A CN113110542A (zh) 2021-07-13
CN113110542B true CN113110542B (zh) 2022-07-15

Family

ID=76717753

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110411015.7A Expired - Fee Related CN113110542B (zh) 2021-04-16 2021-04-16 一种基于单轴云台的四旋翼无人机目标追踪系统

Country Status (1)

Country Link
CN (1) CN113110542B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114489103B (zh) * 2022-01-21 2023-09-08 西南交通大学 一种耦合单轴云台的手机无人机及控制方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101063510A (zh) * 2007-04-02 2007-10-31 北京科技大学 一种电控云台
CN102902884A (zh) * 2012-09-24 2013-01-30 天津市亚安科技股份有限公司 云台摄像机自动定位角度计算方法
WO2021022580A1 (zh) * 2019-08-02 2021-02-11 睿魔智能科技(深圳)有限公司 一种自动跟踪拍摄方法及系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103149939B (zh) * 2013-02-26 2015-10-21 北京航空航天大学 一种基于视觉的无人机动态目标跟踪与定位方法
JP6090786B2 (ja) * 2013-05-09 2017-03-08 国立大学法人 東京大学 背景差分抽出装置及び背景差分抽出方法
CN111596693B (zh) * 2020-06-17 2023-05-26 中国人民解放军国防科技大学 基于云台相机的无人机对地面目标跟踪控制方法及系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101063510A (zh) * 2007-04-02 2007-10-31 北京科技大学 一种电控云台
CN102902884A (zh) * 2012-09-24 2013-01-30 天津市亚安科技股份有限公司 云台摄像机自动定位角度计算方法
WO2021022580A1 (zh) * 2019-08-02 2021-02-11 睿魔智能科技(深圳)有限公司 一种自动跟踪拍摄方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于云台相机的四旋翼无人机跟踪控制系统;车玉涵 等;《吉林大学学报(信息科学版)》;20190531;第37卷(第03期);第278-285页 *

Also Published As

Publication number Publication date
CN113110542A (zh) 2021-07-13

Similar Documents

Publication Publication Date Title
CN108549399B (zh) 飞行器偏航角修正方法、装置及飞行器
García Carrillo et al. Stabilization and trajectory tracking of a quad-rotor using vision
CN108873929B (zh) 一种固定翼飞机自主着舰方法及系统
WO2021098238A1 (zh) 一种车载摄像头云台伺服系统及控制方法
Hançer et al. Robust hovering control of a quad tilt-wing UAV
Walid et al. Modeling and control of a quadrotor UAV
Prabowo et al. Hardware in-the-loop simulation for visual servoing of fixed wing UAV
CN113110542B (zh) 一种基于单轴云台的四旋翼无人机目标追踪系统
CN111661322A (zh) T型双旋翼飞行器及控制方法
CN115686038A (zh) 一种四旋翼吊挂系统的摆角观测与减摆控制系统及方法
Karásek et al. Free flight force estimation of a 23.5 g flapping wing MAV using an on-board IMU
CN114035598A (zh) 一种多旋翼吊挂系统的视觉摆角检测与减摆方法
CN109308074A (zh) 一种无人机重心偏移的补偿方法及系统
CN116149193B (zh) 一种基于视觉的旋翼无人机抗扰动控制方法及系统
Etemadi Mathematical dynamics, kinematics modeling and PID equation controller of quadcopter
Altug et al. Vision-based pose estimation and control of a model helicopter
CN108279562A (zh) 一种基于滑模pid控制的飞行机械臂
Alcan et al. Robust hovering control of a quadrotor using acceleration feedback
Mian et al. Backstepping based PID control strategy for an underactuated aerial robot
CN108475066B (zh) 无人飞行器姿态计算方法、飞行控制器及无人飞行器
Joukhadar et al. Integral LQR-based 6DOF autonomous quadrocopter balancing system control
Abdessameud et al. Dynamic image-based tracking control for VTOL UAVs
CN114510065A (zh) 一种多旋翼无人机地面目标跟踪控制方法
Harkare et al. Design and Control of a Quadcopter
Abrougui et al. Roll control of a Tail-Sitter VTOL UAV

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220715