CN113472464B - 用于实时模拟高动态卫星信道多普勒特性的方法和系统 - Google Patents

用于实时模拟高动态卫星信道多普勒特性的方法和系统 Download PDF

Info

Publication number
CN113472464B
CN113472464B CN202111035878.5A CN202111035878A CN113472464B CN 113472464 B CN113472464 B CN 113472464B CN 202111035878 A CN202111035878 A CN 202111035878A CN 113472464 B CN113472464 B CN 113472464B
Authority
CN
China
Prior art keywords
dynamic
doppler
frequency
delay
signal
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
CN202111035878.5A
Other languages
English (en)
Other versions
CN113472464A (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
Priority to CN202111035878.5A priority Critical patent/CN113472464B/zh
Publication of CN113472464A publication Critical patent/CN113472464A/zh
Application granted granted Critical
Publication of CN113472464B publication Critical patent/CN113472464B/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
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/0009Time-delay networks
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/40Monitoring; Testing of relay systems

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Radio Relay Systems (AREA)

Abstract

本发明提出一种用于实时模拟高动态卫星信道多普勒特性的方法和系统。通过提高中频采样信号的采样率,对中频采样信号进行预处理,以获取射频采样信号;基于射频采样信号的信号参数和多普勒动态参数,确定多普勒模拟动态时延;通过累加控制对多普勒模拟动态时延执行插值处理,以获取多普勒模拟高频度动态时延,多普勒模拟高频度动态时延以动态时延控制字的形式来表征;将动态时延控制字划分为整数部分的第一动态时延和小数部分的第二动态时延;将射频采样信号分别发送至数据缓存区和可变时延滤波器组,以获得经数据缓存区存储转发的第一输出信号和经可变时延滤波器组滤波的第二输出信号;整合第一输出信号与第二输出信号。

Description

用于实时模拟高动态卫星信道多普勒特性的方法和系统
技术领域
本发明属于卫星技术领域,尤其涉及一种用于实时模拟高动态卫星信道多普勒特性的方法和系统。
背景技术
卫星信道模拟器能够实时模拟信号传播特性,是卫星发射前在地面对接测试不可缺少的设备。卫星信号接收系统通过测量信号的传播时延来测量卫星到接收机之间的距离,从而进行定位解算。当接收机与卫星间存在相对速度时,其在信号上体现为信号的多普勒效应,对信号的传播时延测量具有较大的影响。
因此,实现信号高动态多普勒特性模拟对于卫星信道的模拟具有重要意义。如何利用卫星信道模拟器对高动态多普勒特性进行模拟,从而实现实时多普勒模拟是卫星技术领域亟需解决的技术问题之一。
发明内容
针对上述问题,本发明提出一种用于实时模拟高动态卫星信道多普勒特性的方案,以解决上述技术问题。
本发明第一方面公开了一种用于实时模拟高动态卫星信道多普勒特性的方法,所述方法包括:
步骤S1、通过提高中频采样信号的采样率,对所述中频采样信号进行预处理,以获取射频采样信号;
步骤S2、基于所述射频采样信号的信号参数和多普勒动态参数,确定多普勒模拟动态时延,其中:所述信号参数包括多数射频采样信号的初始载波频率和初始相位;所述多普勒动态参数包括多普勒频率和多普勒频率变化率;
步骤S3、利用相位加速度累加器和相位速度累加器,通过累加控制对所述多普勒模拟动态时延执行插值处理,以获取多普勒模拟高频度动态时延,所述多普勒模拟高频度动态时延以动态时延控制字的形式来表征;
步骤S4、将所述动态时延控制字划分为整数部分的第一动态时延和小数部分的第二动态时延,所述第一动态时延通过数据缓存区的存储转发来实现,所述第二动态时延由基于多相滤波结构的可变时延滤波器组来实现;
步骤S5、将所述射频采样信号分别发送至所述数据缓存区和所述可变时延滤波器组,以获得经所述数据缓存区存储转发的第一输出信号和经所述可变时延滤波器组滤波的第二输出信号;
步骤S6、整合所述第一输出信号与所述第二输出信号,经整合后的输出信号表征所述高动态卫星信道多普勒特性。
根据本发明第一方面的方法,所述步骤S2具体包括:
步骤S2-1、基于所述初始载波频率、多普勒频率、多普勒频率变化率,利用公式(1)计算所述射频采样信号的动态载波频率:
Figure GDA0003303934650000021
其中,f表示所述动态载波频率,f0表示所述初始载波频率,fd表示所述多普勒频率,m表示时间,f′d表示所述多普勒频率变化率,*表示点乘;
步骤S2-2、利用公式(2)来表征所述射频采样信号的动态相位:
2πfm+θ0=2πf0(m+τ(m))+θ0 (2)
其中,θ0表示所述初始相位,2πfm+θ0表示基于所述动态载波频率f来表征的所述动态相位,2πf0(m+τ(m))+θ0表征基于所述多普勒模拟动态时延τ(m)来表征的所述动态相位;
步骤S2-3、结合所述公式(1)和所述公式(2),求解所述多普勒模拟动态时延,如公式(3):
Figure GDA0003303934650000031
其中,所述多普勒模拟动态时延的参数频度为毫秒级。
根据本发明第一方面的方法,在所述步骤S3中,获取所述多普勒模拟高频度动态时延具体包括:
步骤S3-1、利用采样脉冲对所述多普勒频率变化率的输入量执行第一累加控制,并将所述第一累加控制的结果发送至所述相位加速度累加器;
步骤S3-2、利用所述采样脉冲对所述相位加速度累加器的输出与所述多普勒频率执行第二累加控制,并将所述第二累加控制的结果发送至所述相位速度累加器;
步骤S3-3、将所述相位速度累加器的输出与所述初始相位相加,以获得所述多普勒模拟高频度动态时延,将所述多普勒模拟高频度动态时延与标称频率控制字相加,以获得所述动态时延控制字;
其中,所述采样脉冲为累加使能脉冲,所述多普勒频率变化率为第一插值参数,所述多普勒频率为第二插值参数,所述初始相位为第三插值参数。
根据本发明第一方面的方法,在所述步骤S4中,所述可变时延滤波器组的多相滤波结构包含若干条具有不同相位的多相支路,每一条支路上均包含带通滤波器和可变时延滤波器,由所述每一条支路上的可变时延滤波器来模拟所述第二动态时延。
根据本发明第一方面的方法,在所述步骤S5中,所述数据缓存区对接收到的所述射频采样信号进行存储,并在转发所述射频采样信号时,生成具有第一动态延时的第一输出信号。
根据本发明第一方面的方法,在所述步骤S6中,对所述第一输出信号和所述第二输出信号进行求和,并对求和后的信号进行下变频处理,以获取所述输出信号,所述输出信号具有实时多普勒特性。
根据本发明第一方面的方法,在所述步骤S6中,对所述第一输出信号和所述第二输出信号进行求和,并对求和后的信号进行下变频处理,以获取所述输出信号,所述输出信号具有实时多普勒特性。
本发明第二方面公开了一种用于实时模拟高动态卫星信道多普勒特性的系统,所述系统包括:
第一处理单元,被配置为,通过提高中频采样信号的采样率,对所述中频采样信号进行预处理,以获取射频采样信号;
第二处理单元,被配置为,基于所述射频采样信号的信号参数和多普勒动态参数,确定多普勒模拟动态时延,其中:
所述信号参数包括多数射频采样信号的初始载波频率和初始相位;所述多普勒动态参数包括多普勒频率和多普勒频率变化率;
第三处理单元,被配置为,利用相位加速度累加器和相位速度累加器,通过累加控制对所述多普勒模拟动态时延执行插值处理,以获取多普勒模拟高频度动态时延,所述多普勒模拟高频度动态时延以动态时延控制字的形式来表征;
第四处理单元,被配置为,将所述动态时延控制字划分为整数部分的第一动态时延和小数部分的第二动态时延,所述第一动态时延通过数据缓存区的存储转发来实现,所述第二动态时延由基于多相滤波结构的可变时延滤波器组来实现;
第五处理单元,被配置为,将所述射频采样信号分别发送至所述数据缓存区和所述可变时延滤波器组,以获得经所述数据缓存区存储转发的第一输出信号和经所述可变时延滤波器组滤波的第二输出信号;
第六处理单元,被配置为,整合所述第一输出信号与所述第二输出信号,经整合后的输出信号表征所述高动态卫星信道多普勒特性。
根据本发明第二方面的系统,所述第二处理单元具体被配置为,采用如下方式确定多普勒模拟动态时延:
基于所述初始载波频率、多普勒频率、多普勒频率变化率,利用公式(1)计算所述射频采样信号的动态载波频率:
Figure GDA0003303934650000051
其中,f表示所述动态载波频率,f0表示所述初始载波频率,fd表示所述多普勒频率,m表示时间,f′d表示所述多普勒频率变化率,*表示点乘;
利用公式(2)来表征所述射频采样信号的动态相位:
2πfm+θ0=2πf0(m+τ(m))+θ0 (2)
其中,θ0表示所述初始相位,2πfm+θ0表示基于所述动态载波频率f来表征的所述动态相位,2πf0(m+τ(m))+θ0表征基于所述多普勒模拟动态时延τ(m)来表征的所述动态相位;
结合所述公式(1)和所述公式(2),求解所述多普勒模拟动态时延,如公式(3):
Figure GDA0003303934650000061
其中,所述多普勒模拟动态时延的参数频度为毫秒级。
根据本发明第二方面的系统,所述第三处理单元具体被配置为,采用如下方式获取所述多普勒模拟高频度动态时延:
利用采样脉冲对所述多普勒频率变化率的输入量执行第一累加控制,并将所述第一累加控制的结果发送至所述相位加速度累加器;
利用所述采样脉冲对所述相位加速度累加器的输出与所述多普勒频率执行第二累加控制,并将所述第二累加控制的结果发送至所述相位速度累加器;
将所述相位速度累加器的输出与所述初始相位相加,以获得所述多普勒模拟高频度动态时延,将所述多普勒模拟高频度动态时延与标称频率控制字相加,以获得所述动态时延控制字;
其中,所述采样脉冲为累加使能脉冲,所述多普勒频率变化率为第一插值参数,所述多普勒频率为第二插值参数,所述初始相位为第三插值参数。
根据本发明第二方面的系统,所述可变时延滤波器组的多相滤波结构包含若干条具有不同相位的多相支路,每一条支路上均包含带通滤波器和可变时延滤波器,由所述每一条支路上的可变时延滤波器来模拟所述第二动态时延。
根据本发明第二方面的系统,所述数据缓存区对接收到的所述射频采样信号进行存储,并在转发所述射频采样信号时,生成具有第一动态延时的第一输出信号。
根据本发明第二方面的系统,所述第五处理单元具体被配置为:对所述可变时延滤波器组的所述每一条支路上的滤波结果进行加权求和,以获取具有第二动态延时的第二输出信号。
根据本发明第二方面的系统,所述第六处理单元具体被配置为:对所述第一输出信号和所述第二输出信号进行求和,并对求和后的信号进行下变频处理,以获取所述输出信号,所述输出信号具有实时多普勒特性。
本发明第三方面公开了一种电子设备。所述电子设备包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时,实现本发明第一方面所述的一种用于实时模拟高动态卫星信道多普勒特性的方法中的步骤。
本发明第四方面公开了一种计算机可读存储介质。所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时,实现本发明第一方面所述的一种用于实时模拟高动态卫星信道多普勒特性的方法中的步骤。
综上,本发明的技术方案通过高速插值获得刷新率较高的多普勒动态参数,以满足卫星信道的近实时动态模拟需求;同时利用基于多相滤波结构的可变时延滤波器,以较低的资源消耗和极少的实时参数更新量,从而实现高动态卫星信道多普勒实时模拟。该技术方案将多普勒特性等效为信号时延的动态变化,通过对信号相位的实时动态控制来实现实时多普勒模拟。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为根据本发明实施例的一种用于实时模拟高动态卫星信道多普勒特性的方法的流程图;
图2为根据本发明实施例的对多普勒模拟动态时延执行插值处理的示意图;
图3为根据本发明实施例的基于多相滤波结构的可变时延滤波器组的结构示意图;
图4为根据本发明实施例的各个支路上的可变时延滤波器的结构示意图;
图5为根据本发明实施例的一种用于实时模拟高动态卫星信道多普勒特性的系统的结构图;
图6为根据本发明实施例的一种电子设备的结构图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例只是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明第一方面公开了一种用于实时模拟高动态卫星信道多普勒特性的方法。图1为根据本发明实施例的一种用于实时模拟高动态卫星信道多普勒特性的方法的流程图;如图1所示,所述方法包括:
步骤S1、通过提高中频采样信号的采样率,对所述中频采样信号进行预处理,以获取射频采样信号;
步骤S2、基于所述射频采样信号的信号参数和多普勒动态参数,确定多普勒模拟动态时延,其中:所述信号参数包括多数射频采样信号的初始载波频率和初始相位;所述多普勒动态参数包括多普勒频率和多普勒频率变化率;
步骤S3、利用相位加速度累加器和相位速度累加器,通过累加控制对所述多普勒模拟动态时延执行插值处理,以获取多普勒模拟高频度动态时延,所述多普勒模拟高频度动态时延以动态时延控制字的形式来表征;
步骤S4、将所述动态时延控制字划分为整数部分的第一动态时延和小数部分的第二动态时延,所述第一动态时延通过数据缓存区的存储转发来实现,所述第二动态时延由基于多相滤波结构的可变时延滤波器组来实现;
步骤S5、将所述射频采样信号分别发送至所述数据缓存区和所述可变时延滤波器组,以获得经所述数据缓存区存储转发的第一输出信号和经所述可变时延滤波器组滤波的第二输出信号;
步骤S6、整合所述第一输出信号与所述第二输出信号,经整合后的输出信号表征所述高动态卫星信道多普勒特性。
在步骤S1中,通过提高中频采样信号的采样率,对所述中频采样信号进行预处理,以获取射频采样信号。
在一些实施例中,对中频采样信号S0(n)进行上变频处理,将中频采样信号的所有频点的采样率提高一倍,从而得到射频采样信号S(n)。
在步骤S2中,基于所述射频采样信号的信号参数和多普勒动态参数,确定多普勒模拟动态时延,其中:所述信号参数包括多数射频采样信号的初始载波频率和初始相位;所述多普勒动态参数包括多普勒频率和多普勒频率变化率。
在一些实施例中,所述步骤S2具体包括:
步骤S2-1、基于所述初始载波频率、多普勒频率、多普勒频率变化率,利用公式(1)计算所述射频采样信号的动态载波频率:
Figure GDA0003303934650000091
其中,f表示所述动态载波频率,f0表示所述初始载波频率,fd表示所述多普勒频率,m表示时间,f′d表示所述多普勒频率变化率,*表示点乘;
步骤S2-2、利用公式(2)来表征所述射频采样信号的动态相位:
2πfm+θ0=2πf0(m+τ(m))+θ0 (2)
其中,θ0表示所述初始相位,2πfm+θ0表示基于所述动态载波频率f来表征的所述动态相位,2πf0(m+τ(m))+θ0表征基于所述多普勒模拟动态时延τ(m)来表征的所述动态相位;
步骤S2-3、结合所述公式(1)和所述公式(2),求解所述多普勒模拟动态时延,如公式(3):
Figure GDA0003303934650000101
其中,所述多普勒模拟动态时延的参数频度为毫秒级。
在步骤S3,利用相位加速度累加器和相位速度累加器,通过累加控制对所述多普勒模拟动态时延执行插值处理,以获取多普勒模拟高频度动态时延,所述多普勒模拟高频度动态时延以动态时延控制字的形式来表征。
在一些实施例中,在所述步骤S3中,获取所述多普勒模拟高频度动态时延具体包括:
步骤S3-1、利用采样脉冲对所述多普勒频率变化率的输入量执行第一累加控制,并将所述第一累加控制的结果发送至所述相位加速度累加器;
步骤S3-2、利用所述采样脉冲对所述相位加速度累加器的输出与所述多普勒频率执行第二累加控制,并将所述第二累加控制的结果发送至所述相位速度累加器;
步骤S3-3、将所述相位速度累加器的输出与所述初始相位相加,以获得所述多普勒模拟高频度动态时延,将所述多普勒模拟高频度动态时延与标称频率控制字相加,以获得所述动态时延控制字;
其中,所述采样脉冲为累加使能脉冲,所述多普勒频率变化率为第一插值参数,所述多普勒频率为第二插值参数,所述初始相位为第三插值参数。
具体如图2所示(图2为根据本发明实施例的对多普勒模拟动态时延执行插值处理的示意图),en_vel_ch为累加使能脉冲(采样脉冲),用于控制累加的速度,累加使能脉冲(采样脉冲)的频率越高,累加速度越快。各个插值参数根据采样脉冲的频率实时更新,插值参数1(第一插值参数)为整个插值间隔内的多普勒频率变化率,插值参数2(第二插值参数)为插值时刻处的多普勒频率,插值参数3(第三插值参数)为插值时刻处的初始相位。通过各个插值参数更新多普勒模拟高频度动态时延TW,当卫星和信号接收机之间不存在相对动态时,多普勒模拟高频度动态时延TW为0。将给定的标称频率控制字TW0与多普勒模拟高频度动态时延TW相加得到动态时延控制字T_ctrl(n)。
在步骤S4,将所述动态时延控制字划分为整数部分的第一动态时延和小数部分的第二动态时延,所述第一动态时延通过数据缓存区的存储转发来实现,所述第二动态时延由基于多相滤波结构的可变时延滤波器组来实现。
具体地,所述动态时延控制字以公式(4)来表示:
T_ctrl(n)=τd(n)+τs(n) (4)
其中,τd(n)为整数部分的第一动态时延,τs(n)为小数部分的第二动态时延。
在一些实施例中,所述可变时延滤波器组的多相滤波结构包含若干条具有不同相位的多相支路,每一条支路上均包含带通滤波器和可变时延滤波器,由所述每一条支路上的可变时延滤波器来模拟所述第二动态时延。
具体如图3所示(图3为根据本发明实施例的基于多相滤波结构的可变时延滤波器组的结构示意图),多相滤波可以按照相位(通过均匀划分)将数字滤波器的系统函数分解成若干个具有不同相位的多相支路(支路0和支路1),在各个支路上对射频采样信号进行滤波,由此可以利用多个阶数较低的滤波器(BPF为带通滤波器,VFD为可变时延录波器)来实现原本阶数较高的滤波器,在而且在资源消耗较低的情况下保证了处理过程中信号的数据率和效果保持不变。
在步骤S5,将所述射频采样信号分别发送至所述数据缓存区和所述可变时延滤波器组,以获得经所述数据缓存区存储转发的第一输出信号和经所述可变时延滤波器组滤波的第二输出信号。
在一些实施例中,在所述步骤S5中,所述数据缓存区对接收到的所述射频采样信号进行存储,并在转发所述射频采样信号时,生成具有第一动态延时的第一输出信号。
具体地,对于整数部分的第一动态时延的时延设计,采用将射频采样信号存储在数据缓存区中,经过整数时延后转发的方式实现;第一输出信号Sd(n)的表达式如公式(5)所示:
sd(n)=S(n+τd(n)) (5)
在一些实施例中,在所述步骤S5中,对所述可变时延滤波器组的所述每一条支路上的滤波结果进行加权求和,以获取具有第二动态延时的第二输出信号。
具体如图4所示(图4为根据本发明实施例的各个支路上的可变时延滤波器的结构示意图),可变时延滤波器组的输出信号y(n)可以用小数部分的第二动态时延τs(n)对N+1条支路上的滤波器加权系数Ci(z),i=O,1,2,...,N进行实时更新,对输入信号x佃)通过加权求和来表示,其中,当i取0时,该条支路上的滤波系数为常数,即,非动态变化量。由此获得第二输出信号Ss(n)的表达式如公式(6)所示:
ss(n)=s(n+τs(n)) (6)
由此,通过更新第一动态时延和第二动态时延能够实现可变时延的控制,可变时延滤波器组的系数根据所需延时量实时更新,从而保证在高动态多普勒变化的条件下,仍然可以对多普勒特性进行高精度模拟。
在步骤S6,整合所述第一输出信号与所述第二输出信号,经整合后的输出信号表征所述高动态卫星信道多普勒特性。
在一些实施例中,对所述第一输出信号和所述第二输出信号进行求和,并对求和后的信号进行下变频处理,以获取所述输出信号,所述输出信号具有实时多普勒特性。
求和后的信号s′(n)的表达式如公式(7)所示:
S′(n)=sd(n)+ss(n) (7)
对求和后的信号(射频信号)s′(n)进行数字下变频处理获得中频信号作为具有实时多普勒特性的输出信号。
综上,本发明的技术方案通过高速插值获得刷新率较高的多普勒动态参数,以满足卫星信道的近实时动态模拟需求;同时利用基于多相滤波结构的可变时延滤波器,以较低的资源消耗和极少的实时参数更新量,从而实现高动态卫星信道多普勒实时模拟。该技术方案将多普勒特性等效为信号时延的动态变化,通过对信号相位的实时动态控制来实现实时多普勒模拟。
本发明第二方面公开了一种用于实时模拟高动态卫星信道多普勒特性的系统。图5为根据本发明实施例的一种用于实时模拟高动态卫星信道多普勒特性的系统的结构图;如图5所示,所述系统500包括:
第一处理单元501,被配置为,通过提高中频采样信号的采样率,对所述中频采样信号进行预处理,以获取射频采样信号;
第二处理单元502,被配置为,基于所述射频采样信号的信号参数和多普勒动态参数,确定多普勒模拟动态时延,其中:
所述信号参数包括多数射频采样信号的初始载波频率和初始相位;所述多普勒动态参数包括多普勒频率和多普勒频率变化率;
第三处理单元503,被配置为,利用相位加速度累加器和相位速度累加器,通过累加控制对所述多普勒模拟动态时延执行插值处理,以获取多普勒模拟高频度动态时延,所述多普勒模拟高频度动态时延以动态时延控制字的形式来表征;
第四处理单元504,被配置为,将所述动态时延控制字划分为整数部分的第一动态时延和小数部分的第二动态时延,所述第一动态时延通过数据缓存区的存储转发来实现,所述第二动态时延由基于多相滤波结构的可变时延滤波器组来实现;
第五处理单元505,被配置为,将所述射频采样信号分别发送至所述数据缓存区和所述可变时延滤波器组,以获得经所述数据缓存区存储转发的第一输出信号和经所述可变时延滤波器组滤波的第二输出信号;
第六处理单元506,被配置为,整合所述第一输出信号与所述第二输出信号,经整合后的输出信号表征所述高动态卫星信道多普勒特性。
根据本发明第二方面的系统,所述第二处理单元502具体被配置为,采用如下方式确定多普勒模拟动态时延:
基于所述初始载波频率、多普勒频率、多普勒频率变化率,利用公式(1)计算所述射频采样信号的动态载波频率:
Figure GDA0003303934650000151
其中,f表示所述动态载波频率,f0表示所述初始载波频率,fd表示所述多普勒频率,m表示时间,f′d表示所述多普勒频率变化率,*表示点乘;
利用公式(2)来表征所述射频采样信号的动态相位:
2πfm+θ0=2πf0(m+τ(m))+θ0 (2)
其中,θ0表示所述初始相位,2πfm+θ0表示基于所述动态载波频率f来表征的所述动态相位,2πf0(m+τ(m))+θ0表征基于所述多普勒模拟动态时延τ(m)来表征的所述动态相位;
结合所述公式(1)和所述公式(2),求解所述多普勒模拟动态时延,如公式(3):
Figure GDA0003303934650000152
其中,所述多普勒模拟动态时延的参数频度为毫秒级。
根据本发明第二方面的系统,所述第三处理单元503具体被配置为,采用如下方式获取所述多普勒模拟高频度动态时延:
利用采样脉冲对所述多普勒频率变化率的输入量执行第一累加控制,并将所述第一累加控制的结果发送至所述相位加速度累加器;
利用所述采样脉冲对所述相位加速度累加器的输出与所述多普勒频率执行第二累加控制,并将所述第二累加控制的结果发送至所述相位速度累加器;
将所述相位速度累加器的输出与所述初始相位相加,以获得所述多普勒模拟高频度动态时延,将所述多普勒模拟高频度动态时延与标称频率控制字相加,以获得所述动态时延控制字;
其中,所述采样脉冲为累加使能脉冲,所述多普勒频率变化率为第一插值参数,所述多普勒频率为第二插值参数,所述初始相位为第三插值参数。
根据本发明第二方面的系统,所述可变时延滤波器组的多相滤波结构包含若干条具有不同相位的多相支路,每一条支路上均包含带通滤波器和可变时延滤波器,由所述每一条支路上的可变时延滤波器来模拟所述第二动态时延。
根据本发明第二方面的系统,所述数据缓存区对接收到的所述射频采样信号进行存储,并在转发所述射频采样信号时,生成具有第一动态延时的第一输出信号。
根据本发明第二方面的系统,所述第五处理单元505具体被配置为:对所述可变时延滤波器组的所述每一条支路上的滤波结果进行加权求和,以获取具有第二动态延时的第二输出信号。
根据本发明第二方面的系统,所述第六处理单元506具体被配置为:对所述第一输出信号和所述第二输出信号进行求和,并对求和后的信号进行下变频处理,以获取所述输出信号,所述输出信号具有实时多普勒特性。
本发明第三方面公开了一种电子设备。所述电子设备包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时,实现本发明第一方面所述的一种用于实时模拟高动态卫星信道多普勒特性的方法中的步骤。
图6为根据本发明实施例的一种电子设备的结构图;如图6所示,电子设备包括通过系统总线连接的处理器、存储器、通信接口、显示屏和输入装置。其中,该电子设备的处理器用于提供计算和控制能力。该电子设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该电子设备的通信接口用于与外部的终端进行有线或无线方式的通信,无线方式可通过WIFI、运营商网络、近场通信(NFC)或其他技术实现。该电子设备的显示屏可以是液晶显示屏或者电子墨水显示屏,该电子设备的输入装置可以是显示屏上覆盖的触摸层,也可以是电子设备外壳上设置的按键、轨迹球或触控板,还可以是外接的键盘、触控板或鼠标等。
本领域技术人员可以理解,图6中示出的结构,仅仅是与本公开的技术方案相关的部分的结构图,并不构成对本申请方案所应用于其上的电子设备的限定,具体的电子设备可以包括比图中所示更多或更少的部件,或者组合某些部件,或者具有不同的部件布置。
本发明第四方面公开了一种计算机可读存储介质。所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时,实现本发明第一方面所述的一种用于实时模拟高动态卫星信道多普勒特性的方法中的步骤。
综上,本发明的技术方案通过高速插值获得刷新率较高的多普勒动态参数,以满足卫星信道的近实时动态模拟需求;同时利用基于多相滤波结构的可变时延滤波器,以较低的资源消耗和极少的实时参数更新量,从而实现高动态卫星信道多普勒实时模拟。该技术方案将多普勒特性等效为信号时延的动态变化,通过对信号相位的实时动态控制来实现实时多普勒模拟。
请注意,以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。

Claims (9)

1.一种用于实时模拟高动态卫星信道多普勒特性的方法,其特征在于,所述方法包括:
步骤S1、通过提高中频采样信号的采样率,对所述中频采样信号进行预处理,以获取射频采样信号;
步骤S2、基于所述射频采样信号的信号参数和多普勒动态参数,确定多普勒模拟动态时延,其中:
所述信号参数包括多数射频采样信号的初始载波频率和初始相位;所述多普勒动态参数包括多普勒频率和多普勒频率变化率;
步骤S3、利用相位加速度累加器和相位速度累加器,通过累加控制对所述多普勒模拟动态时延执行插值处理,以获取多普勒模拟高频度动态时延,所述多普勒模拟高频度动态时延以动态时延控制字的形式来表征;
步骤S4、将所述动态时延控制字划分为整数部分的第一动态时延和小数部分的第二动态时延,所述第一动态时延通过数据缓存区的存储转发来实现,所述第二动态时延由基于多相滤波结构的可变时延滤波器组来实现;
步骤S5、将所述射频采样信号分别发送至所述数据缓存区和所述可变时延滤波器组,以获得经所述数据缓存区存储转发的第一输出信号和经所述可变时延滤波器组滤波的第二输出信号;
步骤S6、整合所述第一输出信号与所述第二输出信号,经整合后的输出信号表征所述高动态卫星信道多普勒特性;
其中,在所述步骤S3中,获取所述多普勒模拟高频度动态时延具体包括:
步骤S3-1、利用采样脉冲对所述多普勒频率变化率的输入量执行第一累加控制,并将所述第一累加控制的结果发送至所述相位加速度累加器;
步骤S3-2、利用所述采样脉冲对所述相位加速度累加器的输出与所述多普勒频率执行第二累加控制,并将所述第二累加控制的结果发送至所述相位速度累加器;
步骤S3-3、将所述相位速度累加器的输出与所述初始相位相加,以获得所述多普勒模拟高频度动态时延,将所述多普勒模拟高频度动态时延与标称频率控制字相加,以获得所述动态时延控制字;
其中,所述采样脉冲为累加使能脉冲,所述多普勒频率变化率为第一插值参数,所述多普勒频率为第二插值参数,所述初始相位为第三插值参数。
2.根据权利要求1所述的一种用于实时模拟高动态卫星信道多普勒特性的方法,其特征在于,所述步骤S2具体包括:
步骤S2-1、基于所述初始载波频率、多普勒频率、多普勒频率变化率,利用公式(1)计算所述射频采样信号的动态载波频率:
Figure FDA0003303934640000021
其中,f表示所述动态载波频率,f0表示所述初始载波频率,fd表示所述多普勒频率,m表示时间,fd表示所述多普勒频率变化率,*表示点乘;
步骤S2-2、利用公式(2)来表征所述射频采样信号的动态相位:
2πfm+θ0=2πf0(m+τ(m))+θ0 (2)
其中,θ0表示所述初始相位,2πfm+θ0表示基于所述动态载波频率f来表征的所述动态相位,2πf0(m+τ(m))+θ0表征基于所述多普勒模拟动态时延τ(m)来表征的所述动态相位;
步骤S2-3、结合所述公式(1)和所述公式(2),求解所述多普勒模拟动态时延,如公式(3):
Figure FDA0003303934640000031
其中,所述多普勒模拟动态时延的参数频度为毫秒级。
3.根据权利要求2所述的一种用于实时模拟高动态卫星信道多普勒特性的方法,其特征在于,在所述步骤S4中,所述可变时延滤波器组的多相滤波结构包含若干条具有不同相位的多相支路,每一条支路上均包含带通滤波器和可变时延滤波器,由所述每一条支路上的可变时延滤波器来模拟所述第二动态时延。
4.根据权利要求3所述的一种用于实时模拟高动态卫星信道多普勒特性的方法,其特征在于,在所述步骤S5中,所述数据缓存区对接收到的所述射频采样信号进行存储,并在转发所述射频采样信号时,生成具有第一动态延时的第一输出信号。
5.根据权利要求4所述的一种用于实时模拟高动态卫星信道多普勒特性的方法,其特征在于,在所述步骤S5中,对所述可变时延滤波器组的所述每一条支路上的滤波结果进行加权求和,以获取具有第二动态延时的第二输出信号。
6.根据权利要求5所述的一种用于实时模拟高动态卫星信道多普勒特性的方法,其特征在于,在所述步骤S6中,对所述第一输出信号和所述第二输出信号进行求和,并对求和后的信号进行下变频处理,以获取所述输出信号,所述输出信号具有实时多普勒特性。
7.一种用于实时模拟高动态卫星信道多普勒特性的系统,其特征在于,所述系统包括:
第一处理单元,被配置为,通过提高中频采样信号的采样率,对所述中频采样信号进行预处理,以获取射频采样信号;
第二处理单元,被配置为,基于所述射频采样信号的信号参数和多普勒动态参数,确定多普勒模拟动态时延,其中:
所述信号参数包括多数射频采样信号的初始载波频率和初始相位;所述多普勒动态参数包括多普勒频率和多普勒频率变化率;
第三处理单元,被配置为,利用相位加速度累加器和相位速度累加器,通过累加控制对所述多普勒模拟动态时延执行插值处理,以获取多普勒模拟高频度动态时延,所述多普勒模拟高频度动态时延以动态时延控制字的形式来表征;
第四处理单元,被配置为,将所述动态时延控制字划分为整数部分的第一动态时延和小数部分的第二动态时延,所述第一动态时延通过数据缓存区的存储转发来实现,所述第二动态时延由基于多相滤波结构的可变时延滤波器组来实现;
第五处理单元,被配置为,将所述射频采样信号分别发送至所述数据缓存区和所述可变时延滤波器组,以获得经所述数据缓存区存储转发的第一输出信号和经所述可变时延滤波器组滤波的第二输出信号;
第六处理单元,被配置为,整合所述第一输出信号与所述第二输出信号,经整合后的输出信号表征所述高动态卫星信道多普勒特性;
其中,所述第三处理单元具体被配置为,以如下方式获取所述多普勒模拟高频度动态时延:
利用采样脉冲对所述多普勒频率变化率的输入量执行第一累加控制,并将所述第一累加控制的结果发送至所述相位加速度累加器;
利用所述采样脉冲对所述相位加速度累加器的输出与所述多普勒频率执行第二累加控制,并将所述第二累加控制的结果发送至所述相位速度累加器;
将所述相位速度累加器的输出与所述初始相位相加,以获得所述多普勒模拟高频度动态时延,将所述多普勒模拟高频度动态时延与标称频率控制字相加,以获得所述动态时延控制字;
其中,所述采样脉冲为累加使能脉冲,所述多普勒频率变化率为第一插值参数,所述多普勒频率为第二插值参数,所述初始相位为第三插值参数。
8.一种电子设备,其特征在于,所述电子设备包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时,实现权利要求1至6中任一项所述的一种用于实时模拟高动态卫星信道多普勒特性的方法中的步骤。
9.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时,实现权利要求1至6中任一项所述的一种用于实时模拟高动态卫星信道多普勒特性的方法中的步骤。
CN202111035878.5A 2021-09-06 2021-09-06 用于实时模拟高动态卫星信道多普勒特性的方法和系统 Active CN113472464B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111035878.5A CN113472464B (zh) 2021-09-06 2021-09-06 用于实时模拟高动态卫星信道多普勒特性的方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111035878.5A CN113472464B (zh) 2021-09-06 2021-09-06 用于实时模拟高动态卫星信道多普勒特性的方法和系统

Publications (2)

Publication Number Publication Date
CN113472464A CN113472464A (zh) 2021-10-01
CN113472464B true CN113472464B (zh) 2021-11-12

Family

ID=77867522

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111035878.5A Active CN113472464B (zh) 2021-09-06 2021-09-06 用于实时模拟高动态卫星信道多普勒特性的方法和系统

Country Status (1)

Country Link
CN (1) CN113472464B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114006668B (zh) * 2021-10-29 2024-02-20 中国人民解放军国防科技大学 卫星信道免系数更新的高精确时延滤波方法和装置
CN114374451B (zh) * 2022-01-10 2023-08-22 中国人民解放军国防科技大学 基于光矩阵交换的大规模多入多出信道模拟方法和装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101261317A (zh) * 2008-04-25 2008-09-10 浙江大学 高动态多模式卫星导航信号源模拟方法及其装置
CN101509968A (zh) * 2009-03-19 2009-08-19 北京理工大学 高动态高精度中频模拟卫星信号产生方法
CN102932072A (zh) * 2012-08-08 2013-02-13 中国传媒大学 高动态卫星信道模拟器
CN108449126A (zh) * 2018-02-08 2018-08-24 中国航天时代电子有限公司 一种多星地面组网远程对接测试系统
CN111064503A (zh) * 2019-12-30 2020-04-24 长光卫星技术有限公司 一种卫星信道高动态时延多普勒模拟系统
CN112311449A (zh) * 2020-11-16 2021-02-02 东南大学 一种低轨卫星通信动态时延和多普勒模拟方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3483632A1 (en) * 2017-11-09 2019-05-15 Septentrio N.V. Method for correcting a pseudorange in a receiver for satellite navigation
CN108111450B (zh) * 2017-11-22 2020-09-22 西南电子技术研究所(中国电子科技集团公司第十研究所) 采用非相干解调模拟目标动态特性的方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101261317A (zh) * 2008-04-25 2008-09-10 浙江大学 高动态多模式卫星导航信号源模拟方法及其装置
CN101509968A (zh) * 2009-03-19 2009-08-19 北京理工大学 高动态高精度中频模拟卫星信号产生方法
CN102932072A (zh) * 2012-08-08 2013-02-13 中国传媒大学 高动态卫星信道模拟器
CN108449126A (zh) * 2018-02-08 2018-08-24 中国航天时代电子有限公司 一种多星地面组网远程对接测试系统
CN111064503A (zh) * 2019-12-30 2020-04-24 长光卫星技术有限公司 一种卫星信道高动态时延多普勒模拟系统
CN112311449A (zh) * 2020-11-16 2021-02-02 东南大学 一种低轨卫星通信动态时延和多普勒模拟方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
GNSS信号软件处理关键技术研究;牟卫华;《中国优秀博士学位论文全文数据库 信息科技辑》;20200215;全文 *
卫星信道模拟器延时技术研究;岳宇航;《中国优秀硕士学位论文全文数据库 信息科技辑》;20180415;全文 *
卫星导航信号模拟器研究与实现;刘志国;《中国优秀硕士学位论文全文数据库 信息科技辑》;20140715;全文 *
基于三状态Markov链的卫星信道模拟器设计与实现;马上等;《系统仿真学报》;20070905(第17期);全文 *
高动态GPS卫星信号模拟器设计与实现;贾科军等;《宇航计测技术》;20061025(第05期);全文 *

