CN107864106A - 一种适用于非数据辅助的mpsk载波同步方法 - Google Patents

一种适用于非数据辅助的mpsk载波同步方法 Download PDF

Info

Publication number
CN107864106A
CN107864106A CN201711192653.4A CN201711192653A CN107864106A CN 107864106 A CN107864106 A CN 107864106A CN 201711192653 A CN201711192653 A CN 201711192653A CN 107864106 A CN107864106 A CN 107864106A
Authority
CN
China
Prior art keywords
mrow
centerdot
mfrac
msub
frequency
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
Application number
CN201711192653.4A
Other languages
English (en)
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.)
Chengdu Jiu Jin Technology Co Ltd
Original Assignee
Chengdu Jiu Jin Technology Co Ltd
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 Chengdu Jiu Jin Technology Co Ltd filed Critical Chengdu Jiu Jin Technology Co Ltd
Priority to CN201711192653.4A priority Critical patent/CN107864106A/zh
Publication of CN107864106A publication Critical patent/CN107864106A/zh
Pending legal-status Critical Current

Links

Classifications

    • 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/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2657—Carrier synchronisation
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明涉及数字通信领域。目的是提供一种适用于非数据辅助的MPSK载波同步方法。采用的技术方案是：一种适用于非数据辅助的MPSK载波同步方法，步骤包括：根据输入信号的调制阶数M，将输入信号进行M次方运算；计算瞬时相位和瞬时频率；利用瞬时频率计算频偏、利用瞬时相位计算相偏；并进行频偏及相偏补偿，去掉频率及相位偏差。本发明提出的方法不存在鉴相、跟踪以及环路调整过程，能适应高速率、实时性的信号处理过程，非常适合于数字信号处理器或者软件加以实现；无需前导数据辅助、无反馈，可进行连续解调，计算量小，易于工程实现。

Description

