一种高动态环境下双M元编码扩频方法及系统与流程

一种高动态环境下双m元编码扩频方法及系统
技术领域
1.本发明属于扩频通信技术领域，特别涉及高动态环境下扩频通信方法及系统。


背景技术：

2.早在20世纪40年代，科学研究人员就开始对扩展频谱通信技术进行研究。扩展频谱通信技术由于其良好的抗干扰、抗衰落和抗多径性能最初主要被应用于军事领域。随着扩展频谱通信技术的发展，该技术现已渗透到民用通信的各个领域，被广泛应用于测控、卫星导航定位等诸多领域。直接序列扩频技术是最常见的扩频技术，但其主要针对传输数据率低、频带利用率有限的情况。随着现代科技水平的发展，人们对通信系统性能的要求越来越高，高数据率、多用户数据传输的业务越来越多。因此，提高扩频通信系统的传输效率和抗干扰能力，降低系统的设计复杂度成为一种迫切需求。此外，扩频技术中建立一个快速精准的同步机制尤为关键，而同步的过程大致分为捕获和跟踪。在捕获方法中，串行搜索法电路设计和运算较为简单，但是当扩频码长度增加时，搜索时间也随之增加；并行搜索法搜索时间短，电路设计和运算较为复杂。此外，在高动态环境下，常见算法的捕获概率较低，捕获性能欠佳。因此，如何有效提高扩频系统中扩频码的相位跟踪精度和频偏估计精度显得尤为重要。
3.有的学者将经典的捕获算法应用于直接序列扩频通信系统，有的学者结合部分匹配滤波器与快速傅里叶变换来改善捕获效果，但直接序列扩频通信系统的性能欠佳，且现有的捕获方法存在参数估计精度低等问题，使得高动态坏境下扩频通信系统的抗干扰性能差，捕获概率低，缺乏实际应用价值。


技术实现要素：

