CN109001774B - 一种针对发射机通道非理想性的导航信号校正方法 - Google Patents
一种针对发射机通道非理想性的导航信号校正方法 Download PDFInfo
- Publication number
- CN109001774B CN109001774B CN201810523233.8A CN201810523233A CN109001774B CN 109001774 B CN109001774 B CN 109001774B CN 201810523233 A CN201810523233 A CN 201810523233A CN 109001774 B CN109001774 B CN 109001774B
- Authority
- CN
- China
- Prior art keywords
- signal
- filter
- transmitter
- linear
- transmitter channels
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/40—Correcting position, velocity or attitude
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Signal Processing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Transmitters (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
本发明公开了一种针对发射机通道非理想性的导航信号校正方法,解决现有技术不能有效地解决工作在饱和点的HPA的失真补偿、计算复杂度高、资源消耗大、以及不适合星载产品使用的问题。本发明包括以下步骤:步骤1:获取发射机通道的线性和非线性传输特性;步骤2:建立发射机通道非理想性失真补偿仿真模型;步骤3:设计复系数FIR数字基带滤波器;步骤4:在发射机的导航信号基带产生单元中实现步骤3所设计的复系数FIR数字基带滤波器,并评估其效果。本发明方法通过较少阶数的复系数FIR数字滤波器即可实现导航信号的失真校正效果,计算复杂度低,资源消耗小,能适用于发射机非理想性造成的现有的全部导航信号的信号质量校正。
Description
技术领域
本发明涉及发射机非线性特性的校正领域,具体涉及一种针对发射机通道非理想性的导航信号校正方法。
背景技术
全球导航卫星系统(GNSS)已经成为使用范围最广的导航工具,能够为陆地、航空、航海用户提供精确的、可靠的定位导航服务,GNSS信号质量是影响系统关键性能和指标实现的重要因素。随着BOC、MBOC、AltBOC调制等新型调制技术以及Interplex、CASM、POCET等多种恒包络复用技术的使用,单一频点上往往调制着多个分量,占据着更大的带宽。发射机是实现导航信号放大的关键业务单元,受限于目前的工艺水平,发射机通道非理想性——通道幅频响应不平坦、相频响应非线性、高功率放大器(HPA)AM-AM、AM-PM失真——对宽带信号的信号质量影响尤为显著,具体表现为信号带内功率谱失真,相关损失和S曲线过零点偏差偏离正常值,最终影响导航定位性能。
在数字基带信号生成单元中实现导航信号预失真是解决发射机通道非理想性失真最简单易行的手段。但目前大量针对HPA非线性失真的研究均不能有效地解决工作在饱和点的HPA的失真补偿问题,并且这些算法计算复杂度高,资源消耗大,不适合星载产品使用。同时,在HPA非线性失真的作用下,线性补偿的效果也将大大减弱。因此,在不回避HPA非线性失真的情况下,如何改善发射机输出导航信号的质量是一个亟需解决的问题。
发明内容
本发明要解决的技术问题是:提供一种针对发射机通道非理想性的导航信号校正方法,解决现有技术不能有效地解决工作在饱和点的HPA的失真补偿、计算复杂度高、资源消耗大、以及不适合星载产品使用的问题。
为实现上述目的,本发明采用的技术方案如下:
一种针对发射机通道非理想性的导航信号校正方法,包括以下步骤:
步骤1:获取发射机通道的线性和非线性传输特性;
步骤2:根据步骤1获取的发射机通道线性和非线性传输特性建立发射机通道非理想性失真补偿仿真模型;
步骤3:根据步骤2建立的发射机通道非理想性失真补偿仿真模型设计复系数FIR数字基带滤波器;
步骤4:在发射机的导航信号基带产生单元中实现步骤3所设计的复系数FIR数字基带滤波器,并评估其效果。
具体地说,在所述步骤1中,因HPA前置滤波器、HPA、以及多路复用器(OMUX)是发射机通道线性失真和非线性失真的主要来源,因此建立采用矢量网络分析仪对HPA前置滤波器、HPA、以及OMUX这三者的级联链路进行测量的发射机通道传输特性测量系统,以获取发射机通道的线性和非线性传输特性。
更具体地说,在所述步骤1中,通过发射机通道传输特性测量系统来获取发射机通道的线性和非线性传输特性曲线的具体方法为:
获取发射机通道线性传输特性曲线:将矢量网络分析仪置于频率扫描模式,输出激励大小以能使HPA工作在饱和区为准,扫描的频带范围需覆盖信号的发射带宽,获取发射机通道的S21参数,即可得到发射机通道线性传输特性曲线;
获取发射机通道非线性传输特性曲线:矢量网络分析仪置于功率扫描模式,扫描的连续波频率为导航信号中心频点,功率扫描范围选取在工作点附近,获取发射机通道的S21参数,即可得到发射机通道非线性传输特性曲线。
进一步地,在所述步骤2中,所述发射机通道非理想性失真补偿仿真模型包括基带信号生成单元、与所述基带信号生成单元的输出端信号连接的HPA非线性模拟单元、与所述HPA非线性模拟单元的输出端信号连接的发射机通道线性响应模拟单元、分别与所述发射机通道线性响应模拟单元的输出端信号连接的数字滤波器设计单元和信号质量评估单元,所述发射机通道线性响应模拟单元包括与所述HPA非线性模拟单元的输出端信号连接的HPA前置滤波器、以及与所述HPA前置滤波器的输出端信号连接的OMUX,所述信号质量评估单元与所述OMUX的输出端信号连接,用于通过导航信号功率谱、相关损失和S曲线过零点偏差对由所述OMUX输出的最终信号质量进行评估,同时将评估结果反馈至所述数字滤波器设计单元,所述数字滤波器设计单元以所述发射机通道线性响应模拟单元为目标,同时结合所述信号质量评估单元所反馈的评估结果来设计数字滤波器,所述基带信号生成单元包括数字基带信号生成器、由所述数字滤波器设计单元设计并与所述数字基带信号生成器的输出端信号连接的数字滤波器、以及与所述数字滤波器的输出端信号连接并且将卫星发射机数字基带低采样信号转换成信号质量评估所需的高采样信号的采样率转换器,所述HPA非线性模拟单元与所述采样率转换器的输出端信号连接。
具体地说,在所述步骤3中,设计复系数FIR数字基带滤波器的具体步骤为:
步骤3-1:所述矢量网络分析仪在频率扫描模式扫描得到的发射机通道线性传输特性曲线为一个带通滤波器传输特性曲线H0(jf),可表示为其中A0(f)表示幅频响应,表示相频响应,e表示自然底数,j表示虚数;在频域上将带通滤波器的幅频响应曲线和相频响应曲线由导航信号的中心频点线性搬移至0Hz,分别进行平滑处理得到其中,H1(jf)、A1(f)和分别表示频谱线性搬移之后再平滑后得到的通道线性传输特性曲线、幅频响应和相频响应;
步骤3-2:由H1(jf)得到目标滤波器的响应,在信号发射带宽以内,目标滤波器应与H1(jf)幅度互为倒数,并且群时延之和为常数,才能补偿H1(jf)带来的线性失真,因此目标滤波器的响应H2(jf)为:
其中,g为可调固定群时延,Btran为信号发射带宽,选取适当的g值使目标滤波器的群时延在滤波器阶数一半附近;
步骤3-3:设计复系数FIR数字基带滤波器,用数字基带信号的采样率对H2(jf)的频率进行归一化处理,H2(jf)共有K点数据,其频率矩阵为W=[f1;f2;...;fK],则H2(jf)的数据P为P=[H2(jf1);H2(jf2);…;H2(jfK)];N阶FIR滤波器系数A为A=[a(1);a(2);…;a(N)],其系统函数HF(z)可以表示为:
HF(z)=a(1)+a(2)z-1+a(3)z-2+...+a(N)z-(N-1)
HF(z)的傅里叶变换为HF(jf)=a(1)+a(2)e-j2πf+a(3)e-2j2πf+...+a(N)e-(N-1)j2πf,HF(z)在W上的值Q分别为Q=[HF(jf1);HF(jf2);...;HF(jfK)];定义一个矩阵R,使Q=R×A,R可表示为:
则滤波器系数矩阵A应满足:
min||(R×A-P)||
上式是一个经典的最小二乘优化问题,利用MATLAB优化工具箱即可得到最优解;
步骤3-4:将步骤3-3所设计的复系数FIR数字基带滤波器加入步骤2所述仿真模型中,评估经过OMUX之后的导航信号质量,若不满足指标要求,则利用信号质量评估单元通过导航信号功率谱、相关损失和S曲线过零点偏差对由所述OMUX输出的最终信号质量进行评估的评估结果来对目标滤波器的响应进行修正;对幅频响应的修正可通过在幅度上加修正量解决;在相位上乘一个放大因子对相频响应进行有效的修正,得到新的目标滤波器响应;
步骤3-5:重复步骤3-3和步骤3-4的过程,直至设计最后的信号质量满足要求。
进一步地,在所述步骤4中,在发射机的导航信号基带产生单元中实现步骤3所设计的复系数FIR数字基带滤波器,并评估其效果的具体步骤为:
步骤4-1:在发射机的导航信号基带产生单元中实现步骤3所设计的复系数FIR数字基带滤波器;
步骤4-2:采集发射机输出的真实导航信号,评估其信号质量,若不满足指标要求,则利用评估结果对目标滤波器的响应进行修正,得到新的目标滤波器响应;
步骤4-3:重复步骤3-3、步骤3-4、步骤4-1和步骤4-2,直至最后的信号质量满足要求。
与现有技术相比,本发明具有以下有益效果:
(1)本发明方法通过实际提取发射机通道线性和非线性传输特性为依据,通过建立发射机链路模型,设计的基带预失真滤波器能有效改善发射的导航信号质量。
(2)本发明方法通过较少阶数的复系数FIR数字滤波器即可实现导航信号的失真校正效果,计算复杂度低,资源消耗小。
(3)本发明方法不受信号调制方式的约束,能适用于发射机非理想性造成的现有的全部导航信号的信号质量校正。
附图说明
图1为本发明方法的流程图。
图2为本发明实施例中所建立的发射机通道传输特性测量系统框图。
图3为本发明发射机通道非理想性失真补偿仿真模型框图。
图4(1)为本发明实例中矢量网络分析仪测得的发射机通道线性特性曲线图。
图4(2)为本发明实例中矢量网络分析仪测得的发射机通道非线性特性曲线图。
图5为本发明实例仿真模型中得到的某信号S曲线过零点偏差校正效果图。
图6为本发明实例中设计的复系数FIR滤波器的幅频及相频响应图。
图7为本发明实例中某信号实际S曲线过零点偏差校正效果图。
具体实施方式
下面结合附图说明和实施例对本发明作进一步说明,本发明的方式包括但不仅限于以下实施例。
本发明提供的一种针对发射机通道非理想性的导航信号校正方法,该方法的流程如图1所示,其包括以下步骤:
步骤1:获取发射机通道的线性和非线性传输特性;
步骤2:根据步骤1获取的发射机通道线性和非线性传输特性建立发射机通道非理想性失真补偿仿真模型;
步骤3:根据步骤2建立的发射机通道非理想性失真补偿仿真模型设计复系数FIR数字基带滤波器;
步骤4:在发射机的导航信号基带产生单元中实现步骤3所设计的复系数FIR数字基带滤波器,并评估其效果。
其中,在所述步骤1中,因HPA前置滤波器、HPA、以及多路复用器(OMUX)是发射机通道线性失真和非线性失真的主要来源,因此建立如图2所示的发射机通道传输特性测量系统,该系统采用矢量网络分析仪对HPA前置滤波器、HPA、以及OMUX这三者的级联链路进行测量,以获取发射机通道的线性和非线性传输特性。
本发明通过发射机通道传输特性测量系统来获取发射机通道的线性和非线性传输特性曲线的具体方法为:
获取发射机通道线性传输特性曲线:将矢量网络分析仪置于频率扫描模式,输出激励大小以能使HPA工作在饱和区为准,扫描的频带范围需覆盖信号的发射带宽,获取发射机通道的S21参数,即可得到发射机通道线性传输特性曲线;获取发射机通道非线性传输特性曲线:矢量网络分析仪置于功率扫描模式,扫描的连续波频率为导航信号中心频点,功率扫描范围选取在工作点附近,获取发射机通道的S21参数,即可得到发射机通道非线性传输特性曲线。
在所述步骤2中,建立发射机通道非理想性失真补偿仿真模型,模型框图如图3所示。所述发射机通道非理想性失真补偿仿真模型包括基带信号生成单元、与所述基带信号生成单元的输出端信号连接的HPA非线性模拟单元、与所述HPA非线性模拟单元的输出端信号连接的发射机通道线性响应模拟单元、分别与所述发射机通道线性响应模拟单元的输出端信号连接的数字滤波器设计单元和信号质量评估单元,所述发射机通道线性响应模拟单元包括与所述HPA非线性模拟单元的输出端信号连接的HPA前置滤波器、以及与所述HPA前置滤波器的输出端信号连接的OMUX,所述信号质量评估单元与所述OMUX的输出端信号连接,用于通过导航信号功率谱、相关损失、S曲线过零点偏差、信号分量间载波相位关系、信号分量间码相位关系、信号分量功率比、S曲线斜率偏差对由所述OMUX输出的最终信号质量进行评估,同时将评估结果反馈至所述数字滤波器设计单元,所述数字滤波器设计单元以所述发射机通道线性响应模拟单元为目标,同时结合所述信号质量评估单元所反馈的评估结果来设计数字滤波器,所述基带信号生成单元包括数字基带信号生成器、由所述数字滤波器设计单元设计并与所述数字基带信号生成器的输出端信号连接的数字滤波器、以及与所述数字滤波器的输出端信号连接并且将卫星发射机数字基带低采样信号转换成信号质量评估所需的高采样信号的采样率转换器,所述HPA非线性模拟单元与所述采样率转换器的输出端信号连接。
建立发射机通道非理想性失真补偿仿真模型的具体步骤为:
a、建立数字基带信号产生单元。根据信号处理理论,射频信号的线性化处理过程可由基带等效,同时通道非线性失真中HPA的AM-AM、AM-PM失真作用于信号时域波形包络,而时域波形包络由基带波形决定,因此可以用基带信号处理过程来模拟导航射频信号在发射机的实际工作过程。此外,数字基带信号生成器中的采样频率应和卫星发射机基带信号产生单元的工作时钟一致。
b、建立HPA非线性模拟单元。HPA非线性由矢量网络分析仪扫描得到的AM-AM、AM-PM特性曲线模拟,具体过程是将基带信号生成单元输出的信号以其功率为索引对AM-AM、AM-PM特性曲线查表得到其输出信号。
c、建立通道线性响应模拟单元。矢量网络分析仪扫描得到的通道线性传输特性包括HPA前置滤波器和OMUX二者作用之后,其中又以后者的贡献为主,因此将通道线性响应模拟单元置于HPA非线性模拟单元之后。
d、建立信号质量评估单元。此模块用于评估最终输出的信号质量,具体评估的指标包括导航信号功率谱,相关损失,S曲线过零点偏差等。
e、建立数字滤波器设计单元。此模块以通道线性响应为目标,同时结合信号质量评估单元得到的反馈结果设计数字滤波器。
在所述步骤3中,设计复系数FIR数字基带滤波器的具体步骤为:
步骤3-1:所述矢量网络分析仪在频率扫描模式扫描得到的发射机通道线性传输特性曲线为一个带通滤波器传输特性曲线H0(jf),可表示为其中A0(f)表示幅频响应,表示相频响应,e表示自然底数,j表示虚数;在频域上将带通滤波器的幅频响应曲线和相频响应曲线由导航信号的中心频点线性搬移至0Hz,分别进行平滑处理得到其中,H1(jf)、A1(f)和分别表示频谱线性搬移之后再平滑后得到的通道线性传输特性曲线、幅频响应和相频响应;
步骤3-2:由H1(jf)得到目标滤波器的响应,在信号发射带宽以内,目标滤波器应与H1(jf)幅度互为倒数,并且群时延之和为常数,才能补偿H1(jf)带来的线性失真,因此目标滤波器的响应H2(jf)为:
其中,g为可调固定群时延,Btran为信号发射带宽,选取适当的g值使目标滤波器的群时延在滤波器阶数一半附近;
步骤3-3:设计复系数FIR数字基带滤波器,用数字基带信号的采样率对H2(jf)的频率进行归一化处理,H2(jf)共有K点数据,其频率矩阵为W=[f1;f2;...;fK],则H2(jf)的数据P为P=[H2(jf1);H2(jf2);…;H2(jfK)];N阶FIR滤波器系数A为A=[a(1);a(2);…;a(N)],其系统函数HF(z)可以表示为:
HF(z)=a(1)+a(2)z-1+a(3)z-2+...+a(N)z-(N-1)
HF(z)的傅里叶变换为HF(jf)=a(1)+a(2)e-j2πf+a(3)e-2j2πf+...+a(N)e-(N-1)j2πf,HF(z)在W上的值Q分别为Q=[HF(jf1);HF(jf2);...;HF(jfK)];定义一个矩阵R,使Q=R×A,R可表示为:
则滤波器系数矩阵A应满足:
min||(R×A-P)||
上式是一个经典的最小二乘优化问题,利用MATLAB优化工具箱即可得到最优解;
步骤3-4:将步骤3-3所设计的复系数FIR数字基带滤波器加入步骤2所述仿真模型中,评估经过OMUX之后的导航信号质量,若不满足指标要求,则利用信号质量评估单元通过导航信号功率谱、相关损失、S曲线过零点偏差、信号分量间载波相位关系、信号分量间码相位关系、信号分量功率比、S曲线斜率偏差对由所述OMUX输出的最终信号质量进行评估的评估结果来对目标滤波器的响应进行修正;对幅频响应的修正可通过在幅度上加修正量解决;在相位上乘一个放大因子对相频响应进行有效的修正,得到新的目标滤波器响应;
步骤3-5:重复步骤3-3和步骤3-4的过程,直至设计最后的信号质量满足要求。
在所述步骤4中,在发射机的导航信号基带产生单元中实现步骤3所设计的复系数FIR数字基带滤波器,并评估其效果的具体步骤为:
步骤4-1:在发射机的导航信号基带产生单元中实现步骤3所设计的复系数FIR数字基带滤波器;
步骤4-2:采集发射机输出的真实导航信号,评估其信号质量,若不满足指标要求,则利用评估结果对目标滤波器的响应进行修正,得到新的目标滤波器响应;
步骤4-3:重复步骤3-3、步骤3-4、步骤4-1和步骤4-2,直至最后的信号质量满足要求。
为了使本领域技术人员能够更加清楚地理解本发明,特提供以下实例加以说明。
某信号分量在经过某一发射机通道后,其S曲线过零点偏差在1个码片相关器间距下高达1.6ns,此失真将会直接导致此相关器间距下的接收机产生1.6ns的伪距偏差。利用本方法对其进行校正的步骤如下:
第一步,获取发射机通道的线性和非线性传输特性。采用如图2所示的连接关系,利用矢量网络分析仪对发射机通道的测量结果如图4(1)和图4(2)所示。
第二步,建立发射机通道非理想性失真补偿仿真模型。在MATLAB平台上建立如图3所示的发射机通道非理想性失真补偿仿真模型。然后将发射机通道的非线性响应曲线和线性响应曲线分别导入HPA非线性模拟单元和通道线性响应模拟单元。
第三步,设计复系数FIR数字基带滤波器。由实测发射机通道线性响应曲线得到目标滤波器的响应。利用MATLAB优化工具箱得到数字滤波器参数,然后将预失真滤波器加入发射机通道非理想性失真补偿仿真模型中,评估其信号校正效果。经过三次迭代设计,如图5所示,最终将S曲线过零点偏差由1.5ns以上校正至0.1ns以下。三次迭代后设计的滤波器响应如图6所示:
第四步,在发射机中实现数字滤波器,评估其效果。在发射机中实现第三步所设计的数字滤波器,最终的效果如图7所示,S曲线偏差由1.6ns减小至0.2ns以下。
以上仅以S曲线过零点偏差的校正过程为例说明了本发明的实施过程,但是本发明对其它信号体制和导航信号质量指标同样适用。
本发明通过对发射机进行建模,准确表征了发射通道失真情况;通过仿真建立发射通道失真补偿模型,对发射通道带来的导航信号失真进行补偿;通过对真实发射机的实现效果迭代,最终完成对导航信号失真的校正。
上述实施例仅为本发明的优选实施方式之一,不应当用于限制本发明的保护范围,但凡在本发明的主体设计思想和精神上作出的毫无实质意义的改动或润色,其所解决的技术问题仍然与本发明一致的,均应当包含在本发明的保护范围之内。
Claims (4)
1.一种针对发射机通道非理想性的导航信号校正方法,其特征在于,包括以下步骤:
步骤1:获取发射机通道的线性和非线性传输特性;
步骤2:根据步骤1获取的发射机通道线性和非线性传输特性建立发射机通道非理想性失真补偿仿真模型;
步骤3:根据步骤2建立的发射机通道非理想性失真补偿仿真模型设计复系数FIR数字基带滤波器;
步骤4:在发射机的导航信号基带产生单元中实现步骤3所设计的复系数FIR数字基带滤波器,并评估其效果;
在所述步骤1中,因HPA前置滤波器、HPA、以及OMUX是发射机通道线性失真和非线性失真的主要来源,因此建立采用矢量网络分析仪对HPA前置滤波器、HPA、以及OMUX这三者的级联链路进行测量的发射机通道传输特性测量系统,以获取发射机通道的线性和非线性传输特性;通过发射机通道传输特性测量系统来获取发射机通道的线性和非线性传输特性曲线的具体方法为:
获取发射机通道线性传输特性曲线:将矢量网络分析仪置于频率扫描模式,输出激励大小以能使HPA工作在饱和区为准,扫描的频带范围需覆盖信号的发射带宽,获取发射机通道的S21参数,即可得到发射机通道线性传输特性曲线;
获取发射机通道非线性传输特性曲线:矢量网络分析仪置于功率扫描模式,扫描的连续波频率为导航信号中心频点,功率扫描范围选取在工作点附近,获取发射机通道的S21参数,即可得到发射机通道非线性传输特性曲线。
2.根据权利要求1所述的一种针对发射机通道非理想性的导航信号校正方法,其特征在于,在所述步骤2中,所述发射机通道非理想性失真补偿仿真模型包括基带信号生成单元、与所述基带信号生成单元的输出端信号连接的HPA非线性模拟单元、与所述HPA非线性模拟单元的输出端信号连接的发射机通道线性响应模拟单元、分别与所述发射机通道线性响应模拟单元的输出端信号连接的数字滤波器设计单元和信号质量评估单元,所述发射机通道线性响应模拟单元包括与所述HPA非线性模拟单元的输出端信号连接的HPA前置滤波器、以及与所述HPA前置滤波器的输出端信号连接的OMUX,所述信号质量评估单元与所述OMUX的输出端信号连接,用于通过导航信号功率谱、相关损失和S曲线过零点偏差对由所述OMUX输出的最终信号质量进行评估,同时将评估结果反馈至所述数字滤波器设计单元,所述数字滤波器设计单元以所述发射机通道线性响应模拟单元为目标,同时结合所述信号质量评估单元所反馈的评估结果来设计数字滤波器,所述基带信号生成单元包括数字基带信号生成器、由所述数字滤波器设计单元设计并与所述数字基带信号生成器的输出端信号连接的数字滤波器、以及与所述数字滤波器的输出端信号连接并且将卫星发射机数字基带低采样信号转换成信号质量评估所需的高采样信号的采样率转换器,所述HPA非线性模拟单元与所述采样率转换器的输出端信号连接。
3.根据权利要求2所述的一种针对发射机通道非理想性的导航信号校正方法,其特征在于,在所述步骤3中,设计复系数FIR数字基带滤波器的具体步骤为:
步骤3-1:所述矢量网络分析仪在频率扫描模式扫描得到的发射机通道线性传输特性曲线为一个带通滤波器传输特性曲线H0(jf),可表示为其中A0(f)表示幅频响应,表示相频响应,e表示自然底数,j表示虚数;在频域上将带通滤波器的幅频响应曲线和相频响应曲线由导航信号的中心频点线性搬移至0Hz,分别进行平滑处理得到其中,H1(jf)、A1(f)和分别表示频谱线性搬移之后再平滑后得到的通道线性传输特性曲线、幅频响应和相频响应;
步骤3-2:由H1(jf)得到目标滤波器的响应,在信号发射带宽以内,目标滤波器应与H1(jf)幅度互为倒数,并且群时延之和为常数,才能补偿H1(jf)带来的线性失真,因此目标滤波器的响应H2(jf)为:
其中,g为可调固定群时延,Btran为信号发射带宽,选取适当的g值使目标滤波器的群时延在滤波器阶数一半附近;
步骤3-3:设计复系数FIR数字基带滤波器,用数字基带信号的采样率对H2(jf)的频率进行归一化处理,H2(jf)共有K点数据,其频率矩阵为W=[f1;f2;...;fK],则H2(jf)的数据P为P=[H2(jf1);H2(jf2);…;H2(jfK)];N阶FIR滤波器系数A为A=[a(1);a(2);…;a(N)],其系统函数HF(z)可以表示为:
HF(z)=a(1)+a(2)z-1+a(3)z-2+...+a(N)z-(N-1)
HF(z)的傅里叶变换为HF(jf)=a(1)+a(2)e-j2πf+a(3)e-2j2πf+...+a(N)e-(N-1)j2πf,HF(z)在W上的值Q分别为Q=[HF(jf1);HF(jf2);...;HF(jfK)];定义一个矩阵R,使Q=R×A,R可表示为:
则滤波器系数矩阵A应满足:
min||(R×A-P)||
上式是一个经典的最小二乘优化问题,利用MATLAB优化工具箱即可得到最优解;
步骤3-4:将步骤3-3所设计的复系数FIR数字基带滤波器加入步骤2所述仿真模型中,评估经过OMUX之后的导航信号质量,若不满足指标要求,则利用信号质量评估单元通过导航信号功率谱、相关损失和S曲线过零点偏差对由所述OMUX输出的最终信号质量进行评估的评估结果来对目标滤波器的响应进行修正;对幅频响应的修正可通过在幅度上加修正量解决;在相位上乘一个放大因子对相频响应进行有效的修正,得到新的目标滤波器响应;
步骤3-5:重复步骤3-3和步骤3-4的过程,直至设计最后的信号质量满足要求。
4.根据权利要求3所述的一种针对发射机通道非理想性的导航信号校正方法,其特征在于,在所述步骤4中,在发射机的导航信号基带产生单元中实现步骤3所设计的复系数FIR数字基带滤波器,并评估其效果的具体步骤为:
步骤4-1:在发射机的导航信号基带产生单元中实现步骤3所设计的复系数FIR数字基带滤波器;
步骤4-2:采集发射机输出的真实导航信号,评估其信号质量,若不满足指标要求,则利用评估结果对目标滤波器的响应进行修正,得到新的目标滤波器响应;
步骤4-3:重复步骤3-3、步骤3-4、步骤4-1和步骤4-2,直至最后的信号质量满足要求。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810523233.8A CN109001774B (zh) | 2018-05-28 | 2018-05-28 | 一种针对发射机通道非理想性的导航信号校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810523233.8A CN109001774B (zh) | 2018-05-28 | 2018-05-28 | 一种针对发射机通道非理想性的导航信号校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109001774A CN109001774A (zh) | 2018-12-14 |
CN109001774B true CN109001774B (zh) | 2019-07-19 |
Family
ID=64574198
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810523233.8A Active CN109001774B (zh) | 2018-05-28 | 2018-05-28 | 一种针对发射机通道非理想性的导航信号校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109001774B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110417702B (zh) * | 2019-07-23 | 2021-06-15 | 三维通信股份有限公司 | 滤波器系数生成方法、系统和降低信号峰均比的系统 |
CN111190197B (zh) * | 2020-01-16 | 2022-10-28 | 中国人民解放军63921部队 | 一种导航卫星信号质量在轨优化与保持方法 |
CN111504551B (zh) * | 2020-03-10 | 2022-05-20 | 天津大学 | 一种基于最小二乘复指数法的应变式力矩仪带宽扩展方法 |
CN111934694B (zh) * | 2020-07-08 | 2022-03-15 | 电子科技大学 | 一种宽带零中频收发系统失真补偿装置 |
CN112068159B (zh) * | 2020-08-25 | 2024-03-26 | 西安空间无线电技术研究所 | 一种导航卫星信号质量在轨优化方法 |
CN115021836A (zh) * | 2022-05-31 | 2022-09-06 | 哲库科技(北京)有限公司 | 信号补偿方法及装置、频域补偿数据确定方法及装置 |
CN115390104B (zh) * | 2022-10-26 | 2023-04-07 | 中国人民解放军国防科技大学 | 一种导航卫星时延偏差建模方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6118811A (en) * | 1997-07-31 | 2000-09-12 | Raytheon Company | Self-calibrating, self-correcting transceivers and methods |
CN101521651B (zh) * | 2009-04-09 | 2012-01-11 | 哈尔滨工业大学 | 下一代通信系统中宽带卫星通信链路多址信号处理方法 |
CN102893399B (zh) * | 2012-05-24 | 2015-09-09 | 华为技术有限公司 | 预失真校正方法、预失真校正装置、发射机及基站 |
CN103281268B (zh) * | 2013-04-25 | 2016-01-13 | 西安空间无线电技术研究所 | 一种卫星导航信号生成系统带内群时延波动的预补偿方法 |
-
2018
- 2018-05-28 CN CN201810523233.8A patent/CN109001774B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109001774A (zh) | 2018-12-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109001774B (zh) | 一种针对发射机通道非理想性的导航信号校正方法 | |
CN101080911B (zh) | 对高频发射机中模拟i/q调制器幅度和相位失衡以及dc偏移的校准方法及装置 | |
Remley et al. | Millimeter-wave modulated-signal and error-vector-magnitude measurement with uncertainty | |
US6636816B1 (en) | Vector signal analysis method and apparatus therefor | |
CN101868053B (zh) | 零中频信号修正的方法和装置 | |
CN105393452B (zh) | 用于多频带功率放大器中数字预失真器的基带等效伏尔泰拉级数 | |
CN102291154B (zh) | 极坐标发射机 | |
US7593459B1 (en) | Wireless link simulator | |
US8953711B2 (en) | Configurable pre-emphasis component for transmission circuitry | |
CN105591656A (zh) | 一种收发信机的增益平坦度补偿方法 | |
Sombrin | On the formal identity of EVM and NPR measurement methods: Conditions for identity of error vector magnitude and noise power ratio | |
CN109188479A (zh) | 一种高精度卫星导航信号预失真方法 | |
CN103166906A (zh) | 基于最小二乘估计法的极发射器中的非线性频率失真的适应补偿 | |
CN110336572A (zh) | 一种收发信机的增益平坦度补偿方法 | |
CN104954294B (zh) | 一种发射机的支路相位失配检测和校正系统 | |
CN105656434A (zh) | 基于修改分段线性函数的功放数字预失真装置及方法 | |
CN108988962A (zh) | 一种宽带矢量调制信号误差矢量幅度的修正装置与方法 | |
CN105721078B (zh) | 一种通用群时延纹波自动化测试方法 | |
CN105897350A (zh) | 一种发射机芯片的测试方法和装置 | |
CN102904652B (zh) | 一种直放站功放非线性检测方法 | |
CN109583575B (zh) | 基于深度学习提高仪器矢量信号分析性能的处理方法 | |
CN106130573B (zh) | NB-IoT导频信号发生装置的导频信号发生方法 | |
Myoung et al. | Mismatch detection and compensation method for the LINC system using a closed-form expression | |
Prakoso et al. | Antenna representation in two‐port network scattering parameter | |
CN101518014A (zh) | 发射机的模拟i/q调制器的i/q不均衡和dc偏移校准 |
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 |