首页 / 专利库 / 视听技术与设备 / 均值漂移过程 / 一种基于蚁群Unscented粒子滤波算法的组合定姿方法

一种基于蚁群Unscented粒子滤波算法的组合定姿方法

阅读:287发布:2020-11-18

专利汇可以提供一种基于蚁群Unscented粒子滤波算法的组合定姿方法专利检索,专利查询,专利分析的服务。并且一种基于蚁群Unscented粒子滤波 算法 的组合定姿方法,本 发明 涉及一种惯性/天文组合定姿方法。该方法首先利用惯性量测信息进行补偿陀螺输出数据,通过 姿态 解算,得到载体姿态信息;然后利用天文量测信息,通过确定性算法,获得所需天文姿态信息;最后利用蚁群Unscented粒子滤波(Unscented Particle Filter)算法将天文姿态信息和载体姿态信息相融合,解决系统非线性和噪声非高斯问题,求解高 精度 载体姿态信息,估计陀螺漂移,并反馈校正载体姿态和补偿陀螺漂移;最终实现基于天文量测信息实时消除惯性/天文组合 导航系统 陀螺随机误差的在线修正,完成对 航天器 的长时间、高精度组合定姿。,下面是一种基于蚁群Unscented粒子滤波算法的组合定姿方法专利的具体信息内容。

