CN112255612A - 一种雷达航迹随机抖动性评估方法 - Google Patents
一种雷达航迹随机抖动性评估方法 Download PDFInfo
- Publication number
- CN112255612A CN112255612A CN202011209612.3A CN202011209612A CN112255612A CN 112255612 A CN112255612 A CN 112255612A CN 202011209612 A CN202011209612 A CN 202011209612A CN 112255612 A CN112255612 A CN 112255612A
- Authority
- CN
- China
- Prior art keywords
- matrix
- track
- radar
- point
- error
- 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
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000004458 analytical method Methods 0.000 claims abstract description 25
- 238000012360 testing method Methods 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 88
- 230000001427 coherent effect Effects 0.000 claims description 22
- 230000008859 change Effects 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 10
- 238000005259 measurement Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 5
- 238000001514 detection method Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 description 6
- 238000011156 evaluation Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000009897 systematic effect Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供了一种雷达航迹随机抖动性评估方法,该方法将雷达航迹点位置映射到极坐标系中,从角度和距离两个维度上计算测量雷达航迹点相对于参考雷达航迹点的位置偏差情况,采用多尺度相干分析方法,量化和评估多个周期内的测试航迹相对于标准航迹的随机抖动情况。
Description
技术领域
本发明属于空中交通管理领域,尤其涉及一种雷达航迹随机抖动性评估方法。
背景技术
雷达误差一般有***误差和随机误差。***误差是对同一被测量物体进行了多次测量,将所得结果的平均值与被测物体的真实值之差。***误差一般来源于***本身,可以通过精确校准来减少误差,或者通过数学补偿的方法可以基本去除***误差。随机误差又称为偶然误差或者不定误差,是在测量过程中一系列有关因素微小的无序波动而形成的具有相互抵偿性的误差,其主要来源为接收机的热噪声和杂波干扰,雷达随机误差没有有效办法能够彻底消除,因此需要对雷达随机误差进行分析。
民航管制自动化***接入雷达信号,经过格式解析获取目标信息,然后将目标信息显示在管制***显示界面上。管制员通过观察管制界面上的飞机位置、高度、速度等信息进行日常管制活动。为了实现精确管制,保障飞机飞行安全,管制员需要知道飞机的精确位置信息,通过雷达航迹的随机抖动性评估,管制员能够准确区分雷达源的数据好坏,避免管制判断出错。
之前的雷达航迹抖动性评估方式都是计算单个点的误差,通过计算误差的标准差或者方差来表示航迹的随机抖动性。方差或标准差体现的是数据整体离散程度,没有考虑到雷达航迹在时间上连续、相关的特性。
发明内容
发明目的:本发明的目的在于提供一种实用的一种雷达航迹随机抖动性评估方法,通过前期的数据处理,得到完整可用的航迹数据,然后从距离和角度两个方面计算误差,运用相干分析方法对误差数据进行量化分析得到航迹抖动性变化曲线。本发明具体包括:
步骤1、剔除野值。航空器在飞行过程中,其高度、速度、位置数据会随着时间变化,但是其高度、速度、位置在一定雷达扫描周期内的变化是有限的。雷达扫描到的目标信息可能会因为某些不确定因素影响出现明显偏离实际的情况,因此要剔除这部分数据;
步骤2、数据补偿;
步骤3、进行时间对齐;被测雷达A和标准雷达B,对同一目标的位置测量值分别为(xA,yA)、(xB,yB),速度、角度、时间分别为(vA,cA,tA)、(vB,cB,tB)。由于两部雷达的报出时间不一定统一,因此需要对相近时间内的数据点进行时间对齐操作。因为两部雷达的时间差很短,小于等于两部雷达中最小扫描周期,可以认为是一个很短的时间,因此时间对齐操作采用线性插值的方法实现;
步骤4、进行坐标转换;被测雷达A和标准雷达B,对同一目标的位置测量值分别为(xA,yA)、(xB,yB)。雷达A的坐标为(xRA,yRA),需要将(xA,yA)、(xB,yB)转换到以(xRA,yRA)为极点的极坐标系中,得到(ρA,θA),(ρB,θB);
步骤5、计算距离误差;经过坐标转换后,被测雷达A和标准雷达B探测到的目标位置分别为(ρA,θA),(ρB,θB),计算两个目标到极点距离的差值。假设目标航迹一共有K个点,则得到K个距离误差值(Δρ1,Δρ2,Δρ3…ΔρK);
步骤6、计算角度误差;经过坐标转换后,被测雷达A和标准雷达B探测到的目标位置分别为(ρA,θA),(ρB,θB),计算两个目标极角的差值。假设目标航迹一共有K个点,则得到K个角度误差值(Δθ1,Δθ2,Δθ3…ΔθK);
步骤7、多尺度相干分析;根据步骤4和步骤5计算得到两个点的距离误差和角度误差。分别对距离误差和角度误差其进行相干分析计算。例如,利用相干分析法计算距离抖动性:选取尺度N(N为大于1的正整数),如N取5,将5个航迹点的距离误差作为一组得到一个1*5维向量即X1=(Δρ1,Δρ2,Δρ3,Δρ4,Δρ5),选取连续的五组得到一个5阶矩阵即求矩阵的协方差矩阵E,求协方差矩阵E的特征值,根据相干系数计算公式得到相干系数。
步骤8、绘制抖动性曲线。如果尺度N=5,将窗口矩阵右移一个航迹点,得到一个重叠了五分之四的新矩阵,按照步骤7的方法重新求得相干系数。窗口继续右移,就会得到相干系数,将相邻的相干系数连接起来,就得到相干系数曲线。
步骤1中,设定剔除门限,航空器在飞行过程中,其高度、速度、位置数据会随着时间变化,但是其高度、速度、位置在一定雷达扫描周期内的变化是有限的。雷达扫描到的目标信息可能会因为某些不确定因素影响出现明显偏离实际的情况,因此要剔除这部分数据。设定一架民航飞机的最大速度为v(m/s),雷达扫描周期为T(s),则连续两个探测点的距离应当不大于v*T(m),设定一个随机误差固定容错率20%,则相邻两个点的距离应当小于v*T*1.2(m),将不满足上述条件的点剔除。
一般民航飞机在飞行过程中其飞行速度是有限的,在直线平飞过程中,民航飞机的速度一般为900km/h。假设飞机的速度为900km/h且处于直线飞行,若雷达扫描周期为4s,在不计入随机误差的情况下,相邻两个航迹点的距离大约为250m/s*4s=1000m,考虑随机误差,其相邻两个航迹点间的距离正常情况下也不会超过2000m。因此,在进行计算时,需要将部分明显异常的数据剔除掉。设定剔除门限,剔除门限可以在理论飞行距离的基础上加上一个固定值来得到,根据实际飞行过程中的大量航迹点数据,可以将这个固定值设置为1000m。
步骤2包括:经过步骤1之后,可能出现如下情况,被测雷达现有三个航迹点A1、A2、A3,经过步骤1的剔除野值后,剩下航迹点A1、A3,已知被测雷达扫描周期为T,此时航迹点A1和点A3之间的时间差值变为2T,为了避免不同的时间尺度影响,因此需要采用线性差分方法,补偿一个数据点代替点A2;
因为经过野值剔除后,可能缺陷缺点的情况,或者由于雷达扫描和报出缺陷造成航迹点未能报出的情况。此时需要利用线性差值方法,将缺失的数据补齐。已知雷达扫描周期为T,设定航迹点A1位置为(x1,y1,v1,ω1),航迹点A3位置为(x3,y3,v2,ω2),其中(x1,y1)代表航迹点A1二维坐标数据,v1代表航迹点A1的速度,ω1代表航迹点A1的航向,(x3,y3)代表航迹点A3二维坐标数据,v2代表航迹点A1的速度,ω2代表航迹点A3的航向,采用线性差分方法,补偿一个数据点A2代替航迹点A2,则A2的位置(x2,y2)为:
步骤3包括:需要进行时间对齐:补齐后的航迹点数据分为被测雷达数据和标准雷达数据,由于两部雷达的报出时间不一定统一,因此需要对相近时间内的数据点进行时间对齐操作。因为两部雷达的时间差很短,小于等于两部雷达中最小扫描周期,可以认为是一个很短的时间,因此时间对齐操作采用线性插值的方法实现。设定被测雷达A探测到的目标点位置为A1,坐标为(xA,yA),速度、角度、时间分别为vA,cA,tA;标准雷达B探测到的目标位置为B1,坐标为(xB,yB),速度、角度、时间分别为vB,cB,tB,经过线性插分后得到与点B1在同一时间的点A’1,点A’1的位置为(x’A,y’A),计算公式如下:
步骤4包括:需要进行坐标转换。经过数据处理后得到完整的航迹点对应信息,由于雷达航迹信息此时为直角坐标系数据。为了得到相应的距离误差和角度误差需要进行坐标转换,将被测雷达A的位置(xRA,yRA)作为极坐标的极点。被测雷达A和标准雷达B,对同一目标的位置测量值分别为(xA,yA)、(xB,yB),需要将(xA,yA)、(xB,yB)转换到以(xRA,yRA)为极点的极坐标系中,得到(ρA,θA),(ρB,θB),其中ρA代表被测雷达A极坐标的距离,θA代表被测雷达A极坐标的角度,ρB代表标准雷达B极坐标的距离,θB代表标准雷达B极坐标的角度:
步骤5包括:需要进行距离误差计算:距离误差为测试雷达目标位置到极点的距离与标准雷达目标位置到极点的距离差值。距离误差存在正负值。设定目标航迹一共有K个点,则得到K个距离误差值(Δρ1,Δρ2,Δρ3…ΔρK),ΔρK表示第K个距离误差值;为了简化后续计算,需要对(Δρ1,Δρ2,Δρ3…ΔρK)进行归一化处理得到K个值(Δξ1,Δξ2,Δξ3…ΔξK):
其中,i=1,2,3…K,ΔξK表示ΔρK归一化处理后的值。
步骤6包括:需要进行角度误差计算:角度误差为测试雷达目标与极轴之间的夹角与标准雷达目标与极轴之间的夹角差值。角度误差存在正负值。设定目标航迹一共有K个点,则得到K个角度误差值(Δθ1,Δθ2,Δθ3…ΔθK),ΔθK表示第K个角度误差值,为了简化后续计算,需要对(Δθ1,Δθ2,Δθ3…ΔθK)进行归一化处理得到K个值(Δτ1,Δτ2,Δτ3…ΔτK),ΔτK表示ΔθK根据公式(4)归一化处理后的值。
步骤7包括:根据步骤4和步骤5计算得到两个点的距离误差和角度误差,分别对距离误差和角度误差其进行相干分析计算。例如,利用相干分析法计算距离抖动性:选取尺度N(N为大于1的正整数),如N取5,将5个航迹点的距离误差作为一组得到一个1*5维向量即X1=(Δρ1,Δρ2,Δρ3,Δρ4,Δρ5),选取连续的5组得到一个5阶矩阵即求矩阵的协方差矩阵E,求协方差矩阵E的特征值,根据相干系数计算公式得到相干系数。
获取到每个点的误差矩阵D,根据矩阵D求得相干系数,具体包括:
步骤7-1:对于m×n维矩阵D,矩阵D由连续航迹点的距离误差或者角度误差组成。求其n维协方差矩阵B;
步骤7-2:求协方差矩阵B的特征值λ1、λ2…λn,λn表示第n个特征值;
步骤7-3:根据如下公式计算相干系数k:
其中λmax为λ1、λ2…λn中的最大值;k值越大,表示数据越不相似;
步骤7-4,分别对航迹距离误差和航迹角度误差进行相干分析。
步骤7-4中,对航迹距离误差进行相干分析具体包括如下步骤:
步骤a1,计算航迹距离误差变化的多阶矩阵F1:
其中,ΔξN表示测试雷达和标准雷达的距离误差经过归一化后的值,K代表选定航迹点的总数,N表示尺度;
步骤a2,选择N×N大小的窗口进行相干分析,选取的窗口矩阵为D1,D1是从矩阵F1中截取的连续N列的矩阵:
其中XN代表D1中每一列的矩阵数据。
根据如下公式求矩阵D1的协方差矩阵E1:
其中cov(Xi,Xj)=E[(Xi-E(Xi))(Xj-E(Xj))],1≤i,j≤n;cov(Xi,Xj)代表向量Xi,Xj的协方差值。求矩阵E的特征值,根据公式(5)计算相干系数λ;
步骤a3,将窗口矩阵右移一格,求得当前窗口的相干系数,依次计算右移后新的窗口矩阵的相干系数,直到窗口到达矩阵最右端;
步骤a4,将相邻的相干系数连接起来得到距离误差变化相干系数曲线。
步骤7-4中,对航迹角度误差进行相干分析具体包括如下步骤:
步骤b1,计算航迹角度误差变化的多阶矩阵F2:
步骤b2,选择N×N大小的窗口进行相干分析,选取的窗口矩阵为D2,D2是从矩阵F2中截取的连续N列的矩阵:
其中XN代表D2中每一列的矩阵数据。
根据公式(6)求矩阵D2的协方差矩阵E2,求矩阵E2的特征值,根据公式(5)计算相干系数λ;
步骤b3,将窗口矩阵右移一格,求得当前窗口的相干系数,依次计算右移后新的窗口矩阵的相干系数,直到窗口到达矩阵最右端;
步骤b4,将相邻的相干系数连接起来得到角度误差变化相干系数曲线。
有益效果:本发明具有如下技术效果:
1、考虑了数据补偿和时间对齐准则,增强了数据的完整性和可用性,并且增加了最终相干系数曲线的平滑性和连续性。
2、从距离和角度两个方向较为全面的评估航迹抖动性。
3、在评估抖动性时使用了相干分析方法,并且设计了多尺度窗口矩阵的模式,增强了结果的鲁棒性和抗噪性。
4、采用了窗口矩阵的方法,借鉴了信号处理的分帧思想,前一个窗口矩阵和后一个窗口矩阵存在部分重叠区域,目的是为了使窗口矩阵在移动过程中能够平滑过渡。
附图说明
下面结合附图和具体实施方式对本发明做更进一步的具体说明,本发明的上述和/或其他方面的优点将会变得更加清楚。
图1是距离、角度误差计算示意图。
图2是多尺度窗口矩阵相干分析示意图。
图3是误差相干系数变化曲线示意图。
图4是航迹随机抖动性评估流程图示意图。
具体实施方式
如图4所示,本发明提供了一种雷达航迹随机抖动性评估方法,包括以下步骤:
步骤1、野值剔除。剔除雷达扫描到的出现明显偏离实际情况部分数据。例如现有雷达扫描到的同一目标的两个周期内的航迹点A1、A2。航迹点A1的速度为250m/s,航迹点A2的速度为800m/s,对于民航客机而言,A2点的速度明显超出正常范围,予以剔除。
步骤2、数据补偿。经过步骤一之后,可能出现如下情况,被测雷达现有航迹点A1、A2、A3,经过野值剔除后,剩下航迹点A1、A3,为了避免不同的时间尺度影响,因此需要采用线性差分方法,补偿一个数据点代替点A2。
假设雷达扫描周期为4s,现有航迹点A1(-4800m,2000m,250m/s,π/2)、A3(-3000,2000,240m/s,π/2),则航迹点A2的位置(x2,y2)为:
步骤3、时间对齐。被测雷达A和标准雷达B,对同一目标的位置测量值分别为(xA,yA)、(xB,yB),速度、角度、时间分别为(vA,cA,tA)、(vB,cB,tB)。由于两部雷达的报出时间不一定统一,因此需要对相近时间内的数据点进行时间对齐操作。因为两部雷达的时间差很短,小于等于两部雷达中最小扫描周期,可以认为是一个很短的时间,因此时间对齐操作采用线性插值的方法实现。
假设被测雷达A探测到的目标点位置为A1(-4800m,2000m),速度、角度、时间为(250m/s,π/2,123456(10ms),);标准雷达B探测到的目标位置为B1,时间为tB=123445(10ms),经过线性插分后得到与点B1在同一时间的点A’1,点A’1的位置为(X’A,y’A):
步骤4、坐标转换。被测雷达A和标准雷达B,对同一目标的位置测量值分别为(xA,yA)、(xB,yB)。雷达A的坐标为(xRA,yRA),需要将(xA,yA)、(xB,yB)转换到以(xRA,yRA)为极点的极坐标系中,得到(ρA,θA),(ρB,θB)。
假设被测雷达A探测到的目标点位置为A1(-4827.5m,2000m),标准雷达B探测到的目标位置为B1(-4900m,1900m),被测雷达A的坐标为(0,0)。经过极坐标转换公式(3)转换后得到A1极坐标(5225.4,2.749),B1点的极坐标为(5255.5,2.772)。
步骤5、距离误差计算。经过坐标转换后,被测雷达A和标准雷达B探测到的目标位置分别为(ρA,θA),(ρB,θB),计算两个目标到极点距离的差值。假设目标航迹一共有K个点,则得到K个距离误差值(Δρ1,Δρ2,Δρ3…ΔρK)。
步骤6、角度误差计算。经过坐标转换后,被测雷达A和标准雷达B探测到的目标位置分别为(ρA,θA),(ρB,θB),计算两个目标极角的差值。假设目标航迹一共有K个点,则得到K个角度误差值(Δθ1,Δθ2,Δθ3…ΔθK)。
步骤5和步骤6计算距离和角度误差数据,如图1所示,可以看出对同一目标被测雷达A和标准雷达B在两个周期内测得的目标位置分别为A1(ρA1,θA1)、A2(ρA2,θA2)和B1(ρB1,θB1)、B2(ρB2,θB2),A1、A2为经过了前期数据处理的位置。
例如现有A1,A2,A3,A4,A5五个测试航迹点,经过误差计算后得到(-30.1,-0.023),(-42.1,-0.02),(-26.5,-0.025),(-32,-0.02),(-40,-0.022)五组误差数据。根据公式4分别对距离误差和角度误差进行归一化处理得到:(0.769,0.4),(0,1),(1,0),(0.647,1),(0.135,0.4)。
步骤7、多尺度相干分析。根据步骤4和步骤5计算得到两个点的距离误差和角度误差。分别对距离误差和角度误差其进行相干分析计算。
假设窗口长度N=3,距离误差为(a1,a2,a3,…aK),则将(a1,a2,a3)的转置作为矩阵的第一列,然后右移一个数据点将(a2,a3,a4)的转置作为第二列,以此类推,直到移动到第K-2个点,得到(aK-2,aK-1,aK),将其转置作为矩阵的第K-2列,最终形成的N×(K-N+1)阶矩阵,如图2所示。
进行相干分析时,选取N×N大小的窗口,从矩阵的第一列开始根据公式5和公式6进行相干分析计算,得到第一个相干系数,然后窗口右移得到第二个窗口矩阵,直至窗口右移到矩阵的最右端,在图2上表现为窗口矩阵1、窗口矩阵2直到窗口矩阵6,对窗口矩阵分别求协方差矩阵1、协方差矩阵2等,然后求协方差矩阵的特征值,根据公式5求得相干系数。
求得EA的特征值为:1,1,1
根据公式5有k=1-1/(1+1+1)=0.67
本发明中,窗口矩阵选择有两个尺度,N=3或5。根据雷达扫描周期的不同和对不同窗口尺度的大小进行试验,发现当窗口大小设置为3或者5时,对航迹随机抖动性评估效果最好。
步骤8、抖动性曲线绘制。将窗口矩阵右移一个航迹点,得到一个重叠了((N-1)/N)的新矩阵,按照步骤7的方法重新求得相干系数。窗口继续右移,就会得到相干系数,将相邻的相干系数连接起来,就得到相干系数曲线。如图3所示,分别采用N=3和N=5两个尺度的窗口对测试航迹进行距离误差相干分析,得到两条误差相干系数曲线,从两条误差相干系数曲线可以看出,两个尺度的窗口对测试航迹的随机抖动性分析基本符合测试航迹和标准航迹的随机偏离度。
本发明提供了一种雷达航迹随机抖动性评估方法,具体实现该技术方案的方法和途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
Claims (10)
1.一种雷达航迹随机抖动性评估方法,其特征在于,包括如下步骤:
步骤1、剔除野值;
步骤2、数据补偿;
步骤3、进行时间对齐;
步骤4、进行坐标转换;
步骤5、计算距离误差;
步骤6、计算角度误差;
步骤7、进行多尺度相干分析;
步骤8、绘制抖动性曲线。
2.根据权利要求1所述的方法,其特征在于,步骤1中,设定一架民航飞机的最大速度为v(m/s),雷达扫描周期为T(s),则连续两个探测点的距离应当不大于v*T(m),设定一个随机误差固定容错率20%,则相邻两个点的距离应当小于v*T*1.2(m),将不满足上述条件的点剔除。
7.根据权利要求6所述的方法,其特征在于,步骤6包括:设定目标航迹一共有K个点,则得到K个角度误差值(Δθ1,Δθ2,Δθ3…ΔθK),ΔθK表示第K个角度误差值,对(Δθ1,Δθ2,Δθ3…ΔθK)进行归一化处理得到K个值(Δτ1,Δτ2,Δτ3…ΔτK),ΔτK表示ΔθK根据公式(4)归一化处理后的值。
9.根据权利要求8所述的方法,其特征在于,步骤7-4中,对航迹距离误差进行相干分析具体包括如下步骤:
步骤a1,计算航迹距离误差变化的多阶矩阵F1:
其中,ΔξN表示测试雷达和标准雷达的距离误差经过归一化后的值,K代表选定航迹点的总数,N表示尺度;
步骤a2,选择N×N大小的窗口进行相干分析,选取的窗口矩阵为D1,D1是从矩阵F1中截取的连续N列的矩阵:
其中XN代表D1中每一列的矩阵数据;
根据如下公式求矩阵D1的协方差矩阵E1:
其中cov(Xi,Xj)=E[(Xi-E(Xi))(Xj-E(Xj))],1≤i,j≤n;cov(Xi,Xj)代表向量Xi,Xj的协方差值;求矩阵E的特征值,根据公式(5)计算相干系数λ;
步骤a3,将窗口矩阵右移一格,求得当前窗口的相干系数,依次计算右移后新的窗口矩阵的相干系数,直到窗口到达矩阵最右端;
步骤a4,将相邻的相干系数连接起来得到距离误差变化相干系数曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011209612.3A CN112255612B (zh) | 2020-11-03 | 2020-11-03 | 一种雷达航迹随机抖动性评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011209612.3A CN112255612B (zh) | 2020-11-03 | 2020-11-03 | 一种雷达航迹随机抖动性评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112255612A true CN112255612A (zh) | 2021-01-22 |
CN112255612B CN112255612B (zh) | 2024-02-13 |
Family
ID=74268583
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011209612.3A Active CN112255612B (zh) | 2020-11-03 | 2020-11-03 | 一种雷达航迹随机抖动性评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112255612B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113838313A (zh) * | 2021-11-29 | 2021-12-24 | 中国民用航空总局第二研究所 | 一种航向信标航道余隙抖动的障碍物识别方法 |
CN113836123A (zh) * | 2021-07-22 | 2021-12-24 | 南京沃旭通讯科技有限公司 | 一种基于距离、角度的轨迹清洗方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000284049A (ja) * | 1999-03-31 | 2000-10-13 | Mitsubishi Electric Corp | レーダ統合追尾方式 |
CN108896973A (zh) * | 2018-07-13 | 2018-11-27 | 北京无线电测量研究所 | 雷达数据的标定方法、点迹质量评估方法、存储介质 |
CN110990756A (zh) * | 2019-12-03 | 2020-04-10 | 中国电子科技集团公司第二十八研究所 | 基于可拓型逐层分析的空中目标航迹质量评估方法 |
CN111366900A (zh) * | 2020-02-18 | 2020-07-03 | 上海机电工程研究所 | 基于残差统计的跟踪雷达航迹质量评估方法、***及介质 |
CN111477036A (zh) * | 2020-04-08 | 2020-07-31 | 中国电子科技集团公司第二十八研究所 | 一种空管自动化***航空器高度异常检测方法 |
CN111680870A (zh) * | 2020-04-29 | 2020-09-18 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 目标运动轨迹质量综合评估方法 |
CN111736172A (zh) * | 2020-08-24 | 2020-10-02 | 中国人民解放军国防科技大学 | 一种基于大气扰动相干激光探测的空中目标探测方法 |
-
2020
- 2020-11-03 CN CN202011209612.3A patent/CN112255612B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000284049A (ja) * | 1999-03-31 | 2000-10-13 | Mitsubishi Electric Corp | レーダ統合追尾方式 |
CN108896973A (zh) * | 2018-07-13 | 2018-11-27 | 北京无线电测量研究所 | 雷达数据的标定方法、点迹质量评估方法、存储介质 |
CN110990756A (zh) * | 2019-12-03 | 2020-04-10 | 中国电子科技集团公司第二十八研究所 | 基于可拓型逐层分析的空中目标航迹质量评估方法 |
CN111366900A (zh) * | 2020-02-18 | 2020-07-03 | 上海机电工程研究所 | 基于残差统计的跟踪雷达航迹质量评估方法、***及介质 |
CN111477036A (zh) * | 2020-04-08 | 2020-07-31 | 中国电子科技集团公司第二十八研究所 | 一种空管自动化***航空器高度异常检测方法 |
CN111680870A (zh) * | 2020-04-29 | 2020-09-18 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 目标运动轨迹质量综合评估方法 |
CN111736172A (zh) * | 2020-08-24 | 2020-10-02 | 中国人民解放军国防科技大学 | 一种基于大气扰动相干激光探测的空中目标探测方法 |
Non-Patent Citations (1)
Title |
---|
李民权: "飞行目标的抖动及雷达散射截面计算", 《合肥工业大学(自然科学版)》, vol. 26, no. 5, pages 1021 - 1024 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113836123A (zh) * | 2021-07-22 | 2021-12-24 | 南京沃旭通讯科技有限公司 | 一种基于距离、角度的轨迹清洗方法 |
CN113838313A (zh) * | 2021-11-29 | 2021-12-24 | 中国民用航空总局第二研究所 | 一种航向信标航道余隙抖动的障碍物识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112255612B (zh) | 2024-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106646450B (zh) | 基于距离分步聚类的雷达航迹抗差关联方法 | |
JP4989331B2 (ja) | 航跡統合装置及びプログラム及び航跡統合方法 | |
CN110058205B (zh) | 一种基于迭代最近点算法的警戒雷达***误差校正方法 | |
CN112255612A (zh) | 一种雷达航迹随机抖动性评估方法 | |
RU2524208C1 (ru) | Способ радиолокационного обнаружения маневра баллистической цели на пассивном участке траектории | |
CN106896363B (zh) | 一种水下目标主动跟踪航迹起始方法 | |
CN110045342B (zh) | 雷达相对***误差估值有效性评价方法 | |
CN103792522B (zh) | 基于可信关联对的多雷达海上目标抗差关联算法 | |
CN110738275A (zh) | 基于ut-phd的多传感器序贯融合跟踪方法 | |
CN112229406A (zh) | 一种多余度引导全自动着陆信息融合方法及*** | |
CN109856623B (zh) | 一种针对多雷达直线航迹线的目标状态估计方法 | |
CN113238208B (zh) | 非规则轨迹水中运动目标前向声散射多普勒频移计算方法 | |
EP2929366B1 (en) | Improvements in and relating to course and/or speed data | |
RU2721623C1 (ru) | Способ определения мгновенного положения точки промаха беспилотного летательного аппарата по информации угломерного канала | |
CN110196409B (zh) | 一种基于区域集合相对距离的抗差异步航迹关联方法 | |
CN112697340A (zh) | 一种固定翼飞机大气数据***及其故障检测方法 | |
CN110646792A (zh) | 一种基于观察哨数字望远镜的雷达搜索窗口设置方法 | |
CN107271974B (zh) | 一种基于稳定角点的时空误差求取方法 | |
CN109856624B (zh) | 一种针对单雷达直线航迹线的目标状态估计方法 | |
RU2714884C1 (ru) | Способ определения курса объекта на линейной траектории с использованием измерений его радиальной скорости | |
CN113705701A (zh) | 一种多源探测信号整合处理方法 | |
CN110780153B (zh) | 输电线路故障定位方法 | |
CN102645646B (zh) | 一种多信源不确定融合定位方法 | |
CN112198504B (zh) | 一种主被动观测特征交织的融合滤波方法 | |
CN111948638A (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 |