CN104243070B - 面向深空通信的无线信道实时模拟器及其模拟方法 - Google Patents
面向深空通信的无线信道实时模拟器及其模拟方法 Download PDFInfo
- Publication number
- CN104243070B CN104243070B CN201410519736.XA CN201410519736A CN104243070B CN 104243070 B CN104243070 B CN 104243070B CN 201410519736 A CN201410519736 A CN 201410519736A CN 104243070 B CN104243070 B CN 104243070B
- Authority
- CN
- China
- Prior art keywords
- formula
- deep space
- unit
- noise
- channel
- 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.)
- Active
Links
- 238000004891 communication Methods 0.000 title claims abstract description 69
- 238000000034 method Methods 0.000 title claims abstract description 22
- 230000033001 locomotion Effects 0.000 claims description 18
- 238000005562 fading Methods 0.000 claims description 14
- 241000209094 Oryza Species 0.000 claims description 12
- 235000007164 Oryza sativa Nutrition 0.000 claims description 12
- 235000009566 rice Nutrition 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 9
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims description 6
- 239000002245 particle Substances 0.000 claims description 5
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- MYMOFIZGZYHOMD-UHFFFAOYSA-N Dioxygen Chemical compound O=O MYMOFIZGZYHOMD-UHFFFAOYSA-N 0.000 claims description 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 3
- 238000000205 computational method Methods 0.000 claims description 3
- 229910052760 oxygen Inorganic materials 0.000 claims description 3
- 239000001301 oxygen Substances 0.000 claims description 3
- 230000001902 propagating effect Effects 0.000 claims description 3
- 230000005855 radiation Effects 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 239000000523 sample Substances 0.000 claims description 3
- 230000009897 systematic effect Effects 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 230000000644 propagated effect Effects 0.000 claims description 2
- 238000004088 simulation Methods 0.000 abstract description 16
- 230000005540 biological transmission Effects 0.000 abstract description 7
- 230000008054 signal transmission Effects 0.000 abstract description 5
- 238000012360 testing method Methods 0.000 abstract description 5
- 238000002474 experimental method Methods 0.000 abstract description 4
- 238000011156 evaluation Methods 0.000 abstract description 3
- 230000007423 decrease Effects 0.000 description 12
- 238000006243 chemical reaction Methods 0.000 description 7
- 239000000203 mixture Substances 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Landscapes
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
一种面向深空通信的无线信道实时模拟器及其模拟方法,支持上行链路和下行链路实时双工通信信道模拟,包括主控单元、基带信道模拟单元、变频器单元以及衰减器单元,所述主控单元通过CPCI总线与基带信道模拟单元、变频器单元以及衰减器单元进行数据及命令传输;所述基带信道模拟单元、变频器单元和衰减器单元间通过SMA接口进行模拟信号传输,并结合其模拟方法可有效实时模拟深空背景下无线信号的传播情况,从而在地面实验室实现深空通信设备的调试、测试及性能评估。
Description
技术领域
本发明属于面向深空通信技术领域,具体涉及一种面向深空通信的无线信道实时模拟器及其模拟方法。
背景技术
当前,国内外航天事业蓬勃发展,发送了大量的航天器对地外空间及地外天体开展科学研究。在这些航天工程中,确保航天器与地面站、航天器与航天器之间的可靠通信是航天任务成败与否的关键。深空通信环境完全不同于地面通信环境,比较复杂,故对航天工程中深空背景下的无线信道研究,特别是深空信道建模及实验室模拟对深空通信系统的研制及测试,乃至航天事业的发展意义重大。
深空通信距离遥远,导致延时长,传播损耗大,同时还受到太阳活动(如太阳风、日冕)、月球及其它天体噪声等影响,使得接收信号功率极其微弱,且信噪比极低。当太阳处在两个通信航天器之间时(即上合时期),无线传播信号将受到太阳风和日冕的影响,这种不规则影响会造成信号幅度的闪烁,使得深空信道具有时变衰落特性。同时,深空信道还具有窗口特性,上下行链路不对称,上行链路主要用于传输遥控指令,码速率较低;下行链路主要用于传输图像及其它探测到的数据信息,码速率较高。此外,深空通信频段高、带宽大,加上航天器之间的相对运动速度较大,将导致大多普勒频移和扩展。
对深空信道的研究包括通信实测、软件仿真和硬件模拟等模拟方法。然而,真实的深空通信链路测试需要借助卫星或者航天器进行,会消耗巨大的人力和物力;软件仿真由于运算量巨大,且无法提供真实信号的输入输出接口而颇受局限。最有效的模拟方法是采用硬件方式实现深空信道的地面实验室可重复模拟,DSP及FPGA强大的运算能力,以及硬件板卡丰富的数模/模数转换接口,使得实验室环境下实时模拟深空信道成为可能。
另外深空信道有别于与地面无线移动信道,深空信号传输距离极其遥远,必然导致巨大的路径损耗;航天器和地球的高速运动造成收发信机巨大的相对速度,从而引入巨大的多普勒频移;深空链路受到太阳活动抛射的非均匀带电粒子的干扰,将引发多径效应,造成信号幅度的闪烁(也称太阳闪烁);深空通信的噪声干扰复杂,包括宇宙背景噪声、太阳噪声、月球噪声、行星噪声、地球站噪声和天线噪声等。
发明内容
本发明的目的提供一种面向深空通信的无线信道实时模拟器及其模拟方法,支持上行链路和下行链路实时双工通信信道模拟,包括主控单元、基带信道模拟单元、变频器单元以及衰减器单元,所述主控单元通过CPCI总线与基带信道模拟单元、变频器单元以及衰减器单元进行数据及命令传输;所述基带信道模拟单元、变频器单元和衰减器单元间通过SMA接口进行模拟信号传输,并结合其模拟方法可有效实时模拟深空背景下无线信号的传播情况,从而在地面实验室实现深空通信设备的调试、测试及性能评估。
为了克服现有技术中的不足,本发明提供了一种面向深空通信的无线信道实时模拟器及其模拟方法的解决方案,具体如下:
一种面向深空通信的无线信道实时模拟器,支持上行链路和下行链路实时双工通信信道模拟,包括主控单元1-0、基带信道模拟单元2-0、变频器单元3-0以及衰减器单元4-0,所述主控单元1-0通过CPCI总线与基带信道模拟单元2-0、变频器单元3-0以及衰减器单元4-0进行数据及命令传输;所述基带信道模拟单元2-0、变频器单元3-0和衰减器单元4-0间通过SMA接口进行模拟信号传输:
所述上行链路包括主控单元中的深空通信场景输入模块1-1,上行链路信道参数估计模块1-2,基带信道模拟单元中的上行链路多径衰落模拟模块2-1,变频器单元中的下变频器模块3-1和上变频器模块3-2,衰减器单元中的衰减器模块4-1,上行链路发射信号,通过上行链路输入5-1端口经下变频器模块3-1下变频为中频信号,该中频信号输入到上行链路多径衰落模拟模块2-1叠加多径衰落信道,之后经上变频器模块3-2上变频到深空通信频段,输出射频信号经衰减器单元4-1衰减后从上行链路输出5-2端口输出的信号即为模拟经过深空链路后的上行链路接收端信号,所述的下行链路的组成和实现方法与所述的上行链路一致。
所述主控单元1-0为带有CPU计算单元的主板,并连接有与之兼容的键盘、鼠标、显示器这样的输入输出设备,带有CPU计算单元的主板还包括主控软件,主控软件实现的功能模块包括深空通信场景输入模块1-1,上行链路信道参数估计模块1-2和下行链路信道参数估计模块1-3,所述的深空通信场景输入模块1-1接受用户输入的通信场景参数,通信场景参数包括深空探测器运动轨迹、通信频段、星体半径、电导率以及介电常数,所述的上行链路信道参数估计模块1-2和下行链路信道参数估计模块1-3从该场景参数中提取出信道参数,包括多径时延、多普勒频移、多普勒扩展、莱斯因子以及噪声功率/信噪比,并将该信道参数经CPCI总线配置给基带信道模拟单元2-0、变频器单元3-0及衰减器单元4-0。
所述基带信道模拟单元2-0为由DSP和FPGA组成的信号处理板卡,包含模数转换模块和数模转换模块,实现上行链路多径衰落模拟和下行链路多径衰落模拟的模块分别为上行链路多径衰落模拟模块2-1和下行链路多径衰落模拟模块2-2,基带信道模拟单元2-0接收主控单元1-0经CPCI总线配置的信道参数,包括多径时延、多普勒频移、多普勒扩展、莱斯因子以及噪声功率/信噪比,基于该硬件信道参数实时产生与之相符的信道衰落,并叠加至输入信号后输出。
所述变频器单元3-0为由混频器和滤波器组成的信号变频单元板卡,包括两个下变频器模块,两个下变频器模块分别为第一下变频器模块3-1和第二下变频器模块3-4,两个上变频器模块,两个上变频器模块分别为第一上变频器模块3-2和第二上变频器模块3-3,变频器单元3-0接收主控单元1-0经CPCI总线配置的工作频段数据,频段包含深空通信的常用频段,VHF频段、L频段、S频段、C频段、X频段、Ku频段以及Ka频段,其中第一下变频器模块3-1和第二下变频器模块3-4的功能为将输入射频深空信号下变频为中频信号,第一上变频器模块3-2和第二上变频器模块3-3的功能为将叠加信道后的中频信号上变频到深空通信所用的射频段。
所述衰减器单元4-0为由程控衰减电路构成的信号衰减单元板卡,内含两个衰减器模块,两个衰减器模块分别为第一衰减器模块4-1和第二衰减器模块4-2,第一衰减器模块4-1和第二衰减器模块4-2都分别由一个10dB每步进十五大步进衰减器和一个1dB每步进的小步进衰减器组成,总衰减量程为大步进衰减器量程和小步进衰减器量程之和,衰减步进与小步进衰减器一致,衰减器单元4-0接收主控单元1-0经CPCI总线配置的衰减值,实时控制深空通信信号的总衰减量。
所述的面向深空通信的无线信道实时模拟器的模拟方法,步骤如下:
步骤1:用户在主控单元1-0的主控软件中输入深空通信场景参数,包括深空探测器飞行轨迹文件、通信频段、星体半径、电导率以及介电常数,主控软件基于收发端的运动轨迹计算传播场景的几何参量,包括收发端运动速度、直线距离和反射点坐标;
步骤2:通过传播场景的几何参量及系统参数计算得到不同时刻的信道状态参数,包括深空链路传播损耗L1(t),L2(t)、多普勒频移fd,1(t),fd,2(t)、各径时延τ1(t),τ2(t)、莱斯因子K和噪声功率PN这样的信道参数,步骤2的具体方法如下:
1)深空传播链路损耗计算
所述的深空传播链路损耗计算为传播链路损耗LLoss考虑自由空间传播损耗Lf,大气传播损耗La,降雨损耗Lr,收发天线增益Gr、Gt这样的因素,计算方法如式(1)所示:
LLoss=Lf+La+Lr-Gr-Gt(dB) (1)
其中如式(2)所示:
Lf=92.44+20lg(d/103)+20lg(f/109)(dB) (2)
式中,f为频率,单位Hz;d为通信距离,单位m;
如式(3)所示:
式中,r0和rw分别为氧气分子和水蒸气分子对电波的衰减率;hs为地面站海拔高度;θ为地面站天线波束仰角;g(h0)、g(hw)是与氧气层等效厚度h0、水蒸气层等效厚度hw相关的函数;
如式(4)所示:
Lr=γrLE(dB) (4)式中,γr为降雨损耗率;LE为等效路径长度;
2)多普勒频移计算
所述的多普勒频移计算为凭借深空通信中,通信双方本身的运动速度很大,巨大的径向运动速度以及超高的通信频率产生了多普勒效应,这样直射路径的最大多普勒频移能够表示为如式(5)所示:
式中,分别为地面站和航天器的移动速度矢量,θB,θA则表示速度方向和直射路径的夹角,f0为信号频率,c为光速;
3)莱斯因子计算
所述的莱斯因子计算根据莱斯因子K与太阳活动造成的太阳闪烁指数具有如式(6)所示的关系而得:
其中,太阳闪烁指数m与太阳-地球-航天器夹角相关,当该夹角大于4°时m可忽略;当该夹角小于4°时,m值随着夹角的减小而增大;夹角趋近0°时m为1,此时K→∞,该情况下无直射路径;
4)信道噪声计算
所述的信道噪声计算是根据式(7)所示的深空信道噪声用等效噪声温度来计算,
PN=k(TRF+TC+TSUN+TMOON+TPL+TG+TAN)BN (7)
式中,PN为噪声功率;k为玻尔兹曼常量并且k=1.3806505×10-23J/K·Hz;BN为接收机等效噪声带宽;TRF为低噪声放大器的噪声温度,Tc为宇宙背景噪声温度,取2.66±0.77K;TG为地球站噪声,取值为0.1K;TAN为天线噪声温度;TSUN为太阳噪声温度,如(8)所示为
式中,S为单位太阳粒子辐射流密度;A为天线有效孔径面积;η为天线效率;r为考虑深空站系统接收机非线性因素后的减少因子;
TMOON为月球噪声温度,S波段,X波段和Ka波段噪声温度分别由式(9)、式(10)和式(11)求得,
TMOON(φ)=136±1.4(K)(9)
式中,φ为天线主波束宽度与月球圆盘中心之间的偏移角度;
TPL为行星噪声温度,如式(12)所示:
式中,Tk为行星黑体温度;d为行星直径;R为与行星的距离;θ为行星中心到天线波束中心的夹角;θ0为天线半功率波束宽度;G为天线增益;
步骤3:根据主控单元1-0经CPCI总线配置的信道状态参数,基带信道模拟单元2-0的DSP首先将其转为硬件参数并进行定点化,实时配置给对应的FPGA工作模块,步骤3的方法具体如下:
1)计算各径的相对时延τ=τ1(t)-τ2(t),其中τ1(t)、τ2(t)分别为直射径和散射径时延,将时延τ转化为定点的FPGA硬件参数,并实时配置给时延模块;
2)计算各径的相对损耗因子α1,α2,α1,α2应当满足(α2/α1)2=K的条件,将其转化为定点的FPGA硬件参数,并实时配置给损耗模块;
3)计算各径的多普勒频移fd,1(t),fd,2(t),将其转化为定点的FPGA硬件参数,并实时配置给多普勒频移模块;
4)计算信道衰落的定点化硬件参数,并实时配置给信道衰落模块;
5)计算深空信道噪声功率PN的定点化参数,并实时配置给高斯噪声模块;
步骤4:根据主控单元1-0经CPCI总线配置的通信频段,变频器单元设置上变频器模块及下变频器模块的工作频段;
步骤5:根据主控单元1-0经CPCI总线配置的衰减参数,衰减器单元实时设置衰减值。
所述的步骤1中的运动轨迹能够进行内插,内插倍数用户能够设置。
由这些技术特征,本发明可有效实时模拟深空背景下无线信号的传播情况,从而在地面实验室实现深空通信设备的调试、测试及性能评估。
附图说明
图1为本发明系统结构示意图。
图2为本发明深空基带信道硬件模拟流程图。
具体实施方式
本发明的目的是研制一种面向深空通信的无线信道实时模拟器及其模拟方法。本发明综合考虑上述深空通信链路的特点,将深空信道建模为
式中,为复基带输入信号;为经过深空环境传播后的复基带接收信号;fd为最大多普勒频率,e为自然对数底数;L(t)表示传播链路损耗,与通信距离和频率有关;为信道干扰噪声。为太阳闪烁导致的闪烁系数,与太阳的闪烁因子有关,且可表示为非零均值的复高斯随机过程,即
式中,是由带电粒子散射导致非视距传输的复高斯随机过程,ρsc(t)与视距接收信号功率成正比。的包络的统计特性服从莱斯分布,其概率密度函数为
式中,为莱斯因子,是零阶第一类修正贝塞尔(Bessel)函数。
下面结合附图和实施例对本发明做进一步说明:
如图1所示,面向深空通信的无线信道实时模拟器,支持上行链路和下行链路实时双工通信信道模拟,包括主控单元1-0、基带信道模拟单元2-0、变频器单元3-0以及衰减器单元4-0,所述主控单元1-0通过CPCI总线与基带信道模拟单元2-0、变频器单元3-0以及衰减器单元4-0进行数据及命令传输;所述基带信道模拟单元2-0、变频器单元3-0和衰减器单元4-0间通过SMA接口进行模拟信号传输;所述上行链路包括主控单元中的深空通信场景输入模块1-1,上行链路信道参数估计模块1-2,基带信道模拟单元中的上行链路多径衰落模拟模块2-1,变频器单元中的下变频器3-1和上变频器3-2,衰减器单元中的衰减器模块4-1,上行链路发射信号,通过上行链路输入5-1端口经下变频器3-1下变频为中频信号,该中频信号输入到上行链路多径衰落模拟模块2-1叠加多径衰落信道,之后经上变频器3-2上变频到深空通信频段,输出射频信号经衰减器模块4-1衰减后从上行链路输出5-2端口输出的信号即为模拟经过深空链路后的上行链路接收端信号,所述的下行链路的组成和实现方法与所述的上行链路一致。所述主控单元1-0为带有CPU计算单元的主板,并连接有与之兼容的键盘、鼠标、显示器这样的输入输出设备,带有CPU计算单元的主板还包括主控软件,主控软件实现的功能模块包括深空通信场景输入模块1-1,上行链路信道参数估计模块1-2和下行链路信道参数估计模块1-3,所述的深空通信场景输入模块1-1接受用户输入的通信场景参数,通信场景参数包括深空探测器运动轨迹、通信频段、星体半径、电导率以及介电常数,所述的上行链路信道参数估计模块1-2和下行链路信道参数估计模块1-3从该场景参数中提取出信道参数,包括多径时延、多普勒频移、多普勒扩展、莱斯因子以及噪声功率/信噪比,并将该信道参数经CPCI总线配置给基带信道模拟单元2-0、变频器单元3-0及衰减器单元4-0。所述基带信道模拟单元2-0为由DSP和FPGA组成的信号处理板卡,包含模数转换模块和数模转换模块,实现上行链路多径衰落模拟和下行链路多径衰落模拟的模块分别为上行链路多径衰落模拟模块2-1和下行链路多径衰落模拟模块2-2,基带信道模拟单元2-0接收主控单元1-0经CPCI总线配置的硬件信道参数,包括多径时延、多普勒频移、多普勒扩展、莱斯因子以及噪声功率/信噪比,基于该硬件信道参数实时产生与之相符的信道衰落,并叠加至输入信号后输出。所述变频器单元3-0为由混频器和滤波器组成的信号变频单元板卡,包括两个下变频器模块,两个下变频器模块分别为第一下变频器模块3-1和第二下变频器模块3-4,两个上变频器模块,两个上变频器模块分别为第一上变频器模块3-2和第二上变频器模块3-3,变频器单元3-0接收主控单元1-0经CPCI总线配置的工作频段数据,频段包含深空通信的常用频段,VHF频段、L频段、S频段、C频段、X频段、Ku频段以及Ka频段,其中第一下变频器模块3-1和第二下变频器模块3-4的功能为将输入射频深空信号下变频为中频信号,第一上变频器模块3-2和第二上变频器模块3-3的功能为将叠加信道后的中频信号上变频到深空通信所用的射频段。所述衰减器单元4-0为由程控衰减电路构成的信号衰减单元板卡,内含两个衰减器模块,两个衰减器模块分别为第一衰减器模块4-1和第二衰减器模块4-2,第一衰减器模块4-1和第二衰减器模块4-2都分别由一个10dB每步进的大步进衰减器和一个1dB每步进的小步进衰减器组成,总衰减量程为大步进衰减器量程和小步进衰减器量程之和,衰减步进与小步进衰减器一致,衰减器单元4-0接收主控单元1-0经CPCI总线配置的衰减值,实时控制深空通信信号的总衰减量。
所述的面向深空通信的无线信道实时模拟器的模拟方法,步骤如下:
步骤1:用户在主控单元1-0的主控软件中输入深空通信场景参数,包括深空探测器飞行轨迹文件、通信频段、星体半径、电导率以及介电常数,主控软件基于收发端的运动轨迹计算传播场景的几何参量,包括收发端运动速度、直线距离和反射点坐标;
步骤2:通过传播场景的几何参量及系统参数计算得到不同时刻的信道状态参数,包括深空链路传播损耗L1(t),L2(t)、多普勒频移fd,1(t),fd,2(t)、各径时延τ1(t),τ2(t)、莱斯因子K和噪声功率PN这样的信道参数,步骤2的具体方法如下:
1)深空传播链路损耗计算
所述的深空传播链路损耗计算为传播链路损耗LLoss考虑自由空间传播损耗Lf,大气传播损耗La,降雨损耗Lr,收发天线增益Gr、Gt这样的因素,计算方法如式(4)所示:
LLoss=Lf+La+Lr-Gr-Gt(dB) (4)
其中如式(5)所示:
Lf=92.44+20lg(d/103)+20lg(f/109)(dB) (5)
式中,f为频率,单位Hz;d为通信距离,单位m;
如式(6)所示:
式中,r0和rw分别为氧气分子和水蒸气分子对电波的衰减率;hs为地面站海拔高度;θ为地面站天线波束仰角;g(h0)、g(hw)是与氧气层等效厚度h0、水蒸气层等效厚度hw相关的函数;
如式(7)所示:
Lr=γrLE(dB) (7)
式中,γr为降雨损耗率;LE为等效路径长度;
2)多普勒频移计算
所述的多普勒频移计算为凭借深空通信中,通信双方本身的运动速度很大,巨大的径向运动速度以及超高的通信频率产生了多普勒效应,这样直射路径的最大多普勒频移能够表示为如式(8)所示:
式中,分别为地面站和航天器的移动速度矢量,θB,θA则表示速度方向和直射路径的夹角,f0为信号频率,c为光速;
3)莱斯因子计算
所述的莱斯因子计算根据莱斯因子K与太阳活动造成的太阳闪烁指数具有如式(9)所示的关系而得:
其中,太阳闪烁指数m与太阳-地球-航天器夹角相关,当该夹角大于4°时m可忽略;当该夹角小于4°时,m值随着夹角的减小而增大;夹角趋近0°时m为1,此时K→∞,该情况下无直射路径。
4)信道噪声计算
所述的信道噪声计算是根据式(10)所示的深空信道噪声用等效噪声温度来计算,
PN=k(TRF+TC+TSUN+TMOON+TPL+TG+TAN)BN (10)
式中,PN为噪声功率;k为玻尔兹曼常量并且k=1.3806505×10-23J/K·Hz;BN为接收机等效噪声带宽;TRF为低噪声放大器的噪声温度,Tc为宇宙背景噪声温度(本实施案例取2.66±0.77K);TG为地球站噪声,本实施例TG取值为0.1K;TAN为天线噪声温度;TSUN为太阳噪声温度,如(11)所示为
式中,S为单位太阳粒子辐射流密度;A为天线有效孔径面积;η为天线效率;r为考虑深空站系统接收机非线性因素后的减少因子;
TMOON为月球噪声温度,S波段,X波段和Ka波段噪声温度分别由式(12)、式(13)和式(14)求得,
TMOON(φ)=136±1.4(K) (12)
式中,φ为天线主波束宽度与月球圆盘中心之间的偏移角度;
TPL为行星噪声温度,如式(15)所示:
式中,Tk为行星黑体温度;d为行星直径;R为与行星的距离;θ为行星中心到天线波束中心的夹角;θ0为天线半功率波束宽度;G为天线增益;
步骤3:根据主控单元1-0经CPCI总线配置的信道状态参数,基带信道模拟单元2-0的DSP首先将其转为硬件参数并进行定点化,实时配置给对应的FPGA工作模块(如图2所示),步骤3的方法具体如下:
1)计算各径的相对时延τ=τ1(t)-τ2(t),其中τ1(t)、τ2(t)分别为直射径和散射径时延,将时延τ转化为定点的FPGA硬件参数,并实时配置给时延模块;
2)计算各径的相对损耗因子α1,α2,α1,α2应当满足(α2/α1)2=K的条件,将其转化为定点的FPGA硬件参数,并实时配置给损耗模块;
3)计算各径的多普勒频移fd,1(t),fd,2(t),将其转化为定点的FPGA硬件参数,并实时配置给多普勒频移模块;
4)计算信道衰落的定点化硬件参数,并实时配置给信道衰落模块;
5)计算深空信道噪声功率PN的定点化参数,并实时配置给高斯噪声模块;
步骤4:根据主控单元1-0经CPCI总线配置的通信频段,变频器设置上变频及下变频器的工作频段;本实施例中,变频器有VHF、L、S、C、X、Ku和Ka七个可选工作频段,选择VHF频段400MHz信号作为输入1输入下变频器3-1,本实施例设置下变频器的工作频段为400MHz,则下变频器3-1输出为70MHz中频信号。该中频信号经基带信道模拟单元叠加信道后输入到上变频器3-2中,本实施例设置上变频器3-2的工作频段为400MHz,则该叠加信道的中频信号经上变频器3-2上变频到400MHz。
步骤5:根据主控单元1-0经CPCI总线配置的衰减参数,衰减器单元实时设置衰减值。鉴于深空通信的远距离大衰落特性,本实施案例采用衰减器模块的方式实时设置路径损耗量。一个衰减器模块由一个大步进衰减器和一个小步进衰减器组成,本实施例中大步进衰减器的衰减范围为0~250dB,步进为10dB,小步进衰减器的衰减范围为0~10dB,步进为1dB,则大小步进衰减器组成的衰减器组的衰减范围为0~260dB,步进为1dB。衰减器模块实时设置的衰减值为步骤二中算得的总路径损耗减去基带信道模拟单元定点化导致的固定增益。
为实现信道状态更加平滑的过渡,可对运动轨迹进行内插,即所述的步骤1中的运动轨迹能够进行内插,内插倍数用户能够设置,本实施例内插200倍,即采用0.005s分辨率。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容做出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案内容,依据本发明的技术实质,在本发明的精神和原则之内,对以上实施例所作的任何简单的修改、等同替换与改进等,均仍属于本发明技术方案的保护范围之内。
Claims (2)
1.一种面向深空通信的无线信道实时模拟器的模拟方法,其特征在于,步骤如下:
步骤1:用户在主控单元的主控软件中输入深空通信场景参数,包括深空探测器飞行轨迹文件、通信频段、星体半径、电导率以及介电常数,主控软件基于收发端的运动轨迹计算传播场景的几何参量,包括收发端运动速度、直线距离和反射点坐标;
步骤2:通过传播场景的几何参量及系统参数计算得到不同时刻的信道状态参数,包括深空链路传播损耗L1(t),L2(t)、多普勒频移fd,1(t),fd,2(t)、各径时延τ1(t),τ2(t)、莱斯因子K和噪声功率PN这样的信道参数,步骤2的具体方法如下:
1)深空传播链路损耗计算
所述的深空传播链路损耗计算为传播链路损耗LLoss考虑自由空间传播损耗Lf,大气传播损耗La,降雨损耗Lr,收发天线增益Gr、Gt这样的因素,计算方法如式(1)所示:
LLoss=Lf+La+Lr-Gr-Gt(dB) (1)
其中如式(2)所示:
Lf=92.44+20lg(d/103)+20lg(f/109)(dB) (2)
式中,f为频率,单位Hz;d为通信距离,单位m;
如式(3)所示:
式中,r0和rw分别为氧气分子和水蒸气分子对电波的衰减率;hs为地面站海拔高度;θ为地面站天线波束仰角;g(h0)、g(hw)是与氧气层等效厚度h0、水蒸气层等效厚度hw相关的函数;
如式(4)所示:
Lr=γrLE(dB) (4)
式中,γr为降雨损耗率;LE为等效路径长度;
2)多普勒频移计算
所述的多普勒频移计算为凭借深空通信中,通信双方本身的运动速度很大,巨大的径向运动速度以及超高的通信频率产生了多普勒效应,这样直射路径的最大多普勒频移能够表示为如式(5)所示:
式中,分别为地面站和航天器的移动速度矢量,θB,θA则表示速度方向和直射路径的夹角,f0为信号频率,c为光速;
3)莱斯因子计算
所述的莱斯因子计算根据莱斯因子K与太阳活动造成的太阳闪烁指数具有如式(6)所示的关系而得:
其中,太阳闪烁指数m与太阳-地球-航天器夹角相关,当该夹角大于4°时m可忽略;当该夹角小于4°时,m值随着夹角的减小而增大;夹角趋近0°时m为1,此时K→∞,该情况下无直射路径;
4)信道噪声计算
所述的信道噪声计算是根据式(7)所示的深空信道噪声用等效噪声温度来计算,
PN=k(TRF+TC+TSUN+TMOON+TPL+TG+TAN)BN (7)
式中,PN为噪声功率;k为玻尔兹曼常量并且k=1.3806505×10-23J/K·Hz;BN为接收机等效噪声带宽;TRF为低噪声放大器的噪声温度,Tc为宇宙背景噪声温度,取2.66±0.77K;TG为地球站噪声,取值为0.1K;TAN为天线噪声温度;TSUN为太阳噪声温度,如(8)所示为
式中,S为单位太阳粒子辐射流密度;A为天线有效孔径面积;η为天线效率;r为考虑深空站系统接收机非线性因素后的减少因子;
TMOON为月球噪声温度,S波段,X波段和Ka波段噪声温度分别由式(9)、式(10)和式(11)求得,
TMOON(φ)=136±1.4(K) (9)
式中,φ为天线主波束宽度与月球圆盘中心之间的偏移角度;
TPL为行星噪声温度,如式(12)所示:
式中,Tk为行星黑体温度;d为行星直径;R为与行星的距离;θ为行星中心到天线波束中心的夹角;θ0为天线半功率波束宽度;G为天线增益;
步骤3:根据主控单元经CPCI总线配置的信道状态参数,基带信道模拟单元的DSP首先将其转为硬件参数并进行定点化,实时配置给对应的FPGA工作模块,步骤3的方法具体如下:
1)计算各径的相对时延τ=τ1(t)-τ2(t),其中τ1(t)、τ2(t)分别为直射径和散射径时延,将时延τ转化为定点的FPGA硬件参数,并实时配置给时延模块;
2)计算各径的相对损耗因子α1,α2,α1,α2应当满足(α2/α1)2=K的条件,将其转化为定点的FPGA硬件参数,并实时配置给损耗模块;
3)计算各径的多普勒频移fd,1(t),fd,2(t),将其转化为定点的FPGA硬件参数,并实时配置给多普勒频移模块;
4)计算信道衰落的定点化硬件参数,并实时配置给信道衰落模块;
5)计算深空信道噪声功率PN的定点化参数,并实时配置给高斯噪声模块;
步骤4:根据主控单元经CPCI总线配置的通信频段,变频器单元设置上变频器模块及下变频器模块的工作频段;
步骤5:根据主控单元经CPCI总线配置的衰减参数,衰减器单元实时设置衰减值。
2.根据权利要求1所述的面向深空通信的无线信道实时模拟器的模拟方法,其特征在于所述的步骤1中的运动轨迹能够进行内插,内插倍数用户能够设置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410519736.XA CN104243070B (zh) | 2014-09-29 | 2014-09-29 | 面向深空通信的无线信道实时模拟器及其模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410519736.XA CN104243070B (zh) | 2014-09-29 | 2014-09-29 | 面向深空通信的无线信道实时模拟器及其模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104243070A CN104243070A (zh) | 2014-12-24 |
CN104243070B true CN104243070B (zh) | 2017-03-15 |
Family
ID=52230505
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410519736.XA Active CN104243070B (zh) | 2014-09-29 | 2014-09-29 | 面向深空通信的无线信道实时模拟器及其模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104243070B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104753740B (zh) * | 2015-03-30 | 2018-04-27 | 中国空间技术研究院 | 一种深空通信仿真演示验证方法及仿真演示验证系统 |
CN105160153B (zh) * | 2015-08-04 | 2018-01-12 | 中国科学院上海天文台 | 一种基于先验模型的深空探测器三向多普勒测量方法 |
CN106026929B (zh) * | 2016-05-11 | 2019-04-02 | 中国科学院遥感与数字地球研究所 | 宽带Ka频段的下变频装置及下变频方法 |
CN106067772B (zh) * | 2016-06-12 | 2018-05-18 | 广州杰赛科技股份有限公司 | 一种射频模块的增益补偿方法与装置 |
CN106452628A (zh) * | 2016-09-04 | 2017-02-22 | 航天恒星科技有限公司 | 卫星信号模拟器及卫星信号模拟方法 |
CN108039931B (zh) | 2017-12-11 | 2020-10-23 | 上海埃威航空电子有限公司 | 一种基于地理信息的中高频通信信道的模拟测试系统 |
CN108270498B (zh) * | 2018-01-25 | 2020-07-17 | 清华大学 | 基于轨道驱动的深空信道模拟方法及装置 |
CN109495199B (zh) * | 2018-10-30 | 2021-07-13 | 航天恒星科技有限公司 | 一种基于半实物仿真的fpga参数调试系统 |
CN111163106B (zh) * | 2020-01-02 | 2022-03-25 | 上海航天测控通信研究所 | 一种基于天地网关的高速网络预处理装置 |
CN112087268B (zh) * | 2020-08-17 | 2021-06-29 | 北京航空航天大学 | 面向台风探测应用的临空无线信道模拟机 |
CN112671488A (zh) * | 2020-12-15 | 2021-04-16 | 上海卫星工程研究所 | 一种航天器信道模拟装置及其应用方法 |
CN113691335B (zh) * | 2021-08-23 | 2022-06-07 | 北京航空航天大学 | 一种覆盖多类损耗因素的通用电磁信号数据集构建方法 |
CN114040487B (zh) * | 2021-12-27 | 2022-08-12 | 扬州宇安电子科技有限公司 | 一种卫星上下行链路模拟设备的实时衰减控制方法 |
CN114675348B (zh) * | 2022-02-21 | 2024-03-08 | 北京航天飞行控制中心 | 一种深空数据传输窗口确定方法及装置 |
CN115996101B (zh) * | 2022-11-28 | 2024-05-24 | 西安电子科技大学 | 面向月球表面多场景通信的无线信道建模方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100386985C (zh) * | 2005-03-04 | 2008-05-07 | 华为技术有限公司 | 一种可扩展信道模拟范围的信道模拟装置和方法 |
CN201491016U (zh) * | 2009-06-11 | 2010-05-26 | 中国人民解放军92728部队 | 一种用于航空无线通信系统测试的信道模拟装置 |
CN102932072B (zh) * | 2012-08-08 | 2015-04-15 | 中国传媒大学 | 高动态卫星信道模拟器 |
WO2014049194A1 (en) * | 2012-09-27 | 2014-04-03 | Elektrobit System Test Ltd. | Radio channel emulation |
CN103873132A (zh) * | 2012-12-13 | 2014-06-18 | 许雪梅 | 基于PXIe总线的卫星通信系统模拟器 |
-
2014
- 2014-09-29 CN CN201410519736.XA patent/CN104243070B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104243070A (zh) | 2014-12-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104243070B (zh) | 面向深空通信的无线信道实时模拟器及其模拟方法 | |
Wang et al. | 6G wireless channel measurements and models: Trends and challenges | |
CN103532644B (zh) | 多径阴影复合衰落信道模拟装置及其工作方法 | |
CN106452628A (zh) | 卫星信号模拟器及卫星信号模拟方法 | |
Zhu et al. | 3D non‐stationary geometry‐based multi‐input multi‐output channel model for UAV‐ground communication systems | |
CN105049142A (zh) | 一种双通道静态基带信道模拟装置及方法 | |
WO2023169590A1 (zh) | 一种适用于全频段全场景的6g普适信道建模方法 | |
Bai et al. | Channel modeling for satellite communication channels at Q-band in high latitude | |
CN102447518A (zh) | 一种用于临近空间高超声速条件下的信道综合处理方法 | |
CN105099588A (zh) | 航空通信系统干扰传播信道模拟装置及方法 | |
CN114499724B (zh) | 一种面向低轨卫星通信的空时频非平稳传输特性分析方法 | |
Brown | A physical meteor-burst propagation model and some significant results for communication system design | |
CN108270498B (zh) | 基于轨道驱动的深空信道模拟方法及装置 | |
Cameron et al. | Evaluation of the effect of sporadic-E on high frequency radio wave propagation in the Arctic | |
Wang et al. | 6G THz propagation channel characteristics and modeling: Recent developments and future challenges | |
CN103457678A (zh) | 多云天气条件下的气象卫星复合信道仿真模型的建立方法 | |
Yang et al. | Statistical modeling of the high altitude platform dual-polarized MIMO propagation channel | |
Latif et al. | Simulation of radio signals from cosmic-ray cascades in air and ice as observed by in-ice Askaryan radio detectors. | |
Osborne et al. | Propagation characterization of LEO/MEO satellite systems at 900-2100 MHz | |
Wu et al. | Research on Mars Communication Link Budget | |
Polo | Model to investigate the interference induced by NGSO constellations on GSO system ground stations | |
Parks et al. | Wideband characterisation and modelling of the mobile satellite propagation channel at L-and S-bands | |
CN106788672A (zh) | 一种提高卫星系统间干扰仿真步长精度的方法 | |
Wahl | Using Systems Tool Kit to Model Possible CubeSat Orbits | |
Zhang et al. | Real-time aeronautical channel simulator |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |