CN109188479B - 一种高精度卫星导航信号预失真方法 - Google Patents

一种高精度卫星导航信号预失真方法 Download PDF

Info

Publication number
CN109188479B
CN109188479B CN201810987419.9A CN201810987419A CN109188479B CN 109188479 B CN109188479 B CN 109188479B CN 201810987419 A CN201810987419 A CN 201810987419A CN 109188479 B CN109188479 B CN 109188479B
Authority
CN
China
Prior art keywords
signal
filter
inverse characteristic
characteristic curve
predistortion
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
CN201810987419.9A
Other languages
English (en)
Other versions
CN109188479A (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.)
Xian Institute of Space Radio Technology
Original Assignee
Xian Institute of Space Radio 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 Xian Institute of Space Radio Technology filed Critical Xian Institute of Space Radio Technology
Priority to CN201810987419.9A priority Critical patent/CN109188479B/zh
Publication of CN109188479A publication Critical patent/CN109188479A/zh
Application granted granted Critical
Publication of CN109188479B publication Critical patent/CN109188479B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/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/42—Determining position
    • 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/20—Integrity monitoring, fault detection or fault isolation of space segment

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Security & Cryptography (AREA)
  • Amplifiers (AREA)
  • Transmitters (AREA)

Abstract

本发明涉及一种高精度卫星导航信号预失真方法。导航信号生成处理硬件通常由数字部件、DAC(数字模拟转换器)、滤波器、变频器、功率放大器、天线等组成，无线信道中还有大气衰减、干扰等影响，造成接收端导航测量信号失真，传统通信预失真在数字端的反特性结构会破坏信号恒包络特性，导致在功率放大器处附加了新的失真，新的失真无法补偿，这是一个矛盾。本方法通过对通道幅频、群时延特性测量，对导航信号质量失真特性综合评估获取失真特性，综合进行预失真器的迭代设计，完成对导航信号质量的自动校正。

Description

