一种测试信号分离性能的星载ADS-B信号形成方法

摘要

本发明公开了一种测试信号分离性能的星载ADS‑B信号形成方法，通过在设置飞行器的数量，卫星天线的数量，综合考虑不同卫星天线接收到的同一个报文信号功率不同，卫星天线的接收方向不同，按照每个卫星天线的接收报文信号时间将W个报文信号叠加，得到测试接收机多通道信号分离性能的星载ADS‑B信号，最终在多个通道输出功率不同的相关ADS‑B报文信号，测试多通道接收机的信号分离的性能。

说明书

技术领域

本发明涉及信号处理领域，具体涉及一种测试信号分离性能的星载ADS-B信号形成方法。

背景技术

星载ADS-B(广播式自动相关监视)系统中，接收机在卫星上接收飞行器发送的ADS-B报文信号。由于卫星的高度优势，星载ADS-B接收机能够接收到更多飞行器发送的报文信号，因此在单位时间内，信号密度更大，信号发生交织的概率更高，信号交叠成为一个很严重的问题。在这种情况下，星载接收机必须具有良好的解信号分离能力。多通道分离算法是利用多天线接收机对ADS-B报文信号进行接收，是一种具有很大优势的信号分离算法。在现有的ADS-B信号源中，比如KLJ Instruments公司的一款ADS-B信号源SQTR-2M，ADSBator信号发生器，国内的三航公司的ADS-B和A/C信号模拟发生器等，都无法测试这类采用多通道分离算法解信号交织的接收机性能。

发明内容

针对现有技术中的上述不足，本发明提供了一种测试信号分离性能的星载ADS-B信号形成方法解决了现有技术无法测试采用多通道分离算法解信号交织的接收机性能的问题。

为了达到上述发明目的，本发明采用的技术方案为：一种测试信号分离性能的星载ADS-B信号形成方法，包括以下步骤：

S1、预设卫星的天线个数N、卫星的初始位置信息和速度信息；

S2、预设飞行器的个数M、初始位置信息、速度信息、ID和飞行器的发射功率为T_power_j，其中，1≤j≤M；

S3、根据DO-260标准中的发报规则，确定每架飞行器的发报时刻点和报文类型，得到报文信号总数W；

S4、根据W个报文信号与M个飞行器的对应关系以及飞行器的发射功率T_power_j，得到报文信号的发射功率power_k，其中，1≤k≤W；

S5、根据卫星的初始位置信息、M架飞行器的初始位置信息和卫星及飞行器的运动方程，得到M个飞行器和卫星在W个发报时刻点的发报飞行器和卫星的位置信息；

S6、根据W个发报时刻点的发报飞行器的位置信息、速度信息和ID，得到PPM基带信号F(τ)_k，0≤τ≤T_a，T_a为信号长度；

S7、将PPM基带信号F(τ)_k分别加载在频率为f的射频载波信号上，得到一般ADS-B射频信号S_a(k)(τ)；

S8、根据报文信号的发射功率power_k和ADS-B射频信号S_a(k)(τ)，得到W个报文信号到达卫星的N个天线的星载ADS-B信号0≤i≤N；

S9、根据W个发报时刻点的卫星和发报飞行器的位置信息，得到W个报文信号被N个卫星天线接收的接收时间；

S10、根据W个报文信号被N个卫星天线接收的接收时间，预设W个报文信号的接收时间与第一个报文信号的接收时间差为Δτ_k；

S11、按照接收时间差将每个卫星天线接收到的W个星载ADS-B信号叠加，得到测试信号分离性能的星载ADS-B信号S_i(τ)。

进一步地：所述步骤S5中的运动方程表达式为：

其中，t₀为卫星或飞行器的初始时刻，v_l为卫星或飞行器的水平速度，t为卫星或飞行器的飞行时间，v_l(t)为t时刻卫星或飞行器的水平飞行速度，v_c为卫星或飞行器的垂直飞行速度，r为地球半径，θ为卫星或飞行器的航向角，h(t)为t时刻卫星或飞行器的飞行高度，lat(t)为t时刻卫星或飞行器的纬度，lon(t)为t时刻卫星或飞行器的经度。

进一步地：所述步骤S6中得到PPM基带信号F(τ)_k包括以下步骤：

S61、对W个发报时刻点的发报飞行器的位置信息、速度信息和ID进行编码，得到W个报文信号的位置编码、速度编码和ID编码；

S62、对W个报文信号的位置编码、速度编码和ID编码进行CRC校验，得到W个ADS-B报文数据链比特位；

S63、根据W个ADS-B报文数据链比特位，得到PPM基带信号F(τ)_k。

进一步地：所述步骤S7中ADS-B射频信号S_a(k)(τ)的表达式为：

其中，为初始相位，

进一步地：所述步骤S8得到W个报文信号到达N个卫星天线的星载ADS-B信号包括以下步骤：

