异构模型统一化方法

异构模型统一化方法(heterogeneous multi-agent)


搜索关键词:异构多智能体编队控制 UAV与USV、UGV heterogeneous multi-agent fotmation control

引言

​ 编队控制中一大分支是异构系统的编队控制。由于单个类型的多智能体能够执行的任务有限,又各自有各自的缺点。因此,异构系统的编队开始受到研究学者的关注。如无人机(UAV)和无人车(UGV)的陆空组合编队,又或者如无人艇(USV)和无人机的海空编队。

​ 但在进行编队控制之前,面临着一个相当重要需要解决的问题。异构说明系统的模型不相同,不相同的话如何进行控制器的设计那?这显然是不切实际的。

​ 在经过文献检索和阅读后,我发现尽管对系统的控制方法、具体系统不一致,但进行控制器设计之前,都将系统转换为统一的一阶积分器模型或者二阶积分器模型。甚至有的大牛论文上来直接就展示转换好的模型,到具体实验时,才用UAV、USV等具体的模型仿真实验[1]。同时我在经典的著作《FORMATION CONTROL OF Multi-Agent Systems:A Graph Rigidity Approach》中也看到作者专门针对一阶积分器和二阶积分器模型的编队控制进行了两章的说明[2]。不过我还没看,如有相关的知识,在我看后回来补充。

​ 下面举个例子说明系统是如何通过转换为积分器模型从而实现异构模型的统一化的。

在文献[3]中针对UGV、UAV的异构编队控制进行控制器的设计。

UGV模型

p ˙ x i ( t ) = V i ( t ) cos ⁡ θ i ( t ) , p ˙ y i ( t ) = V i ( t ) sin ⁡ θ i ( t ) , θ ˙ i ( t ) = ω i ( t ) , V ˙ i ( t ) = F i ( t ) M r i , ω ˙ i ( t ) = τ i ( t ) J r i , (1) \begin{array}{l} \dot{p}_{x i}(t)=V_{i}(t) \cos \theta_{i}(t), \\ \dot{p}_{y i}(t)=V_{i}(t) \sin \theta_{i}(t), \\ \dot{\theta}_{i}(t)=\omega_{i}(t), \\ \dot{V}_{i}(t)=\frac{F_{i}(t)}{M_{r i}}, \\ \dot{\omega}_{i}(t)=\frac{\tau_{i}(t)}{J_{r i}}, \end{array}\tag{1} p˙xi(t)=Vi(t)cosθi(t),p˙yi(t)=Vi(t)sinθi(t),θ˙i(t)=ωi(t),V˙i(t)=MriFi(t),ω˙i(t)=Jriτi(t),(1)

​ 其中 p i ( t ) = [ p x i ( t ) , p y i ( t ) ] T p_i(t)=[p_{xi}(t) ,p_{yi}(t)]^T pi(t)=[pxi(t),pyi(t)]T是第 i i i个UGV的位置, V i ( t ) V_i(t) Vi(t)是线速度, w i ( t ) w_i(t) wi(t)是角速度, θ i ( t ) \theta_i(t) θi(t)是第 i i i个UGV的方位角, M r i M_{ri} Mri是UGV质量, J r i J_{ri} Jri是转动惯量, τ i \tau_i τi是UGV的输入扭矩, F i ( t ) F_i(t) Fi(t)是输入力。

​ 由于非完整约束,系统状态矩阵的行列式为零。为了解决这个问题,人们认为机器人的前端必须在队形中。机器人的前点可定义为手点
p x i h ( t ) = p x i ( t ) + L r c o s θ i ( t ) , p y i h ( t ) = p y i ( t ) + L r s i n θ i ( t ) ​ (2) p_{xi}^h(t)=p_{xi}(t)+L_rcos\theta_i(t),\\ p_{yi}^h(t)=p_{yi}(t)+L_rsin\theta_i(t)\\​\tag{2} pxih(t)=pxi(t)+Lrcosθi(t),pyih(t)=pyi(t)+Lrsinθi(t)(2)

UAV模型