1.一种基于蚁群Unscented粒子滤波算法的组合定姿方法,其实现步骤如下:
(1)利用惯性量测信息补偿陀螺输出数据,通过姿态解算,得到载体姿态信息;
(2)利用天文量测信息,通过确定性算法,获得所需天文姿态信息;
(3)利用蚁群Unscented粒子滤波(Unscented Particle Filter)算法将天文姿态信息和载体姿态信息相融合,求解高精度的载体姿态信息,估计陀螺漂移,并反馈校正载体姿态和补偿陀螺漂移补偿;最终实现基于天文量测信息实时消除惯性/天文组合导航系统陀螺随机误差的在线修正,完成对航天器的长时间、高精度组合定姿;
所述步骤(3)中利用蚁群UPF算法步骤为:
(3.1)采样时间t=0时,初始化:
对初始的先验概率密度p(x0)进行采样,生成N个服从p(x0)分布的粒子 i=1,…,N,生成的粒子 的均值和方差满足:
其中,为 的均值,P0为 的方差,E[·]为求取[]内元素的期望,将p(x0)分布取为均值为 方差为P0正态分布;
(3.2)采样时间t≥1时,步骤如下:
①采样
利用(3.1)中生成的服从p(x0)分布的粒子 进行下一时刻的采样,用
Unscented卡尔曼滤波对粒子 进行估计,得到 采样
得到更新的粒子 i=1,…,
N,
其中, 别为k-1时刻状态对应的第i个粒子, 为k-1时刻的粒子的误差方差阵,和Pk_UKF分别为利用UKF算法根据k-1时刻的粒子估计出的第k时刻状态估计值和估计误差方差阵,x0:k-1为第0~k-1时刻的状态估计值,y1:k为第1~k时刻的状态观测值,q(xk|x0:k-1,y1:k)为重要性概率密度,此处选为 为均值为
方差为Pk_UKF的正态分布;
② 利 用 ① 中UKF 更 新 的 粒 子 计 算 粒 子 的 权 重:
归一化权重:
其中, 为k时刻第i个粒子的权重, 为归一化后的权重, 为所有粒子的权重的和, 为对应于观测模型的系统状态的观测似然概率密度, 对应于系
统的模型的系统状态转移概率密度, 为重要性概率密度,初始时刻粒子的权值 i=1,…,N;
③利用①中得出的粒子和②中得出的粒子的权重使用蚁群算法进行重采样,选取优等粒子(权值较大的粒子),剔出低等(权值较小的粒子)的粒子,以解决粒子枯竭问题,蚁群算法优化后的粒子为
④输出
按照最小方差准则,载体姿态的最优估计值就是条件分布的均值,即:
其中,为k时刻载体姿态的最优估计, 为蚁群算法优化后的k时刻第i个粒子的权值, 为蚁群算法优化后的第k时刻第i个粒子的值,为蚁群算法优化后的第k时刻粒子的估计值, 为从i=1到N求和,pk为蚁群算法优化后的第k时刻载体姿态的方差。
2.根据要求1所述的基于蚁群Unscented粒子滤波算法的组合定姿方法,其特征在于:
所述利用蚁群算法实现步骤如下:
首先引入如下符号:
m——蚁群中蚂蚁的数量;
dij——两城市i和j之间的距离;
ηij(t)——边(i,j)的能见度,反映由城市i转移到城市j的启发程度,这个量在蚂蚁系统的运行中不改变;
τij(t)——t时刻边(i,j)上的信息素轨迹强度;
Δτij——蚂蚁k在边(i,j)上的留下的单位长度轨迹信息素量;
——蚂蚁k的转移概率,j为未访问的城市;
每只蚂蚁都是具有如下特征的简单主体:
I从城市i到城市j的运动过程中或是在完成一次循环后,蚂蚁在边(i,j)上释放的一种物质,称为信息素轨迹;
II蚂蚁概率的选择下一个将要访问的城市,这个概率是两城市间距离和连接两城市的路径上存有轨迹量的函数;
III为了满足问题的约束条件,在完成一次循环之前,不允许蚂蚁选择已经访问过的城市;
具体实现过程为:
a.初始化
令时间t=0,迭代次数Nc=0,信息素τij(0)=C,C为正常数,根据具体应用进行设置,此处随意设置为C=1,τij(0)为t=0时边(i,j)上的信息素轨迹强度,(i,j)为某时刻蚂蚁所处的位置
b.对N个粒子的权值进行一次排序选择权值最大的点作为起点,将m只蚂蚁置于起点,各只蚂蚁,按照下列转移概率公式,采用赌轮选择方式移动,
其中,allowedk表示蚂蚁k下一步允许走过的路径点的集合;α为启发式因子,β为期望启发式因子,分别反映了蚂蚁在运动过程中所积累的信息和启发信息在蚂蚁选择路径中的相对重要性,可设置α=1,β=2,城市i转移到城市j的能见度ηij(t)=1/dij(t),令dij(t)为第i个粒子和第j个粒子的权重的差值;
c.按照各只蚂蚁的目标函数值Fk,并记录该次循环的最优解;选择下一个粒子(下一个目标城市)的权重 作为目标函数值Fk;
d.按照以下公式修正信息素强度:
τij(t+n)=ρτij(t)+(1-ρ)Δτij,
式中,参数ρ(0≤ρ≤1)为信息素残留因子,1-ρ表示信息素衰减度; 表示第k只蚂蚁在本次循环中留在节点(i,j)上的信息量;常数Q是信息素强度,取Q=100;
e.令t=t+n,Nc=Nc+1,经过n个时刻,完成一次循环时间t加n,循环次数Nc加1;
f.若N<NCmax,则转步骤b,否则转步骤f,其中NCmax为循环次数;
g.输出最优解。
3.根据要求1所述的基于蚁群Unscented粒子滤波算法的组合定姿方法,其特征在于:
所述确定性算法实现如下:
A.定义3×3的矩阵w,v,B和S,3×1的列向量z,a,标量σ,4×1的列向量q;
其中,w=[w1 w2 w3]为k时刻观测的三颗星的星光在星敏感器坐标系中的坐标矢量,v=[v1 v2 v3]为k时刻该三颗星的星光在地心惯性坐标系中的参考矢量+’T
S=B+BT, a=[a1 a2 a3] 为非负的加权系数,σ=tr(B)为矩阵B的秩,T
q=[q1 q2 q3 q4] 为待求解的姿态四元数,标量在前形式,
定义姿态矩阵K阵如下:
I为单位阵,
姿态矩阵K阵的最大特征值所对应的特征矢量是最小均方差意义下的最优估计,即Kq=λmaxq,q为求解所得姿态四元数,λmax为最大特征值;
T
B.由q=[q1 q2 q3 q4],计算姿态余弦阵C′为:
C.由姿态余弦阵C′即可求解载体的实时高精度天文姿态信息,步骤如下:
俯仰、航向和横滚三姿态 的解算公式如下:
-1
俯仰角θ值为:θ=sin (C23);
航向角 值的计算如下表1所示:
横滚角γ值的计算如下表2所示:
C13值判断 C33值判断 横滚角γ值
=0 <0 -π
>0 <0 atan-1(-C13/C33)-π
>0 =0 -π/2
任意值 >0 atan-1(-C13/C33)
<0 =0 π/2
<0 <0 atan-1(-C13/C33)+π

