CN112577489A - 一种基于交互多模型滤波的导引头视线转率提取方法 - Google Patents
一种基于交互多模型滤波的导引头视线转率提取方法 Download PDFInfo
- Publication number
- CN112577489A CN112577489A CN202011422730.2A CN202011422730A CN112577489A CN 112577489 A CN112577489 A CN 112577489A CN 202011422730 A CN202011422730 A CN 202011422730A CN 112577489 A CN112577489 A CN 112577489A
- Authority
- CN
- China
- Prior art keywords
- sight
- equation
- missile
- line
- target
- 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.)
- Granted
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 24
- 238000000605 extraction Methods 0.000 title claims abstract description 23
- 230000002452 interceptive effect Effects 0.000 title claims abstract description 22
- 238000005259 measurement Methods 0.000 claims abstract description 25
- 238000006243 chemical reaction Methods 0.000 claims abstract description 23
- 238000000034 method Methods 0.000 claims abstract description 23
- 230000001133 acceleration Effects 0.000 claims description 42
- 239000011159 matrix material Substances 0.000 claims description 15
- 230000008859 change Effects 0.000 claims description 6
- 230000007704 transition Effects 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 8
- 230000003287 optical effect Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000012937 correction Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
Abstract
本发明的一个实施例公开一种基于交互多模型滤波的导引头视线转率提取方法,包括:S10、根据导引头视线空间运动方程得到视线转率方程;S20、根据视线转率方程得到系统状态方程;S30、根据系统状态方程得到系统量测方程;S40、基于IMM视线转率滤波得到导引头视线转率。本发明采用一种基于交互多模型的导引头视线转率提取方法,可对末制导视线转率进行高精度滤波估计,基于IMM方法可以实现高精度视线转率提取,从而进行高精度制导控制。
Description
技术领域
本发明涉及导引头视线转率提取领域,更具体地,涉及一种基于交互多模型滤波的导引头视线转率提取方法、系统、计算机设备和存储介质。
背景技术
当前目标隐身机动突防性能越来越强,如何在拦截过程中提高拦截概率对导引头视线转率提取精度提出了更高的要求。弹上捷联导引头具有测量精度高、高可靠性、低成本的特点,但是由于导引头测量的视线转率包含了弹体姿态信息,需要解耦算法去除弹体姿态带来的干扰信息,得到真实的导引头信息,且导引头在探测目标过程中引入较大的量测噪声,直接影响制导控制系统稳定性。
目前视线转率提取方法主要有三种思路:第一种采用基于扩维模型的容积卡尔曼滤波算法,该算法虽然考虑扩维模型提高了视线转率提取精度,但是针对高机动目标或者组合机动目标则该算法具有一定的局限性;第二种方法先提取出惯性视线角,然后对惯性视线角进行微分或滤波估计获取惯性视线角速度信息,这类处理方法一方面涉及多次坐标变换,该方法计算复杂,需要多次转换坐标系,实现复杂,并且导航坐标系选择受到限制;第三种方法采用导引头测量的视线角信息构造出惯性视线角速度信息,这类处理方法由于采用了简化的模型计算惯性视线角速度,只适用于低动态的制导武器。
发明内容
本发明的一个目的在于提供一种基于交互多模型滤波的导引头视线转率提取方法。
本发明的一个目的在于提供一种基于交互多模型滤波的导引头视线转率提取系统。
本发明的一个目的在于提供一种计算机设备。
本发明的另一个目的在于提供一种存储介质。
为达到上述目的,本发明采用下述技术方案:
本发明第一方面提供一种基于交互多模型滤波的导引头视线转率提取方法,包括:
S10、根据导引头视线空间运动方程得到视线转率方程;
S20、根据视线转率方程得到系统状态方程;
S30、根据系统状态方程得到系统量测方程;
S40、基于IMM视线转率滤波得到导引头视线转率。
在一个具体实施例中,在步骤S10中,目标位于于距离原点O的xp轴方向,O点位于导引头旋转中心,目标在视线系x、y和z三个方向下的加速度分别为atx、aty和atz,导弹在视线系x、y和z三个方向下的加速度分别为asx、asy和asz,目标加速度为,导弹加速度为
目标和导弹之间的相对关系为:
由公式(2)和(3)可以得到目标加速度与导弹加速度的相对关系为:
在一个具体实施例中,在步骤S20中,由惯性系到视线系转换矩阵L(qγ,qλ)为:
其中,qγ为视线高低角;qλ为视线方位角。
由视线系角速度与视线角之间的转换关系为:
视线运动方程为:
由公式(4)、公式(7)-公式(9)可以得出三方向加速度分量为:
其中,atsx代表目标相对导弹的加速度在视线坐标系Oxs上的分量;atsy代表目标相对导弹的加速度在视线坐标系Oys上的分量;atsz代表目标相对导弹的加速度在视线坐标系Ozs上的分量;为R的二阶导。
设定三个姿态角为俯仰角θ,偏航角ψ,滚转角γ
建立状态方程如下:
其中,jtsx为目标相对导弹的加速度变化率在视线坐标系Oxs上的分量。
式(12)中状态方程为:
将三个姿态角变化率向机体轴上投影,得欧拉角和机体系下的角速度关系为:
根据公式(14)中的状态变量得到状态方程:
从而得到总状态方程为:
在一个具体实施例中,在步骤S30中,根据导弹脱靶量和弹目相对距离关系,确定系统量测信息包括导引头测量视线高低角qγ,视线方位角qλ及导弹脱靶量ρ,系统量测量为z=[qγ,qγ,ρ]
量测方程如下:
其中,Rij代表地面系到弹体系的转换矩阵R中第i行第j列元素;v1代表视线高低角的量测噪声;v2代表视线方位角的量测噪声;v3代表脱靶量的量测噪声。
在一个具体实施例中,所述S40还包括:
S400、状态估计的交互式作用
设从模型i转移到模型j的转移概率矩阵为:
式(19)中,可以得到
S402、模型修正
S404、模型似然性计算
模型j滤波残差为vj(k),相应的协方差为Sj(k),并服从高斯分布,模型j的似然性Λj(k)为
其中,
S406、模型概率更新
模型j的概率更新为:
S408、模型输出
本发明第二方面提供一种基于交互多模型滤波的导引头视线转率提取系统,包括:
第一模块:用于根据导引头视线空间运动方程得到视线转率方程;
第二模块:用于根据视线转率方程得到系统状态方程;
第三模块:用于根据系统状态方程得到系统量测方程;
第四模块:用于基于IMM视线转率滤波得到导引头视线转率。
本发明第三方面提供一种计算机设备,包括:
一个或多个处理器;
存储装置,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现如本发明第一方面所述的方法。
本发明第四方面提供一种存储介质,其上存储有计算机程序,该程序被处理器执行时实现如本发明第一方面所述的方法。
本发明的有益效果如下:
本发明采用一种基于交互多模型的导引头视线转率提取方法,可对末制导视线转率进行高精度滤波估计,基于IMM方法可以实现高精度视线转率提取,从而进行高精度制导控制。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明。
图1示出根据本申请的实施例的基于交互多模型滤波的视线转率提取方法的流程图。
图2示出根据本申请的实施例的导弹与目标相对位置估计误差示意图。
图3示出根据本申请的实施例的导弹与目标相对速度估计误差示意图。
图4示出根据本申请的实施例的视线高低角误差示意图。
图5示出根据本申请的实施例的视线高低角速率误差示意图。
图6示出根据本申请的实施例的视线方位角误差示意图。
图7示出根据本申请的实施例的视线方位角速率误差示意图。
图8示出根据本申请的实施例的实施本发明的视线转率提取方法的计算机设备的结构示意图。
具体实施方式
为了更清楚地说明本发明,下面结合优选实施例和附图对本发明做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解,下面所具体描述的内容是说明性的而非限制性的,不应以此限制本发明的保护范围。
实施例一
图1是本发明实施例提供的基于交互多模型滤波的视线转率提取方法的流程图。如图1所示,该方法包括如下步骤:
S10、根据导引头视线空间运动方程得到视线转率方程。
在一个具体实施例中,目标位于于距离原点O的xp轴方向,O点位于导引头旋转中心,目标在视线系x、y和z三个方向下的加速度分别为atx、aty和atz,导弹在视线系x、y和z三个方向下的加速度分别为asx、asy和asz,目标加速度为,导弹加速度为
目标和导弹之间的相对关系为:
由公式(2)和(3)可以得到目标加速度与导弹加速度的相对关系为:
S20、根据视线转率方程得到系统状态方程。
在一个具体实施例中,由惯性系到视线系转换矩阵L(qγ,qλ)为:
其中,qγ为视线高低角;qλ为视线方位角。
由视线系角速度与视线角之间的转换关系为:
视线运动方程为:
由公式(4)、公式(7)-公式(9)可以得出三方向加速度分量为:
其中,atsx代表目标相对导弹的加速度在视线坐标系Oxs上的分量;atsy代表目标相对导弹的加速度在视线坐标系Oys上的分量;atsz代表目标相对导弹的加速度在视线坐标系Ozs上的分量;为R的二阶导。
设定三个姿态角为俯仰角θ,偏航角ψ,滚转角γ
建立状态方程如下:
其中,jtsx为目标相对导弹的加速度变化率在视线坐标系Oxs上的分量。
式(12)中状态方程为:
将三个姿态角变化率向机体轴上投影,得欧拉角和机体系下的角速度关系为:
根据公式(14)中的状态变量得到状态方程:
从而得到总状态方程为:
S30、根据系统状态方程得到系统量测方程。
在一个具体实施例中,根据导弹脱靶量和弹目相对距离关系,确定系统量测信息包括导引头测量视线高低角qγ,视线方位角qλ及导弹脱靶量ρ,系统量测量为z=[qγ,qγ,ρ]
量测方程如下:
其中,Rij代表地面系到弹体系的转换矩阵R中第i行第j列元素;v1代表视线高低角的量测噪声;v2代表视线方位角的量测噪声;v3代表脱靶量的量测噪声。
S40、基于IMM视线转率滤波得到导引头视线转率。
在一个具体实施例中,所述S40还包括:
S400、状态估计的交互式作用
设从模型i转移到模型j的转移概率矩阵为:
式(19)中,可以得到
S402、模型修正
S404、模型似然性计算
模型j滤波残差为vj(k),相应的协方差为Sj(k),并服从高斯分布,模型j的似然性Λj(k)为
其中,
S406、模型概率更新
模型j的概率更新为:
S408、模型输出
曲线见附图说明图2-图7,分别示出了导弹与目标相对位置估计误差图、导弹与目标相对速度估计误差图、视线高低角误差图、视线高低角速率误差图、视线方位角误差图和视线方位角速率误差图。从图2-7中可以看出,采用基于交互多模型的导引头视线转率提取方法,可对末制导视线转率进行高精度滤波估计。
实施例二
本实施例提供了一种基于交互多模型滤波的导引头视线转率提取系统,包括:
第一模块:用于根据导引头视线空间运动方程得到视线转率方程;
第二模块:用于根据视线转率方程得到系统状态方程;
第三模块:用于根据系统状态方程得到系统量测方程;
第四模块:用于基于IMM视线转率滤波得到导引头视线转率。
实施例三
如图8所示,本发明的一个实施例提供了一种计算机设备的结构示意图,图8显示的计算机设备12仅仅是一个示例,不应对本发明实施例的功能和使用范围带来任何限制。
如图8所示,计算机设备12以通用计算设备的形式表现。计算机设备12的组件可以包括但不限于:一个或者多个处理器或者处理单元16,系统存储器28,连接不同系统组件(包括系统存储器28和处理单元16)的总线18。
总线18表示几类总线结构中的一种或多种,包括存储器总线或者存储器控制器,外围总线,图形加速端口,处理器或者使用多种总线结构中的任意总线结构的局域总线。举例来说,这些体系结构包括但不限于工业标准体系结构(ISA)总线,微通道体系结构(MAC)总线,增强型ISA总线、视频电子标准协会(VESA)局域总线以及外围组件互连(PCI)总线。
计算机设备12典型地包括多种计算机系统可读介质。这些介质可以是任何能够被计算机设备12访问的可用介质,包括易失性和非易失性介质,可移动的和不可移动的介质。
系统存储器28可以包括易失性存储器形式的计算机系统可读介质,例如随机存取存储器(RAM)30和/或高速缓存存储器32。计算机设备12可以进一步包括其它可移动/不可移动的、易失性/非易失性计算机系统存储介质。仅作为举例,存储系统34可以用于读写不可移动的、非易失性磁介质(图8未显示,通常称为“硬盘驱动器”)。尽管图8中未示出,可以提供用于对可移动非易失性磁盘(例如“软盘”)读写的磁盘驱动器,以及对可移动非易失性光盘(例如CD-ROM,DVD-ROM或者其它光介质)读写的光盘驱动器。在这些情况下,每个驱动器可以通过一个或者多个数据介质接口与总线18相连。存储器28可以包括至少一个程序产品,该程序产品具有一组(例如至少一个)程序模块,这些程序模块被配置以执行本发明各实施例的功能。
具有一组(至少一个)程序模块42的程序/实用工具40,可以存储在例如存储器28中,这样的程序模块42包括但不限于操作系统、一个或者多个应用程序、其它程序模块以及程序数据,这些示例中的每一个或某种组合中可能包括网络环境的实现。程序模块42通常执行本发明所描述的实施例中的功能和/或方法。
计算机设备12也可以与一个或多个外部设备14(例如键盘、指向设备、显示器24等)通信,还可与一个或者多个使得用户能与该计算机设备12交互的设备通信,和/或与使得该计算机设备12能与一个或多个其它计算设备进行通信的任何设备(例如网卡,调制解调器等等)通信。这种通信可以通过输入/输出(I/O)接口22进行。并且,计算机设备12还可以通过网络适配器20与一个或者多个网络(例如局域网(LAN),广域网(WAN)和/或公共网络,例如因特网)通信。如图8所示,网络适配器20通过总线18与计算机设备12的其它模块通信。应当明白,尽管图8中未示出,可以结合计算机设备12使用其它硬件和/或软件模块,包括但不限于:微代码、设备驱动器、冗余处理单元、外部磁盘驱动阵列、RAID系统、磁带驱动器以及数据备份存储系统等。
处理器单元16通过运行存储在系统存储器28中的程序,从而执行各种功能应用以及数据处理,例如实现实施例一所提供的方法。
实施例四
本发明的另一个实施例提供了一种存储介质,其上存储有计算机程序,该程序被处理器执行时实现如上述实施例一所提供的方法。
在实际应用中,所述计算机可读存储介质可以采用一个或多个计算机可读的介质的任意组合。计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质。计算机可读存储介质例如可以是但不限于电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括:具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本实施例中,计算机可读存储介质可以是任何包含或存储程序的有形介质,该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。
计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号,其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式,包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质,该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。
计算机可读介质上包含的程序代码可以用任何适当的介质传输,包括但不限于无线、电线、光缆、RF等等,或者上述的任意合适的组合。
可以以一种或多种程序设计语言或其组合来编写用于执行本发明操作的计算机程序代码,所述程序设计语言包括面向对象的程序设计语言—诸如Java、Smalltalk、C++,还包括常规的过程式程序设计语言—诸如“C”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中,远程计算机可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)—连接到用户计算机,或者,可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。
Claims (8)
1.一种基于交互多模型滤波的导引头视线转率提取方法,其特征在于,包括:
S10、根据导引头视线空间运动方程得到视线转率方程;
S20、根据视线转率方程得到系统状态方程;
S30、根据系统状态方程得到系统量测方程;
S40、基于IMM视线转率滤波得到导引头视线转率。
2.根据权利要求1所述的方法,其特征在于,在步骤S10中,目标位于于距离原点O的xp轴方向,O点位于导引头旋转中心,目标在视线系x、y和z三个方向下的加速度分别为atx、aty和atz,导弹在视线系x、y和z三个方向下的加速度分别为asx、asy和asz,目标加速度为,导弹加速度为
目标和导弹之间的相对关系为:
由公式(2)和(3)可以得到目标加速度与导弹加速度的相对关系为:
3.根据权利要求1所述的方法,其特征在于,在步骤S20中,由惯性系到视线系转换矩阵L(qγ,qλ)为:
其中,qγ为视线高低角;qλ为视线方位角;
由视线系角速度与视线角之间的转换关系为:
视线运动方程为:
由公式(4)、公式(7)-公式(9)可以得出三方向加速度分量为:
其中,atsx代表目标相对导弹的加速度在视线坐标系Oxs上的分量;atsy代表目标相对导弹的加速度在视线坐标系Oys上的分量;atsz代表目标相对导弹的加速度在视线坐标系Ozs上的分量;为R的二阶导;
设定三个姿态角为俯仰角θ,偏航角ψ,滚转角γ
建立状态方程如下:
其中,jtsx为目标相对导弹的加速度变化率在视线坐标系Oxs上的分量;
式(12)中状态方程为:
将三个姿态角变化率向机体轴上投影,得欧拉角和机体系下的角速度关系为:
根据公式(14)中的状态变量得到状态方程:
5.根据权利要求1所述的方法,其特征在于,所述S40还包括:
S400、状态估计的交互式作用;
设从模型i转移到模型j的转移概率矩阵为:
式(19)中,可以得到
S402、模型修正;
S404、模型似然性计算;
模型j滤波残差为vj(k),相应的协方差为Sj(k),并服从高斯分布,模型j的似然性Λj(k)为
其中,
S406、模型概率更新
模型j的概率更新为:
S408、模型输出
6.一种基于交互多模型滤波的导引头视线转率提取系统,其特征在于,包括:
第一模块:用于根据导引头视线空间运动方程得到视线转率方程;
第二模块:用于根据视线转率方程得到系统状态方程;
第三模块:用于根据系统状态方程得到系统量测方程;
第四模块:用于基于IMM视线转率滤波得到导引头视线转率。
7.一种计算机设备,其特征在于,包括:
一个或多个处理器;
存储装置,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现如权利要求1-5中任一项所述的方法。
8.一种存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现如权利要求1-5中任一项所述方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011422730.2A CN112577489B (zh) | 2020-12-08 | 2020-12-08 | 一种基于交互多模型滤波的导引头视线转率提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011422730.2A CN112577489B (zh) | 2020-12-08 | 2020-12-08 | 一种基于交互多模型滤波的导引头视线转率提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112577489A true CN112577489A (zh) | 2021-03-30 |
CN112577489B CN112577489B (zh) | 2024-05-07 |
Family
ID=75127698
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011422730.2A Active CN112577489B (zh) | 2020-12-08 | 2020-12-08 | 一种基于交互多模型滤波的导引头视线转率提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112577489B (zh) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4312262A (en) * | 1979-02-22 | 1982-01-26 | General Electric Company | Relative velocity gunsight system and method |
GB8502369D0 (en) * | 1985-01-30 | 1996-07-24 | Secr Defence | Missile guidance system |
ITRM970588A1 (it) * | 1995-11-06 | 1999-03-30 | Secr Defence Brit | Procedimento di guida d'arma mediante stima dello stato del bersaglio |
EP1094292A1 (de) * | 1999-10-21 | 2001-04-25 | Bodenseewerk Gerätetechnik GmbH | Verfahren zum Führen eines Flugkörpers auf ein Ziel bei Zielverlust |
US20080246653A1 (en) * | 2007-04-05 | 2008-10-09 | Sirf Technology, Inc. | Gps navigation using interacting multiple model (imm) estimator and probabilistic data association filter (pdaf) |
US20130092785A1 (en) * | 2008-07-11 | 2013-04-18 | Davidson Technologies, Inc. | System and method for guiding and controlling a missile using high order sliding mode control |
CN104111078A (zh) * | 2014-04-29 | 2014-10-22 | 北京理工大学 | 消除全捷联导引头制导回路刻度尺系数误差的装置和方法 |
CN107255924A (zh) * | 2017-06-14 | 2017-10-17 | 哈尔滨工业大学 | 基于扩维模型的容积卡尔曼滤波提取捷联导引头制导信息的方法 |
CN108009358A (zh) * | 2017-12-01 | 2018-05-08 | 哈尔滨工业大学 | 基于imm_ukf的三维制导律辨识滤波方法 |
CN110095116A (zh) * | 2019-04-29 | 2019-08-06 | 桂林电子科技大学 | 一种基于lift的视觉定位和惯性导航组合的定位方法 |
CN111121770A (zh) * | 2019-12-31 | 2020-05-08 | 西北工业大学 | 一种交互式多弹多模型航迹融合方法 |
CN111797478A (zh) * | 2020-07-27 | 2020-10-20 | 北京电子工程总体研究所 | 一种基于变结构多模型的强机动目标跟踪方法 |
-
2020
- 2020-12-08 CN CN202011422730.2A patent/CN112577489B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4312262A (en) * | 1979-02-22 | 1982-01-26 | General Electric Company | Relative velocity gunsight system and method |
GB8502369D0 (en) * | 1985-01-30 | 1996-07-24 | Secr Defence | Missile guidance system |
ITRM970588A1 (it) * | 1995-11-06 | 1999-03-30 | Secr Defence Brit | Procedimento di guida d'arma mediante stima dello stato del bersaglio |
EP1094292A1 (de) * | 1999-10-21 | 2001-04-25 | Bodenseewerk Gerätetechnik GmbH | Verfahren zum Führen eines Flugkörpers auf ein Ziel bei Zielverlust |
US20080246653A1 (en) * | 2007-04-05 | 2008-10-09 | Sirf Technology, Inc. | Gps navigation using interacting multiple model (imm) estimator and probabilistic data association filter (pdaf) |
US20130092785A1 (en) * | 2008-07-11 | 2013-04-18 | Davidson Technologies, Inc. | System and method for guiding and controlling a missile using high order sliding mode control |
CN104111078A (zh) * | 2014-04-29 | 2014-10-22 | 北京理工大学 | 消除全捷联导引头制导回路刻度尺系数误差的装置和方法 |
CN107255924A (zh) * | 2017-06-14 | 2017-10-17 | 哈尔滨工业大学 | 基于扩维模型的容积卡尔曼滤波提取捷联导引头制导信息的方法 |
CN108009358A (zh) * | 2017-12-01 | 2018-05-08 | 哈尔滨工业大学 | 基于imm_ukf的三维制导律辨识滤波方法 |
CN110095116A (zh) * | 2019-04-29 | 2019-08-06 | 桂林电子科技大学 | 一种基于lift的视觉定位和惯性导航组合的定位方法 |
CN111121770A (zh) * | 2019-12-31 | 2020-05-08 | 西北工业大学 | 一种交互式多弹多模型航迹融合方法 |
CN111797478A (zh) * | 2020-07-27 | 2020-10-20 | 北京电子工程总体研究所 | 一种基于变结构多模型的强机动目标跟踪方法 |
Non-Patent Citations (4)
Title |
---|
张家斌: "捷联导引头视线角速率提取研究", 中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑, no. 6, 15 June 2018 (2018-06-15), pages 44 - 57 * |
梁彦刚, 郝道亮, 唐国金: "捷联导引头视线转率估计的交互式多模型样条滤波算法", 国防科技大学学报, vol. 36, no. 5, 31 October 2014 (2014-10-31), pages 70 - 74 * |
秦雷, 李君龙: "临近空间目标非弹道式机动模式跟踪滤波技术", 系统仿真学报, vol. 29, no. 6, 31 December 2017 (2017-12-31), pages 1380 - 1385 * |
臧路尧;林德福;杨哲;王武刚;: "红外制导弹药协同定位跟踪方法研究", 红外与激光工程, no. 04, 25 April 2018 (2018-04-25), pages 124 - 130 * |
Also Published As
Publication number | Publication date |
---|---|
CN112577489B (zh) | 2024-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109297510B (zh) | 相对位姿标定方法、装置、设备及介质 | |
CN110780285B (zh) | 激光雷达与组合惯导的位姿标定方法、系统及介质 | |
JP7299261B2 (ja) | 車両のデッドレコニング方法、装置、デバイス、記憶媒体、及びプログラム | |
CN106767767A (zh) | 一种微纳多模星敏感器系统及其数据融合方法 | |
CN116067370B (zh) | 一种imu姿态解算方法及设备、存储介质 | |
CN109781117B (zh) | 一种组合定位方法和系统 | |
CN103512584A (zh) | 导航姿态信息输出方法、装置及捷联航姿参考系统 | |
US20230366680A1 (en) | Initialization method, device, medium and electronic equipment of integrated navigation system | |
CN114018274A (zh) | 车辆定位方法、装置及电子设备 | |
CN113759384B (zh) | 一种传感器位姿转换关系确定方法、装置、设备和介质 | |
CN109631939B (zh) | 一种基于磁强计和加速度计的快速对准方法 | |
CN111596329A (zh) | 车辆定位方法、装置、设备及车辆 | |
CN114179825A (zh) | 多传感器融合获取量测值置信度方法及自动驾驶车辆 | |
CN111238474A (zh) | 基于倾斜坐标系的捷联导引头非奇异视线角速度提取方法 | |
CN111121755A (zh) | 一种多传感器的融合定位方法、装置、设备及存储介质 | |
CN110096134B (zh) | 一种vr手柄射线抖动矫正方法、装置、终端和介质 | |
CN117213480B (zh) | 一种传递对准方法、系统、设备及存储介质 | |
CN116608853B (zh) | 载体动态姿态估计方法、设备、存储介质 | |
CN112577489B (zh) | 一种基于交互多模型滤波的导引头视线转率提取方法 | |
CN110285811A (zh) | 卫星定位和惯性导航的融合定位方法和装置 | |
CN113218380B (zh) | 一种电子罗盘的校正方法、装置、电子设备及存储介质 | |
CN112270066B (zh) | 一种计算卫星刚性耦合系数的优化方法和计算机设备 | |
CN112649823A (zh) | 无人机导航定位方法及装置 | |
CN115326106B (zh) | 一种惯性传感器数据仿真的改进方法 | |
CN114279396B (zh) | 运动姿态确定方法及装置、计算机可读介质和电子设备 |
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 |