一种适用于非数据辅助的MPSK载波同步方法
技术领域
本发明属于数字通信领域，具体涉及一种适用于非数据辅助的MPSK载波同步方法。
背景技术
载波同步的概念起源于相干解调同相***。在数字通信***中，信号的处理与传输都必须确保在规定的时隙内完成。因此，为了使得信号接收机端能够准确地解调信号，接收机端和发送端必须保持严格意义上的同步。载波同步是PSK信号相干解调中的一项关键技术，它直接关系到解调的指标。
载波同步技术按照采用的辅助手段，分为数据辅助、非数据辅助与编码辅助三个大类。而非数据辅助有锁相环和盲估计两种工作方式。锁相环利用鉴相器提取出同步参数信息，盲估计则利用发送信号调制类型的统计特性，对接收信号进行非线性运算提取出同步参数的估计值。
目前实现PSK信号载波同步的方法主要有：非线性变换--滤波法、同相正交环、逆调制环和判决反馈环，这些方法均是基于反馈的闭环结构；还有基于数据辅助的最大似然参数估计法。
发明内容
本发明的目的是提供一种适用于非数据辅助的MPSK载波同步方法。为实现上述发明目的，本发明所采用的技术方案是：一种适用于非数据辅助的MPSK载波同步方法。
优选的：一种适用于非数据辅助的MPSK载波同步方法，其步骤如下：
(1)根据输入信号的调制阶数M，将输入信号进行M次方运算；
(2)计算瞬时相位和瞬时频率；
(3)利用瞬时频率计算频偏、利用瞬时相位计算相偏；并进行频偏及相偏补偿，去掉频率及相位偏差。
优选的，上述适用于非数据辅助的MPSK载波同步方法，其根据输入信号的调制阶数M，将输入信号进行M次方运算的具体步骤为：
假设接收到的MPSK信号，经过符号同步之后其最佳采样点的表达式为：
其中，f为频差、p为相差、M为调制阶数，fs为采样率，symb(n)为星座映射符号，其与通信原理中MPSK信号对应的星座图对应；
因此，MPSK信号可以将其表示为：
其中，s(n)∈(0，1，...，M-1)
其中，s(n)为信号信息；
通信***中，当发射信号为BPSK时，s(n)为[0,1]伪随机数代表的信息；当为QPSK时，symb(n)为[0,1,2,3]伪随机数代表的信息；当为8PSK时，symb(n)为[0,1,2,3,4,5,6,7]伪随机数代表的信息，以此类推其它调制阶数的MPSK信号；
因此，接收信号完整表达式为：
优选的，上述适用于非数据辅助的MPSK载波同步方法，其计算瞬时相位和瞬时频率的具体步骤为：
对上述接收到的R(n)进行M次方运算：
根据欧拉公式：
ej·6＝cos(θ)+j·sin(θ)
可知：
ej·2πs(n)＝cos(2·π·s(n))+j·sin(2·π·s(n))
由于s(n)为整数，则cos(2·π·s(n))＝1，sin(2·π·s(n))＝0，即得到ej·2·π·s(n)＝1；
因此可以去掉不同符号带来的影响，因此：
对上式进行瞬时相位提取：
其中，angle(RM(n))指取复数RM(n)的相位，如其中actan( )指取反正切；
再进行瞬时频率提取：
因此频差f可由上式推导出：
优选的，上述适用于非数据辅助的MPSK载波同步方法，其利用瞬时频率计算频偏、利用瞬时相位计算相偏；并进行频偏及相偏补偿，去掉频率及相位偏差的具体步骤为：
(1)单次瞬时频率的估计，容易受到噪声影响，误差较大，这里采用多点平均取出其直流分量，可以得到频率偏差的准确估计值：
(2)通过每个瞬时相位计算初始相位：
同样，单次初始相位估计也会噪声影响，同时还会受到频率影响，因此这里频率采用准确估计值，同时也采用多次平均的方式去掉噪声的影响，初始相位估计值为：
(3)根据频偏和相偏，可以对原始接收符号进行补偿，得到：
相应的，上述适用于非数据辅助的MPSK载波同步方法在数字通信***中的应用。
本发明具有以下有益效果：
(1)本发明提出的方法不存在鉴相、跟踪以及环路调整过程，能适应高速率、实时性的信号处理过程，非常适合于数字信号处理器或者软件加以实现。
(2)本发明无需前导数据辅助、无反馈可进行连续解调，计算量小，易于工程实现。
附图说明
图1为PSK类信号载波同步实现框图；
图2为QPSK信号的标准星座图；
图3为8PSK信号的标准星座图；
图4为QPSK信号载波同步前星座图；
图5为QPSK信号载波频偏补偿后星座图；
图6为QPSK信号载波频偏、相偏补偿后星座图；
图7为8PSK信号载波同步前星座图；
图8为8PSK信号载波频偏补偿后星座图；
图9为8PSK信号载波频偏、相偏补偿后星座图。
具体实施方式
下面结合附图，对本发明具体实施步骤进行描述。
实施例一
(1)根据输入信号的调制阶数M，将输入信号进行M次方运算：
假设接收到的MPSK信号，经过符号同步之后其最佳采样点的表达式为：
其中，f为频差、p为相差、M为调制阶数，fs为采样率，symb(n)为星座映射符号，其与通信原理中MPSK信号对应的星座图对应，如图2、图3所示。
因此，MPSK信号可以将其表示为：
其中，s(n)∈(0，1，…，M-1)
其中，s(n)为信号信息。
通信***中，当发射信号为BPSK时，s(n)为[0,1]伪随机数代表的信息；当为QPSK时，symb(n)为[0,1,2,3]伪随机数代表的信息；当为8PSK时，symb(n)为[0,1,2,3,4,5,6,7]伪随机数代表的信息，以此类推其它调制阶数的MPSK信号。
因此，接收信号完整表达式为：
(2)计算瞬时相位和瞬时频率：
对(1)接收到的R(n)进行M次方运算可以得到：
根据欧拉公式：
ej·θ＝cos(θ)+j·sin(θ)
可知：
ej·2πs(n)＝cos(2·π·s(n))+j·sin(2·π·s(n))
由于s(n)为整数，则cos(2·π·s(n))＝1，sin(2·π·s(n))＝0，即得到ej·2·π·s(n)＝1。
因此可以去掉不同符号带来的影响，因此：
对上式进行瞬时相位提取：
其中，angle(RM(n))指取复数RM(n)的相位，如其中actan( )指取反正切；
再进行瞬时频率提取：
因此频差f可由上式推导出：
(3)利用瞬时频率计算频偏、利用瞬时相位计算相偏，并进行频偏及相偏补偿，去掉频率及相位偏差：
1)单次瞬时频率的估计，容易受到噪声影响，误差较大，这里采用多点平均取出其直流分量，可以得到频率偏差的准确估计值：
2)通过每个瞬时相位计算初始相位：
同样，单次初始相位估计也会噪声影响，同时还会受到频率影响，因此这里频率采用准确估计值，同时也采用多次平均的方式去掉噪声的影响，初始相位估计值为：
3)根据频偏和相偏，可以对原始接收符号进行补偿，得到：
实施例二：仿真操作1
(1)仿真设置：采样率8GHz；待载波同步信号1：QPSK，信号载波频偏：3MHz，信号载波相偏：0.2·π，符号率：500MHz，发射滤波器为：升余弦滤波器、滚降因子为0.35，加入高斯加性噪声，信噪比(SNR)为30dB，信号符号数：1500(24000个采样点)。
(2)仿真结果：
1)频偏估计值：3.00014788MHz；相偏估计值：-0.299744·π。
从估计值可以看出，频偏估计值与真实值间相差极小，相对误差为4.93e-5；而相偏估计值-0.299744·π与实际相偏0.2·π之间的差为0.499744·π，极为接近0.5·π，即代表反向补偿相偏使星座图接近QPSK的标准星座图(图2)。
2)仿真图
仿真图如图4、图5、图6所示。
从图4可以看出，频偏和相偏补偿前的星座图为一圆环；如图5所示，补偿频偏后，信号最佳采样点聚合为4个点，但与QPSK信号的标准星座图(图2)相比，矢量点的位置有偏差，表明信号还存在相偏；当将相偏补偿后，图6与QPSK信号的标准星座图(图2)基本一致。
(3)结果分析
从仿真结果可以看出，在未补偿频偏和相偏时，由于PSK信号是恒模的，频偏与相偏使星座图发生了相同幅度的旋转，从图中可以看出QPSK的星座图(图4)为一个圆环；当补偿了频偏之后，只有相偏影响信号的星座图，故QPSK星座图(图5)上的矢量点有一个固定的旋转偏移；最后，当去除了相偏之后就是同步后的QPSK星座图(图6)。
实施例三：仿真操作2
(1)仿真设置：采样率8GHz；待载波同步信号2：8PSK，信号载波频偏：3MHz，信号载波相偏：，符号率：200MHz，发射滤波器为：根升余弦滤波器、滚降因子为0.3，加入高斯加性噪声，信噪比(SNR)为30dB，信号符号数：1500(60000个采样点)。
(2)仿真结果
1)频偏估计值：3.0001931MHz；相偏估计值：-0.201998·π。
从估计值可以看出，频偏估计值与真实值间相差极小，相对误差为6.437e-5；而相偏估计值-0.201998·π与实际相偏0.3·π之间的差为0.501998·π，极为接近0.5·π，即代表反向补偿相偏使星座图接近8PSK的标准星座图(图3)。
2)仿真图
从图7可以看出，频偏和相偏补偿前的星座图为一圆环；如图8所示，补偿频偏后，信号最佳采样点聚合为8个点，但与8PSK信号的标准星座图(图3)相比，矢量点的位置有偏差，表明信号还存在相偏；当将相偏补偿后，图9与8PSK信号的标准星座图(图3)基本一致。
(3)结果分析
从仿真结果可以看出，在未补偿频偏和相偏时，由于PSK信号是恒模的，频偏与相偏使星座图发生了相同幅度的旋转，从图7可以看出8PSK的星座图为一个圆环；当补偿了频偏之后，只有相偏影响信号的星座图，故8PSK星座图(图8)上的矢量点有一个固定的旋转偏移；最后，当去除了相偏之后就是同步后的8PSK星座图(图9)。