4.本发明的目的是解决在高动态环境下，现有扩频系统扩频码产生方式复杂，抗干扰性能差，且对扩频码的捕获速度慢、捕获精度低等问题，采用双m元编码算法和高动态环境下pmf
‑
h
‑
fft
‑
lf捕获算法，能够更好的降低扩频码产生的复杂度，实现扩频码相位的正确搜索和频率偏差的精准评估，提高扩频通信系统工作效率，并在抗衰落、抗多径干扰等方面呈现出良好的性能。
5.为了实现上述目的，本发明的技术方案如下：
6.一种高动态环境下双m元编码扩频系统，包括：
7.发射端，用于产生原始信号，并对所述原始信号进行双调制处理，获得发送信号；
8.接收端，用于接收发送信号，并对接收到的信号进行处理，恢复原始信号。
9.作为本发明的一种优选方式，所述的双调制处理过程为：基于双m元编码扩频方式，生成扩频码，将生成的扩频码与原始信号进行扩频调制得到扩频信号，再将扩频信号经bpsk二次调制，获得发送信号。
10.进一步地，所述双m元编码扩频方式具体为：
11.随机选择一条扩频码作为原型扩频码，此扩频码的长度为l，用c
0～l
‑1表示扩频码
中的码片，其中l≥m，m＝2
k
；k表示一个传输符号对应的信息序列的比特数；
12.将原型扩频码顺着同一方向依次移动，每次移动d个单位长度，d为正整数，经次循环移位后，对应生成条扩频码；
13.对已生成的条扩频码取反，获得剩余的条扩频码。
14.进一步优选地，所述的接收端包括扩频码捕获单元，所述扩频码捕获单元用于对接收的信号进行扩频码同步处理，具体为：
15.将采样后的信号与本地nco产生的本振信号混频后；
16.将混频信号输入部分匹配滤波器进行相关运算；
17.对部分匹配滤波器输出的数据进行加海明窗处理；
18.采用fft算法进行频域分析；
19.进行阈值逼近判决，完成扩频码相位和频率偏差的二维搜索。
20.进一步优选地，所述接收端还包括解扩解调单元，用于将捕获的信号的频带恢复为原始信息序列的带宽，然后进行解调处理，恢复出所传输的信息。
21.本发明还提出一种高动态环境下双m元编码扩频方法，该方法包括：
22.发射端产生原始信号，并对所述原始信号进行双调制处理，获得发送信号；
23.接收端接收发送信号，并对接收的信号进行解扩解调处理，恢复原始信号。
24.进一步地，对所述原始信号进行双调制处理包括：基于双m元编码扩频方式，生成扩频码，将生成的扩频码与原始信号进行扩频调制得到扩频信号，再将扩频信号经bpsk二次调制，获得发送信号。
25.进一步地，所述双m元编码扩频方式具体为：
26.随机选择一条扩频码作为原型扩频码，此扩频码的长度为l，用c
0～l
‑1表示扩频码中的码片，其中l≥m，m＝2
k
；k表示一个传输符号对应的信息序列的比特数；
27.将原型扩频码在同一方向顺次移动，每次移动d个单位长度，d为正整数，经次循环移位后，对应生成条扩频码；
28.对已生成的条扩频码取反，获得剩余的条扩频码。
29.进一步优选地，所述接收端通过以下方法对接收的信号进行扩频码同步处理：
30.将采样后的信号与本地nco产生的本振信号混频后；
31.将混频信号输入部分匹配滤波器进行相关运算；
32.对部分匹配滤波器输出的数据进行加海明窗处理；
33.采用fft算法进行频域分析；
34.进行阈值逼近判决，完成扩频码相位和频率偏差的二维搜索。
35.本发明与现有技术相比，具有如下有益效果：
36.(1)提出了采用双m元编码方式生成扩频码，降低生成扩频码的复杂度。
37.(2)提出了pmf
‑
h
‑
fft
‑
lf捕获算法，即通过部分匹配滤波器进行相关处理、对其输
出的数据加海明窗、接着进行傅里叶变换、最后利用似然函数估计所需参数，从而完成扩频码相位的正确搜索和频率偏差的精准评估，确保扩频码的同步捕获，提高了整个扩频通信系统的工作效率，并在抗衰落、抗多径干扰等方面呈现出良好的性能。
附图说明
38.图1是本发明实施例提供的一种高动态环境下双m元编码扩频系统框图；
39.图2是本发明实施例提供的一种高动态环境下双m元编码扩频方法流程图；
40.图3是本发明实施例提供的扩频码生成器示意图；
41.图4是本发明实施例提供的利用双m元扩频编码方式产生扩频码的示意图；
42.图5是本发明实施例提供的海明窗的仿真示意图；
43.图6是本发明实施例提供的pmf
‑
h
‑
fft
‑
lf捕获算法流程示意图；
具体实施方式
44.为了便于理解本发明，下面结合附图和具体实施例，对本发明进行更详细的说明。附图中给出了本发明的较佳的实施例。但是，本发明可以以许多不同的形式来实现，并不限于本说明书所描述的实施例。相反地，提供这些实施例的目的是使对本发明的公开内容的理解更加透彻全面。
45.本实施例提供的高动态环境下双m元编码扩频方法，是针对高动态环境下扩频通信系统的精准发送与同步捕获接收方法，如图1所示，具体包括：
46.步骤s1、数字信号的产生与处理
47.发射端产生的信号s
k
为信息流，则发送信号s(t)可表示为：
[0048][0049]
其中，r
s
为码元速度，t
s
为码元宽度，s
n
为信息码，为信息码，t
s
＝1/r
s
。
[0050]
此外，需要对信号流进行串并转换处理，并依据扩频系统的设计要求进行初始化参数设置。
[0051]
步骤s2、对发送信号进行双调制处理。
[0052]
对发送信号进行双调制处理，分别为扩频调制和bpsk调制。其中，扩频调制基于双m元编码扩频方式进行，此种扩频方式产生的扩频码，不仅保证了良好的自相关特性，且可降低选码难度，降低实际系统开发的复杂度。
[0053]
1.扩频调制
[0054]
本发明中，基于m序列来构造性能良好的扩频码。一般用生成函数来表示一个序列，如图3所示，若{b
m
}＝{b0,b1,b2,
…
}表示移位寄存器的输出序列，其下标的数字表示时间，则输出序列的生成函数g(x)为：
[0055][0056]
若
[0057][0058]
在本发明中，待发送的信息序列按每k比特划分为一个符号，那么所需m条扩频码对应传输每个输入符号的信息内容。双m元编码的具体步骤如下：
[0059]
(1)首先随机选择一条扩频码作为原型扩频码，此扩频码的长度为l，用c
0～l
‑1表示扩频码中的码片，其中l≥m，m＝2
k
；
[0060]
(2)将原型扩频码在同一方向顺次移动，每次移动d个单位长度，d为正整数，d的取值需保证循环移位后的扩频码与已产生的扩频码不重复，经次循环移位后，可对应生成条扩频码；
[0061]
(3)剩余的条扩频码可对已生成的扩频码取反获得；取反的含义为：对已生成的条扩频码中
‑
1对应变成1,1对应变成
‑
1。
[0062]
循环移位后的扩频码可用q
i
进行表示，具体产生方式见图4。
[0063]
由于采用双m元编码扩频方式，因此只需要产生条扩频码即可，极大降低了生成扩频码的复杂度。
[0064]
将生成的m条扩频码与原始信号进行扩频调制得到扩频信号。
[0065]
2.bpsk调制
[0066]
扩频后的信号，需要通过bpsk进行二次调制。信号经过双调制后，可表示为：
[0067][0068]
其中，e代表平均功率，p(t)为扩频码波形，θ为随机载波相位，f
c
为载波频率，τ为接收的扩频码相位偏差，f
d
为频率偏差。
[0069]
对于未知的τ和f
d
来说，两者的平均似然函数是r的条件概率密度函数中关于θ期望的极限。其表示为：
[0070][0071][0072]
步骤s3、输入信道进行传输
[0073]
在信道传输过程中，假想干扰是不存在的，仅仅存在awgn信道的噪声对信号产生干扰，噪声信号用n(t)表示，e[n(t)]＝0。由于噪声是零均值的，故e[r
i
]＝s
i
，
且高斯白噪声的自相关函数为
[0074]
步骤s4、获取接收信号
[0075]
通过awgn信道后，接收信号可以表示为：
[0076]
r(t)＝s(t) n(t)0≤t≤t
[0077]
其中，s(t)为发送信号，n(t)为加性高斯白噪声，t为观察持续的时间。
[0078]
由于在时间间隔t内保持正交的基函数完备集满足：
[0079][0080]
那么发送信号和接收信号可表示为：
[0081][0082][0083][0084]
利用基函数完备集和求得上述式子的展开系数，分别表示为：
[0085][0086][0087][0088]
步骤s5、扩频码捕获处理
[0089]
如图6所示，扩频码捕获过程主要是：将采样的信号与本地nco产生的本振信号混频后，输入到部分匹配滤波器进行相关运算；对部分匹配滤波器输出的数据进行加海明窗处理，再采用fft算法进行频域分析，最后进行阈值逼近判决，从而完成扩频码相位和频率偏差的二维搜索。
[0090]
其中，采用海明窗进行处理，可使旁瓣显著减小，其序列为：
[0091][0092]
如图5所示，经过海明窗进行处理，能够减小傅里叶变换输出频谱峰值衰减程度并减小频谱泄漏。
[0093]
若fft计算出的频谱峰值高于捕获判决预设的阈值，则判定当前接收信号的扩频码与本地扩频码相位一致，从而完成扩频码相位和频率偏差的二维搜索。此时接收信号的扩频码长度为l，频率偏差为f
d
，长度为w的部分匹配滤波器的数量为q，经归一化处理，并对q个部分匹配滤波器输出的相关值进行n点fft运算后，fft结果中的第k点的归一化幅频响应为：
[0094][0095]
若q≠2
n
，那么fft结果中的第k点的归一化幅频响应为：
[0096][0097]
步骤s6、参数估计及优化。
[0098]
在r(t)展开式中，系数是高斯的、统计独立的，方差为因此，
[0099][0100]
带入似然函数的求解公式当中，将最大似然估计无关因子消去，可得：
[0101][0102]
通过积分运算，交换求极限和积分运算的次序，可得：
[0103][0104][0105]
那么平均似然函数为：
[0106][0107]
其中，s为信号波形在观察时间间隔内的能量。因为s在所选参数范围内变化不大，因此可以将包含其的因子删去，可得：
[0108][0109]
假设接收载波的相位θ在[0,2π]内均匀分布，关于θ的三角展开式及积分可写为：
[0110][0111]
其中，
[0112]
[0113][0114]
通过遍历τ和f
d
的所有可能值，并选择r(τ,f
d
)的最大值来确定估计参数的值，则有：
[0115][0116]
步骤s7、解扩解调恢复信号
[0117]
为在接收端获得有效的数据信息，在接收端需要对接收的信号进行正确解扩，而解扩的前提是扩频码的同步，本发明的捕获算法有效的保证了扩频码的同步。在接收端进行解扩处理，将信号的频带恢复为原始信号序列的带宽，即为中频调制信号。然后紧接着进行解调处理，最后恢复出所传输的信息s(t)，完成信息的传输。
[0118]
本发明提供的另一个实施例是一种高动态环境下双m元编码扩频系统，如图2所示，该系统包括发射端101、接收端102和传输信道103。
[0119]
其中，发射端101用于产生原始数字信号，并对原始信号串并转换后，基于双m元编码扩频方式，生成扩频码，将生成的扩频码与原始信号进行扩频调制得到扩频信号，再将扩频信号经bpsk二次调制，获得发送信号。
[0120]
发送信号经传输信道103传输至接收端102。
[0121]
接收端102，首先对接收的信号采用pmf
‑
h
‑
fft
‑
lf捕获算法，实现扩频码的同步，具体为：
[0122]
捕获单元将采样的发送信号与本地nco产生的本振信号混频后，输入到部分匹配滤波器进行相关运算；对部分匹配滤波器输出的数据进行加海明窗处理，再采用fft算法进行频域分析，最后进行阈值逼近判决，从而完成扩频码相位和频率偏差的二维搜索。
[0123]
然后，解扩解调单元对接收的有效信号进行解扩处理，将信号的频带恢复为原始信号序列的带宽，即为中频调制信号。然后紧接着进行解调处理，最后恢复出所传输的信息，完成信息的传输。
[0124]
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。