专利名称：一种伪码序列的捕获方法
技术领域：
本发明涉及直接序列扩频通信技术领域，特别涉及一种伪码序列 的捕获方法。
背景技术：
软件无线电的概念是指用可编程微处理器代替模拟或者数字器 件实现信号处理的大部分功能，软件无线电把接收机的仿真与实现紧 密联系起来，具有更大的灵活性。自从软件无线电的概念出现以后，
从80年代以硬件为主的数字无线电技术开始向着今后更加灵活、更加 先进的方向发展。
软件捕获信号算法的实现可以使接收机更加灵活的适应不同码 长的不同伪随机噪声序列、载波频率、相位变化等，也为接收机的广 泛使用提供了保障。1)接收机的灵活配置为接收机随时对所需要的 信号进行实时模拟处理提供了可能性。可以针对不同的情况采用不同 的算法，对数据进行采集、分析和测试。2)接收机的软件模块化设 计也可以对各模块的输入和输出信号进行可视化监测，因此在不同算 法后可以对各个模块的性能进行评估，以便择优使用。3)接收机的 软件化能够大大提高接收机的利用效率。接收机的软件化不仅可以灵 活接收处理来自不同伪随机噪声码，也可以根据接收机所处的环境灵 活使用不同的算法来处理数据、分析数据等，这样可以提高接收机的
使用效率。
软件接收机包括硬件和软件设计。硬件设计主要是模数转换器
(ADC)和数字信号处理器(DSP)的设计，软件设计包括各种处理算法，
如相关算法、定位算法等。
对同步捕获方法的研究是随着扩频技术的出现而开始的，早期曾 经采用发送参考序'列、发送特定同步码、统一定时以及序列状态估计 等方法实现码同歩，但是这些都不是很理想。比较成熟的技术就是相 关搜索技术，即采用不同相位的本地序列与接收序列做相关运算，通 过相关值的大小判断序列是否同步。而这些方法目前都是靠使用硬件 来实现的。
就目前而言，关于接收机伪随机噪声码的同步捕获方法的研究主 要集中在捕获方法、捕获时间、检测判决变量在各种信道下的统计分 布三个方面。通常在一般的接收机中，对于码相位的同步，需要对所 有的码相位进行搜索，计算量大，搜索时间长。
发明内容
本发明所要解决的技术问题是提供一种计算量小、搜索时间短的 伪码序列的捕获方法。
本发明的原理是这样的根据伪随机码序列的相关性知道，对于 任意一个伪随机码序列的自相关函数具有如下性质
其中/ 为伪随机码序列的周期，即码长；71.为序列的码元持续时间， 即码元宽度。上式也就是说明了如果两个伪随机码序列的相位差在一 个码元内，两个序列归一化的相关结果比较大，否则，两个伪随机码
序列归一化的相关结果只是-丄。
为了方便表示，这里假设接收到两个伪随机噪声码的混合序列， 并且只用其伪码形式表示，则接收码为
xOO^/W^ + rJ + P^G + q),其中/W,表示第/个伪随机噪声序列码；「 表示接收到的相应码所对应的码相位。
再假设本地复现了第l个伪随机噪声序列码，本地复现码可以表
示为
少(,),("r;)
2)
其中/W,表示第l个伪随机噪声序列码；z";表示复现的第l个伪随机噪
声序列码序列所对应的码相位。
如果把接收的码序列和复现的码序列做相关运算，则会有
柳-x(,).:K,)
=[/W,(, + r,) +所2(〖+ ^)] ZW,(, + r;)
+ r,) + r;) + iW2(/ + r2) + r;)
(1-3)
上式中的相关结果可以看成有两部分，其中第一部分v ,")是同
一码族中的两个码序列之间的相关部分，第二部分/^(t)是不同码族 中的两个码序列之间的相关部分。从伪随机码序列的性质可以知道，
'两个不同码族中的码序列之间的相关值为O (理想情况下)。从对第 二部分的说明能够看出，即使接收到的码中有多个伪随机噪声序列码 存在，都属于这一部分，其值理论上都为0 (实际上也是接近0的一个 值)。
下面再来说明第一部分的情况。关于第一部分也可以分为两种情
况
1、当接收到的序列与本地复现的序列是同族序列中的相位不同
(指相位差大于一个码片)的两个序列时，根据式a-1)，第一部
分可以写成
<formula>formula see original document page 7</formula>(l-4)
其中，<formula>formula see original document page 7</formula>7;为码元宽度。上式说明随着p的增大，这样的两 个序列相关值也是比较小的。
2当接收到的序列与本地复现的序列是同族序列中的相位接近(指
相位差小于一个码片)两个的序列时，同样根据式(i-i)，第一部
分可以写成
<formula>formula see original document page 7</formula>(1-5)
其中，h-r;^7;.， 7;为码元宽度。上式说明随着两个序列的相位差越
小，其归一化的相关值越大，二者呈现一个线性增长关系。
通过上述分析可以得到不论复现的序列是否与接收序列是同 族，只要复现序列与接收序列的相位相差不在一个码片之内，那么其 归一划的相关值就是比较小。
利用上述的分析结果，可以得出如果把本地复现的多个不同相 位的序列叠加在一起，用它们的和来与接收码序列相关，根据相关结 果就可以判断出多个不同相位序列中的一个序列是否出现接收序列。 利用这个分析结果，把前面的算法进行了改进。 为了说明该算法，进一步假设接收的码序列为
<formula>formula see original document page 8</formula>，
本地复现的码序列可以表示为少<formula>formula see original document page 8</formula>，其中 /W/表示第l种伪随机噪声序列码；r;表示复现的第l种伪随机噪声序列 码的第一个码序列所对应的码相位；、表示复现的第1种伪随机噪声 序列码的第二个码序列所对应的码相位，这里还要假设h -r;^7;和 z",-r;^7;，表示复现序列中的第一个序列的码相位与接收序列的码 相位的相位差在一个码片或者码元内，而第二个复现序列中的码相位 与接收序列的码相位的相位差大于一个码元周期。
同样如果把接收的序列码和复现的序列码做相关运算，则会有
<formula>formula see original document page 8</formula>
(1-6)
上式中一共有三项，其中第三项还包括两部分，分别是<formula>formula see original document page 9</formula>
和
<formula>formula see original document page 9</formula>
根据前面分析知道，这两项都是不同族的两个序列进行相关运算，其 归一化的相关值是很小的，可以忽略不计的。
式G-6)中的第二项是<formula>formula see original document page 9</formula>，表示的是同族的 两个相位不同的序列进行相关运算。由于已经假设卜,-r;^z;，因此
该项归一化的相关结果可以直接得到，应该是-1。
式(1-6)中的第一项是/W," + O"A^ + 0，表示的是同族的
两个相位不同的序列进行相关运算。由于已经假设h-巧l^7;，因此
该项归一化的相关结果也可以直接从式(1-5)得到，应该是
1-
到此可以得到接收码序列与复现码序列的相关值为:
<formula>formula see original document page 9</formula>
(1-7)
上式与式(1-5)相比，多了一项-丄，显然，是因为经过求和才使
得相关结果中增加了该项，会导致结果有一定的衰减。
需要补充说明一点，由于上面推导假设了求和中有一个复现码序 列与接收序列是同族但是相位差大于一个码元宽度，如果有仿个这样
的序列出现在和项中，式(卜7)需要重新写成下式<formula>formula see original document page 10</formula>(l-8)
其中参数范表示参与求和的序列中与接收序列同族但是与接收序列 码相位的相位差大于一个码元的序列个数。
根据上述原理，本发明的技术方案是这样实现的， 一种伪码序 列的捕获方法，其特征在于， '第一步，由伪随机码产生器产生位于接收机上空的卫星伪码序列 (卫星伪随机码)的本地复制码，将此本地复制码的所有不同码相位 序列均匀的划分成n分组，每分组包含m个不同码相位序列，将每分 组中的m个码相位序列相加，得到该分组的复合码序列，记为Pi， i 为分组号，i=l，2,3……n;
第二步，将复合码序列P(分别与接收到的卫星伪码序列进行相关 积分，如果复合码序列Pi相关积分值大于阈值Thl，则表示检测的卫 星伪码序列包含在该分组的码相位序列中；如果复合码序列P,相关积 分值小于阈值Thl，则继续进行复合码序列Pw的相关积分运算；如 果所有的复合码序列P'相关积分值都不大于阈值Thl，则改为用每个 单独的码相位序列进行搜索，即采用传统码相位序列搜索方法；
第三步，对于上述复合码序列Pi相关积分值大于阈值Thl的相应 分组，将该分组中的m个码相位序列再分别和接收到的卫星伪码序列 进行相关积分，如果某个码相位序列相关积分值大于阈值Th2，该码 相位序列就和卫星伪码序列同步，即卫星伪码序列被捕获。
所述卫星伪码序列是GPS信号。
所述将此本地复制码的所有不同码相位序列均匀的划分成n分 组，每分组包含m个不同码相位序列，其中m为10 32， n为32 103。
所述将此本地复制码的所有不同码相位序列均匀的划分成n分 组，每分组包含m个不同码相位序列，其中m为32， n为32。 所述卫星伪码序列是伽利略信号。
本发明对要进行搜索的码相位序列先进行分组，根据比较阈值， 找出包含卫星伪码序列的分组，然后再对分组中的码相位序列进行搜 索，找出最终要伺步的码相位序列，这样将会大大降低搜索量和搜索 时间，提高接收机的性能。
当m-l， n二1023时，相当于每分组一个码相位，有1023组， 这种情况等于没有进行分组，和传统的码相位搜索方法一样了 。
当m=1023， n-l时，相当于只有一个分组，但是这个分组中 有1023个码相位，同样和传统的方法是一样的。
当m-10， n=103,每组10个码相位， 一共103组，最后一组 3个码相位，相关积分计算量Cgal03 + 10-113。
因此，最小相关积分计算量应该是min(〃 + 1023/z7),"e[U023]，可以 求出《*32，最小的相关积分计算量C2-32 + 31-63 。而传统的码相位 检测方法需要的相关积分计算量是1023。
以下结合
和具体的实施方式对本发明作进一步的详细 说明。
图l是本发明的流程图。
图2是两个序列自相关函数归一化的图形，横坐标为两个同族码 相差的码相位r，纵坐标为归一化自相关函数值。-l是归一化相关
值，7;为码元宽度。
图3 GPS信号仿真模块图。
图4是信噪比为-10dB时的情况，横坐标为相关时间点，纵坐标 为归一化自相关函数值。
图4a表示12个复现码序列的延迟码之和都与接收码序列的相位 对齐。
图4b表示12个复现码序列的延迟码之和，其中的一个复现码序 列与接收码序列相位相同。
图4c表示本地复现码的12个延迟码中不存在与接收码对齐的码。
图5是信噪比为-20dB的情况下，横坐标为相关时间点，纵坐标 为归一化自相关函数值。
图5a表示12个复现码序列的延迟码之和都与接收码序列的相位 对齐。
图5b表示12个复现码序列的延迟码之和，其中的一个复现码序 列与接收码序列相位相同。
图5c表示本地复现码的12个延迟码中不存在与接收码对齐的码。
图6是信噪比为-30dB的情况，横坐标为相关时间点，纵坐标为 归一化自相关函数值。图6a表示12个复现码序列的延迟码之和都与接收码序列的相位对齐。
图6b表示12个复现码序列的延迟码之和，其中的一个复现码序
列与接收码序列相位相同。
图6c表示本地复现码的12个延迟码中不存在与接收码对齐的码。
具体实施例方式
如图1所示为本发明的流程图1)首先由伪随机码产生器，产 生位于接收机上空的卫星伪码序列的本地复现码。2)将本地复现码 的所有不同码相位序列，均匀的划分成n分组，每分组包含m个码相 位序列。3)将每分组中的m个码相位序列相加得到复合码序列。4) 将相加后的复合码序列，和接收到的卫星伪码序列进行相关积分，如 果大于阈值Thl ，则表示检测的卫星伪码序列包含在此分组中。如 果不大于阈值Thl，则对下一分组的复合码序列进行相关积分，直到 找到大于阈值Thl的复合码序列。5)对找出的复合码序列，将其相 应分组中的每一个码相位序列与卫星伪码序列再进行相关积分，找出 其中大于阈值Th2的那个码相位序列，这样就可以找出卫星伪码序列 的具体码相位。
实施例当m = 32, n = 32，每分组32个码相位序列， 一共32 分组，最后一分组31个码相位序列，因为与32相差不大，可以取相 同阈值作为判定条件或取阈值的31/32作为判定条件。然后对每一分 组的码相位序列求和，再和卫星伪随机码序列相关积分，假如第12
分组的相关积分值大于阈值2dB,说明这个分组中含有卫星伪随机码
序列，那么就停止进行组搜索。随后在第12分组内进行32个码相位 序列的具体搜索，假如第25个码相位序列进行相关积分后，其值大 于3dB，那么就可以判断找到了卫星伪随机码序列。该码相位序列的 排列次序为11X32+25=377，相关积分的计算量是12+25=37。
下面通过对GPS信号和伽利略信号的捕获实例的仿真进一步说. 明本发明的实施方式。仿真的平台是Matlab6.5。
参照图2，是两个序列自相关函数归一化的图形，两个序列为不 同码相位的同族伪随机码，橫坐标为这两个同族码相差的码相位r ， 纵坐标为归一化自相关函数值。可以看出，在一个周期内，相差的码 相位r越小，相关值越大。
参照图3，因为GPS信号中伪码的长度是1023，将每个GPS随机 码分为86分组，每分组12个码相位序列，每分组中的不同码相位序 列进行相加，然后和接收卫星伪随机码序列进行相关积分运算。图中 的加性白高斯噪声模块(AWGN)，即对相关积分运算的输入信号引入 的是高斯白噪声，并且此噪声是加在输入信号端的。Gold Sequence generator模块是格尔德码产生器，产生伪随机码。Matrix multiply 模块是相关积分的模块。Delay模块是延时模块，延吋一个码相位。 还有一个码相位序列相加器，进行码相位序列的相加，即将12个码 相位序列进行相加。
图4、图5、图6，分别给出了在不同信噪比(GPS信号和加性
高斯噪声之比)三种情况的相关结果。从仿真结果可以看出-
图4 a 、图5 a 、图6 a可以说明如果12个本地复现码的 码相位都与接收码的码相位相同的话，相关结果中的峰值都是比较明 显的。
图4b 、图5b 、图6b可以说明如果12个本地复现码序列 中有一个序列的码相位与接收码的码相位相同的话，其相关结果的峰 值虽然较相应图4 a 、图5 a 、图6 a的峰值有明显下降，但是其 峰值也是很明显的。峰值下降的原因是有11个与接收码的码相位不 同的序列参与了求和，实际上就是因为式(1-8)中的第二项在起作 用。
图4c 、图5c 、图6c可以说明如果本地复现的12个码序 列的码相位中没有一个相位码与接收码的码相位相同，则相关结果是 更小的。
将本发明也可用于捕获伽利略系统E1频段开放信号中。对于伽 利略信号的捕获，与GPS信号的主要区别就是码的长度不同，其他步 骤相同，可以得到一致结果。
因此，上述结果说明，完全可以先对同族码的码相位序列进行分 组求和，确定分组后，再确定具体的码相位序列。这样可以减少相关 积分的次数，减少运算量和运算时间，对于接收机来说，在同样的时 钟周期下，可以縮短捕获信号的吋间。
权利要求
1、一种伪码序列的捕获方法，其特征在于，包括以下步骤第一步，由伪随机码产生器产生位于接收机上空的卫星伪码序列的本地复制码，将此本地复制码的所有不同码相位序列均匀的划分成n分组，每分组包含m个不同码相位序列，将每分组中的m个码相位序列相加，得到该分组的复合码序列，记为Pi，i为分组号，i＝1，2，3……n；第二步，将复合码序列Pi分别与接收到的卫星伪码序列进行相关积分，如果复合码序列Pi相关积分值大于阈值Th1，则表示检测的卫星伪码序列包含在该分组的码相位序列中；如果复合码序列Pi相关积分值小于阈值Th1，则继续进行复合码序列Pi+1的相关积分运算；如果所有的复合码序列Pi相关积分值都不大于阈值Th1，则改为用每个单独的码相位序列进行搜索，即采用传统码相位序列搜索方法；第三步，对于上述复合码序列Pi相关积分值大于阈值Th1的相应分组，将该分组中的m个码相位序列再分别和接收到的卫星伪码序列进行相关积分，如果某个码相位序列相关积分值大于阈值Th2，该码相位序列就和卫星伪码序列同步，即卫星伪码序列被捕获。
2、 根据权利要求1所述的一种伪码序列的捕获方法，其特征在 于，所述卫星伪码序列是GPS信号。
3、 根据权利要求2所述的一种伪码序列的捕获方法，其特征在 于，所述将此本地复制码的所有不同码相位序列均匀的划分成n分 组，每分组包含m个不同码相位序列，其中m为10 32， n为32 103。
4、 根据权利要求3所述的一种伪码序列的捕获方法，其特征在 于，所述m为32， n为32。
5、 根据权利要求1所述的一种伪码序列的捕获方法，其特征在 于，所述卫星伪码序列是伽利略信号。
全文摘要
本发明涉及直接序列扩频通信技术领域，特别涉及一种伪码序列的捕获方法。本发明通过将本地复制码的所有不同码相位序列均匀的划分成n分组，每分组包含m个不同码相位序列，将每分组中的m个码相位序列相加，得到该分组的复合码序列。首先利用复合码序列分别与接收到的卫星伪码序列进行相关积分，搜索含有卫星伪码序列的复合码序列；再通过在该复合码序列相应分组的m个码相位序列中搜索卫星伪码序列。这种方法可以广泛应用GPS信号、伽利略信号以及其它通信信号的软件接收机中，可以灵活捕获不同长度的直接序列扩频码。
文档编号G01S19/30GK101207405SQ20071019923
公开日2008年6月25日 申请日期2007年12月14日 优先权日2007年12月14日
发明者何秋生, 周文益, 王江安 申请人:西安华迅微电子有限公司