CN105681121A - 一种1090es ads-b消息报头的检测方法 - Google Patents
一种1090es ads-b消息报头的检测方法 Download PDFInfo
- Publication number
- CN105681121A CN105681121A CN201410855818.1A CN201410855818A CN105681121A CN 105681121 A CN105681121 A CN 105681121A CN 201410855818 A CN201410855818 A CN 201410855818A CN 105681121 A CN105681121 A CN 105681121A
- Authority
- CN
- China
- Prior art keywords
- sampling point
- time period
- sampling
- power
- condition
- 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
Landscapes
- Synchronisation In Digital Transmission Systems (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明提供一种1090ES?ADS-B消息前导脉冲检测方法,包括:将采样点按采样时序存入检测窗口的对应数据位;每当接收到新采样点,并将新采样点存入检测窗口的D1数据位后,判断条件1和条件2是否均满足,如均满足,执行后步;否则,继续采样;输出信号X1,显示位置为t3+t1;在输出信号X1之后,继续采样,并判断新采样点是否同时满足条件1和条件2,如满足,则持续输出X1,同时,继续采样;只要初次判断出后续的某个新采样点不同时满足条件1和条件2时,即输出信号X2,其显示位置即为报头中第4前导脉冲结束时刻的有效定位。为一种能量动态匹配检测方法,具有更高检测概率,减少因信噪比降低而引起漏检率上升问题。
Description
技术领域
本发明属于信号处理技术领域,具体涉及一种1090ESADS-B消息报头的检测方法。
背景技术
ADS-B,即:广播式自动相关监视技术,是现代集成先进的无线数据通信技术、卫星定位导航信息传输技术和自动监视技术等相关技术于一体的新兴空中交通管制技术,是国际民航组织(ICAO)正在积极推广的一种航行监视技术。
ADS-B系统是一个集通信与监视于一体的信息系统,包括:机载设备、数据链和地面接收站;机载设备通过发射机以周期方式向外发射自己的位置、高度、速度等数据信息,该数据信息以ADS-B消息形式,通过空-空、空-地数据链广播式传播。在地面接收站,对接收到的ADS-B消息进行解算,获得机载设备所发送的数据信息。
如图1所示,为ADS-B消息的信号格式图,包括报头和其后的数据块;其中,报头的长度为8us,包含4个前导脉冲,分别记为:第1前导脉冲、第2前导脉冲、第3前导脉冲和第4前导脉冲;每个前导脉冲的标准宽度均为0.5us;第1前导脉冲、第2前导脉冲、第3前导脉冲和第4前导脉冲的标准相对位置依次为:0.0~0.5us、1.0us~1.5us、3.5us~4.0us、4.5us~5.0us;对于地面接收站,在进行消息解算时,首先需要完成的操作是报头检测,这是一切后续处理的前提和基础,
现有的ADS-B地面接收站,主要采用脉冲位置检测方法进行报头检测,即:主要依靠上升沿标志、下降沿标志和有效脉冲位置标志相结合的检测方法,只有当上述三个标志均通过相应标准后,才能判定一个有效的报头。该种检测方法在信噪比较高时,具有较好的检测性能。但是,当信噪比相对较低时,传统检测方法的报头漏检率较高。例如,当4个前导脉冲的信号幅度发生较强波动时,采用脉冲位置检测方法并不能够准确有效地检测每个前导脉冲的上升沿、下降沿和有效脉冲位置。尤其对于微弱信号,因干扰很容易出现判断错误,从而导致报头漏检。
发明内容
针对现有技术存在的缺陷,本发明提供一种1090ESADS-B消息报头的检测方法,可有效解决上述问题。
本发明采用的技术方案如下:
本发明提供一种1090ESADS-B消息报头的检测方法,包括以下步骤:
S1,设置初始参数值;其中,所述初始参数值包括采样频率f1、检测窗口的窗口长度n、移位寄存器的长度m、0.5us时间段提取的采样点数w和采样频率f1下的延迟时间t1;
其中,检测窗口的窗口长度n代表该检测窗口共有n个数据位,分别记为:D1、D2…Dn;并且,n=5us*f1,即:检测窗口用于存储5us时间长度内的采样点;
移位寄存器的长度m代表该移位寄存器共有m个数据位,分别记为:E1、E2…Em;
S2,按照所述采样频率f1,不断采样输入的1090ESADS-B数字基带信号,并将采样点按采样时序存入所述检测窗口的对应数据位;其中,所述采样点指采样信号的功率值;
所述将采样点按采样时序存入所述检测窗口的对应数据位,具体指:
每当接收到新采样点时,新采样点存储于检测窗口的D1数据位,在新采样点之前的n-1个采样点按接收时间由近及远顺序,依次存储于检测窗口的D2、D3…Dn数据位;
S3,每当接收到新采样点,并将所述新采样点存入检测窗口的D1数据位后,以所述检测窗口的Dn数据位为时间基准0,所述检测窗口的D1数据位为5.0us,使所述检测窗口的各个数据位具有时间属性;然后,判断以下的条件1和条件2是否均满足,如果均满足,则执行S4;否则,返回S2,继续采样数字基带信号;
条件1:
(1)分别从所述检测窗口的0~0.5us时间段、1.0~1.5us时间段、3.5~4.0us时间段和4.5~5.0us时间段提取w个采样点,共得到w*4个采样点,计算w*4个采样点的功率值的和,将功率和值按时序存入移位寄存器;具体存储方法为:每当得到新功率和值时,新功率和值存储于移位寄存器的E1数据位,在新功率和值之前的m-1个功率和值按计算时间由近及远顺序,依次存储于移位寄存器的E2、E3…Em数据位;
(2)判断所述移位寄存器中当前m个数据位所存储的功率值是否满足E1>E2>E3…>Em;
条件2:
(1)分别计算所述检测窗口以下时间段的功率平均值:
0~0.5us时间段的所有采样点的功率平均值P1;
0.5~1.0us时间段的所有采样点的功率平均值P2;
1.0~1.5us时间段的所有采样点的功率平均值P3;
1.5~2.0us时间段的所有采样点的功率平均值P4;
3.0~3.5us时间段的所有采样点的功率平均值P5;
3.5~4.0us时间段的所有采样点的功率平均值P6;
4.0~4.5us时间段的所有采样点的功率平均值P7;
4.5~5.0us时间段的所有采样点的功率平均值P8;
还计算:
2.0~2.5us时间段的所有采样点的功率平均值P9;
2.5~3.0us时间段的所有采样点的功率平均值P10;
(2)判断以下关系式是否均成立;
1/2*P1>P2;
1/2*P3>P4;
P5<1/2*P6;
P7<1/2*P8;
P3>P9;
P10<1/2*P6;
1/2*P3<P8;
P9<P8;
S4,将所述新采样点记为C1,采样点C1的绝对采样时间为t3;
则:输出与所述采样点C1对应的第一类型提示信号X1,由于存在延迟时间t1,则:所述第一类型提示信号X1位于所述采样点C1之后t1时间长度的位置,该位置采样点的绝对采样时间为t3+t1;
S5:
S5.1:在输出所述第一类型提示信号X1之后,按照所述采样频率f1,继续采样数字基带信号,并更新检测窗口,将采样点存入所述检测窗口的D1数据位,在采样点之前的n-1个采样点按接收时间由近及远顺序,依次存储于检测窗口的D2、D3…Dn数据位;
S5.2,判断S3中的条件1和条件2是否均满足,如果满足,则继续输出第一类型提示信号X1,然后,返回S5.1;只要初次判断出后续的某个新采样点不同时满足条件1和条件2时,即执行S6;
S6,将初次不同时满足条件1和条件2的新采样点记为采样点C2,其绝对采样时间为t4;
则:输出与所述采样点C2对应的第二类型提示信号X2,由于存在延迟时间t1,则:所述第二类型提示信号X2位于所述采样点C2之后t1时间长度的位置,该位置采样点的绝对采样时间为t4+t1;
S7,查找到绝对时间为t4+t1的位置,该位置即为报头中第4前导脉冲结束时刻的有效定位。
优选的,采样频率f1为25MHz;检测窗口的窗口长度n为125;移位寄存器的长度m为3;0.5us时间段提取的采样点数w为10个;采样频率f1下的延迟时间t1为40ns。
优选的,还包括:
报头中第4前导脉冲结束时刻即为报头第5us时间点,获得第5us时间点所对应的采样点,向前按采样时间由近而远顺序,依次再获取n-1个采样点,由此定位到0~5us部分的报头信息,进而定位到第1前导脉冲、第2前导脉冲、第3前导脉冲和第4前导脉冲。
优选的,还包括:报文参考功率次大值和报文参考功率次小值的计算:
对于所述0~5us部分的报头信息:
(1)按设定规则提取位于0~0.5us的采样点,计算提取到的采样点的平均功率值,由此得到第1前导脉冲的平均功率值P100;
按设定规则提取位于1.0~1.5us的采样点,计算提取到的采样点的平均功率值,由此得到第2前导脉冲的平均功率值P101;
按设定规则提取位于3.5~4.0us的采样点,计算提取到的采样点的平均功率值,由此得到第3前导脉冲的平均功率值P102;
按设定规则提取位于4.5~50us的采样点,计算提取到的采样点的平均功率值,由此得到第4前导脉冲的平均功率值P103;
(2)在P100、P101、P102和P103中,基于比较算法,得到次大值功率和次小值功率,然后,比较次大值功率和次小值功率的大小关系,较大的功率值即作为报文参考功率次大值;较小的功率值即作为报文参考功率次小值。
优选的,所述比较算法具体为:
(1)比较P100和P101,若P100>=P101,则Pmax01=P100,Pmin01=P101;
(2)比较P102和P103,若P102>=P103,则Pmax23=P102,Pmin23=P103;
(3)比较Pmax01和Pmax23,若Pmax01>=Pmax23,则Pmax23为次大值功率;
(4)比较Pmin01和Pmin23,若Pmin01>=Pmin23,则Pmin01为次小值功率。
优选的,还包括:有效报头位置验证步骤,具体包括:
(1)将报文参考功率次大值记为P200,将报文参考功率次小值记为P300;依下式分别计算得到第一功率阈值P400和第二功率阈值P500;
P400=P200*(1+9/16);
P500=P300*(1-9/16);
(2)判断P100、P101、P102和P103是否均属于区间[P500,P400];如果是,则通过功率一致性验证,执行后续步骤(3);否则,丢弃该报头;
(3)后3us能量验证:
从定位到第4前导脉冲结束时刻开始,将该结束时刻记为5.0us,经过3.5us时间到达一帧1090ESADS-B信号的8.5us处;然后,在当前检测窗口内,存储从3.5us至8.5us时间段内的采样点;
分别提取3.5-4.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P800;
分别提取4.0-4.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P801;
分别提取4.5-5.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P802;
分别提取5.0-5.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P803;
分别提取5.5-6.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P804;
分别提取6.0-6.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P805;
分别提取6.5-7.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P806;
分别提取7.0-7.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P807;
分别提取7.5-8.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P808;
分别提取8.0-8.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P809;
然后,判断以下四个条件是否均满足,如果均满足,则通过报头后3us能量验证,执行后续步骤(4);否则,丢弃该报头;
条件1:1/2*P802>P808;
条件2:1/2*P802>P806;
条件3:1/2*P802>P805;
条件4:1/2*P802>P803;
(4)DF字段脉冲功率一致性验证
从定位到第4前导脉冲结束时刻开始,将该结束时刻记为5.0us,经过8.0us时间到达一帧1090ESADS-B信号的13.0us处;然后,在当前检测窗口内,存储从8.0us至13.0us时间段内的采样点;
分别提取8.0-8.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P900;
分别提取8.5-9.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P901;
分别提取9.0-9.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P902;
分别提取9.5-10.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P903;
分别提取10.0-10.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P904;
分别提取10.5-11.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P905;
分别提取11.0-11.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P906;
分别提取11.5-12.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P907;
分别提取12.0-12.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P908;
分别提取12.5-13.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P909;
然后,判断以下五个条件是否均满足,如果均满足,则通过DF字段脉冲功率一致性验证,并最终确定所定位到的位置为报头中第4前导脉冲结束时刻位置;否则,丢弃该报头;
条件1:判断P900和P901中的较大值是否属于区间[P500,P400];
条件2:判断P902和P903中的较大值是否属于区间[P500,P400];
条件3:判断P904和P905中的较大值是否属于区间[P500,P400];
条件4:判断P906和P907中的较大值是否属于区间[P500,P400];
条件5:判断P908和P909中的较大值是否属于区间[P500,P400]。
优选的,还包括:
对检测到的前后两个报头交织的情况进行判定,并对存在交织的两个报头进行取舍的步骤。
优选的,设被检测1090ESADS-B消息的一帧长度为Zus;
对检测到的前后两个报头交织的情况进行判定,并对存在交织的两个报头进行取舍的步骤,具体为:
S100,当检测到已经通过验证的有效报头A后,判断在有效报头A的5.0us位置之后的Zus-5.0us内,是否又检测到可通过验证的有效报头B;如果未检测到,则判定出不存在与有效报头A发生交织的情况,结束流程;如果检测到,则执行S101;
S101,判断有效报头B的参考功率次小值是否大于2.5倍的有效报头A的参考功率次大值,如果是,则得出有效报头B为强信号,交织于弱信号的有效报头A,则舍弃有效报头A,保留有效报头B;否则,舍弃有效报头B,保留有效报头A。
优选的,Z为64us或120us。
本发明的有益效果如下:
与传统检测方法相比,本发明具有更高的检测概率,有效减少因信噪比降低而引起的漏检率上升问题。
附图说明
图1为现有技术提供的ADS-B消息的信号格式图;
图2为本发明提供的有效报头确认过程中关系曲线对比图;
其中,1代表输入的1090ESADS-B数字基带信号;
2代表25MHz采样时钟下采样得到的数字信号;
3代表定位到的第4前导脉冲结束时刻的提示信号;
4代表定位到第4前导脉冲结束时刻的标记信号;
图3为图2中A区域的局部放大图;
5代表100MHz的时钟信号;
6代表25MHz的采样时钟信号;
图4为本发明提供的1090ESADS-B消息报头的检测方法的流程示意图。
具体实施方式
以下结合附图对本发明进行详细说明:
本发明提供一种1090ESADS-B消息报头的检测方法,核心思想为:
参考图1,为ADS-B消息的信号格式图,包括8us报头、56us或112us数据块。其中,8us报头包括0~5us时间段和5us~8us时间段,而DF字段位于8us~13us时间段。本发明主要通过检测0~5us时间段内脉冲能量是否动态匹配1090ESADS-B的标准脉冲能量,以及通过5us~8us时间段和8us~13us时间段内的DF字段的能量验证,验证检测结果是否正确。
本发明提供的1090ESADS-B消息报头的检测方法,核心步骤为:有效报头的确认、报文参考功率次大值和报文参考功率次小值的计算、有效报头的位置验证、交织情况判断及取舍。
以下分别详细介绍:
(一)有效报头的确认
如图4所示,包括以下步骤:
S1,设置初始参数值;其中,所述初始参数值包括采样频率f1、检测窗口的窗口长度n、移位寄存器的长度m、0.5us时间段提取的采样点数w和采样频率f1下的延迟时间t1;
其中,检测窗口的窗口长度n代表该检测窗口共有n个数据位,分别记为:D1、D2…Dn;并且,n=5us*f1,即:检测窗口用于存储5us时间长度内的采样点;
移位寄存器的长度m代表该移位寄存器共有m个数据位,分别记为:E1、E2…Em;
S2,按照所述采样频率f1,不断采样输入的1090ESADS-B数字基带信号,并将采样点按采样时序存入所述检测窗口的对应数据位;其中,所述采样点指采样信号的功率值;
所述将采样点按采样时序存入所述检测窗口的对应数据位,具体指:
每当接收到新采样点时,新采样点存储于检测窗口的D1数据位,在新采样点之前的n-1个采样点按接收时间由近及远顺序,依次存储于检测窗口的D2、D3…Dn数据位;
S3,每当接收到新采样点,并将所述新采样点存入检测窗口的D1数据位后,以所述检测窗口的Dn数据位为时间基准0,所述检测窗口的D1数据位为5.0us,使所述检测窗口的各个数据位具有时间属性;然后,判断以下的条件1和条件2是否均满足,如果均满足,则执行S4;否则,返回S2,继续采样数字基带信号;
条件1:
(1)分别从所述检测窗口的0~0.5us时间段、1.0~1.5us时间段、3.5~4.0us时间段和4.5~5.0us时间段提取w个采样点,共得到w*4个采样点,计算w*4个采样点的功率值的和,将功率和值按时序存入移位寄存器;具体存储方法为:每当得到新功率和值时,新功率和值存储于移位寄存器的E1数据位,在新功率和值之前的m-1个功率和值按计算时间由近及远顺序,依次存储于移位寄存器的E2、E3…Em数据位;
(2)判断所述移位寄存器中当前m个数据位所存储的功率值是否满足E1>E2>E3…>Em;
条件2:
(1)分别计算所述检测窗口以下时间段的功率平均值:
0~0.5us时间段的所有采样点的功率平均值P1;
0.5~1.0us时间段的所有采样点的功率平均值P2;
1.0~1.5us时间段的所有采样点的功率平均值P3;
1.5~2.0us时间段的所有采样点的功率平均值P4;
3.0~3.5us时间段的所有采样点的功率平均值P5;
3.5~4.0us时间段的所有采样点的功率平均值P6;
4.0~4.5us时间段的所有采样点的功率平均值P7;
4.5~5.0us时间段的所有采样点的功率平均值P8;
还计算:
2.0~2.5us时间段的所有采样点的功率平均值P9;
2.5~3.0us时间段的所有采样点的功率平均值P10;
(2)判断以下关系式是否均成立;
1/2*P1>P2;
1/2*P3>P4;
P5<1/2*P6;
P7<1/2*P8;
P3>P9;
P10<1/2*P6;
1/2*P3<P8;
P9<P8;
S4,将所述新采样点记为C1,采样点C1的绝对采样时间为t3;
则:输出与所述采样点C1对应的第一类型提示信号X1,由于存在延迟时间t1,则:所述第一类型提示信号X1位于所述采样点C1之后t1时间长度的位置,该位置采样点的绝对采样时间为t3+t1;
S5,在输出所述第一类型提示信号X1之后,返回S2,继续采样数字基带信号,并判断新采样点是否同时满足条件1和条件2,如果满足,则持续输出第一类型提示信号X1,同时,继续采样数字基带信号,循环执行S2-S5;只要初次判断出后续的某个新采样点不同时满足条件1和条件2时,即执行S6;
S6,将初次不同时满足条件1和条件2的新采样点记为采样点C2,其绝对采样时间为t4;
则:输出与所述采样点C2对应的第二类型提示信号X2,由于存在延迟时间t1,则:所述第二类型提示信号X2位于所述采样点C2之后t1时间长度的位置,该位置采样点的绝对采样时间为t4+t1;
S7,查找到绝对时间为t4+t1的位置,该位置即为报头中第4前导脉冲结束时刻的有效定位。
报头中第4前导脉冲结束时刻即为报头第5us时间点,获得第5us时间点所对应的采样点,向前按采样时间由近而远顺序,依次再获取n-1个采样点,由此定位到0~5us部分的报头信息,进而定位到第1前导脉冲、第2前导脉冲、第3前导脉冲和第4前导脉冲。
实施例:
以下介绍上述有效报头确认的一个实施例:
S1:设置的初始参数值如下:
采样频率f1为25MHz;
检测窗口的窗口长度n为125;检测窗口可采用一个125位的移位寄存器,为避免与用于存储功率和值的移位寄存器混淆,本例中,仍然采用检测窗口这个词进行介绍;
移位寄存器的长度m为3;
0.5us时间段提取的采样点数w为10个;
采样频率f1下的延迟时间t1为40ns。
则:检测窗口共有125个数据位,分别记为D1、D2…D125;
移位寄存器共有3个数据位,分别记为:E1、E2和E3;
S2,按照25MHz的采样频率,不断采样输入基准频率为100MHz的1090ESADS-B数字基带信号,每当接收到新采样点时,新采样点存储于检测窗口的D1数据位,在新采样点之前的124个采样点按接收时间由近及远顺序,依次存储于检测窗口的D2、D3…D125数据位;
S3,每当接收到新采样点,并将新采样点存入检测窗口的D1数据位后,由于检测窗口所存储的125个数据即对应5us长度的采样数据,因此,判断以下的条件1和条件2是否均满足,如果均满足,则执行S4;否则,返回S2,继续采样数字基带信号;
条件1:(1)分别从所述检测窗口的0~0.5us时间段、1.0~1.5us时间段、3.5~4.0us时间段和4.5~5.0us时间段提取10个采样点,共得到40个采样点,计算40个采样点的功率值的和,将功率和值按时序存入移位寄存器;具体存储方法为:每当得到新功率和值时,新功率和值存储于移位寄存器的E1数据位,在新功率和值之前的2个功率和值按计算时间由近及远顺序,依次存储于移位寄存器的E2和E3数据位;
(2)判断所述移位寄存器中当前3个数据位所存储的功率值是否满足E1>E2>E3;
条件2:
(1)分别计算所述检测窗口以下时间段的功率平均值:
0~0.5us时间段的所有采样点的功率平均值P1;
0.5~1.0us时间段的所有采样点的功率平均值P2;
1.0~1.5us时间段的所有采样点的功率平均值P3;
1.5~2.0us时间段的所有采样点的功率平均值P4;
3.0~3.5us时间段的所有采样点的功率平均值P5;
3.5~4.0us时间段的所有采样点的功率平均值P6;
4.0~4.5us时间段的所有采样点的功率平均值P7;
4.5~5.0us时间段的所有采样点的功率平均值P8;
还计算:
2.0~2.5us时间段的所有采样点的功率平均值P9;
2.5~3.0us时间段的所有采样点的功率平均值P10;
(2)判断以下关系式是否均成立;
1/2*P1>P2;
1/2*P3>P4;
P5<1/2*P6;
P7<1/2*P8;
P3>P9;
P10<1/2*P6;
1/2*P3<P8;
P9<P8;
S4,将所述新采样点记为C1,采样点C1的绝对采样时间为t3;
则:输出与所述采样点C1对应的第一类型提示信号X1,由于存在延迟时间t1,则:所述第一类型提示信号X1位于所述采样点C1之后t1时间长度的位置,该位置采样点的绝对采样时间为t3+t1;
S5:
S5.1:在输出所述第一类型提示信号X1之后,按照所述采样频率f1,继续采样数字基带信号,并更新检测窗口,将采样点存入所述检测窗口的D1数据位,在采样点之前的n-1个采样点按接收时间由近及远顺序,依次存储于检测窗口的D2、D3…Dn数据位;
S5.2,判断S3中的条件1和条件2是否均满足,如果满足,则继续输出第一类型提示信号X1,然后,返回S5.1;只要初次判断出后续的某个新采样点不同时满足条件1和条件2时,即执行S6;
S6,将初次不同时满足条件1和条件2的新采样点记为采样点C2,其绝对采样时间为t4;
则:输出与所述采样点C2对应的第二类型提示信号X2,由于存在延迟时间t1,则:所述第二类型提示信号X2位于所述采样点C2之后t1时间长度的位置,该位置采样点的绝对采样时间为t4+t1;
S7,查找到绝对时间为t4+t1的位置,该位置即为报头中第4前导脉冲结束时刻的有效定位。
参考图3和图4:
(1)假设当前为时钟周期T1,寄存器中E1、E2和E3存储的数据分别为:6574、5368和4223,由于为递增顺序,因此,满足条件1;此时,对当前检测窗口进行条件2判断时,条件2不满足;
(2)因此,当下一个时钟周期T2到达时,采样下一个采样点G1并存入检测窗口D1位置,此时,计算求得功率和值为7873,则寄存器中E1、E2和E3更新为:7873、6574和5368,由于为递增顺序,因此,满足条件1;此时,对当前检测窗口进行条件2判断时,条件2不满足;
(3)因此,当下一个时钟周期T3到达时,采样下一个采样点G2并存入检测窗口D1位置,此时,计算求得功率和值为9177,则寄存器中E1、E2和E3更新为:9177、7873和6574,由于为递增顺序,因此,满足条件1;此时,对当前检测窗口进行条件2判断时,条件2满足;
则输出当前采样点G2对应的第一类型提示信号X1,例如,输出高电平信号,由于通常会存在延迟一个时钟周期T的情况,因此,第一类型提示信号X1显示在位于采样点G2之后的G3采样点位置,G3采样点对应的绝对采样时间为T4=T3+T;
(4)当下一个时钟周期T4到达时,采样下一个采样点G3并存入检测窗口D1位置,此时,计算求得功率和值为10430,则寄存器中E1、E2和E3更新为:10430、9177和7873,由于为递增顺序,因此,满足条件1;此时,对当前检测窗口进行条件2判断时,条件2满足;
T5和T6时钟周期也均满足条件1和条件2,因此,持续输出高电平信号;
(5)当时钟周期T7达到时,新接收到的采样点为G6,寄存器中E1、E2和E3更新为:11635、12013和11509,此时,功率和比前一时钟降低,已不满足条件1,则立刻输出第二类型提示信号X2,例如,输出低电平信号,由于通常会存在延迟一个时钟周期T的情况,因此,第二类型提示信号X2显示在位于采样点G6之后的G7采样点位置,G7采样点对应的绝对采样时间为T8=T7+T。T8即为第4前导脉冲结束时刻,此时,即确定了第4前导脉冲结束时刻的位置。
(二)报文参考功率次大值和报文参考功率次小值的计算
对于上述步骤定位到的0~5us部分的报头信息:
(1)按设定规则提取位于0~0.5us的采样点,计算提取到的采样点的平均功率值,由此得到第1前导脉冲的平均功率值P100;
按设定规则提取位于1.0~1.5us的采样点,计算提取到的采样点的平均功率值,由此得到第2前导脉冲的平均功率值P101;
按设定规则提取位于3.5~4.0us的采样点,计算提取到的采样点的平均功率值,由此得到第3前导脉冲的平均功率值P102;
按设定规则提取位于4.5~50us的采样点,计算提取到的采样点的平均功率值,由此得到第4前导脉冲的平均功率值P103;
上次设定规则为:在选取采样点计算各个前导脉冲功率值时,需选择每个前导脉冲中处于高电平时刻附近的采样点,并可根据实际需求,灵活设定每个前导脉冲选择的采样点数量。
(2)在P100、P101、P102和P103中,基于比较算法,得到次大值功率和次小值功率,然后,比较次大值功率和次小值功率的大小关系,较大的功率值即作为报文参考功率次大值;较小的功率值即作为报文参考功率次小值。
下面介绍一种具体的比较算法:
(1)比较P100和P101,若P100>=P101,则Pmax01=P100,Pmin01=P101;
(2)比较P102和P103,若P102>=P103,则Pmax23=P102,Pmin23=P103;
(3)比较Pmax01和Pmax23,若Pmax01>=Pmax23,则Pmax23为次大值功率;
(4)比较Pmin01和Pmin23,若Pmin01>=Pmin23,则Pmin01为次小值功率。
仍以前述有效报头确认的实施例为例:
(1)当定位到0~5us部分的报头信息后,在0~5us时段内,共有125个采样点,按采样时序由近而远顺序,依次记为:第[0]点、第[1]点…第[124]点。根据4个前导脉冲形状的特点:即:0.05us≤脉冲上升沿的时间≤0.1us;0.05us≤脉冲下降沿的时间≤0.2us。在选取采样点计算各个前导脉冲功率值时,需选择每个前导脉冲中处于高电平时刻附近的采样点,本例如,每个前导脉冲取8个采样点,分别对应的采样点序号为:
对于第4前导脉冲,即“4.5-5.0us”处,选取的采样点序号为:[3],[4],[5],[6],[7],[8],[9],[10];
对于第3前导脉冲,即“3.5-4.0us”处,选取的采样点序号为:[28],[29],[30],[31],[32],[33],[34],[35];
对于第2前导脉冲,即:“1.0-1.5us”处,选取的采样点序号为:[90],[91],[92],[93],[94],[95],[96],[97];
对于第1前导脉冲,即:“0-0.5us”处,选取的采样点序号为:[115],[116],[117],[118],[119],[120],[121],[122]。
然后,分别计算第1前导脉冲的平均功率值、第2前导脉冲的平均功率值、第3前导脉冲的平均功率值和第4前导脉冲的平均功率值;
再在4个前导脉冲的平均功率值中,去除平均功率值最大和最小的前导脉冲,将剩余的2个平均功率值作为4个前导脉冲参考功率的次大值和次小值。
计算出报头参考功率的次大值和次小值,用于为后续有效报头位置验证提供能量比较依据。
(三)有效报头的位置验证
(1)将报文参考功率次大值记为P200,将报文参考功率次小值记为P300;依下式分别计算得到第一功率阈值P400和第二功率阈值P500;
P400=P200*(1+9/16);
P500=P300*(1-9/16);
(2)判断P100、P101、P102和P103是否均属于区间[P500,P400];如果是,则通过功率一致性验证,执行后续步骤(3);否则,丢弃该报头;
此处需要强调的是,本例中,对报头参考功率的次大值和报头参考功率的次小值按一定比例放大和缩小作为参考比较门限,此处,放大和缩小的比例,为本发明人经过多次试验,反复尝试而得,当采用这个放大和缩小的比例时,可以有效提高报头位置的验证准确率。
(3)后3us能量验证:
从定位到第4前导脉冲结束时刻开始,将该结束时刻记为5.0us,经过3.5us时间到达一帧1090ESADS-B信号的8.5us处;然后,在当前检测窗口内,存储从3.5us至8.5us时间段内的采样点;
分别提取3.5-4.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P800;
分别提取4.0-4.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P801;
分别提取4.5-5.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P802;
分别提取5.0-5.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P803;
分别提取5.5-6.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P804;
分别提取6.0-6.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P805;
分别提取6.5-7.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P806;
分别提取7.0-7.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P807;
分别提取7.5-8.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P808;
分别提取8.0-8.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P809;
仍以前述有效报头确认的实施例为例,则:设定数量的采样点为8个,采样点选择原则仍然为:将对应时段中位于边界的采样点去除,尽量选择位于时段中间附近的采样点。即:
对于3.5us至8.5us时间段,共5us时间长度,共有125个采样点,按采样时序由近而远顺序,依次记为:第[0]点、第[1]点…第[124]点。每1us内有两个码片,每个码片内分别取8个采样点:
“3.5-4.0us”时段内,选取的采样点序号为:[122],[121],[120],[119],[118],[117],[116],[115];
“4.0-4.5us”时段内,选取的采样点序号为:[109],[108],[107],[106],[105],[104],[103],[102];
“4.5-5.0us”时段内,选取的采样点序号为:[97],[96],[95],[94],[93],[92],[91],[90];
“5.0-5.5us”时段内,选取的采样点序号为:[84],[83],[82],[81],[80],[79],[78],[77];
“5.5-6.0us”时段内,选取的采样点序号为:[72],[71],[70],[69],[68],[67],[66],[65];
“6.0-6.5us”时段内,选取的采样点序号为:[59],[58],[57],[56],[55],[54],[53],[52];
“6.5-7.0us”时段内,选取的采样点序号为:[47],[46],[45],[44],[43],[42],[41],[40];
“7.0-7.5us”时段内,选取的采样点序号为:[34],[33],[32],[31],[30],[29],[28],[27];
“7.5-8.0us”时段内,选取的采样点序号为:[22],[21],[20],[19],[18],[17],[16],[15];
“8.0-8.5us”时段内,选取的采样点序号为:[9],[8],[7],[6],[5],[4],[3],[2]。
然后,判断以下四个条件是否均满足,如果均满足,则通过报头后3us能量验证,执行后续步骤(4);否则,丢弃该报头;
条件1:1/2*P802>P808;
条件2:1/2*P802>P806;
条件3:1/2*P802>P805;
条件4:1/2*P802>P803。
(4)DF字段脉冲功率一致性验证
从定位到第4前导脉冲结束时刻开始,将该结束时刻记为5.0us,经过8.0us时间到达一帧1090ESADS-B信号的13.0us处;然后,在当前检测窗口内,存储从8.0us至13.0us时间段内的采样点;
分别提取8.0-8.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P900;
分别提取8.5-9.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P901;
分别提取9.0-9.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P902;
分别提取9.5-10.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P903;
分别提取10.0-10.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P904;
分别提取10.5-11.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P905;
分别提取11.0-11.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P906;
分别提取11.5-12.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P907;
分别提取12.0-12.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P908;
分别提取12.5-13.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P909;
仍以前述有效报头确认的实施例为例,则:设定数量的采样点为8个,采样点选择原则仍然为:将对应时段中位于边界的采样点去除,尽量选择位于时段中间附近的采样点。即:
对于8.0us至13.0us时间段,共5us时间长度,共有125个采样点,按采样时序由近而远顺序,依次记为:第[0]点、第[1]点…第[124]点。每1us内有两个码片,每个码片内分别取8个采样点:
“8.0-8.5us”时段内,选取的采样点序号为:[122],[121],[120],[119],[118],[117],[116],[115];
“8.5-9.0us”时段内,选取的采样点序号为:[109],[108],[107],[106],[105],[104],[103],[102];
“9.0-9.5us”时段内,选取的采样点序号为:[97],[96],[95],[94],[93],[92],[91],[90];
“9.5-10.0us”时段内,选取的采样点序号为:[84],[83],[82],[81],[80],[79],[78],[77];
“10.0-10.5us”时段内,选取的采样点序号为:[72],[71],[70],[69],[68],[67],[66],[65];
“10.5-11.0us”时段内,选取的采样点序号为:[59],[58],[57],[56],[55],[54],[53],[52];
“11.0-11.5us”时段内,选取的采样点序号为:[47],[46],[45],[44],[43],[42],[41],[40];
“11.5-12.0us”时段内,选取的采样点序号为:[34],[33],[32],[31],[30],[29],[28],[27];
“12.0-12.5us”时段内,选取的采样点序号为:[22],[21],[20],[19],[18],[17],[16],[15];
“12.5-13.0us”时段内,选取的采样点序号为:[9],[8],[7],[6],[5],[4],[3],[2]。
然后,判断以下五个条件是否均满足,如果均满足,则通过DF字段脉冲功率一致性验证,并最终确定所定位到的位置为报头中第4前导脉冲结束时刻位置;否则,丢弃该报头;
条件1:判断P900和P901中的较大值是否属于区间[P500,P400];
条件2:判断P902和P903中的较大值是否属于区间[P500,P400];
条件3:判断P904和P905中的较大值是否属于区间[P500,P400];
条件4:判断P906和P907中的较大值是否属于区间[P500,P400];
条件5:判断P908和P909中的较大值是否属于区间[P500,P400]。
(四)交织情况判断及取舍
设被检测1090ESADS-B消息的一帧长度为Zus,现有技术中,主要存在两种类型的消息,一种为64us长度的消息,即:报头长度为8us,数据块长度为56us;另一种为120us长度的消息,即:报头长度为8us,数据块长度为112us。
对检测到的前后两个报头交织的情况进行判定,并对存在交织的两个报头进行取舍的步骤,具体为:
S100,当检测到已经通过验证的有效报头A后,判断在有效报头A的5.0us位置之后的59us或115us内,是否又检测到可通过验证的有效报头B;如果未检测到,则判定出不存在与有效报头A发生交织的情况,结束流程;如果检测到,则执行S101;
S101,判断有效报头B的参考功率次小值是否大于2.5倍的有效报头A的参考功率次大值,如果是,则得出有效报头B为强信号,交织于弱信号的有效报头A,则舍弃有效报头A,保留有效报头B;否则,舍弃有效报头B,保留有效报头A。
比较例:
传统检测方法为:依据3个重要标志:有效脉冲位置标志VPP、每个比特位信号的上升沿位置标志LEP、每个比特位信号的下降沿位置标志FEP进行检测。
将传统检测方式和本发明提供的检测方法,在不同信噪比下进行对比检测试验,试验结果表明,在低信噪比下,本发明检测方法的漏检率明显低于传统检测方法。
本发明提供的1090ESADS-B消息前导脉冲检测方法,具有以下优点:
(1)实质为一种能量动态匹配检测方法,充分利用了1090ESADS-B消息前导部分前5us内有脉冲能量和后3us内无脉冲能量的特点,定位报头位置,并通过消息数据块部分中DF字段脉冲功率一致性验证,从而确保1090ESADS-B消息报头检测的严谨性,由于为一种能量动态匹配检测方法,不需要进行传统技术中的上升沿或下降沿匹配,因此,具有更高的检测概率,有效减少因信噪比降低而引起的漏检率上升问题。
(2)本发明提供的检测方法,既适用于不存在交织时的报头检测,也适用于存在交织时的报头检测,适用范围广。
(3)本发明检测交织的时间范围大,是从检测到前一帧消息的报头后就开始交织点检测,直到前一帧信号所有数据位结束,而非传统技术中的0~13us范围内检测交织。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视本发明的保护范围。
Claims (9)
1.一种1090ESADS-B消息报头的检测方法,其特征在于,包括以下步骤:
S1,设置初始参数值;其中,所述初始参数值包括采样频率f1、检测窗口的窗口长度n、移位寄存器的长度m、0.5us时间段提取的采样点数w和采样频率f1下的延迟时间t1;
其中,检测窗口的窗口长度n代表该检测窗口共有n个数据位,分别记为:D1、D2…Dn;并且,n=5us*f1,即:检测窗口用于存储5us时间长度内的采样点;
移位寄存器的长度m代表该移位寄存器共有m个数据位,分别记为:E1、E2…Em;
S2,按照所述采样频率f1,不断采样输入的1090ESADS-B数字基带信号,并将采样点按采样时序存入所述检测窗口的对应数据位;其中,所述采样点指采样信号的功率值;
所述将采样点按采样时序存入所述检测窗口的对应数据位,具体指:
每当接收到新采样点时,新采样点存储于检测窗口的D1数据位,在新采样点之前的n-1个采样点按接收时间由近及远顺序,依次存储于检测窗口的D2、D3…Dn数据位;
S3,每当接收到新采样点,并将所述新采样点存入检测窗口的D1数据位后,以所述检测窗口的Dn数据位为时间基准0,所述检测窗口的D1数据位为5.0us,使所述检测窗口的各个数据位具有时间属性;然后,判断以下的条件1和条件2是否均满足,如果均满足,则执行S4;否则,返回S2,继续采样数字基带信号;
条件1:
(1)分别从所述检测窗口的0~0.5us时间段、1.0~1.5us时间段、3.5~4.0us时间段和4.5~5.0us时间段提取w个采样点,共得到w*4个采样点,计算w*4个采样点的功率值的和,将功率和值按时序存入移位寄存器;具体存储方法为:每当得到新功率和值时,新功率和值存储于移位寄存器的E1数据位,在新功率和值之前的m-1个功率和值按计算时间由近及远顺序,依次存储于移位寄存器的E2、E3…Em数据位;
(2)判断所述移位寄存器中当前m个数据位所存储的功率值是否满足E1>E2>E3…>Em;
条件2:
(1)分别计算所述检测窗口以下时间段的功率平均值:
0~0.5us时间段的所有采样点的功率平均值P1;
0.5~1.0us时间段的所有采样点的功率平均值P2;
1.0~1.5us时间段的所有采样点的功率平均值P3;
1.5~2.0us时间段的所有采样点的功率平均值P4;
3.0~3.5us时间段的所有采样点的功率平均值P5;
3.5~4.0us时间段的所有采样点的功率平均值P6;
4.0~4.5us时间段的所有采样点的功率平均值P7;
4.5~5.0us时间段的所有采样点的功率平均值P8;
还计算:
2.0~2.5us时间段的所有采样点的功率平均值P9;
2.5~3.0us时间段的所有采样点的功率平均值P10;
(2)判断以下关系式是否均成立;
1/2*P1>P2;
1/2*P3>P4;
P5<1/2*P6;
P7<1/2*P8;
P3>P9;
P10<1/2*P6;
1/2*P3<P8;
P9<P8;
S4,将所述新采样点记为C1,采样点C1的绝对采样时间为t3;
则:输出与所述采样点C1对应的第一类型提示信号X1,由于存在延迟时间t1,则:所述第一类型提示信号X1位于所述采样点C1之后t1时间长度的位置,该位置采样点的绝对采样时间为t3+t1;
S5:
S5.1:在输出所述第一类型提示信号X1之后,按照所述采样频率f1,继续采样数字基带信号,并更新检测窗口,将采样点存入所述检测窗口的D1数据位,在采样点之前的n-1个采样点按接收时间由近及远顺序,依次存储于检测窗口的D2、D3…Dn数据位;
S5.2,判断S3中的条件1和条件2是否均满足,如果满足,则继续输出第一类型提示信号X1,然后,返回S5.1;只要初次判断出后续的某个新采样点不同时满足条件1和条件2时,即执行S6;
S6,将初次不同时满足条件1和条件2的新采样点记为采样点C2,其绝对采样时间为t4;
则:输出与所述采样点C2对应的第二类型提示信号X2,由于存在延迟时间t1,则:所述第二类型提示信号X2位于所述采样点C2之后t1时间长度的位置,该位置采样点的绝对采样时间为t4+t1;
S7,查找到绝对时间为t4+t1的位置,该位置即为报头中第4前导脉冲结束时刻的有效定位。
2.根据权利要求1所述的1090ESADS-B消息报头的检测方法,其特征在于,采样频率f1为25MHz;检测窗口的窗口长度n为125;移位寄存器的长度m为3;0.5us时间段提取的采样点数w为10个;采样频率f1下的延迟时间t1为40ns。
3.根据权利要求1所述的1090ESADS-B消息报头的检测方法,其特征在于,还包括:
报头中第4前导脉冲结束时刻即为报头第5us时间点,获得第5us时间点所对应的采样点,向前按采样时间由近而远顺序,依次再获取n-1个采样点,由此定位到0~5us部分的报头信息,进而定位到第1前导脉冲、第2前导脉冲、第3前导脉冲和第4前导脉冲。
4.根据权利要求3所述的1090ESADS-B消息报头的检测方法,其特征在于,还包括:报文参考功率次大值和报文参考功率次小值的计算:
对于所述0~5us部分的报头信息:
(1)按设定规则提取位于0~0.5us的采样点,计算提取到的采样点的平均功率值,由此得到第1前导脉冲的平均功率值P100;
按设定规则提取位于1.0~1.5us的采样点,计算提取到的采样点的平均功率值,由此得到第2前导脉冲的平均功率值P101;
按设定规则提取位于3.5~4.0us的采样点,计算提取到的采样点的平均功率值,由此得到第3前导脉冲的平均功率值P102;
按设定规则提取位于4.5~50us的采样点,计算提取到的采样点的平均功率值,由此得到第4前导脉冲的平均功率值P103;
(2)在P100、P101、P102和P103中,基于比较算法,得到次大值功率和次小值功率,然后,比较次大值功率和次小值功率的大小关系,较大的功率值即作为报文参考功率次大值;较小的功率值即作为报文参考功率次小值。
5.根据权利要求4所述的1090ESADS-B消息报头的检测方法,其特征在于,所述比较算法具体为:
(1)比较P100和P101,若P100>=P101,则Pmax01=P100,Pmin01=P101;
(2)比较P102和P103,若P102>=P103,则Pmax23=P102,Pmin23=P103;
(3)比较Pmax01和Pmax23,若Pmax01>=Pmax23,则Pmax23为次大值功率;
(4)比较Pmin01和Pmin23,若Pmin01>=Pmin23,则Pmin01为次小值功率。
6.根据权利要求4所述的1090ESADS-B消息报头的检测方法,其特征在于,还包括:有效报头位置验证步骤,具体包括:
(1)将报文参考功率次大值记为P200,将报文参考功率次小值记为P300;依下式分别计算得到第一功率阈值P400和第二功率阈值P500;
P400=P200*(1+9/16);
P500=P300*(1-9/16);
(2)判断P100、P101、P102和P103是否均属于区间[P500,P400];如果是,则通过功率一致性验证,执行后续步骤(3);否则,丢弃该报头;
(3)后3us能量验证:
从定位到第4前导脉冲结束时刻开始,将该结束时刻记为5.0us,经过3.5us时间到达一帧1090ESADS-B信号的8.5us处;然后,在当前检测窗口内,存储从3.5us至8.5us时间段内的采样点;
分别提取3.5-4.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P800;
分别提取4.0-4.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P801;
分别提取4.5-5.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P802;
分别提取5.0-5.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P803;
分别提取5.5-6.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P804;
分别提取6.0-6.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P805;
分别提取6.5-7.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P806;
分别提取7.0-7.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P807;
分别提取7.5-8.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P808;
分别提取8.0-8.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P809;
然后,判断以下四个条件是否均满足,如果均满足,则通过报头后3us能量验证,执行后续步骤(4);否则,丢弃该报头;
条件1:1/2*P802>P808;
条件2:1/2*P802>P806;
条件3:1/2*P802>P805;
条件4:1/2*P802>P803;
(4)DF字段脉冲功率一致性验证
从定位到第4前导脉冲结束时刻开始,将该结束时刻记为5.0us,经过8.0us时间到达一帧1090ESADS-B信号的13.0us处;然后,在当前检测窗口内,存储从8.0us至13.0us时间段内的采样点;
分别提取8.0-8.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P900;
分别提取8.5-9.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P901;
分别提取9.0-9.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P902;
分别提取9.5-10.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P903;
分别提取10.0-10.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P904;
分别提取10.5-11.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P905;
分别提取11.0-11.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P906;
分别提取11.5-12.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P907;
分别提取12.0-12.5us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P908;
分别提取12.5-13.0us时间段内设定数量的采样点,并计算得到所提取的采样点的平均功率值P909;
然后,判断以下五个条件是否均满足,如果均满足,则通过DF字段脉冲功率一致性验证,并最终确定所定位到的位置为报头中第4前导脉冲结束时刻位置;否则,丢弃该报头;
条件1:判断P900和P901中的较大值是否属于区间[P500,P400];
条件2:判断P902和P903中的较大值是否属于区间[P500,P400];
条件3:判断P904和P905中的较大值是否属于区间[P500,P400];
条件4:判断P906和P907中的较大值是否属于区间[P500,P400];
条件5:判断P908和P909中的较大值是否属于区间[P500,P400]。
7.根据权利要求6所述的1090ESADS-B消息报头的检测方法,其特征在于,还包括:
对检测到的前后两个报头交织的情况进行判定,并对存在交织的两个报头进行取舍的步骤。
8.根据权利要求7所述的1090ESADS-B消息报头的检测方法,其特征在于,设被检测1090ESADS-B消息的一帧长度为Zus;
对检测到的前后两个报头交织的情况进行判定,并对存在交织的两个报头进行取舍的步骤,具体为:
S100,当检测到已经通过验证的有效报头A后,判断在有效报头A的5.0us位置之后的Zus-5.0us内,是否又检测到可通过验证的有效报头B;如果未检测到,则判定出不存在与有效报头A发生交织的情况,结束流程;如果检测到,则执行S101;
S101,判断有效报头B的参考功率次小值是否大于2.5倍的有效报头A的参考功率次大值,如果是,则得出有效报头B为强信号,交织于弱信号的有效报头A,则舍弃有效报头A,保留有效报头B;否则,舍弃有效报头B,保留有效报头A。
9.根据权利要求8所述的1090ESADS-B消息报头的检测方法,其特征在于,Z为64us或120us。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410855818.1A CN105681121B (zh) | 2014-12-31 | 2014-12-31 | 一种1090es ads-b消息报头的检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410855818.1A CN105681121B (zh) | 2014-12-31 | 2014-12-31 | 一种1090es ads-b消息报头的检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105681121A true CN105681121A (zh) | 2016-06-15 |
CN105681121B CN105681121B (zh) | 2018-10-19 |
Family
ID=56958258
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410855818.1A Expired - Fee Related CN105681121B (zh) | 2014-12-31 | 2014-12-31 | 一种1090es ads-b消息报头的检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105681121B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106597122A (zh) * | 2017-01-10 | 2017-04-26 | 烟台大学 | 一种雷达及通信信号的脉冲宽度检测算法 |
CN107102215A (zh) * | 2017-03-17 | 2017-08-29 | 清华大学 | 接收端的前导帧头脉冲的检测方法及接收端 |
CN108134649A (zh) * | 2017-11-29 | 2018-06-08 | 西安空间无线电技术研究所 | 一种星载ads-b信号突发检测方法及系统 |
CN111211866A (zh) * | 2019-12-27 | 2020-05-29 | 北京航空航天大学杭州创新研究院 | 基于ads-b的新型功率一致性验证方法 |
CN111988115A (zh) * | 2020-08-31 | 2020-11-24 | 四川九洲空管科技有限责任公司 | 通用平台上基于并行计算的ads-b分布式处理系统 |
CN112217584A (zh) * | 2020-09-23 | 2021-01-12 | 中国民航大学 | 高动态ads-b电文的低误报率报头检测方法 |
CN113030870A (zh) * | 2021-05-24 | 2021-06-25 | 成都和为时代科技有限公司 | 一种基于时域特征的iff模式5信号盲识别的方法 |
CN116055927A (zh) * | 2023-04-03 | 2023-05-02 | 深圳市紫光同创电子有限公司 | 数据二倍过采样方法、系统、设备及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103354512A (zh) * | 2013-07-09 | 2013-10-16 | 宁波成电泰克电子信息技术发展有限公司 | 基于相关检测技术的s模式ads-b系统的报头检测方法 |
CN103378934A (zh) * | 2012-04-27 | 2013-10-30 | 中国传媒大学 | 一种1090MHz ES ADS-B帧信号解调方法 |
CN103701488A (zh) * | 2013-12-21 | 2014-04-02 | 中电科航空电子有限公司 | 具有1090兆赫扩展电文的广播式自动相关监视能力的s模式应答机 |
EP2833567A2 (en) * | 2012-03-27 | 2015-02-04 | Huawei Technologies Co., Ltd. | Optical fiber testing method, apparatus and passive optical network system |
-
2014
- 2014-12-31 CN CN201410855818.1A patent/CN105681121B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2833567A2 (en) * | 2012-03-27 | 2015-02-04 | Huawei Technologies Co., Ltd. | Optical fiber testing method, apparatus and passive optical network system |
CN103378934A (zh) * | 2012-04-27 | 2013-10-30 | 中国传媒大学 | 一种1090MHz ES ADS-B帧信号解调方法 |
CN103354512A (zh) * | 2013-07-09 | 2013-10-16 | 宁波成电泰克电子信息技术发展有限公司 | 基于相关检测技术的s模式ads-b系统的报头检测方法 |
CN103701488A (zh) * | 2013-12-21 | 2014-04-02 | 中电科航空电子有限公司 | 具有1090兆赫扩展电文的广播式自动相关监视能力的s模式应答机 |
Non-Patent Citations (1)
Title |
---|
王飞: "基于FPGA的ADS-B信号报头检测系统设计", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106597122A (zh) * | 2017-01-10 | 2017-04-26 | 烟台大学 | 一种雷达及通信信号的脉冲宽度检测算法 |
CN106597122B (zh) * | 2017-01-10 | 2019-03-19 | 烟台大学 | 一种雷达及通信信号的脉冲宽度检测算法 |
CN107102215B (zh) * | 2017-03-17 | 2019-09-20 | 清华大学 | 接收端的前导帧头脉冲的检测方法及接收端 |
CN107102215A (zh) * | 2017-03-17 | 2017-08-29 | 清华大学 | 接收端的前导帧头脉冲的检测方法及接收端 |
CN108134649B (zh) * | 2017-11-29 | 2020-08-14 | 西安空间无线电技术研究所 | 一种星载ads-b信号突发检测方法及系统 |
CN108134649A (zh) * | 2017-11-29 | 2018-06-08 | 西安空间无线电技术研究所 | 一种星载ads-b信号突发检测方法及系统 |
CN111211866A (zh) * | 2019-12-27 | 2020-05-29 | 北京航空航天大学杭州创新研究院 | 基于ads-b的新型功率一致性验证方法 |
CN111211866B (zh) * | 2019-12-27 | 2022-06-24 | 北京航空航天大学杭州创新研究院 | 基于ads-b的新型功率一致性验证方法 |
CN111988115A (zh) * | 2020-08-31 | 2020-11-24 | 四川九洲空管科技有限责任公司 | 通用平台上基于并行计算的ads-b分布式处理系统 |
CN111988115B (zh) * | 2020-08-31 | 2023-02-03 | 四川九洲空管科技有限责任公司 | 通用平台上基于并行计算的ads-b分布式处理系统 |
CN112217584A (zh) * | 2020-09-23 | 2021-01-12 | 中国民航大学 | 高动态ads-b电文的低误报率报头检测方法 |
CN112217584B (zh) * | 2020-09-23 | 2022-03-18 | 中国民航大学 | 高动态ads-b电文的低误报率报头检测方法 |
CN113030870A (zh) * | 2021-05-24 | 2021-06-25 | 成都和为时代科技有限公司 | 一种基于时域特征的iff模式5信号盲识别的方法 |
CN116055927A (zh) * | 2023-04-03 | 2023-05-02 | 深圳市紫光同创电子有限公司 | 数据二倍过采样方法、系统、设备及存储介质 |
CN116055927B (zh) * | 2023-04-03 | 2023-06-02 | 深圳市紫光同创电子有限公司 | 数据二倍过采样方法、系统、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN105681121B (zh) | 2018-10-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105681121A (zh) | 一种1090es ads-b消息报头的检测方法 | |
CN102857996B (zh) | 一种小区搜索定时同步的方法 | |
CN101727654B (zh) | 一种并行流水线实现的点目标连通域实时标记与识别方法 | |
CN102831789B (zh) | 一种基于fpga的s模式ads_b系统的纠检错方法 | |
CN102298146A (zh) | 一种二次雷达信号解码方法 | |
CN103199944A (zh) | 广播式自动相关监视信号检测方法及装置 | |
CN106483508A (zh) | 一种用于数字信道化接收机的同时到达信号检测方法 | |
CN104253721B (zh) | 一种ads‑b系统中s模式应答信号的数据位提取方法 | |
CN104064057A (zh) | 图像跟踪测量数据与雷达测量数据互补融合的实现方法 | |
CN106100786A (zh) | 一种基于报头整体相关的星载ads‑b报头检测方法 | |
CN104502905A (zh) | 基于多通道顺序统计的距离扩展目标检测方法 | |
CN105911570A (zh) | 一种卫星通信突发捕获方法 | |
CN105445699A (zh) | 一种非视距误差消除的测距方法及系统 | |
CN109270520B (zh) | 基于幅度信息获取二次雷达应答目标身份代码的处理方法 | |
CN108134649A (zh) | 一种星载ads-b信号突发检测方法及系统 | |
CN104635242A (zh) | 一种基于多源卫星遥感数据的沙尘暴监测方法 | |
CN103354512A (zh) | 基于相关检测技术的s模式ads-b系统的报头检测方法 | |
CN104966102A (zh) | 基于卫星图像的台风检测方法 | |
CN116545470B (zh) | 一种基于pdw的s应答信号解码器解码方法及fpga实现装置 | |
CN104993908A (zh) | 提高s模式应答信号译码可靠性的方法 | |
CN106950527A (zh) | 一种多基线干涉仪测向体制下的脉冲信号分选方法 | |
CN102480455A (zh) | 长期演进系统中主同步信号的检测方法和检测装置 | |
CN111901268A (zh) | 基于频谱重排及Gumbel分布拟合检验的BPSK/QPSK信号调制识别方法 | |
CN105467383A (zh) | 一种tof技术中基于波形匹配的测距方法 | |
CN104616007A (zh) | 一种基于显著性检测和颜色直方图模型的车辆识别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181019 Termination date: 20211231 |
|
CF01 | Termination of patent right due to non-payment of annual fee |