Claims (5)

1.一种适用于非数据辅助的MPSK载波同步方法，其特征在于：步骤如下：
(1)根据输入信号的调制阶数M，将输入信号进行M次方运算；
(2)计算瞬时相位和瞬时频率；
(3)利用瞬时频率计算频偏、利用瞬时相位计算相偏；并进行频偏及相偏补偿，去掉频率及相位偏差。
2.根据权利要求1所述的一种适用于非数据辅助的MPSK载波同步方法，其特征在于：其步骤(1)的具体步骤为：
假设接收到的MPSK信号，经过符号同步之后其最佳采样点的表达式为：
<mrow> <mi>R</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>s</mi> <mi>y</mi> <mi>m</mi> <mi>b</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mo>&amp;CenterDot;</mo> <mn>2</mn> <mo>&amp;CenterDot;</mo> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>n</mi> <mo>&amp;CenterDot;</mo> <mi>f</mi> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>+</mo> <mi>p</mi> <mo>)</mo> </mrow> </mrow> </msup> </mrow>
其中，f为频差、p为相差、M为调制阶数，fs为采样率，symb(n)为星座映射符号，其与通信原理中MPSK信号对应的星座图对应；
因此，MPSK信号可以将其表示为：
其中，s(n)∈(0，1，…M-1)
其中，s(n)为信号信息；
通信***中，当发射信号为BPSK时，s(n)为[0,1]伪随机数代表的信息；当为QPSK时，symb(n)为[0,1,2,3]伪随机数代表的信息；当为8PSK时，symb(n)为[0,1,2,3,4,5,6,7]伪随机数代表的信息，以此类推其它调制阶数的MPSK信号；
因此，接收信号完整表达式为：
<mrow> <mi>R</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mo>&amp;CenterDot;</mo> <mn>2</mn> <mo>&amp;CenterDot;</mo> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <mi>s</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> <mi>M</mi> </mfrac> </mrow> </msup> <mo>&amp;CenterDot;</mo> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mo>&amp;CenterDot;</mo> <mn>2</mn> <mo>&amp;CenterDot;</mo> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>n</mi> <mo>&amp;CenterDot;</mo> <mi>f</mi> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>+</mo> <mi>p</mi> <mo>)</mo> </mrow> </mrow> </msup> <mo>.</mo> </mrow>
3.根据权利要求1所述的一种适用于非数据辅助的MPSK载波同步方法，其特征在于：步骤(2)的具体步骤为：
对权利要求2接收到的R(n)进行M次方运算：
<mrow> <mi>R</mi> <mi>M</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mo>&amp;CenterDot;</mo> <mn>2</mn> <mo>&amp;CenterDot;</mo> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <mi>s</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </msup> <mo>&amp;CenterDot;</mo> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mn>2</mn> <mo>&amp;CenterDot;</mo> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <mi>M</mi> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>n</mi> <mo>&amp;CenterDot;</mo> <mi>f</mi> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>+</mo> <mi>p</mi> <mo>)</mo> </mrow> </mrow> </msup> </mrow>
根据欧拉公式：
ej·θ＝cos(θ)+j·sin(θ)
可知：
ej·2·π·s(n)＝cos(2·π·s(n))+j·sin(2·π·s(n))
由于s(n)为整数，则cos(2·π·s(n)＝1，sin(2·π·s(n))＝0，即得到ej·2·π·s(n)＝1；
因此可以去掉不同符号带来的影响，因此：
<mrow> <mi>R</mi> <mi>M</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mn>2</mn> <mo>&amp;CenterDot;</mo> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <mi>M</mi> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>n</mi> <mo>&amp;CenterDot;</mo> <mi>f</mi> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>+</mo> <mi>p</mi> <mo>)</mo> </mrow> </mrow> </msup> </mrow>
对上式进行瞬时相位提取：
<mrow> <mi>P</mi> <mi>h</mi> <mi>a</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>a</mi> <mi>n</mi> <mi>g</mi> <mi>l</mi> <mi>e</mi> <mrow> <mo>(</mo> <mi>R</mi> <mi>M</mi> <mo>(</mo> <mi>n</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>=</mo> <mn>2</mn> <mo>&amp;CenterDot;</mo> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <mi>M</mi> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>n</mi> <mo>&amp;CenterDot;</mo> <mi>f</mi> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>+</mo> <mi>p</mi> <mo>)</mo> </mrow> </mrow>
其中，angle(RM(n))指取复数RM(n)的相位，如其中actan()指取反正切；
再进行瞬时频率提取：
<mrow> <mi>f</mi> <mi>r</mi> <mi>q</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>p</mi> <mi>h</mi> <mi>a</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>p</mi> <mi>h</mi> <mi>a</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>=</mo> <mn>2</mn> <mo>&amp;CenterDot;</mo> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <mi>M</mi> <mo>&amp;CenterDot;</mo> <mfrac> <mi>f</mi> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> </mrow>
因此频差f可由上式推导出：
<mrow> <mi>f</mi> <mo>=</mo> <mfrac> <mrow> <mi>f</mi> <mi>r</mi> <mi>q</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msub> <mi>f</mi> <mi>s</mi> </msub> </mrow> <mrow> <mn>2</mn> <mo>&amp;CenterDot;</mo> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <mi>M</mi> </mrow> </mfrac> <mo>.</mo> </mrow>
4.根据权利要求1所述的一种适用于非数据辅助的MPSK载波同步方法，其特征在于：步骤(3)的具体步骤为：
(1)单次瞬时频率的估计，容易受到噪声影响，误差较大，这里采用多点平均取出其直流分量，可以得到频率偏差的准确估计值：
<mrow> <mover> <mi>f</mi> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <msubsup> <mi>&amp;Sigma;</mi> <mn>0</mn> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mfrac> <mrow> <mi>f</mi> <mi>r</mi> <mi>q</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msub> <mi>f</mi> <mi>s</mi> </msub> </mrow> <mrow> <mn>2</mn> <mo>&amp;CenterDot;</mo> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <mi>M</mi> </mrow> </mfrac> </mrow>
(2)通过每个瞬时相位计算初始相位：
<mrow> <msub> <mi>p</mi> <mrow> <mi>i</mi> <mi>n</mi> <mi>i</mi> <mi>t</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>p</mi> <mi>h</mi> <mi>a</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mn>2</mn> <mo>&amp;CenterDot;</mo> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <mi>M</mi> </mrow> </mfrac> <mo>-</mo> <mfrac> <mrow> <mi>n</mi> <mo>&amp;CenterDot;</mo> <mover> <mi>f</mi> <mo>&amp;OverBar;</mo> </mover> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> </mrow>
同样，单次初始相位估计也会噪声影响，同时还会受到频率影响，因此这里频率采用准确估计值，同时也采用多次平均的方式去掉噪声的影响，初始相位估计值为：
<mrow> <mover> <mi>p</mi> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <msubsup> <mi>&amp;Sigma;</mi> <mn>0</mn> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <msub> <mi>p</mi> <mrow> <mi>i</mi> <mi>n</mi> <mi>i</mi> <mi>t</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow>
(3)根据频偏和相偏，可以对原始接收符号进行补偿，得到：
<mrow> <mover> <mrow> <mi>R</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mi>s</mi> <mi>y</mi> <mi>m</mi> <mi>b</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msup> <mi>e</mi> <mrow> <mi>j</mi> <mo>&amp;CenterDot;</mo> <mn>2</mn> <mo>&amp;CenterDot;</mo> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>n</mi> <mo>&amp;CenterDot;</mo> <mi>f</mi> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>+</mo> <mi>p</mi> <mo>)</mo> </mrow> </mrow> </msup> <mo>&amp;CenterDot;</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>j</mi> <mo>&amp;CenterDot;</mo> <mn>2</mn> <mo>&amp;CenterDot;</mo> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>n</mi> <mo>&amp;CenterDot;</mo> <mover> <mi>f</mi> <mo>&amp;OverBar;</mo> </mover> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>+</mo> <mover> <mi>p</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> </mrow> </msup> <mo>&amp;ap;</mo> <mi>s</mi> <mi>y</mi> <mi>m</mi> <mi>b</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
5.权利要求1所述的适用于非数据辅助的MPSK载波同步方法在数字通信***中的应用。
CN201711192653.4A 2017-11-24 2017-11-24 一种适用于非数据辅助的mpsk载波同步方法 Pending CN107864106A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711192653.4A CN107864106A (zh) 2017-11-24 2017-11-24 一种适用于非数据辅助的mpsk载波同步方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711192653.4A CN107864106A (zh) 2017-11-24 2017-11-24 一种适用于非数据辅助的mpsk载波同步方法

