首页 > 专利 > 杭州电子科技大学 > 基于数据驱动学习的滚动优化控制方法专利详情

基于数据驱动学习的滚动优化控制方法   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2020-07-01
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2020-11-24
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2022-05-03
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2040-07-01
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN202010627592.5 申请日 2020-07-01
公开/公告号 CN111897215B 公开/公告日 2022-05-03
授权日 2022-05-03 预估到期日 2040-07-01
申请年 2020年 公开/公告年 2022年
缴费截止日
分类号 G05B13/04 主分类号 G05B13/04
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 0
权利要求数量 1 非专利引证数量 1
引用专利数量 0 被引证专利数量 0
非专利引证 1、CN 105182998 A,2015.12.23CN 107563044 A,2018.01.09王寅 等.在线滚动优化下的无人自转旋翼机自动起飞控制《.控制理论与应用》.2015,第32卷(第11期),第1526-1533页.;
引用专利 被引证专利
专利权维持 2 专利申请国编码 CN
专利事件 事务标签 公开、实质审查、授权
申请人信息
申请人 第一申请人
专利权人 杭州电子科技大学 当前专利权人 杭州电子科技大学
发明人 吕强、胡晨、仲朝亮、石厅、林伟杰 第一发明人 吕强
地址 浙江省杭州市下沙高教园区2号大街 邮编 310018
申请人数量 1 发明人数量 5
申请人所在省 浙江省 申请人所在市 浙江省杭州市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
杭州君度专利代理事务所 代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
杨舟涛
摘要
本发明公开了基于数据驱动学习的滚动优化控制方法。本发明首先将连续的搜寻环境离散化,根据收集到的大量信号强度数据,采用高斯过程回归预测信号源的可能位置;然后,基于信号源可能的位置,设计成本函数,采用滚动优化的方法,产生四旋翼无人机最优的控制序列;最后,根据最优控制序列,进而获得四旋翼无人机的最优运动轨迹,并将第一个位置输入到四旋翼无人机的控制器中,实现对信号源的搜索和定位。本发明弥补了现有技术的不足,通过基于数据驱动学习的滚动优化控制方法,不但使四旋翼无人机利用检测到的信号强度信息预测出信号源的可能位置,而且产生的最优参考轨迹,能够很好的管理不确定事件的发生,快速定位信号源。
  • 摘要附图
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:其中:
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0008]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:5
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:矩系数;
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0017]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0019]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0022]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0024]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0028]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0030]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0031]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0036]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0039]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0040]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0049]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0052]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:和
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0058]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0060]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0063]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0065]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0069]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0071]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0072]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0077]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0080]
    基于数据驱动学习的滚动优化控制方法
  • 说明书附图:[0081]
    基于数据驱动学习的滚动优化控制方法
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2022-05-03 授权
2 2020-11-24 实质审查的生效 IPC(主分类): G05B 13/04 专利申请号: 202010627592.5 申请日: 2020.07.01
3 2020-11-06 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
1.基于数据驱动学习的滚动优化控制方法,该方法具体步骤如下:
第一步:建立四旋翼无人机动力学模型,具体步骤如下:
a.四旋翼无人机动力学模型如公式(1)所示;
其中:X、Y、Z分别表示在惯性坐标系下四旋翼无人机的位置;φ、θ、ψ分别表示在惯性坐标系下四旋翼无人机的横滚角,俯仰角和偏航角; 分别为X、Y、Z的二阶导数;g是重力加速度;m是四旋翼无人机的机体质量; 表示四旋翼无人机
的总推力; 和 代表四个电机的速度;kF表示转子推力系数;sin(·)和cos(·)分别为正弦函数和余弦函数;
b.四旋翼无人机惯性坐标系下姿态动力学模型如(2)所示:
其中:Ix,Iy,和Iz表示三个轴上的惯性矩; 分别为φ、θ、ψ的一阶导数;
分别为φ、θ、ψ的二阶导数; 和
表示四旋翼的扭矩;l表示从转子到四旋翼无人机质量中心的距离;kM表示转子转矩系数;
第二步:基于收集到的信号强度数据,采用高斯过程回归预测出环境内信号源的位置,具体步骤如下:
a.在四旋翼无人机位置数据集Ω上定义理想的信号强度值f(x),对位置数据x1,x2,…,xn∈Ω,n是样本数据的个数,f(x)满足条件(3);
f(x)~GP(m(x),K(x,x))   (3)
T
其中:f(x)=[f(x1),f(x2),...,f(xn)]表示由信号强度组成的向量;m(x)=[m(x1),mT
(x2),...,m(xn)] 表示由均值函数定义的位置向量;GP(·,·)表示高斯过程回归模型,K(x,x)表示协方差矩阵,如公式(4)所示;
其中:
其中:xi和xj是四旋翼无人机的位置,i,j=1,…,n;ζ,σf,和σy是三个调节参数; 是克罗内克函数,如果xi=xj,则为1,否则为0;exp(·)是指数函数;||·||2表示2范数;
b.给定一个样本数据集 xi和yi分别表示已经收集到的四旋翼
无人机的位置和信号强度;M是样本数据的数量;R是实数集;样本数据集中的信号强度满足公式(6);
其中:ε表示独立噪声;N表示正态分布;
为超参数集,根据样本数据,采用公式(7)计算;
T
其中: 表示样本条件概率的负对数似然函数;y=[y1,...,yM] ;
函数 表示在 取最小值时,Θ的取值;
*
d.给定预测集 M是预测样本的数量;x*和f(x*)分别表示待预
测的四旋翼无人机的位置和信号强度,后验概率如公式(8)所示:
其中:
其中:cov(·)为方差;I是单位矩阵;
e.使用公式(8),可以获得环境内每一位置信号强度的概率分布,即建立信号强度分布地图,进一步把信号强度均值最大的位置确定为信号源的可能位置;
第三步:根据信号源的可能位置,对四旋翼无人机路径进行滚动优化,步骤如下:
a.成本函数设计如公式(11)所示:
其中:Jk是成本函数;Nc是预测长度;xk+n|k表示第k时刻预测的第k+Nc步的位置变量;
xk+i|k是在第k时刻预测的第k+i步的位置变量;uk+i|k是在第k时刻预测的第k+i步控制变量;
xd是信号源的可能位置;α,β,γ表示非负加权因子;
b.采用滚动优化控制方法产生四旋翼无人机最优控制序列,如公式(12)所示:
其中:min表示取最小值; 表示最优控制序列 s.t.表示约
束条件;Ω和U是位置变量x和控制变量u的可行区域;因此,通过求解优化问题(12)获得了最优控制序列
c.根据最优控制序列,获得最优轨迹 并将第一个位置
输入到控制器中;
第四步:如果信号源没有定位,则直接返回第二步;如果信号源已经定位,则先标记该信号源,然后返回第二步,检测其它信号源;如果给定的最大搜索时间已经达到,则无人机停止运行,并输出环境内所有信号源的位置。
说明书

