主权项 |
一种直扩信号截获检测方法,其特征在于,该方法包括以下步骤:步骤1:获取离散数据序列x(n),n=0,1,2,…,N‑1:从传感器接收N个采样点的实时采集数据或提取N个采样点存储的现成数据作为待处理的数据序列x(n),n=0,1,2,…,N‑1,其中N为采样点数;步骤2:对所述离散数据序列x(n)进行自相关估计,得到自相关序列r(m),m=0,1,2,…,M‑1,其中M为自相关序列长度,m为自相关序列的离散时间索引;步骤3:对所述自相关序列r(m)的绝对值进行平滑处理,即根据下式计算得到平滑自相关序列R(k):<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>R</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>k</mi><mo>+</mo><mi>K</mi></mrow></mfrac><munderover><mo>Σ</mo><mrow><mi>m</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>k</mi><mo>+</mo><mi>K</mi></mrow></munderover><mo>|</mo><mi>r</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>|</mo><mo>,</mo><mi>k</mi><mo>=</mo><mn>0</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>K</mi><mo>-</mo><mn>1</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>R</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>K</mi><mo>+</mo><mn>1</mn></mrow></mfrac><munderover><mo>Σ</mo><mrow><mi>m</mi><mo>=</mo><mi>k</mi><mo>-</mo><mi>K</mi></mrow><mrow><mi>k</mi><mo>+</mo><mi>K</mi></mrow></munderover><mo>|</mo><mi>r</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>|</mo><mo>,</mo><mi>k</mi><mo>=</mo><mi>K</mi><mo>,</mo><mn>...</mn><mo>,</mo><mi>M</mi><mo>-</mo><mn>1</mn><mo>-</mo><mi>K</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>R</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>M</mi><mo>-</mo><mi>k</mi><mo>+</mo><mi>K</mi></mrow></mfrac><munderover><mo>Σ</mo><mrow><mi>m</mi><mo>=</mo><mi>k</mi><mo>-</mo><mi>K</mi></mrow><mrow><mi>M</mi><mo>-</mo><mn>1</mn></mrow></munderover><mo>|</mo><mi>r</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>|</mo><mo>,</mo><mi>k</mi><mo>=</mo><mi>M</mi><mo>-</mo><mi>K</mi><mo>,</mo><mn>...</mn><mo>,</mo><mi>M</mi><mo>-</mo><mn>1</mn></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001015581990000011.GIF" wi="1062" he="431" /></maths>其中,K为平滑处理的半窗长,且K<M,k为平滑自相关序列的离散时间索引;步骤4:设定判定伪码周期具有等间隔规律的最少自相关峰个数L,L>1,再根据自相关序列长度M设定搜索步进j,j=J<sub>s</sub>,J<sub>s</sub>+1,J<sub>s</sub>+2,…,J<sub>e</sub>,其中J<sub>s</sub>为搜索起始点,J<sub>e</sub>为搜索终点,J<sub>s</sub><J<sub>e</sub>,J<sub>e</sub><M/L;步骤5:按照所述搜索步进j对所述平滑自相关序列R(k)进行累积平均,即从j开始,将L个间隔为j的平滑自相关序列R(k)相加后求平均,得到累积量C(j),计算公式如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>C</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>L</mi></mfrac><mo>(</mo><mrow><mi>R</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><mi>R</mi><mrow><mo>(</mo><mrow><mn>2</mn><mi>j</mi></mrow><mo>)</mo></mrow><mo>+</mo><mn>...</mn><mo>+</mo><mi>R</mi><mrow><mo>(</mo><mrow><mi>L</mi><mi>j</mi></mrow><mo>)</mo></mrow></mrow><mo>)</mo><mo>=</mo><mfrac><mn>1</mn><mi>L</mi></mfrac><munderover><mo>Σ</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><mi>R</mi><mo>(</mo><mrow><mi>i</mi><mi>j</mi></mrow><mo>)</mo><mo>,</mo><mi>j</mi><mo>=</mo><msub><mi>J</mi><mi>s</mi></msub><mo>,</mo><msub><mi>J</mi><mi>s</mi></msub><mo>+</mo><mn>1</mn><mo>,</mo><mn>...</mn><mo>,</mo><msub><mi>J</mi><mi>e</mi></msub></mrow>]]></math><img file="FDA0001015581990000012.GIF" wi="1470" he="134" /></maths>其中,R(ij)为从j开始,以j为间隔,从平滑自相关序列R(k)中抽取的第ij个自相关估计值,ij为上述累积平均中的平滑自相关序列离散时间索引,i为所述L个间隔j的间隔倍数,i=1,2,…,L;步骤6:搜索累积量C(j)的最大值C(j<sub>max</sub>)作为检验统计量,其中j<sub>max</sub>为最大值对应的离散时间索引;步骤7:将所述最大值C(j<sub>max</sub>)与判决门限比较,如大于判决门限,则认为存在直扩信号,并根据下式估计伪码周期T<sub>p</sub>:T<sub>p</sub>=j<sub>max</sub>/f<sub>s</sub>,其中f<sub>s</sub>为离散信号的采样频率;否则,认为离散数据序列x(n)中不含直扩信号。 |