Publications (1)

Publication Number Publication Date
CN107864106A true CN107864106A (zh) 2018-03-30

Family

ID=61703394

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711192653.4A Pending CN107864106A (zh) 2017-11-24 2017-11-24 一种适用于非数据辅助的mpsk载波同步方法

Country Status (1)

Country Link
CN (1) CN107864106A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110505171A (zh) * 2019-07-12 2019-11-26 四川安迪科技实业有限公司 一种无数据辅助频偏估计方法及装置
CN112073351A (zh) * 2020-11-16 2020-12-11 湖南国科锐承电子科技有限公司 一种针对mpsk***的新型载波频偏估计方法
CN114978827A (zh) * 2022-04-22 2022-08-30 深圳市人工智能与机器人研究院 一种基于星座图相位异众比率校正频偏的调制识别方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101626357A (zh) * 2009-09-22 2010-01-13 北京理工大学 一种基于最大似然估计的mpsk***载波同步方法
CN102469060A (zh) * 2011-06-30 2012-05-23 重庆重邮信科通信技术有限公司 一种ofdm***同步估计方法及装置
US20140037023A1 (en) * 2012-07-31 2014-02-06 Wuhan Research Institute Of Posts And Telecommunications Transmitting system and receiving system, and transmitting method and receiving method
CN103856421A (zh) * 2014-03-04 2014-06-11 中国人民解放军理工大学 开闭环结合的mpsk信号非数据辅助信噪比估计方法
CN106534010A (zh) * 2016-11-29 2017-03-22 重庆邮电大学 一种oqpsk信号解调的相位偏差估计与补偿方法
CN106992950A (zh) * 2016-01-20 2017-07-28 晨星半导体股份有限公司 载波频偏估测装置与载波频偏估测方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101626357A (zh) * 2009-09-22 2010-01-13 北京理工大学 一种基于最大似然估计的mpsk***载波同步方法
CN102469060A (zh) * 2011-06-30 2012-05-23 重庆重邮信科通信技术有限公司 一种ofdm***同步估计方法及装置
US20140037023A1 (en) * 2012-07-31 2014-02-06 Wuhan Research Institute Of Posts And Telecommunications Transmitting system and receiving system, and transmitting method and receiving method
CN103856421A (zh) * 2014-03-04 2014-06-11 中国人民解放军理工大学 开闭环结合的mpsk信号非数据辅助信噪比估计方法
CN106992950A (zh) * 2016-01-20 2017-07-28 晨星半导体股份有限公司 载波频偏估测装置与载波频偏估测方法
CN106534010A (zh) * 2016-11-29 2017-03-22 重庆邮电大学 一种oqpsk信号解调的相位偏差估计与补偿方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110505171A (zh) * 2019-07-12 2019-11-26 四川安迪科技实业有限公司 一种无数据辅助频偏估计方法及装置
CN110505171B (zh) * 2019-07-12 2022-03-25 四川安迪科技实业有限公司 一种无数据辅助频偏估计方法及装置
CN112073351A (zh) * 2020-11-16 2020-12-11 湖南国科锐承电子科技有限公司 一种针对mpsk***的新型载波频偏估计方法
CN114978827A (zh) * 2022-04-22 2022-08-30 深圳市人工智能与机器人研究院 一种基于星座图相位异众比率校正频偏的调制识别方法
CN114978827B (zh) * 2022-04-22 2024-05-10 深圳市人工智能与机器人研究院 一种基于星座图相位异众比率校正频偏的调制识别方法

