一种面向星际通信链路的频移对称啁啾扩频调制与解调方法

1.本发明属于通信技术领域，具体涉及一种面向星际通信链路的频移对称啁啾扩频调制与解调方法。


背景技术：

2.由于现阶段仍有大部分区域地面移动蜂窝网络无法全覆盖，需要卫星通信网络作为补充，为这些区域提供基础的通信服务。相比于地面移动蜂窝网络，卫星通信网络的距离相对遥远，从而带来很大的传输链路损耗，受到雨衰的影响，接收链路的噪声性能将进一步恶化；同时，由于卫星为维持运行轨道所需的对地相对速度较大，受到多普效应的影响，发射与接收回路存在较大的频率偏移，且随着通信频段的提高，多普勒效应带来的频移将进一步扩大。
3.频移啁啾扩频调制技术(fscm)是啁啾扩频(css)的一种，由于啁啾扩频对频偏不敏感，频移啁啾扩频技术具有抗多普勒频移能力强、抗干扰能力强和接收灵敏度高的特点，被应用于lora通信协议的调制和解调过程中。频移啁啾扩频将信息调制为啁啾信号的频移量，不同频移量的啁啾信号构成了相互正交的啁啾码片，接收端通过乘以极性相反的啁啾信号得到频移分量，即以频移量为频率的单音信号，通过傅里叶变换后峰值检测得到频移量。由于信息被调制为相对频移量，接收和发射间多普勒效应引入的固定频移偏差将很容易检测到并消除。
4.虽然频移啁啾扩频调制技术具有良好的抗多普勒效应，但相比于直接序列扩频等扩频技术，其啁啾码片间的相关性并不强，当数据吞吐量和数据速率比较高时，其抗噪声性能远不及其他的扩频方式。近些年来，有学者提出将频移啁啾与相位调制结合的方法，命名为psk-lora，通过将部分数据调制为频移啁啾信号的初始相位，从而获得更高的数据吞吐量，但这种方式在接收时初始相位需要通过相干解调方式实现，增加了设备的复杂度。考虑到接收解调算法的复杂度，也有学者提出了使用同相和正交分量同时传输啁啾扩频调制信号的iqcss方法，然而该方法提高了频谱利用率的同时，牺牲了系统的抗噪声性能。
5.针对卫星通信的特点，探寻对多普勒频移不敏感和抗噪声性能强的通信方式已经成为最近研究的热点。


技术实现要素：

