CN116346209A - 一种mpsk卫星通信系统载波同步方法 - Google Patents
一种mpsk卫星通信系统载波同步方法 Download PDFInfo
- Publication number
- CN116346209A CN116346209A CN202310338192.6A CN202310338192A CN116346209A CN 116346209 A CN116346209 A CN 116346209A CN 202310338192 A CN202310338192 A CN 202310338192A CN 116346209 A CN116346209 A CN 116346209A
- Authority
- CN
- China
- Prior art keywords
- phase
- mpsk
- carrier
- unsigned
- value
- 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.)
- Pending
Links
- 102100026758 Serine/threonine-protein kinase 16 Human genes 0.000 title claims abstract description 59
- 101710184778 Serine/threonine-protein kinase 16 Proteins 0.000 title claims abstract description 59
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000004891 communication Methods 0.000 title claims abstract description 25
- 238000001228 spectrum Methods 0.000 claims description 17
- 238000005070 sampling Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000011084 recovery Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000005311 autocorrelation function Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910002056 binary alloy Inorganic materials 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/14—Relay systems
- H04B7/15—Active relay systems
- H04B7/185—Space-based or airborne stations; Stations for satellite systems
- H04B7/1851—Systems using a satellite or space-based relay
- H04B7/18513—Transmission in a satellite or space-based system
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/18—Phase-modulated carrier systems, i.e. using phase-shift keying
- H04L27/20—Modulator circuits; Transmitter circuits
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
- H04L2027/0024—Carrier regulation at the receiver end
- H04L2027/0026—Correction of carrier offset
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Aviation & Aerospace Engineering (AREA)
- General Physics & Mathematics (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
Description
技术领域
本发明属于卫星通信技术领域,具体涉及一种MPSK卫星通信系统载波同步方法,特别是应用于多普勒或干扰等引起的大频偏信道下的高精度频偏校正及相位补偿方法。
背景技术
MPSK调制信号可以有效地利用有限的频带资源,在卫星无线通信中得到广泛应用。在低地球轨道卫星通信系统中或接收机载体高速运动条件下,由于多普勒效应的影响,接收到的信号存在较大的多普勒频偏,甚至达到符号速率的十几倍。由于多普勒频移和收发双方振荡器不稳定等因素的影响,使得接收机接收信号的载波与本地载波存在较大动态范围的频偏,而载波恢复是通信信号接收的一项关键技术,其恢复效果将直接影响接收机信号接收的性能。载波恢复的首要问题是载波频率估计,若载波频率估计不准确,就会产生频率漂移,影响信号恢复效果。在载波频偏的影响被消除后,接下来进行载波相位估计。载波相位估计的目的就是对残留频偏估计误差以及相偏进行纠正,保证系统以较低的误码率恢复出原始信号。
目前,根据是否使用辅助序列可以将载波频偏估计算法分为数据辅助类和非数据辅助类两大类。数据辅助类方法主要是基于一种前导序码或者训练序列来帮助进行频偏估计。该方法性能较好,但是由于在数据辅助估计算法中训练序列占用了一定的资源,从而导致系统的传输效率降低。非数据辅助类方法是利用信道传输信号本身比如循环前缀等来进行盲估计,或者也可以进行简单的直接判决法。后者是利用解调后传输符号的速率来检测相位或者载波的频偏误差,以此进行对系统同步的估计。
非数据辅助频偏估计方法又可以分为两类:基于时域的估计方法和基于频域的估计方法。基于时域估计方法的频偏估计器硬件消耗相对较小,它有几种不同的实现结构,可以归结为两种形式:基于相邻接收信号差分和基于自相关函数的方法。主要有延时共轭相乘算法、Fitz算法等。这些算法均能在一定的信噪比下接近克拉美罗界,缺点是难以同时兼顾频偏估计精度和频偏的估计范围。基于频域估计方法则以寻找周期图的最大值为基础,主要有基于FFT的频偏估计算法、基于CZT的频偏估计算法等。以上这几类方法的频偏估计精度和频偏估计范围不高。此外,这几类方法均忽略了对载波相位估计处理,无法做到载波频偏和载波相位联合估计处理。
发明内容
本发明解决的技术问题是:克服现有技术的不足,本发明提出一种MPSK卫星通信系统载波同步方法,用于实现高动态环境下卫星通信的高精度载波同步。
本发明解决技术的方案是:一种MPSK卫星通信系统载波同步方法,该方法包括如下步骤:
获取两段MPSK数字基带信号,对每一段MPSK数字基带信号分别做M次方运算,得到两段无符号调制数字信号,M为MPSK数字基带信号的调制阶数;
计算两段无符号调制数字信号的无模糊相位估计值,分别记为第一无模糊相位估计值ψ1和第二无模糊相位估计值ψ2;
根据第一无模糊相位估计值ψ1和第二无模糊相位估计值ψ2与MPSK数字基带信号载波频偏、载波相位的关系,建立二元一次方程组,解算二元一次方程组得到MPSK数字基带信号的载波频偏精估计值和载波相位精估计值/>
优选地,所述二元一次方程组为:
其中,N1为第一段MPSK数字基带信号的采样点数,N2为第二段MPSK数字基带信号的采样点数,N1≠N2;T为采样周期。
优选地,无符号调制数字信号的无模糊相位估计值通过如下步骤计算得到:
对无符号调制数字信号进行快速傅里叶变换,得到无符号调制数字信号频谱;
取无符号调制数字信号频谱中最大值所对应的载波频率,作为无符号调制数字信号的载波频偏粗估计值;
对无符号调制数字信号的载波频偏的粗估计值进行累积,得到无符号调制数字信号的载波相位粗估计值;
对无符号调制数字信号频谱中最大值所对应的相位进行补偿,得到无符号调制数字信号的模糊相位值;
根据无符号调制数字信号的载波频偏粗估计值、模糊相位值和载波相位粗估计值,进行解相位模糊处理,获得无符号调制数字信号的无模糊相位估计值。
其中,n0为频谱中最大值所对应的索引值,NFFT为FFT处理点数。
其中,N为MPSK数字基带信号的采样点数。
优选地,无符号调制数字信号的模糊相位值ψ为:
ψ=πM(N-1)ΔfT+Mθ
其中,θ为MPSK数字基带信号的载波相位。
优选地,所述无模糊相位估计值的具体计算步骤如下:
根据载波频偏估计值和载波相位粗估计值,计算出无模糊相位参考值ψ';
优选地,所述无模糊相位参考值ψ'的计算公式如下:
优选地,所述相位模糊数kam的计算公式如下:
其中,[]表示向近取整运算。
优选地,所述解模糊后的无模糊相位值ψ”的计算公式如下:
本发明与现有技术相比的有益效果是:
(1)、本发明通过综合分析频谱峰值点的幅度信息和相位信息获取载波频偏的粗估计值,再利用该载波频偏粗估计值获得载波相位的粗估计值。
(2)、本发明提出了基于频谱峰值点的载波频偏值辅助解相位模糊的方法,可准确求解模糊数,并获得无模糊的峰值点相位信息。
(3)、本发明可同时获得载波频偏和相位的高精度估计值,且易于实现。
(4)、本发明可实现高动态环境下卫星通信的高精度载波同步。
附图说明
图1为本发明实施例所提方法处理流程图;
图2为本发明实施例解相位模糊方法处理流程图。
具体实施方式
下面结合实施例对本发明作进一步阐述。
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述。
本发明提供了一种MPSK卫星通信系统载波同步方法,该方法包括如下步骤:
获取两段MPSK数字基带信号,对每一段MPSK数字基带信号分别做M次方运算,得到两段无符号调制数字信号,M为MPSK数字基带信号的调制阶数;
计算两段无符号调制数字信号的无模糊相位估计值,分别记为第一无模糊相位估计值ψ1和第二无模糊相位估计值ψ2;
根据第一无模糊相位估计值ψ1和第二无模糊相位估计值ψ2与MPSK数字基带信号载波频偏、载波相位的关系,建立二元一次方程组,解算二元一次方程组得到MPSK数字基带信号的载波频偏精估计值和载波相位精估计值/>
所述二元一次方程组为:
其中,N1为第一段MPSK数字基带信号的采样点数,N2为第二段MPSK数字基带信号的采样点数,N1≠N2;T为采样周期。
无符号调制数字信号的无模糊相位估计值通过如下步骤计算得到:
对无符号调制数字信号进行快速傅里叶变换,得到无符号调制数字信号频谱;
取无符号调制数字信号频谱中最大值所对应的载波频率,作为无符号调制数字信号的载波频偏粗估计值;
对无符号调制数字信号的载波频偏的粗估计值进行累积,得到无符号调制数字信号的载波相位粗估计值;
对无符号调制数字信号频谱中最大值所对应的相位进行补偿,得到无符号调制数字信号的模糊相位值;
根据无符号调制数字信号的载波频偏粗估计值、模糊相位值和载波相位粗估计值,进行解相位模糊处理,获得无符号调制数字信号的无模糊相位估计值。实施例:
本发明实施例提供一种MPSK卫星通信系统载波同步方法。结合图1,本发明所提方法的具体包括如下步骤:
S1、获取第一段MPSK数字基带信号;
不考虑符号同步误差,MPSK数字基带信号模型可以表示为:
x(k)=akexp(j2πΔfkT+jθ) (1)
需要说明的是,本发明所提方法适用于采用相位调制类型的通信系统,其具体特征为可通过对接收信号作M次方使得调制相位恒定为1,进而消除该项的影响。此处统一采用MPSK调制模型描述不同调制阶数的接收信号。具体地,当M=2时,ak={1,-1},为BPSK信号;当M=4时,为QPSK信号。
S2、对MPSK数字基带信号做M次方运算,得到无符号调制数字信号;
例如,对BPSK信号(M=2)做平方运算,对QPSK信号(M=4)做4次方运算。其M次方运算后的结果可以表示为:
y(k)=xM(k)=exp(j2πMΔfkT+jMθ) (2)
S3、对无符号调制数字信号进行快速傅里叶变换(FFT),得到无符号调制数字信号的频谱;
对步骤S2中的结果进行FFT处理,其结果可以表示为(n=0,1,,NFFT-1,NFFT为FFT处理点数,NFFT大于等于N,通常情况下,NFFT取2的整数次幂:
S4、根据无符号调制数字信号的频谱中最大值所对应的载波频率,计算无符号调制数字信号的载波频偏粗估计值;
本步骤对FFT处理结果的频谱取最大值,并记其索引值为n0。根据该索引值可得到载波频偏的粗估计值,其可以表示为:
S5、对无符号调制数字信号的频谱中最大值所对应的相位进行补偿,得到无符号调制数字信号的模糊相位值;
频谱最大值处的相位提取值为:
其中,angle{·}表示取相位运算。
提取的相位值中的第二项为已知量,与待估计参数没有关系,因此可对其补偿,补偿的相位项可以表示为:
补偿后的提取相位可以表示为:
ψ=πM(N-1)ΔfT+Mθ (7)
S6、根据无符号调制数字信号的载波频偏的粗估计值,计算得到无符号调制数字信号的载波相位粗估计值;
在接收信号中由于载波频偏所引起的相位变化可以表示为:
提取接收信号中的相位值,再结合载波频偏所引起的相位变化,可以获得载波相位的粗估计值。实际中为提高估计精度可对N个数据点的估计值取均值。因此,载波相位的粗估计值可以表示为:
其中angle{·}表示取相位运算。
S7、根据无符号调制数字信号的无符号调制数字信号的载波频偏粗估计值、模糊相位值和载波相位粗估计值,进行解相位模糊处理,获得无符号调制数字信号的无模糊相位估计值;
具体地,结合图2,解相位模糊处理的实施步骤如下:
步骤S7.1、根据载波频偏估计值和载波相位粗估计值,计算出无模糊相位参考值ψ';
无模糊相位参考值ψ'可通过下式获得:
相位模糊数kam可按下式求得:
其中[·]表示取整运算。
S8、获取第二段MPSK数字基带信号,重复步骤S2~步骤S7;
S9、将两段MPSK数字基带信号处理得到的第一无模糊相位估计值ψ1和第二无模糊相位估计值ψ2代入无模糊相位与载波频偏、载波相位的关系式中,建立二元一次方程组,解算二元一次方程组,得到MPSK数字基带信号的载波频偏精估计值和载波相位精估计值/>
实际中,可采用两帧的接收数据,或者将一帧数据分成两段,从而得到第一段MPSK数字基带信号和第二段MPSK数字基带信号。假设两段接收数据的长度分别为N1和N2(注意N1≠N2),对这两段数据分别采用上述的处理方法,获得无模糊的相位值,将其记为ψ1和ψ2,则ψ1和ψ2满足以下方程:
本发明虽然已以较佳实施例公开如上,但其并不是用来限定本发明,任何本领域技术人员在不脱离本发明的精神和范围内,都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改,因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰,均属于本发明技术方案的保护范围。
Claims (10)
3.根据权利要求1所述的一种MPSK卫星通信系统载波同步方法,其特征在于无符号调制数字信号的无模糊相位估计值通过如下步骤计算得到:
对无符号调制数字信号进行快速傅里叶变换,得到无符号调制数字信号频谱;
取无符号调制数字信号频谱中最大值所对应的载波频率,作为无符号调制数字信号的载波频偏粗估计值;
对无符号调制数字信号的载波频偏的粗估计值进行累积,得到无符号调制数字信号的载波相位粗估计值;
对无符号调制数字信号频谱中最大值所对应的相位进行补偿,得到无符号调制数字信号的模糊相位值;
根据无符号调制数字信号的载波频偏粗估计值、模糊相位值和载波相位粗估计值,进行解相位模糊处理,获得无符号调制数字信号的无模糊相位估计值。
6.根据权利要求3所述的一种MPSK卫星通信系统载波同步方法,其特征在于无符号调制数字信号的模糊相位值ψ为:
ψ=πM(N-1)ΔfT+Mθ
其中,θ为MPSK数字基带信号的载波相位。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310338192.6A CN116346209A (zh) | 2023-03-31 | 2023-03-31 | 一种mpsk卫星通信系统载波同步方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310338192.6A CN116346209A (zh) | 2023-03-31 | 2023-03-31 | 一种mpsk卫星通信系统载波同步方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116346209A true CN116346209A (zh) | 2023-06-27 |
Family
ID=86878611
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310338192.6A Pending CN116346209A (zh) | 2023-03-31 | 2023-03-31 | 一种mpsk卫星通信系统载波同步方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116346209A (zh) |
-
2023
- 2023-03-31 CN CN202310338192.6A patent/CN116346209A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7466768B2 (en) | IQ imbalance compensation | |
US10778496B2 (en) | OFDM system synchronization tracking method and apparatus | |
CN102546500B (zh) | 基于导频和软信息联合辅助的soqpsk载波同步方法 | |
CN111683034B (zh) | 一种基于ofdm的大多普勒无线通信时频同步方法 | |
CN101626357B (zh) | 一种基于最大似然估计的mpsk系统载波同步方法 | |
US20040190663A1 (en) | Method and system for synchronizing in a frequency shift keying receiver | |
CN110300079B (zh) | 一种msk信号相干解调方法及系统 | |
CN104022981A (zh) | 一种正交幅度调制信号的盲载波频偏估计方法 | |
CN103281280B (zh) | 基于旋转平均周期图和解调软信息的载波同步方法 | |
US12081377B2 (en) | Demodulating modulated signals | |
CN107342960B (zh) | 一种适合幅度相移键控的非数据辅助频偏估计方法 | |
CN113726716B (zh) | 一种基于判决反馈的载波相位误差鉴别方法 | |
CN110278169A (zh) | 一种时域频偏估计算法 | |
US20140146859A1 (en) | Efficient frequency estimation | |
CN111800364B (zh) | 基于波形匹配的编码cpm信号频偏估计和校正方法 | |
CN105610755B (zh) | 一种突发信号的频偏估计方法及装置 | |
CN105681245B (zh) | 一种基于连续分段滑窗估计频偏的频率跟踪方法 | |
US7266162B2 (en) | Carrier frequency offset estimator for OFDM systems | |
CN109714286B (zh) | 一种Pi/8 D8PSK解调的载波频偏估计方法 | |
CN112838998B (zh) | 一种适用于单载波的iq不平衡估计方法和系统 | |
US8675744B1 (en) | Channel tracking in an orthogonal frequency-division multiplexing system | |
CN105933265B (zh) | 一种对qam信号进行非数据辅助的相位噪声盲估计方法 | |
CN105721075B (zh) | 数据辅助同频复用卫星干扰信号强度测量方法及装置 | |
CN108964783B (zh) | 大频偏条件下相干光接收机载波同步方法 | |
CN116545824A (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 |