Similar Documents

Publication Publication Date Title
US6535549B1 (en) Method and apparatus for carrier phase tracking
US6560294B1 (en) Phase estimation in carrier recovery of phase-modulated signals such as QAM signals
CN110300079B (zh) 一种msk信号相干解调方法及***
CN107864106A (zh) 一种适用于非数据辅助的mpsk载波同步方法
CN103023831B (zh) 一种适用于突发波形的载波频偏估计方法
CN109831402A (zh) 16apsk信号载波相位同步及其判锁定方法
CN105871765A (zh) 一种基于fft辅助s-pll的无线通信载波跟踪方法
CN111800364B (zh) 基于波形匹配的编码cpm信号频偏估计和校正方法
CN110278169A (zh) 一种时域频偏估计算法
EP1751880B1 (en) Cellular communications system using baseband carrier injection and related methods
CN110266625B (zh) 一种基于帧同步码隐含于信息码的qpsk调制的动态门限判决的解调方法及***
CN114039829B (zh) 一种基于导频辅助的相位恢复方法及***
CN107579942A (zh) 一种apsk调制信号的解调方法
CN106508104B (zh) 一种扩展遥测相干接收机频偏估计范围的方法
CN104935543B (zh) 用于相位调制信号的同步解调器电子电路
Sun et al. Carrier frequency offset synchronization algorithm for short burst communication system
CN110730149B (zh) 一种平移式特殊qpsk***联合捕获方法
CN111917679A (zh) 在低信噪比条件下载波及符号定时同步方法
CN103391269B (zh) 8qam接收机及其提取本地载波相位误差的方法
CN113411279B (zh) 一种基于dvb-s2***的可复用q次方解映射方法及***
CN105049172A (zh) 一种短突发的载波频偏估计方法
CN106411802B (zh) 一种基于无人机数据链传输的频偏校正方法
CN108023847B (zh) 调制符号间相位交错bpsk调制方法
CN101102301A (zh) 一种应用于通信***的正交幅度调制解调方法和装置
CN109391572A (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20180330

RJ01 Rejection of invention patent application after publication