6.鉴于上述，本发明提供了一种面向星际通信链路的频移对称啁啾扩频调制和解调方法。该方法在保持频移啁啾扩频方式抗多普勒频移等优良特性的基础上，通过对称啁啾方式加强了信号的相关性能，提高了系统的抗噪声性能。
7.一种面向星际通信链路的频移对称啁啾扩频调制和解调方法，其具体为：
8.发射端按照扩频因子sf和支路数r将输入信息比特划分为r个深度为sf的信息块，并以信息块的高m位作为支路的索引码，其中m为log2r。每个信息块经过gray编码后以十进制信息表示，并称第i个信息块得到的十进制信息为di。信息块以(di·
bw/2
sf
)作为相对频
移量对原始对称啁啾信号进行频移啁啾扩频调制，其中bw表示啁啾信号的传输带宽，并称第i个信息块调制得到的频移对称啁啾信号为si。r个信息块对应的r路频移对称啁啾信号经线性叠加后在信道中传输。
9.其中，对称啁啾信号由一对频率变化率极性相反的啁啾信号组成。啁啾信号按照频率变化率极性可分为上啁啾信号和下啁啾信号，其中上啁啾信号的频率变化率为正数，即瞬时频率由小变大，下啁啾信号的频率变化率为负数，即瞬时频率由大变小。所述对称啁啾信号由一个上啁啾信号和一个下啁啾信号拼接构成，按照信号的先后排列顺序可分为两种：上啁啾信号在前、下啁啾信号在后的称为正对称啁啾信号；而下啁啾信号在前、上啁啾信号在后的称为反对称啁啾信号。所述对称啁啾信号采用的上、下啁啾信号的频率变化率互为相反数，频率变化范围都为-bw/2～bw/2。
10.对称啁啾信号的上、下啁啾信号拼接过程保证了瞬时频率连续性，即在一个啁啾信号的周期内，上啁啾信号结束频率与下啁啾信号起始频率相等，上啁啾信号起始频率与下啁啾信号结束频率相等。对称啁啾信号的上、下啁啾信号拼接过程保证了相位连续性。即在一个啁啾信号的周期内，上啁啾信号结束相位与下啁啾信号起始相位相等，上啁啾信号起始相位与下啁啾信号结束相位相等。
11.对称啁啾信号的频移过程是一种瞬时频率的循环移位，这保证了对称啁啾信号频移仍然保持相位和瞬时频率的连续性。其具体为：当瞬时频率在频率变化范围边界将产生跳变，即从bw/2跳变到-bw/2，或者是从-bw/2跳变到bw/2。当频移量为δf，正对称啁啾信号的频率将按照(-bw/2 δf)～(bw/2)、(-bw/2)～(-bw/2 δf)、(-bw/2 δf)～(-bw/2)、(bw/2)～(-bw/2 δf)的顺序线性变化；当频移量为δf，反对称啁啾信号的频率将按照：(-bw/2 δf)～(-bw/2)、(bw/2)～(-bw/2 δf)、(-bw/2 δf)～(bw/2)、(-bw/2)～(-bw/2 δf)的顺序线性变化。
12.接收端对接收到的信号依次进行解对称啁啾处理、快速傅里叶变换、频域峰值检测和信息解码与重组。
13.其中，接收端的解对称啁啾过程，是通过对接收信号乘以极性相反的对称啁啾信号得到相对频移量(di·
bw/2
sf
)：若发射时使用正对称啁啾信号调制，接收时使用反对称啁啾信号解调；若发射时使用反对称啁啾信号调制，接收时使用正对称啁啾信号解调。
14.进一步地，解对称啁啾后的信号是以相对频移量(di·
bw/2
sf
)为频率的三角函数的组合，通过快速傅里叶变换得到频域特征，根据解对称啁啾后的信号频域的峰值位置得到相对频移量，并计算每个信息块的十进制信息di。并根据每个信息块中位于高m位的索引码拼接r个信息块还原出信息比特。
15.收发过程采用的帧结构包括前导码，同步字和用户数据，其中前导码个数为n
pre
，同步字个数为n
sync
，用户数据为n
data
。n
pre
和n
sync
由收发双方相互约定，n
data
则由突发用户数据包长度决定。其中，帧结构中的前导码由连续的正、反对称啁啾信号构成，具体为：当采用正对称啁啾信号作为发射调制、反对称啁啾信号作为接收接调时，前导码由n
pre
个连续不加频移调制的正对称啁啾信号和n
sync
个不加频移调制的反对称啁啾信号构成。反之，所述当采用反对称啁啾信号作为发射调制、正对称啁啾信号作为接收接调时，前导码由n
pre
个连续不加频移调制的反对称啁啾信号和n
sync
个不加频移调制的正对称啁啾信号构成。而用户数据则为n
data
个符号，每个符号为r路频移对称啁啾信号的线性叠加。
16.通过帧结构中通过引入前导码快速实现发射端与接收端间多普勒频偏的消除和时间同步。前导码采用未经调制的对称啁啾信号，用于频偏和时偏的消除，其特征在于：对前导码进行解对称啁啾和快速傅里叶变化，以频谱分量的平均值和差值分别估计频偏和时偏。
17.通过帧结构中通过引入同步字可以实现前导码和用户数据分离由于同步字采用与前导码极性相反的对称啁啾信号，在接收时可利用正反扫描的方法快速定位同步字，并实现前导码与用户数据分离。所述正反扫描的方法，是接收端对收到的信号先用正对称啁啾信号进行解对称啁啾并通过快速傅里叶变化得到正向扫描频谱，后用反对称啁啾信号进行解对称啁啾并通过快速傅里叶变化得到反向扫描频谱，通过分析正、反向扫描频谱的差异定位同步字，并在接收帧中区分前导码和用户数据。
附图说明
18.图1a为本发明方法的发送端调制框图
19.图中各环节为：101—串并转换，102—索引添加，103—gray编码，104—二、十进制转换，105—频移量计算，106—频移对称啁啾信号调制，107—线性叠加，108—发射天线。
20.图1b为本发明方法的接收端解调框图
21.图中各环节为：201—并串转换，202—索引去除，203—gray译码，204—十、二进制转换，205—峰值检索，206—快速傅里叶变换，207—解对称啁啾，208—接收天线。
22.图2为本发明方法的帧结构
23.图3a为对称啁啾信号的瞬时频率随时间变化曲线
24.图3b为对称啁啾信号的相位随时间变化曲线
25.图3c为对称啁啾信号的幅度随时间变化曲线
26.图4a为频移后的对称啁啾信号的瞬时频率随时间变化曲线
27.图4b为频移后的对称啁啾信号的相位随时间变化曲线
28.图4c为频移后的对称啁啾信号的幅度随时间变化曲线
29.图5a为啁啾信号的模糊函数图
30.图5b为对称啁啾信号的模糊函数图
具体实施方式
31.为了更为具体地描述本发明，下面结合附图及具体实施方式对本发明的技术方案进行详细说明。
32.如图1a和1b所示，描述了一种面向星际通信链路的频移对称啁啾扩频调制和解调方法的流程示意图。该方法的调制参数主要有扩频因子sf，传输带宽bw，支路数r，并设m＝log2r,n＝2
sf
。
33.如图1a，发送端流程为：
34.1.串并转换：按照扩频因子和支路数将输入信息比特划分成r个等长度的信息块，每个信息块长度为(sf-m)。如式(1)所示，bi表示第i个信息块，b
i,j
表示第i个信息块内第j位信息比特，i的取值范围1～r，j的取值范围1～(sf-m)。
35.bi＝{b
i,1
,b
i,2
,
…
,b
i,j
,
…
,b
i,sf-m
}#(1)
36.2.索引添加：以信息块序号i的二进制作为索引码ci，与信息块中的信息码拼接。拼接时，索引码ci在前，信息码bi在后，如式(2)所示。
37.《ci,bi》＝{c
i,1
,c
i,2
,
…
,c
i,m
,b
i,1
,b
i,2
,
…
,b
i,j
,
…
,b
i,sf-m
,}#(2)
38.3.gray编码：对信息块《ci,bi》进行gray编码，设编码后的信息块为《c
′i,b
′i》。
39.《c
′i,b
′i》＝{c
′
i,1
,c
′
i,2
,
…
,c
′
i,m
,b
′
i,1
,b
′
i,2
,
…
,b
′
i,j
,
…
,b
′
i,sf-m
}#(3)
40.4.进制转换：对编码后的信息块《c
′i,b
′i》进制转换，设转换后的十进制为di。
41.di＝c
′
i,1
·2sf-1

