发明名称 一种基于压缩感知技术的BOC信号无模糊捕获算法
摘要 本发明公开了一种基于压缩感知技术的GNSS(全球导航卫星系统)中BOC(二进制偏置载波)调制信号的无模糊捕获算法,该方法将导航信号的捕获过程分为两个阶段完成,在第一阶段采用压缩感知技术对一个伪码周期内的接收信号进行处理,大大缩小接收信号初始相位估计值的范围,第二阶段仅对缩小范围后的数据采用GRASS(General Removing Ambiguity via Sidepeak Suppression)算法进行处理,完成信号的捕获。该方法采用压缩感知技术与GRASS算法相结合的方式,完成BOC信号的无模糊捕获,适用于各族BOC信号,且与传统的GRASS算法相比,所需的硬件资源大大减少。
申请公布号 CN104199057A 申请公布日期 2014.12.10
申请号 CN201410460465.5 申请日期 2014.09.11
申请人 重庆大学 发明人 谭晓衡;张颜;胡小楠;王斌;姚超;陈国庆
分类号 G01S19/24(2010.01)I 主分类号 G01S19/24(2010.01)I
代理机构 代理人
主权项 一种基于压缩感知技术的BOC信号无模糊捕获算法,在接收端将射频前端天线接收的导航卫星信号经放大、下变频处理后,得到中频信号,信号捕获过程包括以下步骤:A、根据BOC调制信号的调制指数α确定采样频率,对中频信号进行数字化,同时根据多普勒频率估计值Δf产生本地载波信号,完成接收信号的载波消除,取一个伪码周期长度的接收信号r={r[0],r[1],…,r[αL‑1]}<sup>T</sup>,L=1023;B、首先采用转换矩阵Ψ将接收信号转换到其稀疏域,得到稀疏域信号R,同时将R与第一阶段观测矩阵Φ<sub>1</sub>相乘得到第一阶段观测信号y<sub>1</sub>;对y<sub>1</sub>进行FWHT变换得到信号Z<sub>1</sub>,找出Z<sub>1</sub>中大于门限r<sub>1</sub>的且不多于N<sub>p</sub>个数,并用w<sub>1</sub>存放这些数在Z<sub>1</sub>中的坐标;对||w<sub>1</sub>||<sub>0</sub>>0进行判决,如果||w<sub>1</sub>||<sub>0</sub>>0,计算出W<sub>1</sub>,这里<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>W</mi><mn>1</mn></msub><mo>=</mo><munder><mo>&cup;</mo><mrow><msub><mi>w</mi><mn>1</mn></msub><mo>&Subset;</mo><msub><mi>w</mi><mn>1</mn></msub></mrow></munder><mo>{</mo><msub><mi>&eta;w</mi><mn>1</mn></msub><mo>+</mo><mn>1</mn><mo>,</mo><msub><mi>&eta;w</mi><mn>1</mn></msub><mo>+</mo><mn>2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>min</mi><mo>{</mo><mi>&eta;</mi><mrow><mo>(</mo><msub><mi>w</mi><mn>1</mn></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mi>&alpha;L</mi><mo>}</mo><mo>}</mo><mo>,</mo></mrow>]]></math><img file="FDA0000568595040000011.GIF" wi="1050" he="110" /></maths>如果||w<sub>1</sub>||<sub>0</sub><0,将多普勒频移估计值Δf增加固定步进Δp,返回步骤A,继续进行捕获;其特征在于:C、令<img file="FDA0000568595040000014.GIF" wi="397" he="102" />以及W<sub>1</sub>中各个元素的坐标为ν,ν∈{1,2,…,||W<sub>1</sub>||<sub>0</sub>},生成大小为<img file="FDA0000568595040000015.GIF" wi="443" he="93" />的哈达玛矩阵<img file="FDA0000568595040000016.GIF" wi="146" he="90" />生成第二阶段观测矩阵Φ<sub>2</sub>和第二阶段稀疏域信号R′,再计算第二阶段观测信号y<sub>2</sub>:y<sub>2</sub>=Φ<sub>2</sub>R′;Φ<sub>2</sub>的表达式为<img file="FDA0000568595040000012.GIF" wi="382" he="95" />R′的第m行第l列元素可表示为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>R</mi><mrow><mi>m</mi><mo>.</mo><mi>l</mi></mrow><mo>&prime;</mo></msubsup><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msup><mi>R</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>aR</mi><mrow><mi>B</mi><mo>/</mo><mi>L</mi></mrow><mn>2</mn></msubsup><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow><mo>,</mo></mtd><mtd><mi>l</mi><mo>=</mo><msub><mi>W</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&nu;</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>0</mn><mo>,</mo></mtd><mtd><mi>l</mi><mo>&NotEqual;</mo><msub><mi>W</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&nu;</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000568595040000013.GIF" wi="756" he="169" /></maths>m∈{1,2,…,M<sub>2</sub>},a=2α‑3,R(l),R<sub>B/L</sub>(l)分别为R,R<sub>B/L</sub>中第l个元素,R<sub>B/L</sub>是由GRASS算法计算得到的辅助信号;D、完成y<sub>2</sub>的FWHT变换,得到信号Z<sub>2</sub>,找出Z<sub>2</sub>中最大值w<sub>2</sub>;将w<sub>2</sub>与门限r<sub>2</sub>进行比较,如果w<sub>2</sub><r<sub>2</sub>,将多普勒频移估计值Δf增加固定步进Δp后返回步骤A,继续进行捕获,如果w<sub>2</sub>>r<sub>2</sub>,则得到扩频码相位W<sub>1</sub>(w<sub>2</sub>)和多普勒频移Δf,捕获成功。
地址 400044 重庆市沙坪坝区沙正街174号重庆大学通信工程学院