CN103916341A - 升余弦成型mpsk信号的盲符号速率估计与定时方法 - Google Patents
升余弦成型mpsk信号的盲符号速率估计与定时方法 Download PDFInfo
- Publication number
- CN103916341A CN103916341A CN201410152732.2A CN201410152732A CN103916341A CN 103916341 A CN103916341 A CN 103916341A CN 201410152732 A CN201410152732 A CN 201410152732A CN 103916341 A CN103916341 A CN 103916341A
- Authority
- CN
- China
- Prior art keywords
- symbol rate
- timing
- value
- signal
- timing synchronization
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 18
- 102100026758 Serine/threonine-protein kinase 16 Human genes 0.000 title abstract description 10
- 101710184778 Serine/threonine-protein kinase 16 Proteins 0.000 title abstract description 10
- 230000003595 spectral effect Effects 0.000 claims abstract description 13
- 230000001360 synchronised effect Effects 0.000 claims abstract description 7
- 238000000605 extraction Methods 0.000 claims description 26
- 238000005070 sampling Methods 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 4
- YBJHBAHKTGYVGT-ZKWXMUAHSA-N (+)-Biotin Chemical compound N1C(=O)N[C@@H]2[C@H](CCCCC(=O)O)SC[C@@H]21 YBJHBAHKTGYVGT-ZKWXMUAHSA-N 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- FEPMHVLSLDOMQC-UHFFFAOYSA-N virginiamycin-S1 Natural products CC1OC(=O)C(C=2C=CC=CC=2)NC(=O)C2CC(=O)CCN2C(=O)C(CC=2C=CC=CC=2)N(C)C(=O)C2CCCN2C(=O)C(CC)NC(=O)C1NC(=O)C1=NC=CC=C1O FEPMHVLSLDOMQC-UHFFFAOYSA-N 0.000 claims description 3
- 238000000465 moulding Methods 0.000 claims 5
- 230000010354 integration Effects 0.000 claims 1
- 238000000926 separation method Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000001568 sexual effect Effects 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
- Synchronisation In Digital Transmission Systems (AREA)
Abstract
本发明公开了一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法。该方法包括：计算输入基带信号的功率谱密度；根据所述的功率谱密度确定符号速率范围；根据定时同步环路能同步的最大定时频率偏差和所述符号速率范围确定不同的初始符号速率估计值；依次利用所述初始符号速率估计值进行定时同步，并将能使定时同步环路同步的初始符号速率估计值加上此时的定时频率偏差作为符号速率值。本发明的特点在于：（1）具有较好实时性；（2）符号速率估计与定时同步都不受载波频率偏差的影响；（3）在进行盲符号速率估计的同时完成了符号定时同步。本方法也适用于根升余弦成型的MPSK信号。
Description
技术领域
本发明涉及数字通信中一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法。
背景技术
多进制相位键控（M-ary Phase Shift Keying，MPSK）调制信号具有良好的抗干扰性能和工程上易于实现等优点，在航天航空领域的飞行器遥测、遥控、遥感、数传、通信等领域得到广泛应用。而在非协作的通信信号分析中，符号速率作为一个解调的重要参数，需要准确估计。符号速率的准确估计，对于定时同步、调制方式的识别、特定信号搜索以及解调等均有重要意义。
目前，常见的实现PSK信号盲符号速率估计的几种代表性算法主要有：（1）基于循环谱的载频和码速率估计算法，该算法需要较大的采样数据量，不能满足工程上实时性的要求。（2）基于小波变换的突变点和码速率估计算法，该方法存在最优尺度和母函数的选取问题。（3）基于最大似然估计的PSK信号参数估计算法，该算法需要多维搜索，实时性仍较差。
发明内容
鉴于上述现有技术所存在的问题，本发明的目的是提供一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法。特点在于：（1）具有较好实时性；（2）符号速率估计与定时同步都不受载波频率偏差的影响；（3）在进行符号速率估计的同时完成了符号定时同步。本方法也适用于根升余弦成型的MPSK信号。
本发明的目的通过以下技术方案实现的：
一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法，其特征在于，包括：
A、根据输入的基带信号求信号功率谱密度；
B、根据所述的功率谱密度求出符号速率范围；
C、根据定时同步环路的最大定时频率偏差，在所述符号速率范围内，等间隔的确定n个符号速率值作为初始符号速率值；
D、依次选取所述n个初始符号速率值，计算定时同步环路之前的抽取系数，然后对抽取后的信号进行定时同步。当定时同步后，输出定时同步信号和符号速率。
所述步骤A进一步包括：
A1、设置采样点数Nsample，采样频率fs；
A2、设置FFT点数NFFT；
A3、设置窗函数类型与窗函数长度Nwindow；
A4、设置重叠数据长度Noverlap；
A5、设置抽取模块1的抽取系数为S1=1；
A6、计算基带信号的功率谱密度Sy(f)。
所述步骤B进一步包括：
B1、搜索Sy(f)找到其最大值，其值记为v；
B2、在第一奈奎斯特频率范围(0～0.5fs)内搜索满足Sy(f)≥0.5v的最大频率点，记为fL；
B3、在第二奈奎斯特频率范围(0.5fs～fs)内搜索满足Sy(f)≥0.5v的最小频率点，记为fH；
B4、将(fs-fH+fL)作为基带信号带宽的估计值
B5、根据基带信号带宽的估计值如果满足则设置S1=S1×2，并返回步骤A6；如果满足则设置S1=S1÷2，并返回步骤A6；若都不满足，则执行步骤B6；
B6、根据信号带宽估计值求得信号速率估计范围为
所述步骤C进一步包括：
C1、确定定时同步环路支持的最大定时频率偏差ΔfTMAX；
C2、根据所述符号速率估计范围和最大定时频率偏差ΔfTMAX，确定n个长度为2ΔfTMAX的频率区间，每个区间中心频率值为将其作为定时同步环路的初始符号速率。
所述步骤D进一步包括：
D1、设置i=1，设置每个区间的定时检测时间为tp；
D2、选取初始符号速率估计值计算抽取模块2的抽取因子S2；
D3、等待tp时间后，检测环路滤波器的积分支路，若其值不收敛，令i=i+1，返回步骤D2；若其值收敛，则表示定时已同步，输出符号定时同步信号和相对定时频率偏差ΔfT，进入步骤D4；
D4、将此时的初始符号速率估计值和定时频率偏差ΔfT相加的结果作为基带信号的符号速率估计值，并输出该值。
由上述本发明提供的技术方案可以看出，本发明原理简单，计算量小，具有较好实时性，利于工程实现，在进行符号速率估计的同时实现了符号定时同步，并且符号速率估计与定时同步都不受载波频率偏差干扰；通过对实际采集的信号分析表明，由该方法对符号速率进行估计时，具有优越的性能。
附图说明
图1为本发明的系统框图
图2为本发明的符号速率范围分割图
图3为本发明的处理流程图
具体实施方式
本发明提供了一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法。本发明的核心为：对接收到的基带信号进行功率谱密度估计，通过其功率谱密度估计出符号速率范围，根据定时同步环路能同步的最大定时频率偏差在符号速率范围中确定不同的符号速率估计值，将能使定时同步环路同步的符号速率估计值加上相对定时频率偏差作为符号速率输出，同时输出定时同步信号。
本发明所述针对升余弦成型MPSK信号的盲符号速率估计与定时装置的结构示意图如图1所示，其中包括：
1、抽取模块1：用于对基带信号进行抽取，其抽取系数受到后面信号带宽估计模块的控制。
2、功率谱估计模块：用于估计接收到基带信号的功率谱密度，并且将估计出来的功率谱密度送入后面的信号带宽估计模块；本发明中设置FFT的点数为256点，设置的窗函数类型为海明窗，窗函数长度为255，参与计算的采样点数为20000，重叠数据的长度128。
3、信号带宽估计模块：用于估计信号的带宽，并根据信号的带宽控制抽取模块1的抽取系数。
设输入的功率谱密度为Sy(f)，则搜索其最大值，记为v，然后在第一奈奎斯特频率范围(0～0.5fs)内搜索满足Sy(f)≥0.5v的最大频率点，记为fL；在第二奈奎斯特频率范围(0.5fs～fs)内搜索满足Sy(f)≥0.5v的最小频率点，记为fH，则信号的带宽估计值当输入信号带宽估计值满足则设置抽取模块1的抽取系数S1=S1×2，当输入信号带宽估计值满足则设置抽取模块1的抽取系数S1=S1÷2，反之则输出信号带宽估计值
4、符号速率范围估计模块：用于估计信号的速率范围，根据输入的信号带宽估计值符号的速率范围为
5、符号速率估计值选取模块：用于从符号速率范围中确定符号速率估计值，同时根据选取的符号速率估计值计算抽取模块2的抽取系数。
确定符号速率估计值的标准如图2所示，为从最小符号速率开始，下一个符号速率估计值距离前一个符号速率估计值的频率间隔为定时同步环路的最大定时频率偏差的两倍。
6、抽取模块2：用于对基带信号进行抽取，使采样率满足所选符号速率估计值的16倍过采样，其抽取系数受到符号速率估计值选取模块的控制。
7、定时同步模块：用于对信号进行定时同步，同时输出定时同步信号和矫正后的符号速率估计值作为符号速率输出。
本发明所述方法的处理流程图如图3所示，包括如下步骤：
步骤31：设置采样点数Nsample，采样频率fs。
步骤32：设置FFT点数NFFT。
步骤33：设置窗函数类型为海明窗与窗函数长度Nwindow。
步骤34：设置重叠数据长度Noverlap。
步骤35：设置抽取模块1的抽取系数为S1=1。
步骤36：计算基带信号的功率谱密度Sy(f)。
步骤37：搜索Sy(f)找到其最大值，其值记为v。
步骤38：在第一奈奎斯特频率范围(0～0.5fs)内搜索满足Sy(f)≥0.5v的最大频率点，记为fL。
步骤39：在第二奈奎斯特频率范围(0.5fs～fs)内搜索满足Sy(f)≥0.5v的最小频率点，记为fH。
步骤310：将(fs-fH+fL)作为基带信号带宽的估计值
步骤311：根据基带信号带宽的估计值如果满足则执行步骤312；如果满足则执行步骤313；若都不满足，则执行步骤314。
步骤312：设置抽取模块1抽取系数S1=S1×2，执行步骤36。
步骤313：设置抽取模块1抽取系数S1=S1÷2，执行步骤36。
步骤314：根据信号带宽估计值求得信号速率估计范围为
步骤315：选择定时同步环路的最大定时频率偏差ΔfTMAX。
步骤316：根据所述符号速率估计范围和最大定时频率偏差ΔfTMAX，确定n个长度为2ΔfTMAX的频率区间，此时n的值为每个区间中心频率值为将其作为定时同步环路的初始符号速率。
步骤317：设置i=1，设置每个区间的定时检测时间为tp。
步骤318：选取初始符号速率估计值计算抽取模块2的抽取因子S2。
步骤319：等待tp时间后，检测环路滤波器的积分支路，若其值不收敛，令i=i+1，执行步骤320；若其值收敛，则表示定时已同步，输出符号定时同步信号和相对定时频率偏差ΔfT，进入步骤321。
步骤320：令i=i+1，选取下一个初始符号速率估计值。
步骤321：输出定时同步信号和相对定时频率偏差ΔfT。
步骤322：输出基带信号的速率估计值
Claims (5)
1.一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法，其特征在于，包括：
A、根据输入的基带信号求信号功率谱密度；
B、根据所述的功率谱密度求出符号速率范围；
C、根据定时同步环路的最大定时频率偏差，在所述符号速率范围内，等间隔的确定n个符号速率值作为初始符号速率值；
D、依次选取所述n个初始符号速率值，计算定时同步环路之前的抽取系数，然后对抽取后的信号进行定时同步。当定时同步后，输出定时同步信号和符号速率。
2.根据权利要求1所述的一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法，其特征在于，所述步骤A进一步包括：
A1、设置采样点数Nsample，采样频率fs；
A2、设置FFT点数NFFT；
A3、设置窗函数类型与窗函数长度Nwindow；
A4、设置重叠数据长度Noverlap；
A5、设置抽取模块1的抽取系数为S1=1；
A6、计算基带信号的功率谱密度Sy(f)。
3.根据权利要求1所述的一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法，其特征在于，所述步骤B进一步包括：
B1、搜索Sy(f)找到其最大值，其值记为v；
B2、在第一奈奎斯特频率范围（0～0.5fs）内搜索满足Sy(f)≥0.5v的最大频率点，记为fL；
B3、在第二奈奎斯特频率范围（0.5fs～fs）内搜索满足Sy(f)≥0.5v的最小频率点，记为fH；
B4、将(fs-fH+fL)作为基带信号带宽的估计值
B5、根据基带信号带宽的估计值如果满足则设置S1=S1×2，并返回步骤A6；如果满足则设置S1=S1÷2，并返回步骤A6；若都不满足，则执行步骤B6；
B6、根据信号带宽估计值求得信号速率估计范围为
4.根据权利要求1所述的一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法，其特征在于，所述步骤C进一步包括：
C1、确定定时同步环路支持的最大定时频率偏差ΔfTMAX；
C2、根据所述符号速率估计范围和最大定时频率偏差ΔfTMAX，确定n个长度为2ΔfTMAX的频率区间，每个区间中心频率值为将其作为定时同步环路的初始符号速率。
5.根据权利要求1所述的一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法，其特征在于，所述步骤D进一步包括：
D1、设置i=1，设置每个区间的定时检测时间为tp；
D2、选取初始符号速率估计值计算抽取模块2的抽取因子S2；
D3、等待tp时间后，检测环路滤波器的积分支路，若其值不收敛，令i=i+1，返回步骤D2；若其值收敛，则表示定时已同步，输出符号定时同步信号和相对定时频率偏差ΔfT，进入步骤D4；
D4、将此时的初始符号速率估计值和定时频率偏差ΔfT相加的结果作为基带信号的符号速率估计值，并输出该值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410152732.2A CN103916341B (zh) | 2014-04-16 | 2014-04-16 | 一种针对升余弦成型mpsk信号的盲符号速率估计与定时方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410152732.2A CN103916341B (zh) | 2014-04-16 | 2014-04-16 | 一种针对升余弦成型mpsk信号的盲符号速率估计与定时方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103916341A true CN103916341A (zh) | 2014-07-09 |
| CN103916341B CN103916341B (zh) | 2017-04-26 |
Family
ID=51041750
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410152732.2A Expired - Fee Related CN103916341B (zh) | 2014-04-16 | 2014-04-16 | 一种针对升余弦成型mpsk信号的盲符号速率估计与定时方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103916341B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112702293A (zh) * | 2020-12-23 | 2021-04-23 | 天津光电通信技术有限公司 | 一种基于循环谱的信号参数分析方法及系统 |
| CN113037663A (zh) * | 2021-03-09 | 2021-06-25 | 山东大学 | 一种改进的适用于非恒包络信号的码元速率估计算法 |
| CN113708913A (zh) * | 2021-09-03 | 2021-11-26 | 四川安迪科技实业有限公司 | 一种定时偏差的对称折半查找估计值的有效性的判定方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6377312B1 (en) * | 1998-08-24 | 2002-04-23 | Samsung Electronics Co., Ltd. | Adaptive fractionally spaced equalizer for received radio transmissions with digital content, such as DTV signals |
| CN101601195A (zh) * | 2007-02-01 | 2009-12-09 | 艾利森电话股份有限公司 | 涉及数字订户线的设备和方法 |
| CN102474497A (zh) * | 2009-07-30 | 2012-05-23 | 松下电器产业株式会社 | 符号率检测器和接收装置 |
| CN103701731A (zh) * | 2013-12-31 | 2014-04-02 | 上海高清数字科技产业有限公司 | 符号率估计方法及装置 |
-
2014
- 2014-04-16 CN CN201410152732.2A patent/CN103916341B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6377312B1 (en) * | 1998-08-24 | 2002-04-23 | Samsung Electronics Co., Ltd. | Adaptive fractionally spaced equalizer for received radio transmissions with digital content, such as DTV signals |
| CN101601195A (zh) * | 2007-02-01 | 2009-12-09 | 艾利森电话股份有限公司 | 涉及数字订户线的设备和方法 |
| CN102474497A (zh) * | 2009-07-30 | 2012-05-23 | 松下电器产业株式会社 | 符号率检测器和接收装置 |
| CN103701731A (zh) * | 2013-12-31 | 2014-04-02 | 上海高清数字科技产业有限公司 | 符号率估计方法及装置 |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112702293A (zh) * | 2020-12-23 | 2021-04-23 | 天津光电通信技术有限公司 | 一种基于循环谱的信号参数分析方法及系统 |
| CN113037663A (zh) * | 2021-03-09 | 2021-06-25 | 山东大学 | 一种改进的适用于非恒包络信号的码元速率估计算法 |
| CN113037663B (zh) * | 2021-03-09 | 2022-04-08 | 山东大学 | 一种改进的适用于非恒包络信号的码元速率估计方法 |
| CN113708913A (zh) * | 2021-09-03 | 2021-11-26 | 四川安迪科技实业有限公司 | 一种定时偏差的对称折半查找估计值的有效性的判定方法 |
| CN113708913B (zh) * | 2021-09-03 | 2023-04-11 | 四川安迪科技实业有限公司 | 一种定时偏差的对称折半查找估计值的有效性的判定方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103916341B (zh) | 2017-04-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102497337B (zh) | 一种基于稀疏度自适应的压缩感知无线通信信道估计方法 | |
| CN108989260B (zh) | 基于Gardner的改进型全数字定时同步方法及装置 | |
| CN108768604B (zh) | 一种用于pcm/fm多符号检测的低复杂度位同步方法 | |
| CN1719818B (zh) | 一种正交频分复用系统的采样频率偏差跟踪信号处理装置及方法 | |
| CN103856431B (zh) | 基于快速傅里叶变换的mpsk调制下开环载波同步算法 | |
| CN106059973B (zh) | 频偏估计方法和系统 | |
| CN103368896B (zh) | 一种高阶调制解调中载波恢复的方法 | |
| CN107342960B (zh) | 一种适合幅度相移键控的非数据辅助频偏估计方法 | |
| CN103916341B (zh) | 一种针对升余弦成型mpsk信号的盲符号速率估计与定时方法 | |
| CN102946372B (zh) | Ofdm系统中频偏和时偏的联合估计方法 | |
| CN107623647B (zh) | 一种基于离散导频辅助的载波同步方法 | |
| CN106953823A (zh) | 一种基于同步序列的无线通信高精度帧同步方法 | |
| CN104717051B (zh) | 一种并行解调位同步中的插值估计方法 | |
| CN105407067B (zh) | 一种突发模式gmsk的相干载波和定时的恢复方法 | |
| CN113938368B (zh) | 基于fpga的lte-v收端同步的方法及系统 | |
| CN102088432B (zh) | 一种ofdm系统的采样频差校正方法和装置 | |
| CN107612860B (zh) | 适用于802.11ac接收机的同步与下采样估计方法 | |
| CN103390428B (zh) | 用于同步对正弦曲线输入的采样的方法 | |
| CN103346992B (zh) | 一种ofdm系统中的多普勒变化率估计方法和系统 | |
| CN110061943A (zh) | 一种mpsk信号解调中的符号定时同步方法 | |
| CN111131110B (zh) | 适合突发通信逆序组合的判决引导载波恢复系统及方法 | |
| CN107864017B (zh) | 一种相位校正方法和装置 | |
| CN103001914B (zh) | 基于循环前缀的小数倍载波同步算法 | |
| CN104486060A (zh) | 一种基于滑动积分的比特位同步方法及系统 | |
| CN108737319A (zh) | 一种目标ofdm信号的实时检测方法及装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170426 Termination date: 20210416 |