S81、根据报文信号的发射功率power_k、机载发射天线增益gain1_k、星载接收天线增益gain2_ik和路径损耗loss_k，得到W个报文信号到达N个卫星天线时的信号功率rec_power_ik；

S82、根据W个报文信号达到卫星的信号功率rec_power_ik和ADS-B射频信号S_a(k)(τ)，得到W个报文信号达到N个卫星天线的星载ADS-B信号

进一步地：所述步骤S81具体为：

A1、根据W个发报时刻点卫星与发报飞行器的距离，得到报文信号传输距离为d_k；

A2、根据传播频率f和报文信号传输距离d_k，得到报文信息的路径损耗loss_k：

loss_k＝32.44+20logf+20logd_k

A3、通过报文信号传播方向与飞行器的天线中心轴的夹角在机载天线功率归一化方向图查找得到径向长度P_k；

A4、根据机载天线的最大天线增益gain1_max和径向长度P_k，得到机载发射天线增益gain1_k：

gain1_k＝gain1_max+10log(P_k)

A5、根据N个卫星天线与卫星的天线中心轴的偏移角度，得到卫星天线的中心点A_i；

A6、建立以地球中心O₁为坐标原点的地心坐标系，建立以卫星中心O₂为坐标原点的站心坐标系，得到中心点A_i在站心坐标系的坐标；

A7、根据发报飞行器和卫星在W个发报时刻点的位置信息，得到W个发报时刻点的发报飞行器和卫星在地心坐标系中的坐标；

A8、根据中心点A_i的坐标，得到中心点A_i在地心坐标系中的坐标A_a(i)；

A9、根据中心点A_a(i)、卫星和发报飞行器在地心坐标系中的坐标，得到W个报文信号传播方向与卫星天线中心的夹角

A10、根据报文信号传播方向与卫星天线中心的夹角在星载天线场强归一化方向图查找得到径向长度U_ik；

A11、根据星载天线的最大天线增益gain2_max和径向长度U_ik，得到星载接收天线增益gain2_ik：

gain2_ik＝gain2_max+10log(U_ik²)

A12、根据报文信号的发射功率power_k、机载发射天线增益gain1_k、星载接收天线增益gain2_ik和路径损耗loss_k，得到W个报文信号达到N个卫星天线的信号功率rec_power_ik的计算公式为：

rec_power_ik＝power_k+gain1_k+gain2_ik+loss_k。

进一步地：所述步骤S9中得到W个报文信号被卫星天线接收的接收时间包括以下步骤：

S91、根据发报时刻卫星和发报飞行器的位置信息，得到卫星与发报飞行器的距离；

S92、根据卫星与发报飞行器的距离，得到W个报文信号的传输延迟；

S93、根据W个报文信号的传输延迟，得到W个报文信号被N个卫星天线接收的接收时间。

进一步地：所述步骤S11中测试信号分离性能的星载ADS-B信号S_i(τ)的表达式为：

本发明的有益效果为：一种测试信号分离性能的星载ADS-B信号形成方法，通过在设置飞行器的数量，卫星天线的数量，综合考虑不同卫星天线接收到的同一个报文信号功率不同，卫星天线的接收方向不同，按照每个卫星天线的接收报文信号时间将W个报文信号叠加，得到测试接收机多通道信号分离性能的星载ADS-B信号，最终在多个通道输出功率不同的相关ADS-B报文信号，测试多通道接收机的信号分离的性能。

附图说明

图1为一种测试信号分离性能的星载ADS-B信号形成方法流程图；

图2为机载天线功率归一化方向图；

图3为星载天线场强归一化方向图。

具体实施方式

下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。

如图1所示，一种测试信号分离性能的星载ADS-B信号形成方法，包括以下步骤：

S1、预设卫星的天线个数N、卫星的初始位置信息和速度信息；

S2、预设飞行器的个数M、初始位置信息、速度信息、ID和飞行器的发射功率为T_power_j，其中，1≤j≤M；

S3、根据DO-260标准中的发报规则，确定每架飞行器的发报时刻点和报文类型，得到报文信号总数W；

S4、根据W个报文信号与M个飞行器的对应关系以及飞行器的发射功率T_power_j，得到报文信号的发射功率power_k，其中，1≤k≤W；

S5、根据卫星的初始位置信息、M架飞行器的初始位置信息和卫星及飞行器的运动方程，得到M个飞行器和卫星在W个发报时刻点的发报飞行器和卫星的位置信息；

步骤S5中的运动方程表达式为：

其中，t₀为卫星或飞行器的初始时刻，v_l为卫星或飞行器的水平速度，t为卫星或飞行器的飞行时间，v_l(t)为t时刻卫星或飞行器的水平飞行速度，v_c为卫星或飞行器的垂直飞行速度，r为地球半径，θ为卫星或飞行器的航向角，h(t)为t时刻卫星或飞行器的飞行高度，lat(t)为t时刻卫星或飞行器的纬度，lon(t)为t时刻卫星或飞行器的经度。