p ¨ x i ( t ) = ( cos ⁡ ϕ i ( t ) sin ⁡ θ i ( t ) cos ⁡ ψ i ( t ) + sin ⁡ ϕ i ( t ) sin ⁡ ψ i ( t ) ) U 1 i M q i , p ¨ y i ( t ) = ( cos ⁡ ϕ i ( t ) sin ⁡ θ i ( t ) sin ⁡ ψ i ( t ) − sin ⁡ ϕ i ( t ) cos ⁡ ψ i ( t ) ) U 1 i M q i , p ¨ z i ( t ) = − g + ( cos ⁡ ϕ i ( t ) cos ⁡ θ i ( t ) ) U 1 i M q , ϕ ¨ i ( t ) = θ ˙ i ( t ) ψ ˙ i ( t ) ( I y − I z I x ) − J q I x θ ˙ i ( t ) Ω i + I q I x U 2 i , θ ¨ i ( t ) = ϕ ˙ i ( t ) ψ ˙ i ( t ) ( I z − I x I y ) − J q I y ϕ ˙ i ( t ) Ω i + l q I y U 3 i , ψ ¨ i ( t ) = ϕ ˙ i ( t ) θ ˙ i ( t ) ( I x − I y I z ) + 1 I z U 4 i , (3) \begin{aligned} \ddot{p}_{x i}(t)= & \left(\cos \phi_{i}(t) \sin \theta_{i}(t) \cos \psi_{i}(t)\right. \\ & \left.+\sin \phi_{i}(t) \sin \psi_{i}(t)\right) \frac{U_{1}^{i}}{M_{q i}}, \\ \ddot{p}_{y i}(t)= & \left(\cos \phi_{i}(t) \sin \theta_{i}(t) \sin \psi_{i}(t)\right. \\ & \left.-\sin \phi_{i}(t) \cos \psi_{i}(t)\right) \frac{U_{1}^{i}}{M_{q i}}, \\ \ddot{p}_{z i}(t)= & -g+\left(\cos \phi_{i}(t) \cos \theta_{i}(t)\right) \frac{U_{1}^{i}}{M_{q}}, \\ \ddot{\phi}_{i}(t)= & \dot{\theta}_{i}(t) \dot{\psi}_{i}(t)\left(\frac{I_{y}-I_{z}}{I_{x}}\right)-\frac{J_{q}}{I_{x}} \dot{\theta}_{i}(t) \Omega^{i}+\frac{I_{q}}{I_{x}} U_{2}^{i}, \\ \ddot{\theta}_{i}(t)= & \dot{\phi}_{i}(t) \dot{\psi}_{i}(t)\left(\frac{I_{z}-I_{x}}{I_{y}}\right)-\frac{J_{q}}{I_{y}} \dot{\phi}_{i}(t) \Omega^{i}+\frac{l_{q}}{I_{y}} U_{3}^{i}, \\ \ddot{\psi}_{i}(t)= & \dot{\phi}_{i}(t) \dot{\theta}_{i}(t)\left(\frac{I_{x}-I_{y}}{I_{z}}\right)+\frac{1}{I_{z}} U_{4}^{i}, \end{aligned}\tag{3} p¨xi(t)=p¨yi(t)=p¨zi(t)=ϕ¨i(t)=θ¨i(t)=ψ¨i(t)=(cosϕi(t)sinθi(t)cosψi(t)+sinϕi(t)sinψi(t))MqiU1i,(cosϕi(t)sinθi(t)sinψi(t)sinϕi(t)cosψi(t))MqiU1i,g+(cosϕi(t)cosθi(t))MqU1i,θ˙i(t)ψ˙i(t)(IxIyIz)IxJqθ˙i(t)Ωi+IxIqU2i,ϕ˙i(t)ψ˙i(t)(IyIzIx)IyJqϕ˙i(t)Ωi+IylqU3i,ϕ˙i(t)θ˙i(t)(IzIxIy)+Iz1U4i,(3)

其中, p i ( t ) = [ p x i ( t ) , p y i ( t ) , p z i ( t ) ] T p_i(t)=[p_{xi}(t),p_{yi}(t),p_{zi}(t)]^T pi(t)=[pxi(t),pyi(t),pzi(t)]T是第 i i i个UAV的位置, ϕ i ( t ) \phi_i(t) ϕi(t)是航向角, θ i ( t ) \theta_i(t) θi(t)是俯仰角, ψ i ( t ) \psi_i(t) ψi(t)是翻滚角。 I x , I y , I z I_x,I_y,I_z IxIy,Iz是机体惯性, M q M_q Mq是UAV质量, J q J_q Jq是机身惯性, L q L_q Lq是UAV长度, Ω i \Omega^i Ωi是第 i i i个UAV的电机输入, U U U是控制输入。

上面两种模型显然是不同的,没有办法直接设计控制器(针对具体问题,也不是说完全没有办法),因此需要进行模型的统一化。

模型统一化