技术领域

[0001] 本发明属于自动化技术领域,具体涉及一种面向信号源检测与定位的基于数据驱动学习的滚动优化控制方法。

背景技术

[0002] 信号源检测与定位对人类安全有着非常重要的意义,如环境监测、传感器定位、遇险者营救等等。针对上述问题,通常采用四旋翼无人机完成信号源的检测与定位。传统的信号源检测与定位方法,通常根据信号强度和方向信息,采用径向基函数网络、Kalman滤波器或者粒子滤波器等方法完成简单环境中信号源的位置评估,但这些方法难以完成复杂环境多信号源位置的同时评估。此外,由于复杂环境中不确定事件可能发生,会对传统的控制方法产生较大的干扰,从而使四旋翼无人机飞行定位产生较大影响,在这一背景下,本发明弥补了现有技术的不足。

发明内容

[0003] 针对现有技术的不足之处,本发明提出基于数据驱动学习的滚动优化控制方法,可以控制四旋翼无人机同时对复杂环境中多个信号源的位置进行评估,在保证信号源检测和定位精度的前提下,可以有效得管理不确定事件,提高抗干扰能力。
[0004] 为实现上述发明目的,本发明的技术方案如下:
[0005] 第一步:建立四旋翼无人机动力学模型,具体步骤如下:
[0006] a.四旋翼无人机动力学模型如公式(1)所示;
[0007]
[0008] X、Y、Z分别表示在惯性坐标系下四旋翼无人机的位置;φ、θ、ψ分别表示在惯性坐标系下四旋翼无人机的横滚角,俯仰角和偏航角; 分别为X、Y、Z的二阶导数;g是重力加速度;m是四旋翼无人机的机体质量; 表示四旋翼无人机的总推力; 和 代表四个电机的速度;kF表示转子推力系数;sin(·)和cos(·)
分别为正弦函数和余弦函数;
[0009] b.四旋翼无人机惯性坐标系下姿态动力学模型如(2)所示:
[0010]
[0011] 其中:Ix,Iy,和Iz表示三个轴上的惯性矩; 分别为φ、θ、ψ的一阶导数;分别为φ、θ、ψ的二阶导数; 和
表示四旋翼的扭矩;l表示从转子到四旋翼无人机质量中心的距离;kM表示转子转矩系数;
[0012] 第二步:基于收集到的信号强度数据,采用高斯过程回归预测出环境内信号源的位置,具体步骤如下:
[0013] a.在四旋翼无人机位置数据集Ω上定义理想的信号强度值f(x),对位置数据x1,x2,…,xn∈Ω,n是样本数据的个数,f(x)满足条件(3);
[0014] f(x)~GP(m(x),K(x,x))           (3)
[0015] 其中:f(x)=[f(x1),f(x2),...,f(xn)]T表示由信号强度组成的向量;m(x)=[mT(x1),m(x2),...,m(xn)]表示由均值函数定义的位置向量;GP(·,·)表示高斯过程回归模型,K(x,x)表示协方差矩阵,如公式(4)所示;
[0016]
[0017] 其中:
[0018]
[0019] 其中:xi和xj是四旋翼无人机的位置,i,j=1,…,n;ζ,σf,和σy是三个调节参数;是克罗内克函数,如果xi=xj,则为1,否则为0;exp(·)是指数函数;||·||2表示2范数;
[0020] b.给定一个样本数据集 xi和yi分别表示已经收集到的四旋翼无人机的位置和信号强度;M是样本数据的数量;R是实数集;样本数据集中的信号强度满足公式(6);
[0021]
[0022] 其中:ε表示独立噪声;N表示正态分布;
[0023] 为超参数集,根据样本数据,采用公式(7)计算;
[0024]
[0025] 其中:θ(Θ)=‑log p(y|D,Θ)表示样本条件概率的负对数似然函数;y=T[y1,...,yM];函数arg minΘθ(Θ)表示在θ(Θ)取最小值时,Θ的取值;
*
[0026] d.给定预测集 M是预测样本的数量;x*和f(x*)分别表示待预测的四旋翼无人机的位置和信号强度,后验概率如公式(8)所示:
[0027]
[0028] 其中:
[0029]
[0030]
[0031] 其中:cov(·)为方差;I是单位矩阵。
[0032] e.使用公式(8),可以获得环境内每一位置信号强度的概率分布,即建立信号强度分布地图,进一步把信号强度均值最大的位置确定为信号源的可能位置;
[0033] 第三步:根据信号源的可能位置,对四旋翼无人机路径进行滚动优化,步骤如下:
[0034] a.成本函数设计如公式(11)所示:
[0035]
[0036] 其中:Jk是成本函数;Nc是预测长度;xk+n|k表示第k时刻预测的第k+Nc步的位置变量;xk+i|k是在第k时刻预测的第k+i步的位置变量;uk+i|k是在第k时刻预测的第k+i步控制变量;xd是信号源的可能位置;α,β,γ表示非负加权因子;
[0037] b.采用滚动优化控制方法产生四旋翼无人机最优控制序列,如公式(12)所示:
[0038]
[0039]
[0040] 其中:min表示取最小值; 表示最优控制序列 s.t.表示约束条件;Ω和U是位置变量x和控制变量u的可行区域;因此,通过求解优化问题(12)获得了最优控制序列
[0041] c.根据最优控制序列,获得最优轨迹 并将第一个位置输入到控制器中;
[0042] 第四步:如果信号源没有定位,则直接返回第二步;如果信号源已经定位,则先标记该信号源,然后返回第二步,检测其它信号源;如果给定的最大搜索时间已经达到,则无人机停止运行,并输出环境内所有信号源的位置。
[0043] 本发明具有以下有益效果:本发明采用高斯过程回归,能够同时评估复杂环境中多个信号源的位置。此外,运用的滚动优化方法可以有效管理不确定事件的发生,提高信号源检测和定位精度,以及抗干扰能力更强。