说明书全文

一种基于蚁群Unscented粒子滤波算法的组合定姿方法

技术领域

[0001] 本发明涉及一种惯性/天文组合定姿方法,特别是一种基于蚁群Unscented粒子滤波(无迹粒子滤波)算法的组合定姿方法,可用于各种航天器的高精度组合定姿。

背景技术

[0002] 为满足天基对地观测、武器精确打击以及空间探索开发的迫切需求,各类地球卫星、深空探测器、载人飞船、弹道导弹和运载火箭等航天器必须具备自主运行和管理能,而高精度的自主定姿是航天器自主运行和管理的核心技术瓶颈。目前,航天器的高精度自主定姿,无法依靠任何一种导航手段独立实现。纯惯性导航系统能够自主、实时提供连续、全面的导航信息,短时精度高,但其误差随工作时间积累,难以满足航天器的长时间高精度定姿要求;天文导航能够提供高精度姿态信息,误差不随时间积累,但易受气候条件限制,且输出信息不连续;将这两者相结合、优势互补,构成惯性/天文组合定姿系统,是实现航天器长时间、高精度定姿的最为有效的手段。
[0003] 在惯性/天文组合定姿技术方面,以往都采用扩展卡尔曼滤波EKF(Extended Kalman Filter)方法,但是EKF仅适用于滤波误差和预测误差很小的情况。近年来提出的Unscented卡尔曼滤波UKF(Unscented Kalman Filter)是一种EKF的改进算法,有效的解决了系统的非线性问题,但其不足是不适用于噪声非高斯分布的系统。粒子滤波PF由于采用蒙特卡洛采样(Monte Carlo sampling)结构而在非线性、非高斯系统状态跟踪上体现出越来越大的优越性,但其缺点是存在退化现象,消除退化现象主要依赖于两个关键技术:适当选取重要密度函数和进行重采样。对于前者的改进方法,可使用EKPF(Extented Kalman Particle Filter)、UPF(Unscented Particle Filter)来进行重要密度函数的选择,其中UPF算法是利用UKF来得到粒子重要性概率密度函数的一种粒子滤波方法,由于该重要密度函数中包含了最新量测信息,因此具有更好的性能。对于后者的改进方法,常用的重采样算法有累积分布重采样(Binary search)、系统重采样(Systematic resampling)、剩余重采样(Residual resampling)等,这些算法通过增加粒子的有效性解决了粒子的退化问题,但是在实际应用中会影响系统的鲁棒性,重采样完成后,重要度高的粒子通过重采样被多次选取,这在一定程度上丢失了粒子的多样性,由此造成的后果是一旦目标丢失或跟踪精度不够,系统自动收敛的可能性很小,为此,很多学者提出了遗传粒子滤波(GPF)算法,GPF算法虽然在保证粒子有效性的同时又增加了粒子的多样性,仍然存在滤波速度慢和鲁棒性差的问题。

发明内容