Also Published As

Publication number Publication date
CN113472464A (zh) 2021-10-01

Similar Documents

Publication Publication Date Title
CN113472464B (zh) 用于实时模拟高动态卫星信道多普勒特性的方法和系统
CN110704751A (zh) 数据处理方法、装置、电子设备及存储介质
CN109725198B (zh) 一种频谱仪的余辉生成方法及频谱仪
CN110212863B (zh) 校准晶体振荡器的方法、装置、电子设备及存储介质
Rainwater et al. Mixed-resolution ensemble data assimilation
CN112193068A (zh) 一种车速显示方法、装置、车辆及计算机存储介质
CN111580137A (zh) 一种高精度导航接收机射频通道群时特性的拟合方法
CN115795282A (zh) 激波管动态压力重构方法、装置、电子设备及存储介质
Boore et al. The effect of sampling rate and anti-aliasing filters on high-frequency response spectra
CN113238260B (zh) 一种信号参数采集方法、系统、存储介质及电子设备
CN104076202A (zh) 信号分析装置及信号分析方法
KR101041990B1 (ko) 레이더 모의표적장치의 도플러 주파수 구현방법
Ssessanga et al. Complementing regional ground GNSS-STEC computerized ionospheric tomography (CIT) with ionosonde data assimilation
CN114647178A (zh) 基于北斗和地面基准传递的原子钟自动校准方法和系统
CN110146855A (zh) 雷达间断干扰抑制门限计算方法及装置
US20210384892A1 (en) Parallel implementations of frame filters with recursive transfer functions
CN114035203A (zh) 伪距信号的发送方法、装置、存储介质及电子装置
CN113343458A (zh) 发动机传感器的选型方法、装置、电子设备及存储介质
Dillon et al. Asymptotic properties of an autocorrelation coefficient for coherent Doppler sonar
CN107664499B (zh) 一种船用捷联惯导系统的加速度计在线降噪方法
CN115963467B (zh) 调频波形参数的处理方法、装置和计算机设备
CN112859160A (zh) 薄互层砂体厚度的预测方法及装置
CN102495824B (zh) 在同一模拟输出卡上产生多路不同频率模拟信号的方法
CN112996105B (zh) 一种基于同时定位与标定的目标定位方法及装置
CN113447962B (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