实施方案

[0044] 以下结合一个实施例对本发明的技术方案作进一步的阐述说明。
[0045] 以无线传感器网络为例,设定需监测的环境范围长为10米,宽为10米,建立坐标系统,将该监测环境分为50×50的小区域。其中,每个传感器节点在短时间内发送信号序列,并假定两个发送间隔5秒。具体实施步骤如下:
[0046] 第一步:建立四旋翼无人机动力学模型,具体步骤如下:
[0047] a.四旋翼无人机动力学模型如公式(13)所示。
[0048]
[0049] 其中:X、Y、Z分别表示在惯性坐标系下四旋翼无人机的位置。φ、θ、ψ分别表示在惯性坐标系下四旋翼无人机的横滚角,俯仰角和偏航角; 分别为X、Y、Z的二阶导数;g为2
9.8m/s ,是重力加速度;m为1.2kg,是四旋翼无人机的机体质量;
‑4 2
表示四旋翼无人机的总推力; 和 代表四个电机的速度;kF为10 Ns ,表示转子推力系数;sin(·)和cos(·)分别为正弦函数和余弦函数。
[0050] b.四旋翼无人机惯性坐标系下姿态动力学模型如公式(14)所示:
[0051]
[0052] 其中:Ix、Iy和Iz表示三个轴上的惯性矩,取值分别为0.05kgm2、0.05kgm2和2
0.1kgm; 分别为φ、θ、ψ的一阶导数; 分别为φ、θ、ψ的二阶导数;
和 表示四旋翼的扭矩;
2
l为0.2m,表示从转子到四旋翼无人机质量中心的距离;kM为0.1Nms,表示转子转矩系数。
[0053] 第二步:基于收集到的信号强度数据,采用高斯过程回归预测出环境内传感器的位置,具体步骤如下:
[0054] a.在四旋翼无人机位置数据集Ω上定义理想的信号强度值f(x),对位置数据x1,x2,…,xn∈Ω,n是样本数据的个数,f(x)满足公式(15)。
[0055] f(x)~GP(m(x),K(x,x))        (15)
[0056] 其中:f(x)=[f(x1),f(x2),...,f(xn)]T表示由信号强度组成的向量;m(x)=[mT(x1),m(x2),...,m(xn)]由均值函数定义的位置向量;GP(·,·)表示高斯过程回归模型,K(x,x)表示协方差矩阵,如(4)所示。
[0057]
[0058] 其中:
[0059]
[0060] 其中:xi和xj是四旋翼无人机的位置(i,j=1,…,n);ζ,σf,和σy是三个调节参数;是克罗内克函数,如果xi=xj,则为1,否则为0;exp(·)是指数函数;||·||2表示2范数。
[0061] b.给定一个样本数据集 xi和yi分别表示已经收集到的四旋翼无人机的位置和信号强度。M是样本数据的数量;R是实数集。样本数据集中的信号强度满足(18)。
[0062]
[0063] 其中:ε表示独立噪声;σy是噪声方差;N表示正态分布;
[0064] 为超参数集,采用公式(19)计算。
[0065] Θ=argminΘθ(Θ)          (19)
[0066] 其中:θ(Θ)=‑logp(y|D,Θ)表示样本条件概率的负对数似然函数;y=[y1,...,TyM];函数arg minΘθ(Θ)表示在θ(Θ)取最小值时,Θ的取值。
[0067] d.给定预测集 M*是预测位置的数量;x*和f(x*)分别表示待预测的四旋翼无人机的位置和信号强度,后验概率如公式(20)所示:
[0068]
[0069] 其中:
[0070]
[0071]
[0072] 其中:cov(·)为方差;I是单位矩阵。
[0073] e.使用公式(20),可以获得环境内每一位置信号强度的概率分布,即建立信号强度分布地图,进一步把信号强度均值最大的位置确定为传感器节点的可能位置。
[0074] 第三步:根据传感器节点的可能位置,对四旋翼无人机路径进行滚动优化,步骤如下:
[0075] a.成本函数设计如公式(23)所示:
[0076]
[0077] 其中:Jk是成本函数;Nc是预测长度; 表示第k时刻预测的第k+Nc步的位置变量;xk+i|k是在第k时刻预测的第k+i步的位置变量;uk+i|k是在第k时刻预测的第k+i步控制变量;xd是传感器节点的可能位置;α为1、β为0.01、γ为1。
[0078] b.采用滚动优化控制方法产生四旋翼无人机最优控制序列,如公式(24)所示:
[0079]
[0080]
[0081] 其中:min表示取最小值; 表示最优控制序列 s.t.表示约束条件;Ω和U是位置变量x和控制变量u的可行区域。因此,通过求解优化问题(24)获得了最优控制序列
[0082] c.根据最优控制序列,获得最优轨迹 并将第一个位置输入到控制器中。
[0083] 第四步:如果信号源没有定位,则直接返回第二步;如果信号源已经定位,则先标记该信号源,然后返回第二步,检测其它信号源;如果给定的最大搜索时间已经达到,则无人机停止运行,并输出环境内所有信号源的位置。
版权所有:盲专网 ©2023 zlpt.xyz  蜀ICP备2023003576号