CN115150012B - 电离层信道高精度色散时延特性实时模拟方法与系统 - Google Patents

电离层信道高精度色散时延特性实时模拟方法与系统 Download PDF

Info

Publication number
CN115150012B
CN115150012B CN202210925414.XA CN202210925414A CN115150012B CN 115150012 B CN115150012 B CN 115150012B CN 202210925414 A CN202210925414 A CN 202210925414A CN 115150012 B CN115150012 B CN 115150012B
Authority
CN
China
Prior art keywords
delay
channel
dispersion
signal
simulation
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
Application number
CN202210925414.XA
Other languages
English (en)
Other versions
CN115150012A (zh
Inventor
李蓬蓬
吕志成
肖伟
王飞雪
刘文祥
陈雷
倪少杰
李柏渝
牟卫华
钟水彬
覃曼丽
谢名赞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by National University of Defense Technology filed Critical National University of Defense Technology
Publication of CN115150012A publication Critical patent/CN115150012A/zh
Application granted granted Critical
Publication of CN115150012B publication Critical patent/CN115150012B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/391Modelling the propagation channel
    • H04B17/3912Simulation models, e.g. distribution of spectral power density or received signal strength indicator [RSSI] for a given geographic region
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/03Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/07Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
    • G01S19/072Ionosphere corrections
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/20Integrity monitoring, fault detection or fault isolation of space segment
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Electromagnetism (AREA)
  • Signal Processing (AREA)
  • Computer Security & Cryptography (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明涉及电离层信道高精度色散时延特性实时模拟方法与系统,方法包括:获取模拟射频信号,将模拟射频信号进行采样及预处理,得到数字信号;根据电离层色散信道伪码群时延模型,计算实现信道模拟伪码群时延所需的参数;对伪码时延进行分离,获得伪码的频率选择群时延特性分量和滞后群时延特性分量,采用不同的数字滤波方法分别进行实时模拟;根据电离层色散信道载波相时延模型,计算实现信道模拟载波超前相时延所需参数,进行精准实时模拟;通过驱动信道模拟器启动电离层信道模拟流程,输出叠加了对应电离层信道特性的射频信号。本发明满足了电离层色散信道频率选择性延迟效应、载波和伪码时延“超前‑滞后”效应的高精度实时模拟需求。

Description

电离层信道高精度色散时延特性实时模拟方法与系统
技术领域
本发明涉及通信、导航信道模拟技术领域,特别是涉及一种电离层信道高精度色散时延特性实时模拟方法。
背景技术
电离层通常指高度位于60km以上到磁层顶层之间的整个空间,由于太阳的辐射,电离层中部分气体分子被电离,形成大量的自由电子和正离子。当电磁波穿过电离层时,电离层是一种色散介质,引起信号传输特性的变化与射频信号的频率直接相关。对于宽带调制的导航信号而言,经过电离层传播后信号带内不同频点处产生的延迟量互不相同,是电离层色散效应的直接体现。对于模拟信号空间环境传输带来各种特性变换的信道模拟器而言,电离层信道模拟是其中的重要环节,而一般电离层信道模拟大多关注与电离层闪烁、折射等带来的信号幅度的衰落特性,将宽带信道内电离层传输带来的时延变化,等效为信号中心频点处的电离层延迟进行简单化处理。这种简化处理处理方式,对于以数据传输为基础的通信信号、或带宽较窄的扩频信号是可行的,而对于以高精度测距为基础的宽带导航信道而言,会产生一定程度的与电离层特性相关的相关损耗及载波跟踪偏差,因此不能一概视之简化处理。因此需要针对性研究宽带信号的电离层信道色散时延特性模拟方法,及其在信道模拟器中的设计与实现。
发明内容
本发明提供了一种电离层信道高精度色散时延特性实时模拟方法,以解决上述现有技术中存在的问题。
为实现上述目的,本发明提供了如下方案:
电离层信道高精度色散时延特性实时模拟方法,包括:
获取模拟射频信号,对所述模拟射频信号进行采样及预处理,得到数字信号;
计算实现信道模拟伪码群时延所需的参数,对伪码频率选择和滞后群时延分量进行分离,获得伪码频率选择群时延分量和滞后群时延分量;
采用任意群时延特性数字全通滤波器设计方法,实现伪码频率选择群时延特性分量实时模拟;
采用线性相位数字滤波器设计方法,实现滞后群时延特性分量实时模拟;
构建电离层色散信道载波相时延模型,获取载波超前相色散时延模拟;
采用基于查表和线性插值结合的可控相位载波生成方法,实现载波超前相时延特性模拟;
将数字信号进行数字上变频及DAC信号恢复,滤波后得到模拟射频信号输出;
将完成所述色散时延模拟的数字基带信号通过驱动信道模拟器启动电离层信道模拟控制,输出叠加了对应电离层信道特性的射频信号。
优选地,对所述模拟射频信号进行采样及预处理,包括:将所述模拟射频信号进行模拟滤波处理后,进行ADC采样,并进行数字滤波、数字下变频、二次量化操作,获得所述数字信号;
进行所述ADC采样时,控制所述射频信号滤波带宽,用于保持采样信号带内群时延特性平坦,避免引入其他色散特性群时延变化;所述二次量化操作用于根据所述ADC采样的数据有效位情况,量化操作需要充分保留采样数据有效位。
优选地,计算色散时延参数,并对所述伪码频率选择和所述滞后群时延分量进行分离,获得所述伪码频率选择群时延分量和滞后群时延分量;
包括:导航卫星播发的导航信号通过电离层后,宽带信道电离层时延呈现与频率相关的色散特性,通过电离层色散信道伪码群时延模型,计算所述色散时延参数Δτ(f):
Figure BDA0003779232530000031
式中,f为频率,a2、a3、a4均为与频率无关的函数系数;
确定宽带信号带内频率最大值fH,将所述色散时延参数表达式进行拆分,如下式:
Figure BDA0003779232530000032
其中,TEC为电子总数,c为光速;
确定伪码的滞后群时延分量Δτ2(fH),以及频率选择群时延分量Δτ1(f)分别如下:
Figure BDA0003779232530000041
优选地,实现所述伪码频率选择群时延特性分量实时模拟,包括:
基于数字全通滤波器模拟实现,数字全通滤波器频率响应函数如下式表达:
Figure BDA0003779232530000042
其中,ω是归一化角频率,n=0,1,…N,an是系数,N(ω)和D(ω)分别是H(e)的分子多项式和分母多项式的相频响应;
根据给定的所述伪码频率选择群时延特性分量,求解分母多项式的群时延函数,得到呈现电离层频率选择色散群时延特性的数字全通滤波器表达式,用于实时模拟伪码频率选择群时延特性分量。
优选地,实现所述滞后群时延特性分量实时模拟,包括:
采用线性相位特性的数字滤波器模拟实现,根据所述伪码滞后时延特性分量,求解滤波器系数,得到符合要求的线性相位数字全通滤波器表达式,用于模拟信号伪码的滞后群时延特性分量。
优选地,构建所述电离层色散信道载波相时延模型,包括:
电离层色散时延会对导航信号的载波产生相时延,所述相时延参数Δφ(f)在忽略高阶项以后为:
Figure BDA0003779232530000051
式中,f是载波频率,TEC是电子总数,c为光速;
若扩频信号的载波频率已知,载波信号频率为f0,则所述电离层色散信道载波相时延模型为:
Figure BDA0003779232530000052
基于所述电离层色散信道载波相时延模型,计算出相应的初始相位,采用加入初始载波相位的方法进行相应的相时延模拟。
优选地,所述载波超前相色散时延模拟,包括:
建立可控初始相位数字载波信号模型,对于载波信号S(t)如下:
s(t)=sin(2πf*t+θ)
其中,f为载波频率,θ为相位;
基于查表和线性插值结合的可控相位载波生成方法,模拟所述载波的超前相时延特性,基于FPGA利用查表法,生成所述初始相位可控的数字载波信号,所述数字载波信号对应载波的超前色散相时延模拟所需相应初始化相位值。
优选地,得到所述模拟射频信号输出包括:
对数字信号进行数模转换与射频信号滤波,将所述完成电离层特性模拟的数字信号进行数字上变频、DAC信号恢复,模拟滤波处理后,获得所述模拟射频信号;
进行所述模拟滤波时,控制所述滤波器带宽,用于保持信号带内群时延特性平坦,避免引入其他色散特性群时延变化,同时滤除相应谐波和镜频分量,保证射频信号输出质量。
优选地,输出所述叠加了对应电离层信道特性的射频信号,包括:
将需要模拟电离层信道特性的宽带信号收发装置输入和输出端口通过射频线缆与信道模拟器连接,操作所述信道模拟器的上位机软件,设定进行电离层信道仿真的场景,计算获取电离层信道特性相关时延快变参数,同时设置信道模拟设备运行的系统时间、信号输入输出功率、信道链路映射关系;启动电离层信道模拟控制流程,完成信道群时延特性模拟,输出所述叠加了对应电离层信道特性的射频信号。
为了实现上述目的,本发明还提供了一种电离层信道高精度色散时延特性实时模拟系统,包括:
存储器,用于存储计算机程序;
处理器,用于当执行所述计算机程序时,实现电离层信道高精度色散时延特性实时模拟方法。
本发明的有益效果为:
本发明基于伪码和载波剥离处理的方法,对伪码色散群时延和载波色散相时延进行了高逼真的模拟实现,满足了信道模拟器对宽带导航信号电离层色散信道频率选择性延迟效应、载波和伪码时延“超前-滞后”效应的高逼真高精度模拟需求。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例的方法流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
窄带信号条件下,电离层信道传输引起的频率选择效应,对于进行伪距测量的导航信号接收性能影响可以忽略,引起的载波和伪码时延“超前-滞后”效应,对于不采用载波辅助跟踪的普通导航信号接收性能的影响可以忽略。在导航信道带宽逐步拓宽的背景下,上述效应都会导致导航信号接收存在相关损耗、载波跟踪偏差及测距精度偏差等不可忽略的影响。因此,电离层色散信道频率选择性延迟效应、载波和伪码时延“超前-滞后”效应,是信道模拟器对宽带导航信号进行高精度信道模拟的关键环节,需要优化设计。
为了满足上述需求,本发明构建了电离层信道高精度色散时延特性实时模拟方法,包括射频信号采样及预处理、电离层色散信道伪码群时延模型建立、伪码的频率选择和滞后群时延分量分离分析、伪码的频率选择群时延分量模拟、伪码的滞后群时延分量模拟、电离层色散信道载波相时延模型建立、载波的“超前”相时延模拟、数模转换与射频信号滤波、驱动信道模拟器完成电离层信道特性模拟操作等,实现了电离层信道高精度色散时延特性实时模拟。
电离层信道模拟具体实施环节如图1所示,需要注意图中所包含的环节不代表实施顺序。其中,射频信号采样及预处理,将接收到的模拟射频信号,进行ADC采样,并进行数字滤波、数字下变频、二次量化等数字信号预处理操作;电离层色散信道伪码群时延模型建立,用以计算信道模拟伪码群时延实现所需的各类参数;伪码的频率选择和滞后群时延分量分离,伪码群时延根据特性区分为两种不同类型分量,包括频率选择群时延分量和滞后群时延分量;伪码的频率选择色散群时延分量模拟,给出了详细的任意群时延特性数字全通滤波器设计方法,并以此模拟伪码的频率选择群时延特性模拟;伪码的滞后群时延分量模拟,采用线性相位数字滤波器,实现精准的时延滞后操作;电离层色散信道载波相时延模型建立,实现载波相时延模拟所需的参数计算;载波的“超前”相时延模拟,给出了基于查表和线性插值结合的可控相位载波生成方法,并以此模拟载波的“超前”色散相时延特性;数模转换与射频信号滤波,将信道模拟后的数字信号转变为模拟信号,进行适当的模拟滤波设计;驱动信道模拟器完成电离层信道特性模拟操作,给出了基于电离层信道数学模型仿真计算参数,驱动信道模拟器实现对通过的射频信号实现电离层信道色散特性模拟的全部操作流程。
射频信号采样及预处理。该步骤将接收到的模拟射频信号,进行恰当的模拟滤波处理后,进行ADC采样,并进行数字滤波、数字下变频、二次量化等数字信号预处理操作,便于后续数字信号处理的展开。
射频信号ADC采样时,需要适当控制射频信号滤波带宽,保持采样信号带内群时延特性良好,避免引入其他的色散特性群时延变化。
对ADC采样信号进行数字滤波,需要注意合理的数字滤波器带宽设计。
二次量化操作根据ADC的数据有效位情况,量化操作需要充分保留采样数据有效位。
电离层色散信道伪码群时延模型建立。导航信号通过电离层色散信道,信号伪码的色散群时延模型与参数计算,该步骤建立了电离层色散信道数学模型,计算信道信号伪码色散时延模拟实现所需参数。
宽带信道电离层时延呈现与频率相关的色散特性,其色散时延参数计算Δτ(f)如下所示:
Figure BDA0003779232530000091
式中,f为频率,a2、a3、a4均为与频率无关的函数系数。
参数计算简化。考虑到高阶分量在实际测量中对测距值的影响可以忽略,因此将相关的系数代入后,上述公式(1)可以简化为如下表达:
Figure BDA0003779232530000101
式中,TEC是电子总数,c为光速。
伪码的频率选择和滞后群时延分量分离。进一步分析发现,上述色散时延参数计算Δτ(f)是一个与f2呈反比关系的表达式。为了便利时延模拟时间,需要将伪码的“滞后”和“频率选择”群时延分离。
模型表达式拆分。假设已知宽带信号带内频率最大值为fH,将表达式拆分成为如下形式:
Figure BDA0003779232530000102
因此,可以根据上式(3),将伪码的“频率选择”和“滞后”群时延分离开来,从而便于后续通过两组不同的时延滤波器分开处理实现。
伪码的“频率选择”和“滞后”群时延分量分离模型:
假设伪码的滞后群时延分量为Δτ2(fH),以及频率选择群时延分量Δτ1(f),则上式(3)进一步表达如下式(4):
Δτ(f)=Δτ1(f)+Δτ2(fH)  (4)
其中伪码的频率选择群时延分量和滞后群时延分量分别满足如下表达式(5):
Figure BDA0003779232530000111
伪码的频率选择色散群时延分量模拟,该步骤给出了详细的任意群时延特性数字全通滤波器设计方法,并以此模拟伪码频率选择群时延特性。
任意群时延特性的数字全通滤波器模型建立。数字全通滤波器频率响应函数如下表达式(6):
Figure BDA0003779232530000112
其中,ω是归一化角频率,n=0,1,…N,an是系数,N(ω)和D(ω)分别是H(e)的分子多项式和分母多项式的相频响应。
显然易见,对于全通滤波器,其幅值满足如下表达式(7):
|H(e)|=1  (7)
上述数字滤波器,在幅频相应特性上呈现全通特性,在群时延特性上可以呈现电离层“频率选择”色散群时延特性。
根据给定的宽带信号电离层“频率选择”色散群时延求解分母多项式的群时延函数。
对于给定的期望宽带信号电离层“频率选择”色散群时延特性τ(ω),通过得到分母多项式的群时延函数τD(ω)如下表达式(8):
Figure BDA0003779232530000113
根据分母多项式的群时延函数τD(ω)与复倒谱系数c(k)之间的数学关系,进行逆傅里叶变换,可以求解除复倒谱系数。
计算分母多项式的系数an
对于已知的复倒谱系数c(k),可以计算分母多项式的系数an,参照如下表达式(9)关系:
Figure BDA0003779232530000121
根据系数an,可以代入表达式(6),从而得到呈现电离层“频率选择”色散群时延特性的数字全通滤波器表达式。
伪码的滞后群时延分量模拟方法。基于线性相位数字滤波器设计实现伪码的“滞后”群时延模拟。
伪码的“滞后”群时延特性滤波器设计。根据信号伪码的“滞后”群时延特性表达式可知,其受到已知宽带信号带内频率最大值为fH和电子总数TEC决定,与宽带信号频率范围无关,在信号带宽范围内呈平坦群时延特性,适合采用线性相位特性的数字滤波器模拟实现。
伪码的“滞后”群时延特性模拟。获得了符合要求的线性相位数字全通滤波器表达式后,即可通过伪码群时延处理流程,实现信号伪码的滞后群时延模拟。
电离层色散信道载波相时延模型建立,进行导航信号通过电离层色散信道,信号载波的“超前”相时延模型与参数计算。
电离层色散时延会对导航信号的载波产生相时延,对于已知信号,其载波值一般是固定的,其相时延参数Δφ(f)在忽略高阶项以后计算如下式(10)所示:
Figure BDA0003779232530000122
式中,f是载波频率,TEC是电子总数,c为光速。
载波的“超前”相时延模型简化。考虑到扩频信号的载波频率已知,假设载波信号频率为f0,则信号载波的“超前”相时延模型简化表达如下式(11):
Figure BDA0003779232530000131
根据表达式(11)可算出相应的初始相位,采用加入初始载波相位的方法进行相应的相时延模拟。
载波的“超前”色散相时延模拟方法。该步骤给出了基于查表和线性插值结合的可控相位载波生成方法,并以此模拟载波的“超前”相时延特性。
可控初始相位数字载波信号模型。对于载波信号S(t)如下:
s(t)=sin(2πf*t+θ)  (12)
其中,f为载波频率,θ为相位。
将载波信号在一个周期T内分成N等份,相邻等分点差值Dk,每个时间间隔TN=T/N。假设T0为数字信号处理系统时钟周期,其中f0=1/T0为系统时钟频率。在每一个间隔TN内用一个数字累加器进行线性插值,对于给定的f0、f、N,在第k和k+1个等分点之间线性插值累加当量为:
Figure BDA0003779232530000132
因此在任意一个TN间隔内,以系统时钟T0为参考时钟,对载波信号幅值增量ΔDk进行累加得到数字化的载波信号,对应相应初始化相位值。
数字载波信号生成。基于FPGA利用查表法,即可基于上述方法生成初始相位可控的数字载波信号,对应载波的“超前”色散相时延模拟所需相应初始化相位值。
数模转换与射频信号滤波。对数字信号进行数模转换,生成模拟射频信号。该步骤将完成色散时延模拟后得的数字信号,进行恰当的数字滤波处理后,进行DAC信号恢复,生成模拟射频信号。
完成色散时延模拟的数字基带信号,需要进行适当的数字插值,用以匹配数字DAC信号恢复需求。
生成的模拟射频信号,需要适当控制射频信号滤波带宽,保持采样信号带内群时延特性良好,避免引入其他的色散特性群时延变化。
驱动信道模拟器完成电离层信道特性模拟操作。该步骤给出了基于电离层信道数学模型仿真计算参数,驱动信道模拟器实现对通过的射频信号实现电离层信道特性模拟的操作实施流程。
将需要模拟电离层信道特性的宽带信号收发装置输入和输出端口通过射频线缆与信道模拟器连接,提供上述设备工作必要的电源、时频信号、控制信号等。
操作信道模拟器的上位机软件,进行电离层信道仿真场景的设定,计算获取电离层信道特性相关时延等快变参数,同时设置信道模拟设备运行的系统时间、信号输入输出功率、信道链路映射关系等必要参数。
启动电离层信道模拟控制,完成电离层数学模型计算所得的信道群时延特性模拟实施,输出叠加了对应电离层信道特性的射频信号。
本发明还涉及一种电离层信道高精度色散时延特性实时模拟系统,包括存储器,用于存储计算机程序;处理器,用于当执行所述计算机程序时,实现如上述的电离层信道高精度色散时延特性实时模拟方法。
本发明基于伪码和载波剥离处理的方法,对伪码色散群时延和载波色散相时延进行了高逼真的模拟实现,满足了信道模拟器对宽带导航信号电离层色散信道频率选择性延迟效应、载波和伪码时延“超前-滞后”效应的高逼真高精度模拟需求。
以上所述的实施例仅是对本发明优选方式进行的描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。

Claims (10)

1.电离层信道高精度色散时延特性实时模拟方法,其特征在于,包括:
获取模拟射频信号,对所述模拟射频信号进行采样及预处理,得到数字信号;
计算实现信道模拟伪码群时延所需的参数,对伪码频率选择和滞后群时延分量进行分离,获得伪码频率选择群时延分量和滞后群时延分量;
采用任意群时延特性数字全通滤波器设计方法,实现伪码频率选择群时延特性分量实时模拟;
采用线性相位数字滤波器设计方法,实现滞后群时延特性分量实时模拟;
构建电离层色散信道载波相时延模型,获取载波超前相色散时延模拟;
采用基于查表和线性插值结合的可控相位载波生成方法,实现载波超前相时延特性模拟;
将数字信号进行数字上变频及DAC信号恢复,滤波后得到模拟射频信号输出;
将完成所述色散时延模拟的数字基带信号通过驱动信道模拟器启动电离层信道模拟控制,输出叠加了对应电离层信道特性的射频信号。
2.根据权利要求1所述的电离层信道高精度色散时延特性实时模拟方法,其特征在于,对所述模拟射频信号进行采样及预处理,包括:将所述模拟射频信号进行模拟滤波处理后,进行ADC采样,并进行数字滤波、数字下变频、二次量化操作,获得所述数字信号;
进行所述ADC采样时,控制所述射频信号滤波带宽,用于保持采样信号带内群时延特性平坦,避免引入其他色散特性群时延变化;所述二次量化操作用于根据所述ADC采样的数据有效位情况,量化操作需要充分保留采样数据有效位。
3.根据权利要求1所述的电离层信道高精度色散时延特性实时模拟方法,其特征在于,计算色散时延参数,并对所述伪码频率选择和所述滞后群时延分量进行分离,获得所述伪码频率选择群时延分量和滞后群时延分量;
包括:导航卫星播发的导航信号通过电离层后,宽带信道电离层时延呈现与频率相关的色散特性,通过电离层色散信道伪码群时延模型,计算所述色散时延参数Δτ(f):
式中,f为频率,a2、a3、a4均为与频率无关的函数系数;
确定宽带信号带内频率最大值fH,将所述色散时延参数表达式进行拆分,如下式:
其中,TEC为电子总数,c为光速;
确定伪码的滞后群时延分量Δτ2(fH),以及频率选择群时延分量Δτ1(f)分别如下:
4.根据权利要求1所述的电离层信道高精度色散时延特性实时模拟方法,其特征在于,实现所述伪码频率选择群时延特性分量实时模拟,包括:
基于数字全通滤波器模拟实现,数字全通滤波器频率响应函数如下式表达:
其中,ω是归一化角频率,n=0,1,…N,an是系数,N(ω)和D(ω)分别是H(e)的分子多项式和分母多项式的相频响应;
根据给定的所述伪码频率选择群时延特性分量,求解分母多项式的群时延函数,得到呈现电离层频率选择色散群时延特性的数字全通滤波器表达式,用于实时模拟伪码频率选择群时延特性分量。
5.根据权利要求1所述的电离层信道高精度色散时延特性实时模拟方法,其特征在于,实现所述滞后群时延特性分量实时模拟,包括:
采用线性相位特性的数字滤波器模拟实现,根据所述伪码滞后时延特性分量,求解滤波器系数,得到符合要求的线性相位数字全通滤波器表达式,用于模拟信号伪码的滞后群时延特性分量。
6.根据权利要求1所述的电离层信道高精度色散时延特性实时模拟方法,其特征在于,构建所述电离层色散信道载波相时延模型,包括:
电离层色散时延会对导航信号的载波产生相时延,所述相时延参数Δφ(f)在忽略高阶项以后为:
式中,f是载波频率,TEC是电子总数,c为光速;
若扩频信号的载波频率已知,载波信号频率为f0,则所述电离层色散信道载波相时延模型为:
基于所述电离层色散信道载波相时延模型,计算出相应的初始相位,采用加入初始载波相位的方法进行相应的相时延模拟。
7.根据权利要求1所述的电离层信道高精度色散时延特性实时模拟方法,其特征在于,所述载波超前相色散时延模拟,包括:
建立可控初始相位数字载波信号模型,对于载波信号s(t)如下:
s(t)=sin(2πf*t+θ)
其中,f为载波频率,θ为相位;
基于查表和线性插值结合的可控相位载波生成方法,模拟所述载波的超前相时延特性,基于FPGA利用查表法,生成所述初始相位可控的数字载波信号,所述数字载波信号对应载波的超前色散相时延模拟所需相应初始化相位值。
8.根据权利要求1所述的电离层信道高精度色散时延特性实时模拟方法,其特征在于,得到所述模拟射频信号输出包括:
对数字信号进行数模转换与射频信号滤波,将完成电离层特性模拟的数字信号进行数字上变频、DAC信号恢复,模拟滤波处理后,获得所述模拟射频信号;
进行所述模拟滤波时,控制所述滤波器带宽,用于保持信号带内群时延特性平坦,避免引入其他色散特性群时延变化,同时滤除相应谐波和镜频分量,保证射频信号输出质量。
9.根据权利要求1所述的电离层信道高精度色散时延特性实时模拟方法,其特征在于,输出所述叠加了对应电离层信道特性的射频信号,包括:
将需要模拟电离层信道特性的宽带信号收发装置输入和输出端口通过射频线缆与信道模拟器连接,操作所述信道模拟器的上位机软件,设定进行电离层信道仿真的场景,计算获取电离层信道特性相关时延快变参数,同时设置信道模拟设备运行的系统时间、信号输入输出功率、信道链路映射关系;启动电离层信道模拟控制流程,完成信道群时延特性模拟,输出所述叠加了对应电离层信道特性的射频信号。
10.一种电离层信道高精度色散时延特性实时模拟系统,其特征在于,包括:
存储器,用于存储计算机程序;
处理器,用于当执行所述计算机程序时,实现如权利要求1-9任一项所述的电离层信道高精度色散时延特性实时模拟方法。
CN202210925414.XA 2022-07-07 2022-08-03 电离层信道高精度色散时延特性实时模拟方法与系统 Active CN115150012B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202210802092X 2022-07-07
CN202210802092 2022-07-07

Publications (2)

Publication Number Publication Date
CN115150012A CN115150012A (zh) 2022-10-04
CN115150012B true CN115150012B (zh) 2023-04-28

Family

ID=83413550

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210925414.XA Active CN115150012B (zh) 2022-07-07 2022-08-03 电离层信道高精度色散时延特性实时模拟方法与系统

Country Status (1)

Country Link
CN (1) CN115150012B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117639963B (zh) * 2024-01-26 2024-04-02 中国人民解放军国防科技大学 导航信道色散时延模拟方法、系统、设备及存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103825662A (zh) * 2014-01-29 2014-05-28 中国电子科技集团公司第十研究所 电离层色散效应模拟器
CN104199066A (zh) * 2014-09-19 2014-12-10 中国人民解放军国防科学技术大学 一种boc信号电离层色散效应的高精度模拟方法
CN104407323A (zh) * 2014-12-11 2015-03-11 中国工程物理研究院电子工程研究所 一种高动态低信噪比扩频信号伪码时延测量方法
CN105659901B (zh) * 2013-01-31 2015-04-01 北京理工大学 适用于宽带星地链路无线电系统的电离层色散模拟方法
EP3139199A1 (en) * 2015-09-04 2017-03-08 Airbus DS GmbH Wireless communication unit, integrated circuit, satellite communication system and method for compensating for ionospheric group delay
CN106932792A (zh) * 2017-03-29 2017-07-07 中国人民解放军国防科学技术大学 一种基于软件接收机的导航信号模拟器时延测量方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105659901B (zh) * 2013-01-31 2015-04-01 北京理工大学 适用于宽带星地链路无线电系统的电离层色散模拟方法
CN103825662A (zh) * 2014-01-29 2014-05-28 中国电子科技集团公司第十研究所 电离层色散效应模拟器
CN104199066A (zh) * 2014-09-19 2014-12-10 中国人民解放军国防科学技术大学 一种boc信号电离层色散效应的高精度模拟方法
CN104407323A (zh) * 2014-12-11 2015-03-11 中国工程物理研究院电子工程研究所 一种高动态低信噪比扩频信号伪码时延测量方法
EP3139199A1 (en) * 2015-09-04 2017-03-08 Airbus DS GmbH Wireless communication unit, integrated circuit, satellite communication system and method for compensating for ionospheric group delay
CN106932792A (zh) * 2017-03-29 2017-07-07 中国人民解放军国防科学技术大学 一种基于软件接收机的导航信号模拟器时延测量方法

Also Published As

Publication number Publication date
CN115150012A (zh) 2022-10-04

Similar Documents

Publication Publication Date Title
US9231733B2 (en) Method and system for updating multi-frequency-band pre-distortion coefficient lookup table
US7432856B1 (en) Signal simulation device
CN103281268B (zh) 一种卫星导航信号生成系统带内群时延波动的预补偿方法
US20060126852A1 (en) Method and system for processing a sound field representation
CN115150012B (zh) 电离层信道高精度色散时延特性实时模拟方法与系统
JP3730961B2 (ja) チャンネル・シミュレーションの実行方法及びチャンネル・シミュレータ
CN111190197B (zh) 一种导航卫星信号质量在轨优化与保持方法
CN109039508B (zh) 无线多径衰落信道仿真系统及方法
US9874640B2 (en) Automated methods for suppression of spurious signals
JPH1098325A (ja) 無線通信システムで干渉を最小化し雑音の影響を減らすための信号処理装置及びその方法
CN103312347A (zh) 一种多通道干扰信号产生装置及干扰信号产生方法
CN103116170A (zh) Gnss的基于天线阵列干扰抑制模块的室内测试系统
JP5368572B2 (ja) データ収集及びシミュレーション
CN104539340A (zh) 一种基于稀疏表示和协方差拟合的稳健波达角估计方法
EP3798682A1 (en) Wavefront global navigation satellite system and interference simulator systems and methods of use thereof
CN101682367B (zh) 线缆网络中的信号质量确定
CN111371454A (zh) 基于低轨卫星的多普勒信号模拟方法、装置及终端
CN103716018B (zh) 数字宽带激励源实现装置及方法
CN113625310A (zh) 一种大范围高线性正交信号幅频动态模拟方法及模拟系统
Alimohammad et al. Filter-based fading channel modeling
CN114465684B (zh) 一种Bi-Gaussian信道的仿真方法及装置
RU2767180C1 (ru) Устройство имитации канала спутниковой связи с выпускным буксируемым антенным устройством в дмв диапазоне
Fitting Wideband troposcatter radio channel simulator
RU2660126C1 (ru) Способ формирования группового навигационного сигнала ГЛОНАСС
Werner et al. High-Fidelity Simulation of a Pico Satellite Link

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