[0004] 本发明的技术解决问题是:克服现有技术的不足,提出一种基于蚁群UPF的组合定姿方法,解决系统非线性和噪声非高斯问题,以快速获得高精度的姿态信息,并能够准确地估计陀螺漂移,实现各种类型航天器长时间、高精度的组合定姿。
[0005] 本发明的技术解决方案为:一种基于蚁群UPF组合定姿方法,其特点在于:利用惯性量测信息和天文量测信息,通过蚁群(Ant Colony Algorithm)UPF(无迹粒子滤波)方法,实现航天器长时间、高精度的快速组合定姿,其实现步骤如下:
[0006] (1)利用惯性量测信息补偿陀螺输出数据,通过姿态解算,得到载体姿态信息;
[0007] (2)利用天文量测信息,通过确定性算法,获得所需的天文姿态信息;
[0008] (3)利用蚁群(Ant Colony Algorithm)Unscented粒子滤波(Unscented Particle Filter)算法将天文姿态信息和载体姿态信息相融合,求解高精度的载体姿态信息,估计陀螺漂移,并反馈校正载体姿态和补偿陀螺漂移补偿;最终实现基于天文量测信息实时消除惯性/天文组合导航系统陀螺随机误差的在线修正,完成对航天器的高精度组合定姿;
[0009] 利用蚁群UPF算法进行信息融合的步骤为:
[0010] (3.1)采样时间t=0时,初始化:
[0011] 对初始的先验概率密度p(x0)进行采样,生成N个服从p(x0)分布的粒子 i=1,…,N,生成的粒子 的均值和方差满足:
[0012]
[0013]
[0014] 其中, 为 的均值, 为 的方差,E[·]为求取[]内元素的期望,将p(x0)分布取为均值为 方差为P0正态分布;
[0015] (3.2)采样时间t≥1时,步骤如下:
[0016] ①采样
[0017] 利用(3.1)中生成的服从p(x0)分布的粒子 进行下一时刻的采样,用Unscented卡尔曼滤波对粒子 进行估计,得到 采样
得到更新的粒子 i=1,…,N,
[0018] 其中, 和 分别为k-1时刻状态对应的第i个粒子和粒子的误差方差阵,和Pk_UKF分别为根据k-1时刻的粒子估计的第k时刻状态估计值和估计误差方差阵,x0:k-1为第0~k-1时刻的状态估计值,y1:k为第1~k时刻的状态观测值,q(xk|x0:k-1,y1:k)为重要性概率密度,此处选为 为均值为 方差为Pk_UKF的正态分布;
[0019] ② 利 用 ① 中UKF 更 新 的 粒 子 计 算 粒 子 的 权 重归一化权重:
[0020] 其中, 为k时刻第i个粒子的权值, 为归一化后的权重, 为所有粒子的权值的和, 为对应于观测模型的系统状态的观测似然概率密度, 为对应于系统的模型的系统状态转移概率密度, 为重要性概率密度;
[0021] ③利用①中得出的粒子和②中得出的粒子的权重使用蚁群算法进行重采样,选取优等粒子(权值较大的粒子),剔出低等(权值较小的粒子)的粒子,以解决粒子枯竭问题,利用蚁群算法进行优化的步骤如下:
[0022] 首先引入如下记号:
[0023] m——蚁群中蚂蚁的数量;
[0024] dij——两城市i和j之间的距离;
[0025] ηij(t)——边(i,j)的能见度,反映由城市i转移到城市j的启发程度,这个量在蚂蚁系统的运行中不改变;
[0026] τij(t)——t时刻边(i,j)上的信息素轨迹强度;
[0027] Δτij——蚂蚁k在边(i,j)上的留下的单位长度轨迹信息素量;
[0028] ——蚂蚁k的转移概率,j为未访问的城市。
[0029] 每只蚂蚁都是具有如下特征的简单主体:
[0030] I从城市i到城市j的运动过程中或是在完成一次循环后,蚂蚁在边(i,j)上释放的一种物质,称为信息素轨迹;
[0031] II蚂蚁概率的选择下一个将要访问的城市,这个概率是两城市间距离和连接两城市的路径上存有轨迹量的函数;
[0032] III为了满足问题的约束条件,在完成一次循环之前,不允许蚂蚁选择已经访问过的城市。
[0033] a.初始化
[0034] 令时间t=0,迭代次数Nc=0,信息素τij(0)=C,C为正常数,根据具体应用进行设置,此处随意设置为C=1,τij(0)为t=0时边(i,j)上的信息素轨迹强度,(i,j)为某时刻蚂蚁所处的位置
[0035] b.对N个粒子的权值进行一次排序选择权值最大的点作为起点,将m只蚂蚁置于起点,各只蚂蚁,按照下列转移概率公式,采用赌轮选择方式移动,
[0036]
[0037] 其中,allowedk表示蚂蚁k下一步允许走过的路径点的集合;α为启发式因子,β为期望启发式因子,分别反映了蚂蚁在运动过程中所积累的信息和启发信息在蚂蚁选择路径中的相对重要性,可设置α=1,β=2,城市i转移到城市j的能见度ηij(t)=1/dij(t),令dij(t)为第i个粒子和第j个粒子的权重的差值;
[0038] c.按照各只蚂蚁的目标函数值Fk,并记录该次循环的最优解;选择下一个粒子(下一个目标城市)的权值 作为目标函数值Fk;
[0039] d.按照以下公式修正信息素强度:
[0040] τij(t+n)=ρτij(t)+(1-ρ)Δτij,
[0041]
[0042] 式中,参数ρ(0≤ρ≤1)为信息素残留因子,1-ρ表示信息素衰减度; 表示第k只蚂蚁在本次循环中留在节点(i,j)上的信息量;常数Q是信息素强度,取Q=100;
[0043] e.令t=t+n,Nc=Nc+1,经过n个时刻,完成一次循环时间t加n,循环次数Nc加1;
[0044] f.若N<NCmax,则转步骤b,否则转步骤f,其中NCmax为循环次数;
[0045] g.输出最优解。
[0046] 蚁群算法优化后的粒子为
[0047] ④输出
[0048] 按照最小方差准则,载体姿态的最优估计值就是条件分布的均值,即:
[0049]
[0050]
[0051] 其中,为k时刻载体姿态的最优估计, 为蚁群算法优化后的k时刻第i个粒子的权值, 为蚁群算法优化后的第k时刻第i个粒子的值,为蚁群算法优化后的第k时刻粒子的估计值, 为从i=1到N求和,pk为蚁群算法优化后的第k时刻载体姿态的方差。
[0052] 本发明的原理是:首先利用陀螺输出数据对惯性量测信息进行补偿,通过姿态解算,得到载体姿态信息;其次利用天文量测信息,通过确定性算法获得特定间隔的天文姿态信息;最后利用蚁群UPF算法将天文姿态信息和载体姿态信息相融合,解决系统非线性和噪声非高斯问题,求解高精度载体姿态信息,估计陀螺漂移,并反馈校正载体姿态和补偿陀螺漂移补偿;最终实现航天器长时间、高精度的组合定姿。
[0053] 本发明与现有技术相比的优点在于:本发明克服了传统组合定姿方法在定姿精度和陀螺漂移估计精度低的不足,利用UPF有效解决了系统非线性和噪声非高斯的问题,利用蚁群算法在路径寻优方面的优势对Unscented粒子滤波的粒子进行优化,有效的解决了粒子滤波的粒子退化和粒子匮乏问题,实现了优等粒子选择的快速性和有效性,提高了组合定姿的速度和精度;将惯性量测信息和天文量测信息相融合,进一步提高了组合定姿的精度,实现了对陀螺漂移的精确估计,满足了航天器长时间、高精度组合定姿的要求。附图说明
[0054] 图1为本发明的一种基于蚁群UPF的组合定姿方法原理图;
[0055] 图2为本发明中用蚁群算法优化粒子的流程图

具体实施方式

[0056] 如图1所示,本发明的具体实施步骤如下:
[0057] 1、首先对惯性量测信息进行补偿陀螺输出数据后,通过姿态解算,得到载体姿态信息,流程如下:
[0058] a.设定初始姿态为 计算得出初始姿态四元数阵q(0):
[0059]
[0060] 其中, θ0,γ0分别为俯仰横滚角和偏航角,q(0)为0时刻的姿态四元数,cos[·],sin[·]分别为求余弦和正弦;
[0061] b.由a中给出的初始姿态四元数阵q(0)推导出更新矩阵为:
[0062]
[0063] n为第n时刻,I为单位四元数,Δφ=[ΔφX ΔφY ΔφZ]为安装在在X,Y,Z三个轴上的陀螺输出角增量,定义[ΔΦ]为:
[0064]
[0065] c·由b中得出的姿态四元数更新矩阵q(n+1)=[q1,n+1 q2,n+1 q3,n+1 q4,n+1]T,计算姿态余弦阵C为:
[0066]
[0067] 其中,q(n+1)为k+1时刻的姿态四元数,C11~C33对应公式最右边矩阵中的元素,q(n+1)=C·q(n);
[0068] d.由方向余弦阵C求解载体的实时姿态信息:
[0069] 俯仰、航向和横滚三姿态角 的解算公式如下:
[0070] 俯仰角θ值为:θ=sin-1(C23);
[0071] 航向角 值的计算如下表所示:
[0072]
[0073]
[0074] 横滚角γ值的计算如下表所示:
[0075]C13值判断 C33值判断 横滚角γ值
=0 <0 -π
>0 <0 atan-1(-C13/C33)-π
>0 =0 -π/2
任意值 >0 atan-1(-C13/C33)
<0 =0 π/2
<0 <0 atan-1(-C13/C33)+π
[0076] 2、利用天文量测信息,通过确定性算法,求解天文姿态信息的步骤为:
[0077] A.定义3×3的矩阵w,v,B和S,3×1的列向量z,a和标量σ,4×1的列向量q;
[0078] 其中,w=[w1 w2 w3]为k时刻观测的三颗星的星光在星敏感器坐标系中的坐标矢量,v=[v1 v2 v3]为k时刻该三颗星的星光在地心惯性坐标系中的参考矢量,T T
S=B+B, a=[a1 a2 a3] 为非负的加权系数,σ=tr(B)为矩阵B的秩,
T
q=[q1 q2 q3 q4] 为待求解的姿态四元数,
[0079]
[0080]
[0081]
[0082]
[0083]
[0084] 定义姿态矩阵K阵如下:
[0085] I为单位阵,
[0086] 姿态矩阵K阵的最大特征值所对应的特征矢量是最小均方差意义下的最优估计,即Kq=λmaxq,q为求解所得姿态四元数,λmax为最大特征值;T
[0087] B.由q=[q1 q2 q3 q4],计算姿态余弦阵C′为:
[0088]
[0089] C.由姿态余弦阵C′即可求解载体的实时高精度天文姿态信息,步骤如下:
[0090] 俯仰、航向和横滚三姿态角 的解算公式如下:
[0091] 俯仰角θ值为:θ=sin-1(C23);
[0092] 航向角 值的计算如下表所示:
[0093]
[0094] 横滚角γ值的计算如下表所示:
[0095]C13值判断 C33值判断 横滚角γ值
=0 <0 -π
>0 <0 atan-1(-C13/C33)-π
>0 =0 -π/2
任意值 >0 atan-1(-C13/C33)
<0 =0 π/2
<0 <0 atan-1(-C13/C33)+π
[0096] 3、利用蚁群算法优化的UPF算法将天文姿态信息和载体姿态信息相融合,完成对航天器长时间、高精度的组合定姿步骤为:
[0097] ①采样
[0098] 利用(3.1)中生成的服从p(x0)分布的粒子 进行下一时刻的采样,用Unscented卡尔曼滤波对粒子 进行估计,得到 采样
得到更新的粒子 i=1,…,N,
[0099] 其中, 别为k-1时刻状态对应的第i个粒子, 为k-1时刻的粒子的误差方差阵, 和Pk_UKF分别为根据k-1时刻的粒子估计的第k时刻状态估计值和估计误差方差阵,x0:k-1为第0~k-1时刻的状态估计值,y1:k为第1~k时刻的状态观测值,q(xk|x0:k-1,y1:k)为重要性概率密度,此处选为 为均值为 方差为Pk_UKF的正态分布;
[0100] ② 利 用 ① 中 UKF更 新 的 粒 子 计 算 粒 子 的 权 重归一化权重:
[0101] 其中, 为k时刻第i个粒子的权值, 为归一化后的权重, 为所有粒子的权值的和, 为对应于观测模型的系统状态的观测似然概率密度, 为对应于系统的模型的系统状态转移概率密度, 为重要性概率密度,;
[0102] ③利用①中得出的粒子和②中得出的粒子的权重使用蚁群算法进行重采样,选取优等粒子(权值较大的粒子),剔出低等(权值较小的粒子)的粒子,以解决粒子枯竭问题,利用蚁群算法进行优化的步骤如下:
[0103] 首先引入如下记号:
[0104] m——蚁群中蚂蚁的数量;
[0105] dij——两城市i和j之间的距离;
[0106] ηij(t)——边(i,j)的能见度,反映由城市i转移到城市j的启发程度,这个量在蚂蚁系统的运行中不改变;
[0107] τij(t)——t时刻边(i,j)上的信息素轨迹强度;
[0108] Δτij——蚂蚁k在边(i,j)上的留下的单位长度轨迹信息素量;
[0109] ——蚂蚁k的转移概率,j为未访问的城市。
[0110] 每只蚂蚁都是具有如下特征的简单主体:
[0111] I从城市i到城市j的运动过程中或是在完成一次循环后,蚂蚁在边(i,j)上释放的一种物质,称为信息素轨迹;
[0112] II蚂蚁概率的选择下一个将要访问的城市,这个概率是两城市间距离和连接两城市的路径上存有轨迹量的函数;
[0113] III为了满足问题的约束条件,在完成一次循环之前,不允许蚂蚁选择已经访问过的城市。
[0114] a.初始化
[0115] 令时间t=0,迭代次数Nc=0,信息素τij(0)=C,C为正常数,根据具体应用进行设置,此处随意设置为C=1,τij(0)为t=0时边(i,j)上的信息素轨迹强度,(i,j)为某时刻蚂蚁所处的位置;
[0116] b.对N个粒子的权值进行一次排序选择权值最大的点作为起点,将m只蚂蚁置于起点,各只蚂蚁,按照下列转移概率公式,采用赌轮选择方式移动,
[0117]
[0118] 其中,allowedk表示蚂蚁k下一步允许走过的路径点的集合;α为启发式因子,β为期望启发式因子,分别反映了蚂蚁在运动过程中所积累的信息和启发信息在蚂蚁选择路径中的相对重要性,可设置α=1,β=2,城市i转移到城市j的能见度ηij(t)=1/dij(t),令dij(t)为第i个粒子和第j个粒子的权重的差值;
[0119] c.按照各只蚂蚁的目标函数值Fk,并记录该次循环的最优解;选择下一个粒子(下一个目标城市)的权值 作为目标函数值Fk;
[0120] d.按照以下公式修正信息素强度:
[0121] τij(t+n)=ρτij(t)+(1-ρ)Δτij,
[0122]
[0123] 式中,参数ρ(0≤ρ≤1)为信息素残留因子,1-ρ表示信息素衰减度; 表示第k只蚂蚁在本次循环中留在节点(i,j)上的信息量;常数Q是信息素强度,取Q=100;
[0124] e.令t=t+n,Nc=Nc+1,经过n个时刻,完成一次循环时间t加n,循环次数Nc加1;
[0125] f.若N<NCmax,则转步骤b,否则转步骤f,其中NCmax为循环次数;
[0126] g.输出最优解。
[0127] 蚁群算法优化后的粒子为
[0128] ④输出
[0129] 按照最小方差准则,载体姿态的最优估计值就是条件分布的均值,即:
[0130]
[0131]
[0132] (3)利用蚁群Unscented粒子滤波(Unscented Particle Filter)算法将天文姿态信息和载体姿态信息相融合,解决系统非线性和噪声非高斯问题,求解高精度的载体姿态信息,估计陀螺漂移,并反馈校正载体姿态和补偿陀螺漂移补偿;最终实现基于天文量测信息实时消除惯性/天文组合导航系统陀螺随机误差的在线修正,完成对航天器的长时间、高精度组合定姿;
[0133] 本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