S6、根据W个发报时刻点的发报飞行器的位置信息、速度信息和ID，得到PPM基带信号F(τ)_k，0≤τ≤T_a，T_a为信号长度；

步骤S6中得到PPM基带信号F(τ)_k包括以下步骤：

S61、对W个发报时刻点的发报飞行器的位置信息、速度信息和ID进行编码，得到W个报文信号的位置编码、速度编码和ID编码；

S62、对W个报文信号的位置编码、速度编码和ID编码进行CRC校验，得到W个ADS-B报文数据链比特位；

S63、根据W个ADS-B报文数据链比特位，得到PPM基带信号F(τ)_k。

S7、将PPM基带信号F(τ)_k分别加载在频率为f的射频载波信号上，得到一般ADS-B射频信号S_a(k)(τ)；

步骤S7中ADS-B射频信号S_a(k)(τ)的表达式为：

其中，为初始相位，

S8、根据报文信号的发射功率power_k和ADS-B射频信号S_a(k)(τ)，得到W个报文信号到达卫星的N个天线的星载ADS-B信号0≤i≤N；

步骤S8中得到W个报文信号到达N个卫星天线的星载ADS-B信号包括以下步骤：

S81、根据报文信号的发射功率power_k、机载发射天线增益gain1_k、星载接收天线增益gain2_ik和路径损耗loss_k，得到W个报文信号到达N个卫星天线时的信号功率rec_power_ik；

步骤S81具体为：

A1、根据W个发报时刻点卫星与发报飞行器的距离，得到报文信号传输距离为d_k；

A2、根据传播频率f和报文信号传输距离d_k，得到报文信息的路径损耗loss_k：

loss_k＝32.44+20logf+20logd_k>

A3、通过报文信号传播方向与飞行器的天线中心轴的夹角在机载天线功率归一化方向图查找得到径向长度P_k，如图2所示；

A4、根据机载天线的最大天线增益gain1_max和径向长度P_k，得到机载发射天线增益gain1_k：

gain1_k＝gain1_max+10log(P_k)>

A5、根据N个卫星天线与卫星的天线中心轴的偏移角度，得到卫星天线的中心点A_i；

A6、建立以地球中心O₁为坐标原点的地心坐标系，建立以卫星中心O₂为坐标原点的站心坐标系，得到中心点A_i在站心坐标系的坐标；

A7、根据发报飞行器和卫星在W个发报时刻点的位置信息，得到W个发报时刻点的发报飞行器和卫星在地心坐标系中的坐标；

A8、根据中心点A_i的坐标，得到中心点A_i在地心坐标系中的坐标A_a(i)；

A9、根据中心点A_a(i)、卫星和发报飞行器在地心坐标系中的坐标，得到W个报文信号传播方向与卫星天线中心的夹角

A10、根据报文信号传播方向与卫星天线中心的夹角在星载天线场强归一化方向图查找得到径向长度U_ik，如图3所示；

A11、根据星载天线的最大天线增益gain2_max和径向长度U_ik，得到星载接收天线增益gain2_ik：

gain2_ik＝gain2_max+10log(U_ik²)>

A12、根据报文信号的发射功率power_k、机载发射天线增益gain1_k、星载接收天线增益gain2_ik和路径损耗loss_k，得到W个报文信号达到N个卫星天线的信号功率rec_power_ik的计算公式为：

rec_power_ik＝power_k+gain1_k+gain2_ik+loss_k>

S82、根据W个报文信号达到卫星的信号功率rec_power_ik和ADS-B射频信号S_a(k)(τ)，得到W个报文信号达到N个卫星天线的星载ADS-B信号

S9、根据W个发报时刻点的卫星和发报飞行器的位置信息，得到W个报文信号被N个卫星天线接收的接收时间；

步骤S9中得到W个报文信号被卫星天线接收的接收时间包括以下步骤：

S91、根据发报时刻卫星和发报飞行器的位置信息，得到卫星与发报飞行器的距离；

S92、根据卫星与发报飞行器的距离，得到W个报文信号的传输延迟；

S93、根据W个报文信号的传输延迟，得到W个报文信号被N个卫星天线接收的接收时间。

S10、根据W个报文信号被N个卫星天线接收的接收时间，预设W个报文信号的接收时间与第一个报文信号的接收时间差为Δτ_k；

S11、按照接收时间差将每个卫星天线接收到的W个星载ADS-B信号叠加，得到测试信号分离性能的星载ADS-B信号S_i(τ)。

步骤S11中测试信号分离性能的星载ADS-B信号S_i(τ)的表达式为：

本发明的有益效果为：一种测试信号分离性能的星载ADS-B信号形成方法，通过在设置飞行器的数量，卫星天线的数量，综合考虑不同卫星天线接收到的同一个报文信号功率不同，卫星天线的接收方向不同，按照每个卫星天线的接收报文信号时间将W个报文信号叠加，得到测试接收机多通道信号分离性能的星载ADS-B信号，最终在多个通道输出功率不同的相关ADS-B报文信号，测试多通道接收机的信号分离的性能。