针对UGV模型中的公式(2)进行求导,得
[ p ¨ x i h ( t ) p ¨ y i h ( t ) ] = R ( θ i ( t ) ) T ( t ) (4) \left[\begin{array}{l} \ddot{p}_{x i}^{h}(t) \\ \ddot{p}_{y i}^{h}(t) \end{array}\right]=R\left(\theta_{i}(t)\right) T(t)\tag{4} [p¨xih(t)p¨yih(t)]=R(θi(t))T(t)(4)
其中
R ( θ i ( t ) ) = [ cos ⁡ θ i ( t ) − sin ⁡ θ i ( t ) sin ⁡ θ i ( t ) cos ⁡ θ i ( t ) ] , T ( t ) = [ F i ( t ) M r i − L r ω i 2 ( t ) τ i ( t ) L r J r i + V i ( t ) ω i ( t ) ] (5) \begin{array}{l} R\left(\theta_{i}(t)\right)=\left[\begin{array}{cc} \cos \theta_{i}(t) & -\sin \theta_{i}(t) \\ \sin \theta_{i}(t) & \cos \theta_{i}(t) \end{array}\right], \\ T(t)=\left[\begin{array}{c} \frac{F_{i}(t)}{M_{r i}}-L_{r} \omega_{i}^{2}(t) \\ \frac{\tau_{i}(t) L_{r}}{J_{r i}}+V_{i}(t) \omega_{i}(t) \end{array}\right] \tag{5} \end{array} R(θi(t))=[cosθi(t)sinθi(t)sinθi(t)cosθi(t)],T(t)=[MriFi(t)Lrωi2(t)Jriτi(t)Lr+Vi(t)ωi(t)](5)
u i ( t ) = R ( θ i ( t ) ) T ( t ) u_i(t)=R(\theta_i(t))T(t) ui(t)=R(θi(t))T(t),则:
[ p ¨ x i h ( t ) p ¨ y i h ( t ) ] = u i ( t ) (6) \left[\begin{array}{l} \ddot{p}_{x i}^{h}(t) \\ \ddot{p}_{y i}^{h}(t) \end{array}\right]=u_i(t)\tag{6} [p¨xih(t)p¨yih(t)]=ui(t)(6)
针对UAV模型,直接可令:
[ p ¨ x i ( t ) p ¨ y i ( t ) p ¨ z i ( t ) ] = [ u x i ( t ) u y i ( t ) u z i ( t ) ] (7) \left[\begin{array}{l} \ddot{p}_{x i}(t) \\ \ddot{p}_{y i}(t)\\ \ddot{p}_{z i}(t) \end{array}\right]=\left[\begin{array}{l} u_{x i}(t) \\ u_{y i}(t)\\ u_{z i}(t)\\ \end{array}\right]\tag{7} p¨xi(t)p¨yi(t)p¨zi(t) = uxi(t)uyi(t)uzi(t) (7)
通过公式(6)和公式(7),可以发现UGV和UAV的模型都转化为了二阶积分器的形式,可以进一步地写成如下形式,从而进行进一步的编队控制器设计。
η ˙ i ( t ) = ϵ i ( t ) ϵ ˙ i ( t ) = U i ( t ) (8) \dot{\eta}_{i}(t)=\epsilon_i(t)\\ \dot{\epsilon}_{i}(t)=U_i(t) \tag{8} η˙i(t)=ϵi(t)ϵ˙i(t)=Ui(t)(8)
其中 η i \eta_i ηi可以认为是位置, ϵ i \epsilon_i ϵi认为是速度, U i U_i Ui是控制输入。

除此之外,在文献[4]中,系统也是这样处理的,不过进一步地说明了说明了为什么要改为加速度为控制输入的统一控制模型,其大致过程跟文献[3]类似。这里也贴出供参考。
在这里插入图片描述
在这里插入图片描述

参考文献:

[1] Hua Y , Dong X , Hu G , et al. Distributed Time-Varying Output Formation Tracking for Heterogeneous Linear Multiagent Systems With a Nonautonomous Leader of Unknown Input[J]. IEEE Transactions on Automatic Control, 2019, 64(10):4292-4299.

[2] M.d. Queiroz,et al.Formation Control of Multi-Agent Systems: A Graph Rigidity Approach[M].Willey.2019.

[3] Rahimi R , Abdollahi F , Naqshi K . Time-varying formation control of a collaborative heterogeneous multi agent system[J]. Robotics and Autonomous Systems, 2014(62-12).

[4] 闫茂德, 周小强, 张逊逊,等. 关于异构无人器地空联合编队控制仿真研究[J]. 计算机仿真, 2016(10):5.

如果觉得我写的不错,请给我一个免费的赞,如有错误,也欢迎向我反馈。

  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值