…
c
′
i,m
·2sf-m
b
′
i,1
·2sf-m-1
.. b
′
i,sf-m
#(4)
42.5.频移量计算：根据转换后的十进制计算信号的频移量fi。
[0043][0044]
6.以相对频移量fi调制对称啁啾信号，频移后的信号为si(n)，见式(6)。
[0045][0046]
其中，μ代表对称啁啾信号的极性，μ为 1代表正对称啁啾信号，μ为-1代表反对称啁啾信号。
[0047]
7.线性叠加：将频移后的信号线性叠加后传输，最终得到的信号为si(n)，见式(7)。
[0048][0049]
如图1b，接收端流程为：
[0050]
1.解对称啁啾：接收端通过与一个极性相反的对称啁啾信号相乘，见式(8)。
[0051][0052]
其中，e(n)表示信道噪声。
[0053]
2.快速傅里叶变换：解对称啁啾后的信号可进一步化简为式(9)。
[0054][0055]
3.对快速傅里叶变化的结果进行峰值检索，得到峰值处的频率为fi，根据式(5)计算得到十进制di。
[0056]
4.进制转换，将十进制di转换成《c
′i,b
′i》。
[0057]
5.gray译码得到《ci,bi》。
[0058]
6.根据索引码ci，将多路信息bi进行拼接后还原得到输入信息比特。
[0059]
实际通信过程中，接收端捕获信号时需要消除时间偏移δt和频率偏移δf的影响，本方法中采用对称啁啾信号作为前导码的方式实现信号的捕获。如图2所示，为本方法采用的帧结构，主要包括前导码、同步字和用户数据。前导码为n
pre
个符号，n
pre
的取值范围
为4～65536，具体由用户定义，其采用的对称啁啾信号与调制所采用的极性相同；同步字为n
sync
个符号，n
pre
的取值范围1～2，具体由用户定义，其采用的对称啁啾信号与调制所采用的极性相反。前导码用于解调时对频偏的快速估计，而同步字则用于区分前导码和用户数据区，实现解调对时间的同步。
[0060]
m1＝δt
·fs
#(10)
[0061][0062]
在通信过程中时间偏移会导致采样点的滞后，频率偏移会引入固有频差，如式(10)和式(11)所示，以m1和m2衡量时偏和频偏对于系统的影响。则接收信号的表达式为式(12)。
[0063][0064]
则解对称啁啾后的信号为式(13)。
[0065][0066]
则傅里叶变化后包含的信号分量k1和k2为式(14)和(15)，其中前导码中di为0。
[0067]
k1＝-μ
·
m1 m2·
n#(14)
[0068]
k2＝μ
·
m1 m2·
n#(15)
[0069]
可由k1和k2估计频偏量m1和时偏量m2。
[0070]
如图3a和3b所示，对称啁啾信号由一个上啁啾信号和一个下啁啾信号拼接而成，虚线处表示在拼接时刻瞬时频率和相位的连续性。
[0071]
如图4a和4b所示，对称啁啾信号的频移是对瞬时频率的循环移位，即当瞬时频率到达bw/2，下一时刻将跳变到-bw/2，反之亦然。同时这种方式也保证了频移后的信号在拼接时刻瞬时频率和相位的连续性。
[0072]
如图5a所示为啁啾信号的模糊函数图，如图5b所示为对称啁啾信号的模糊函数图。其中x轴代表时延量，y轴代表多普勒频移量，z轴表示相关系数。如图5a，除了在时延为0，多普勒频移量为0的位置达到最大相关1，仍存在许多时延与多普勒频移量的组合是得相关系数为1。而在图5b中，除了时延为0，多普勒频移量为0的位置达到最大相关1外，其余位置得到的相关系数均小于0.5。
[0073]
根据模糊函数的定义，如式(16)所示为啁啾信号在时延τ
*
＝fd/u时相关值，如式(17)所示为对称啁啾信号在时延τ
*
＝fd/u时的相关值。从上述两种信号的模糊函数分析，对称啁啾信号比啁啾信号的相关特性更好，当在信号捕获时对称啁啾信号对频偏的影响更不敏感。
[0074]
[0075][0076]
其中τ
*
＝fd/u。
[0077]
上述对实施例的描述是为便于本技术领域的普通技术人员能理解和应用本发明。熟悉本领域技术的人员显然可以容易地对上述实施例做出各种修改，并把在此说明的一般原理应用到其他实施例中而不必经过创造性的劳动。因此，本发明不限于上述实施例，本领域技术人员根据本发明的揭示，对于本发明做出的改进和修改都应该在本发明的保护范围之内。