一种高精度卫星导航信号预失真方法
技术领域
本发明涉及一种高精度卫星导航信号预失真方法，主要用于导航信号生成中预失真处理。
背景技术
卫星导航信号的质量直接决定了卫星导航系统的测量和定位精度。而卫星导航信号生成处理的硬件通道特性直接影响卫星导航信号质量。通常卫星导航信号处理过程在卫星上，载荷部分由数字部件、DAC(数字模拟转换器)、滤波器、变频器、功率放大器、天线等组成，无线信道中还有大气衰减、干扰等影响，包括接收端的不理想共同造成了接收端导航测量信号失真。
一般，无线信道和接收端的失真通过一定方法可以消除，还可用于检测大气特性等领域，因此对信道的补偿主要针对卫星载荷端。
传统通信预失真方法主要是在发射端进行信号加重等预失真方法，这些方法会破坏信号恒包络特性，导致在功率放大器处附加了新的非线性失真，这对通信数据的解调是可以忍受的，但是会严重影响导航信号的测量性能，新的非线性失真无法补偿，这成为这种方法的固有缺陷。
发明内容
本发明的技术解决问题是：克服现有技术的不足，提供了一种高精度卫星导航信号预失真方法。通过对通道幅频、群时延特性测量，对导航信号质量失真特性的综合评估获取失真特性，综合进行预失真器的迭代设计，完成对导航信号质量的自动校正。
本发明的技术解决方案是：
一种高精度卫星导航信号预失真方法，步骤如下：
(1)在导航卫星数字端生成扫频数字信号，通过数模转换后接入频谱仪，使用频谱仪采集频谱衰减曲线；所述频谱衰减曲线是指变量为频率，值为幅度的幅频函数曲线。使用频谱仪采集频谱衰减曲线具体为通过使用频谱仪MAX HOLD功能采集频谱衰减曲线。
导航信号生成载荷包括卫星导航数字端和发射通道，在导航卫星数字端生成扫频数字信号，具体为：
设R为数字端寄存器，INC为累加量，每个时钟周期，INC和R都会以下式更新：
INC＝INCold+1
R＝Rold+INC
其中INCold为上个时钟周期的INC，Rold为上个时钟周期的R。
(2)确定预失真滤波器的参数：根据步骤(1)获取的曲线，提取带内部分，获取反特性曲线，并将所述反特性曲线转换到本地数字域，再通过最小二乘法转换为滤波器冲激响应，冲激响应系数即为预失真滤波器的参数；
获取反特性曲线具体为：
设提取带内部分后的曲线为W(ω),则反特性为
MAX(W(ω))-W(ω)
其中，MAX(W(ω))表示求取W(ω)的最大值。
(3)卫星上生成的导航信号经过所述预失真滤波器后，再经过发射通道，对发射通道的输出信号进行功率谱偏差、相关性能、SCB性能指标测试，判断是否满足预设阈值，如果满足，则步骤(2)得到的预失真滤波器的参数作为最终参数；如不满足，进入步骤(4)；
所述发射通道包括DAC、上变频、滤波器、功率放大器以及多工器；
信号经过DAC进行数模转换，生成的模拟信号依次通过上变频以及滤波器滤波之后，再通过功率放大器进行功率放大，最终通过多工器合路输出。
对发射通道的输出信号进行功率谱偏差、相关性能、SCB性能指标测试，具体为：
(A1)按照卫星导航信号空间接口文件生成一个码周期的理想基带采样信号，采样率大于500Msps，按照发射带宽使用矩形滤波器滤波，将滤波完信号数据归一化，归一化公式为：
Figure BDA0001780013000000031
其中xi为归一化前信号数据，yi为归一化后信号数据,N为一个码周期的数据采样点数；
(A2)按照与(A1)相同采样率采样发射通道的输出信号，并正交下变频、按照发射带宽使用矩形滤波器滤波，使用与(A1)相同方法归一化，得到归一化后实际信号数据zi；
(A3)功率谱偏差计算公式为
Figure BDA0001780013000000032
其中FFT(y)为一个码周期的yi的快速傅里叶变换，FFT(z)为一个码周期的zi的快速傅里叶变换，ΔHn(ω)表示第n个信号分量的功率谱偏差，它以频率ω为自变量；
相关性能计算公式为
ΔPn(d)＝J(IFFT(FFT(y)·conj(FFT(z))))
其中IFFT()表示进行反傅里叶变换；J()表示截取出正负一个码片内的部分，并将该部分曲线右半边减去左半边，得到不对称曲线ΔPn(d)，即为相关性能，其中d是采样点；
SCB计算公式为
ΔSn(d)＝SCB(ΔPn(d))
其中，SCB()表示按照非相干鉴相得到S曲线偏差，ΔSn(d)即为SCB，其中d是采样点。
(4)使用网络分析仪测量发射通道内从上变频器到多工器的带内幅相特性，提取带内部分，获取反特性曲线；
(5)将步骤(4)中反特性曲线与步骤(2)中获取的反特性曲线相结合，并将结合后生成的反特性曲线转换到本地数字域，再通过最小二乘法转换为滤波器冲激响应，冲激响应系数即为更新的预失真滤波器的参数；
将步骤(2)中反特性曲线与步骤(4)中获取的反特性曲线相结合，具体为：
(5.1)将两个反特性曲线分别用最小二乘法拟合为两条同样自变量的曲线W1(ω)、W2(ω)；
(5.2)将两者之和W1(ω)+W2(ω)作为结合后的曲线。
(6)卫星上生成的导航信号经过步骤(5)更新后的预失真滤波器后，再经过发射通道，对发射通道的输出信号进行功率谱偏差、相关性能、SCB性能指标测试(与上面步骤计算方法相同)，判断是否满足预设阈值，如果满足，则步骤(5)更新后的预失真滤波器的参数作为最终参数；如不满足，进入步骤(7)；
(7)调整发射通道中的滤波器的带宽，具体为：
(7.1)按照导航信号空间接口文件生成一个码周期的理想基带采样信号，计算信号总功率，公式为
Figure BDA0001780013000000041
其中H(ω)计算公式为
H(ω)＝|FFT(x)|2，
FFT(x)表示对x进行快速傅里叶变换，x为该理想基带采样信号；
(7.2)按如下公式要求求取临界带宽BW：
Figure BDA0001780013000000042
以临界带宽BW/2对信号xi进行低通滤波，然后计算时域峰均比
MAX(xi)/MEAN(xi)
其中MAX(xi)表示xi中的最大值，MEAN(xi)表示xi的平均值；
(7.3)判断时域峰均比是否大于等于0.8，如果大于等于，则BW即为滤波器的带宽；如果小于0.8，则扩大带宽进行低通滤波，再次计算时域峰均比，直到时域峰均比满足要求；
扩大带宽是指：以临界带宽BW/2+△BW对信号xi进行低通滤波，△BW为步进正值，小于1.023MHz，增大△BW从而扩大带宽。
(8)根据步骤(4)中反特性曲线、步骤(2)中获取的反特性曲线、步骤(3)与(6)测得的性能指标测试结果，重新设计预失真滤波器的参数；
具体为：
(8.1)通过公式
ΔHf(ω)＝＝F1(ΔHDAC(ω),ΔHCH(ω),ΔHk,n(ω))F2(ΔPk,n(d,ΔSk,n(d))
计算得到中间变量ΔHf(ω)，其中，
F1(ΔHDAC(ω),ΔHCH(ω),ΔHk,n(ω))＝
ΔHDAC(ω)·ΔHCH(ω)·F11(ΔHCH(ω),ΔHk,n(ω))；
按下式计算F11
Figure BDA0001780013000000051
按下式计算F2
Figure BDA0001780013000000052
其中取Kj和Qi初值为1；ΔHDAC(ω)表示步骤(2)得到反特性曲线；ΔHCH(ω)表示步骤(4)得到的反特性曲线；步骤(3)、(6)、(8)得到功率谱偏差、相关性能、SCB三项指标分别为ΔHk,n(ω)、ΔPk,n(d、ΔSk,n(d)，其中k表示第k次计算得到的结果，n表示第n个信号分量；
(8.2)将ΔHf(ω)依据最小二乘法转换为滤波器冲激响应，冲激响应系数即为重新设计的预失真滤波器的参数。
取Kj和Qi初值为1，迭代计算时，Kj设为功率谱包络偏差最大数值、Qi以SCB的差异最大值取值。
(9)卫星上生成的导航信号经过步骤(8)重新设计后的预失真滤波器后，再经过发射通道，对发射通道的输出信号进行功率谱偏差、相关性能、SCB性能指标测试(与上面步骤计算方法相同)，判断是否满足预设阈值，如果满足，则步骤(8)重新设计后的预失真滤波器的参数作为最终参数；如不满足，返回步骤(8)再次设计预失真滤波器的参数。
本发明与现有技术相比带来的有益效果为：
(1)本发明方法可以补偿功率放大器后多工器失真特性，同时不过多引入功率放大器非线性失真；
(2)本发明以高精度导航信号测量指标为输入参数进行设计，可以满足对通道群时延精度要求0.1ns级别的测量信号传输；
(3)本发明方法为迭代优化设计，步骤(1)～(8)可以在卫星地面测试期间完成，从而使得步骤(9)也可以在轨自主进行，保证卫星长期在轨性能。
附图说明
图1为本发明方法流程图；
图2为导航信号生成结构图；
图3为DAC Sinc特性及其反特性示意图；
图4为功率谱偏差求取示意图。
具体实施方式
导航信号生成处理硬件通常由数字部件、DAC(数字模拟转换器)、滤波器、变频器、功率放大器、天线等组成，无线信道中还有大气衰减、干扰等影响，造成接收端导航测量信号失真，传统通信预失真在数字端的反特性结构会破坏信号恒包络特性，导致在功率放大器处附加了新的失真，新的失真无法补偿，这是一个矛盾。本方法通过对通道幅频、群时延特性测量，对导航信号质量失真特性综合评估获取失真特性，综合进行预失真器的迭代设计，完成对导航信号质量的自动校正。
如图1所示，本发明提出了一种高精度卫星导航信号预失真方法，步骤如下：
(1)进行DAC特性提取。
理想DAC特性符合Sinc衰减曲线，但实际DAC由于线路匹配、非基带模式等原因输出并非完全符合理论曲线。
在导航卫星数字端生成扫频数字信号，通过数模转换后接入频谱仪，使用频谱仪MAX HOLD功能采集频谱衰减曲线；所述频谱衰减曲线是指变量为频率，值为幅度的幅频函数曲线。
导航信号生成载荷包括卫星导航数字端和发射通道，在导航卫星数字端生成扫频数字信号，具体为：
设R为数字端寄存器，INC为累加量，每个时钟周期，INC和R都会以下式更新：
INC＝INCold+1
R＝Rold+INC
其中INCold为上个时钟周期的INC，Rold为上个时钟周期的R。
(2)预失真器初设计1
确定预失真滤波器的参数：根据步骤(1)获取的曲线，提取带内部分，获取反特性曲线，并将所述反特性曲线转换到本地数字域，再通过最小二乘法转换为滤波器冲激响应，冲激响应系数即为预失真滤波器的参数；
获取反特性曲线具体为：
设提取带内部分后的曲线为W(ω),则反特性为
MAX(W(ω))-W(ω)
其中，MAX(W(ω))表示求取W(ω)的最大值。
(3)信号性能指标测试
卫星上生成的导航信号经过所述预失真滤波器后，再经过发射通道，对发射通道的输出信号进行功率谱偏差、相关性能、SCB性能指标测试，判断是否满足预设阈值，如果满足，则步骤(2)得到的预失真滤波器的参数作为最终参数；如不满足，进入步骤(4)；
如图2所示，发射通道包括DAC、上变频、滤波器、功率放大器以及多工器；
信号经过DAC进行数模转换，生成的模拟信号依次通过上变频以及滤波器滤波之后，再通过功率放大器进行功率放大，最终通过多工器合路输出。
对发射通道的输出信号进行功率谱偏差、相关性能、SCB性能指标测试，具体为：
(A1)按照卫星导航信号空间接口文件生成一个码周期的理想基带采样信号，采样率大于500Msps，按照发射带宽使用矩形滤波器滤波，将滤波完信号数据归一化，归一化公式为：
Figure BDA0001780013000000081
其中xi为归一化前信号数据，yi为归一化后信号数据,N为一个码周期的数据采样点数；
(A2)按照与(A1)相同采样率采样发射通道的输出信号，并正交下变频、按照发射带宽使用矩形滤波器滤波，使用与(A1)相同方法归一化，得到归一化后实际信号数据zi；
(A3)功率谱偏差计算公式为
Figure BDA0001780013000000082
其中FFT(y)为一个码周期的yi的快速傅里叶变换，FFT(z)为一个码周期的zi的快速傅里叶变换，ΔHn(ω)表示第n个信号分量的功率谱偏差，它以频率ω为自变量；
相关性能计算公式为
ΔPn(d)＝J(IFFT(FFT(y)·conj(FFT(z))))
其中IFFT()表示进行反傅里叶变换；J()表示截取出正负一个码片内的部分，并将该部分曲线右半边减去左半边，得到不对称曲线ΔPn(d)，即为相关性能，其中d是采样点；
SCB计算公式为
ΔSn(d)＝SCB(ΔPn(d))
其中，SCB()表示按照非相干鉴相得到S曲线偏差，ΔSn(d)即为SCB，其中d是采样点。
(4)使用标准仪器(网络分析仪)测量发射通道内从上变频器到多工器的带内幅相特性，提取带内部分，获取反特性曲线；
(5)预失真器初设计2
将步骤(4)中反特性曲线与步骤(2)中获取的反特性曲线相结合，并将结合后生成的反特性曲线转换到本地数字域，再通过最小二乘法转换为滤波器冲激响应，冲激响应系数即为更新的预失真滤波器的参数；
将步骤(2)中反特性曲线与步骤(4)中获取的反特性曲线相结合，具体为：
(5.1)将两个反特性曲线分别用最小二乘法拟合为两条同样自变量的曲线W1(ω)、W2(ω)；
(5.2)将两者之和W1(ω)+W2(ω)作为结合后的曲线。
(6)再次进行信号性能指标测试
卫星上生成的导航信号经过步骤(5)更新后的预失真滤波器后，再经过发射通道，对发射通道的输出信号进行功率谱偏差、相关性能、SCB性能指标测试(与上面步骤计算方法相同)，判断是否满足预设阈值，如果满足，则步骤(5)更新后的预失真滤波器的参数作为最终参数；如不满足，进入步骤(7)；
(7)调整发射通道中的滤波器的带宽，具体为：
(7.1)按照导航信号空间接口文件生成一个码周期的理想基带采样信号，计算信号总功率，公式为
Figure BDA0001780013000000101
其中H(ω)计算公式为
H(ω)＝|FFT(x)|2，
FFT(x)表示对x进行快速傅里叶变换，x为该理想基带采样信号；
(7.2)按如下公式要求求取临界带宽BW：
Figure BDA0001780013000000102
以临界带宽BW/2对信号xi进行低通滤波，然后计算时域峰均比
MAX(xi)/MEAN(xi)
其中MAX(xi)表示xi中的最大值，MEAN(xi)表示xi的平均值；
(7.3)判断时域峰均比是否大于等于0.8，如果大于等于，则BW即为滤波器的带宽；如果小于0.8，则扩大带宽进行低通滤波，再次计算时域峰均比，直到时域峰均比满足要求；
扩大带宽是指：以临界带宽BW/2+△BW对信号xi进行低通滤波，△BW为步进正值，小于1.023MHz，增大△BW从而扩大带宽。
在上述条件下，限带对导航信号恒包络影响较小，功率放大器非理想特性附加在信号上的影响比较小。
(8)预失真器初设计3
根据步骤(4)中反特性曲线、步骤(2)中获取的反特性曲线、步骤(3)与(6)测得的性能指标测试结果，重新设计预失真滤波器的参数；
具体为：
(8.1)通过公式
ΔHf(ω)＝F1(ΔHDAC(ω),ΔHCH(ω),ΔHk,n(ω))F2(ΔPk,n(d,ΔSk,n(d))
计算得到中间变量ΔHf(ω)，其中，
F1(ΔHDAC(ω),ΔHCH(ω),ΔHk,n(ω))＝ΔHDAC(ω)·ΔHCH(ω)·F11(ΔHCH(ω),ΔHk,n(ω))；
按下式计算F11
Figure BDA0001780013000000111
按下式计算F2
Figure BDA0001780013000000112
其中取Kj和Qi初值为1；ΔHDAC(ω)表示步骤(2)得到反特性曲线；ΔHCH(ω)表示步骤(4)得到的反特性曲线；步骤(3)、(6)、(8)得到功率谱偏差、相关性能、SCB三项指标分别为ΔHk,n(ω)、ΔPk,n(d、ΔSk,n(d)，其中k表示第k次计算得到的结果，n表示第n个信号分量；
(8.2)将ΔHf(ω)依据最小二乘法转换为滤波器冲激响应，冲激响应系数即为重新设计的预失真滤波器的参数。
取Kj和Qi初值为1，迭代计算时，Kj设为功率谱包络偏差最大数值(dB单位)、Qi以SCB的差异最大值取值(ns单位)，直到满足要求。
(9)卫星上生成的导航信号经过步骤(8)重新设计后的预失真滤波器后，再经过发射通道，对发射通道的输出信号进行功率谱偏差、相关性能、SCB性能指标测试(与上面步骤计算方法相同)，判断是否满足预设阈值，如果满足，则步骤(8)重新设计后的预失真滤波器的参数作为最终参数；如不满足，返回步骤(8)再次设计预失真滤波器的参数。
设计举例：
在导航卫星设计中，有一些高阶BOC(二进制偏移副载波)调制信号是由多个分量合成的，如北斗导航卫星的B1Cp信号由B1Cpa和B1Cpb组成，组合信号对通道要求很高，幅频波动要求在0.5dBpp、群时延精度超过1ns。通道通常无法保证这种性能，这时必须使用预失真对通道进行补偿。但是通常的预失真会造成信号恒包络失真，导致后端功率放大器附加非线性失真，继而使得预失真失去作用。
实际导航卫星载荷设计中，可以进行这样设计：如图3，DAC的反特性以其Sinc特性为基础设计(对应步骤(1)、步骤(2))。
如图4，进行功率谱包络偏差求取(对应步骤(3))。相关损失事实上是计算有用信号功率损失情况，在不同相关间距跟踪情况下，使用通用方法获取。SCB曲线是获取不同相关间距下的过零点偏差，使用标准方法获取。如果这三项指标均无法满足指标要求，则进行下面步骤：
如图2，断开数字信号生成、DAC和天线部分的设备，接入网络分析仪，扫描S12参数的幅度和群时延特性，提取带内部分，参照步骤(2)获取反特性。
将DAC反特性和网络分析仪得到的通道反特性结合，设计预失真器，注入并重新测试。如果三项指标仍无法满足指标要求，则进行下面步骤：
ΔHDAC(ω),ΔHCH(ω),ΔHk,n(ω)已经通过测量得到，因此，取Kj和Qi初值为1。按照(8)，F1(ΔHDAC(ω),ΔHCH(ω),ΔHk,n(ω))、F2(ΔPk,n(d,ΔSk,n(d))可以直接求出。将ΔHf(ω)转换为滤波器冲激响应，将冲激响应系数注入预失真器，重新计算三项指标。
这种迭代计算可以将最优预失真滤波器参数计算出来，不会增加过多的非线性失真，同时达到高阶BOC信号高的通道特性要求。实测对通道的补偿，幅度精度可以达到0.1dB，群时延精度可以达到0.1ns。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。

Claims (11)

1.一种高精度卫星导航信号预失真方法，其特征在于步骤如下：
(1)在导航卫星数字端生成扫频数字信号，通过数模转换后接入频谱仪，使用频谱仪采集频谱衰减曲线；
(2)确定预失真滤波器的参数：根据步骤(1)获取的曲线，提取带内部分，获取反特性曲线，并将所述反特性曲线转换到本地数字域，再通过最小二乘法转换为滤波器冲激响应，冲激响应系数即为预失真滤波器的参数；
(3)卫星上生成的导航信号经过所述预失真滤波器后，再经过发射通道，对发射通道的输出信号进行功率谱偏差、相关性能、SCB性能指标测试，判断是否满足预设阈值，如果满足，则步骤(2)得到的预失真滤波器的参数作为最终参数；如不满足，进入步骤(4)；
(4)使用网络分析仪测量发射通道内从上变频器到多工器的带内幅相特性，提取带内部分，获取反特性曲线；
(5)将步骤(4)中反特性曲线与步骤(2)中获取的反特性曲线相结合，并将结合后生成的反特性曲线转换到本地数字域，再通过最小二乘法转换为滤波器冲激响应，冲激响应系数即为更新的预失真滤波器的参数；
(6)卫星上生成的导航信号经过步骤(5)更新后的预失真滤波器后，再经过发射通道，对发射通道的输出信号进行功率谱偏差、相关性能、SCB性能指标测试，判断是否满足预设阈值，如果满足，则步骤(5)更新后的预失真滤波器的参数作为最终参数；如不满足，进入步骤(7)；
(7)调整发射通道中的滤波器的带宽；
(8)根据步骤(4)中反特性曲线、步骤(2)中获取的反特性曲线、步骤(3)与(6)测得的性能指标测试结果，重新设计预失真滤波器的参数；
(9)卫星上生成的导航信号经过步骤(8)重新设计后的预失真滤波器后，再经过发射通道，对发射通道的输出信号进行功率谱偏差、相关性能、SCB性能指标测试，判断是否满足预设阈值，如果满足，则步骤(8)重新设计后的预失真滤波器的参数作为最终参数；如不满足，返回步骤(8)再次设计预失真滤波器的参数。
2.根据权利要求1所述的一种高精度卫星导航信号预失真方法，其特征在于：所述频谱衰减曲线是指变量为频率，值为幅度的幅频函数曲线。
3.根据权利要求1所述的一种高精度卫星导航信号预失真方法，其特征在于：使用频谱仪采集频谱衰减曲线具体为通过使用频谱仪MAX HOLD功能采集频谱衰减曲线。
4.根据权利要求1所述的一种高精度卫星导航信号预失真方法，其特征在于：获取反特性曲线具体为：
设提取带内部分后的曲线为W(ω),则反特性为
MAX(W(ω))-W(ω)
其中，MAX(W(ω))表示求取W(ω)的最大值。
5.根据权利要求1所述的一种高精度卫星导航信号预失真方法，其特征在于：所述发射通道包括DAC、上变频、滤波器、功率放大器以及多工器；
信号经过DAC进行数模转换，生成的模拟信号依次通过上变频以及滤波器滤波之后，再通过功率放大器进行功率放大，最终通过多工器合路输出。
6.根据权利要求1所述的一种高精度卫星导航信号预失真方法，其特征在于：所述将步骤(2)中反特性曲线与步骤(4)中获取的反特性曲线相结合，具体为：
(6.1)将两个反特性曲线分别用最小二乘法拟合为两条同样自变量的曲线W1(ω)、W2(ω)；
(6.2)将两者之和W1(ω)+W2(ω)作为结合后的曲线。
7.根据权利要求1所述的一种高精度卫星导航信号预失真方法，其特征在于：所述步骤(7)调整发射通道中的滤波器的带宽，具体为：
(7.1)按照导航信号空间接口文件生成一个码周期的理想基带采样信号，计算信号总功率，公式为
Figure FDA0001780012990000031
其中H(ω)计算公式为
H(ω)＝|FFT(x)|2，
FFT(x)表示对x进行快速傅里叶变换，x为该理想基带采样信号；
(7.2)按如下公式要求求取临界带宽BW：
Figure FDA0001780012990000032
以临界带宽BW/2对信号xi进行低通滤波，然后计算时域峰均比
MAX(xi)/MEAN(xi)
其中MAX(xi)表示xi中的最大值，MEAN(xi)表示xi的平均值；
(7.3)判断时域峰均比是否大于等于0.8，如果大于等于，则BW即为滤波器的带宽；如果小于0.8，则扩大带宽进行低通滤波，再次计算时域峰均比，直到时域峰均比满足要求；
扩大带宽是指：以临界带宽BW/2+△BW对信号xi进行低通滤波，△BW为步进正值，小于1.023MHz，增大△BW从而扩大带宽。
8.根据权利要求1所述的一种高精度卫星导航信号预失真方法，其特征在于：导航信号生成载荷包括卫星导航数字端和发射通道，在导航卫星数字端生成扫频数字信号，具体为：
设R为数字端寄存器，INC为累加量，每个时钟周期，INC和R都会以下式更新：
INC＝INCold+1
R＝Rold+INC
其中INCold为上个时钟周期的INC，Rold为上个时钟周期的R。
9.根据权利要求1所述的一种高精度卫星导航信号预失真方法，其特征在于：对发射通道的输出信号进行功率谱偏差、相关性能、SCB性能指标测试，具体为：
(A1)按照卫星导航信号空间接口文件生成一个码周期的理想基带采样信号，采样率大于500Msps，按照发射带宽使用矩形滤波器滤波，将滤波完信号数据归一化，归一化公式为：
Figure FDA0001780012990000041
其中xi为归一化前信号数据，yi为归一化后信号数据,N为一个码周期的数据采样点数；
(A2)按照与(A1)相同采样率采样发射通道的输出信号，并正交下变频、按照发射带宽使用矩形滤波器滤波，使用与(A1)相同方法归一化，得到归一化后实际信号数据zi；
(A3)功率谱偏差计算公式为
Figure FDA0001780012990000042
其中FFT(y)为一个码周期的yi的快速傅里叶变换，FFT(z)为一个码周期的zi的快速傅里叶变换，ΔHn(ω)表示第n个信号分量的功率谱偏差，它以频率ω为自变量；
相关性能计算公式为
ΔPn(d)＝J(IFFT(FFT(y)·conj(FFT(z))))
其中IFFT()表示进行反傅里叶变换；J()表示截取出正负一个码片内的部分，并将该部分曲线右半边减去左半边，得到不对称曲线ΔPn(d)，即为相关性能，其中d是采样点；
SCB计算公式为
ΔSn(d)＝SCB(ΔPn(d))
其中，SCB()表示按照非相干鉴相得到S曲线偏差，ΔSn(d)即为SCB，其中d是采样点。
10.根据权利要求1所述的一种高精度卫星导航信号预失真方法，其特征在于：所述根据步骤(4)中反特性曲线、步骤(2)中获取的反特性曲线、步骤(3)与(6)测得的性能指标测试结果，重新设计预失真滤波器的参数，具体为：
(8.1)通过公式
ΔHf(ω)＝F1(ΔHDAC(ω),ΔHCH(ω),ΔHk,n(ω))F2(ΔPk,n(d),ΔSk,n(d))
计算得到中间变量ΔHf(ω)，其中，
F1(ΔHDAC(ω),ΔHGH(ω),ΔHk,n(ω))＝
ΔHDAC(ω)·ΔHCH(ω)·F11(ΔHCH(ω),ΔHk,n(ω))；
按下式计算F11
Figure FDA0001780012990000051
按下式计算F2
Figure FDA0001780012990000052
其中取Kj和Qi初值为1；ΔHDAC(ω)表示步骤(2)得到反特性曲线；ΔHCH(ω)表示步骤(4)得到的反特性曲线；步骤(3)、(6)、(8)得到功率谱偏差、相关性能、SCB三项指标分别为ΔHk,n(ω)、ΔPk,n(d)、ΔSk,n(d)，其中k表示第k次计算得到的结果，n表示第n个信号分量；
(8.2)将ΔHf(ω)依据最小二乘法转换为滤波器冲激响应，冲激响应系数即为重新设计的预失真滤波器的参数。
11.根据权利要求10所述的一种高精度卫星导航信号预失真方法，其特征在于：取Kj和Qi初值为1，迭代计算时，Kj设为功率谱包络偏差最大数值、Qi以SCB的差异最大值取值。
CN201810987419.9A 2018-08-28 2018-08-28 一种高精度卫星导航信号预失真方法 Active CN109188479B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810987419.9A CN109188479B (zh) 2018-08-28 2018-08-28 一种高精度卫星导航信号预失真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810987419.9A CN109188479B (zh) 2018-08-28 2018-08-28 一种高精度卫星导航信号预失真方法

Publications (2)

Publication Number Publication Date
CN109188479A CN109188479A (zh) 2019-01-11
CN109188479B true CN109188479B (zh) 2021-04-13

Family

ID=64916408

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810987419.9A Active CN109188479B (zh) 2018-08-28 2018-08-28 一种高精度卫星导航信号预失真方法

Country Status (1)

Country Link
CN (1) CN109188479B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110208822B (zh) * 2019-05-28 2021-06-11 西安空间无线电技术研究所 一种基于低轨移动通信卫星的通信方法
CN110429996A (zh) * 2019-07-23 2019-11-08 熊军 一种针对群时延和非线性失真的联合多天线校准装置
CN111190197B (zh) * 2020-01-16 2022-10-28 中国人民解放军63921部队 一种导航卫星信号质量在轨优化与保持方法
CN111371722B (zh) * 2020-03-18 2021-12-31 南京创远信息科技有限公司 针对5g nr带内调制信号实现预失真补偿处理的方法
CN112068159B (zh) * 2020-08-25 2024-03-26 西安空间无线电技术研究所 一种导航卫星信号质量在轨优化方法
CN112305569B (zh) * 2020-09-30 2023-08-29 北京空间飞行器总体设计部 一种导航信号时域波形失真度的测试方法
CN113126131B (zh) * 2021-03-25 2022-04-22 中国电子科技集团公司第五十四研究所 一种超低失真导航信号采集及无混叠分离方法
CN113469896B (zh) * 2021-05-21 2022-04-15 贵州师范学院 一种提高地球同步轨道卫星对地观测影像几何校正精度的方法
CN113411140B (zh) * 2021-05-27 2022-06-03 中国人民解放军国防科技大学 信道自适应通道特性学习方法及其模拟器、卫星导航系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001048998A1 (en) * 1999-12-28 2001-07-05 Motorola, Inc. Predistorter controlled by post-amplifier equalization in the transmitter
CN103499812A (zh) * 2013-09-23 2014-01-08 中国科学院电子学研究所 一种宽频多通道相参雷达成像系统基带信号预失真方法
CN105223589A (zh) * 2015-08-31 2016-01-06 无锡伊佩克科技有限公司 一种高精度车载导航信道自适应预失真系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001048998A1 (en) * 1999-12-28 2001-07-05 Motorola, Inc. Predistorter controlled by post-amplifier equalization in the transmitter
CN103499812A (zh) * 2013-09-23 2014-01-08 中国科学院电子学研究所 一种宽频多通道相参雷达成像系统基带信号预失真方法
CN105223589A (zh) * 2015-08-31 2016-01-06 无锡伊佩克科技有限公司 一种高精度车载导航信道自适应预失真系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Dual loop feedback pre-distortion in satellite communication;Chengkai Tang et al.;《Journal of Systems Engineering and Electronics》;20130831;第24卷(第4期);第586-591页 *
基于稀疏表示的北斗导航卫星预失真滤波器设计方法;田野等;《第九届中国卫星导航学术年会》;20180523;第1-5页 *

Also Published As

Publication number Publication date
CN109188479A (zh) 2019-01-11

Similar Documents

Publication Publication Date Title
CN109188479B (zh) 一种高精度卫星导航信号预失真方法
US10998926B2 (en) Method for compensating gain flatness of transceiver
US6211663B1 (en) Baseband time-domain waveform measurement method
US6512788B1 (en) RF output spectrum measurement analyzer and method
US8085175B2 (en) Linearizer
US20020027958A1 (en) Feedback channel signal recovery
US7161511B2 (en) Linearization system and method
CN103281268B (zh) 一种卫星导航信号生成系统带内群时延波动的预补偿方法
US7649930B2 (en) Filter equalization using magnitude measurement data
CN111190197B (zh) 一种导航卫星信号质量在轨优化与保持方法
CN112004263A (zh) 群时延波动校准方法及装置
CN114374446B (zh) 一种基于脉冲信号的幅相特性测量方法
CN109583575B (zh) 基于深度学习提高仪器矢量信号分析性能的处理方法
US20090192738A1 (en) Calibration technique for power amplifiers
CN103346792A (zh) 消除模数转换中时钟抖动的方法、装置及数字预失真方法
Verspecht et al. The vector component analyzer: A new way to characterize distortions of modulated signals in high-frequency active devices
CN111555995B (zh) 一种宽带采集系统的相频响应测量及补偿方法
US10057020B2 (en) Joint estimation of coefficients for skew, gain imbalance and channel response for signal sources
US11646836B2 (en) System and method of coherent averaging of repetitive signals for measurement
CN113328845B (zh) 一种下行载波平坦度补偿方法及装置
CN108366037A (zh) 基于iq通信收发系统的多通道tiadc失真校正方法
CN112130011A (zh) 一种星载功率放大器evm测试方法
CN112653526A (zh) 一种跳频发射机非线性失真测试装置及测试方法
EP1916533A1 (en) Characterization of a frequency response for a frequency translation device
Li et al. Research on Generation and Pre-distortion of high-frequency and wide-band signal based on FPGA/DA

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