发明名称 基于PN序列自相关的自适应DTMB频谱检测方法
摘要 本发明公开了一种基于PN序列自相关的自适应DTMB频谱检测方法,包括以下步骤:获取DTMB频谱,并进行基于伪噪声序列PN序列自相关的粗判决以获取第一判决统计量T<sub>PNAC</sub>和第一判决门限γ<sub>PNAC</sub>;根据第一判决门限γ<sub>PNAC</sub>设置第一门限阈值γ<sub>1</sub>和第二门限阈值γ<sub>2</sub>,以初步判断状态;在初步判断状态之后,进一步进行基于不等间隔PN序列累积自相关的细判决,以获取判决统计量和判决门限从而判定DTMB频谱的状态。本发明实施例的检测方法,通过设置双门限阈值,实现采用基于伪噪声序列PN序列自相关的粗判决和基于不等间隔PN序列累积自相关的细判决,有效提高频谱检测性能,并且降低计算复杂度,很好地满足用户的使用要求。
申请公布号 CN104506262B 申请公布日期 2017.04.05
申请号 CN201410782702.X 申请日期 2014.12.16
申请人 清华大学;北京数字电视国家工程实验室有限公司 发明人 王军;罗子威;张超;阳辉;薛永林
分类号 H04B17/382(2015.01)I 主分类号 H04B17/382(2015.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 张大威
主权项 一种基于PN序列自相关的自适应DTMB频谱检测方法,其特征在于,包括以下步骤:获取地面数字多媒体广播DTMB频谱,并进行基于伪噪声序列PN序列自相关的粗判决以获取第一判决统计量T<sub>PNAC</sub>和第一判决门限γ<sub>PNAC</sub>,其中,所述基于伪噪声序列PN序列自相关的粗判决以获取第一判决统计量T<sub>PNAC</sub>和第一判决门限γ<sub>PNAC</sub>具体包括:对于固定感知时间内的接收信号y(n)进行相关间隔为信号帧长M<sub>i</sub>、相关长度为帧头长度L<sub>i</sub>的自相关,i为不同帧头长度模式,并且把所述固定感知时间内共C个自相关结果累加平均,第一统计量为:<maths num="0001"><math><![CDATA[<mrow><msub><mi>t</mi><mrow><mi>P</mi><mi>N</mi><mi>A</mi><mi>C</mi></mrow></msub><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>CL</mi><mi>i</mi></msub></mrow></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>c</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>C</mi><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><msub><mi>L</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub></munderover><mi>y</mi><mrow><mo>(</mo><mi>m</mi><mo>+</mo><mi>n</mi><mo>+</mo><msub><mi>cM</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><msup><mi>y</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>m</mi><mo>+</mo><mi>n</mi><mo>+</mo><mo>(</mo><mrow><mi>c</mi><mo>+</mo><mn>1</mn></mrow><mo>)</mo><msub><mi>M</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0001109945730000011.GIF" wi="1174" he="151" /></maths>其中,m为自相关起始位置,c为帧序号,n为帧内符号序号;其次所述第一判决统计量T<sub>PNAC</sub>为:<maths num="0002"><math><![CDATA[<mrow><msub><mi>T</mi><mrow><mi>P</mi><mi>N</mi><mi>A</mi><mi>C</mi></mrow></msub><mo>=</mo><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><mn>0</mn><mo>&le;</mo><mi>m</mi><mo>&le;</mo><msub><mi>M</mi><mi>i</mi></msub></mrow></munder><mo>|</mo><msub><mi>t</mi><mrow><mi>P</mi><mi>N</mi><mi>A</mi><mi>C</mi></mrow></msub><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>|</mo><mo>,</mo></mrow>]]></math><img file="FDA0001109945730000012.GIF" wi="486" he="103" /></maths>以及对于指定的虚警概率P<sub>FA</sub>,所述第一判决门限γ<sub>PNAC</sub>为:<maths num="0003"><math><![CDATA[<mrow><msub><mi>&gamma;</mi><mrow><mi>P</mi><mi>N</mi><mi>A</mi><mi>C</mi></mrow></msub><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mn>0</mn><mn>2</mn></msubsup><mi>l</mi><mi>n</mi><mfrac><mn>1</mn><mrow><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>F</mi><mi>A</mi></mrow></msub><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><msub><mi>M</mi><mi>i</mi></msub></mrow></msup></mrow></mfrac><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>,</mo></mrow>]]></math><img file="FDA0001109945730000013.GIF" wi="661" he="143" /></maths>其中,<img file="FDA0001109945730000014.GIF" wi="278" he="64" />为归一化噪声功率,<img file="FDA0001109945730000015.GIF" wi="59" he="67" />为噪声功率;根据所述第一判决门限γ<sub>PNAC</sub>设置第一门限阈值γ<sub>1</sub>和第二门限阈值γ<sub>2</sub>,以根据所述第一判决统计量T<sub>PNAC</sub>、所述第一门限阈值γ<sub>1</sub>和第二门限阈值γ<sub>2</sub>初步判断所述DTMB频谱的状态;以及在初步判断所述DTMB频谱的状态之后,对所述DTMB频谱进一步进行基于不等间隔PN序列累积自相关的细判决,以获取判决统计量和判决门限从而判定所述DTMB频谱的状态,其中,所述对所述DTMB频谱进一步进行基于不等间隔PN序列累积自相关的细判决,以获取判决统计量和判决门限从而判定所述DTMB频谱的状态具体包括:判断DTMB帧头的模式;如果所述DTMB帧头为固定相位模式,则对于所述固定感知时间内的接收信号y(n)进行相关间隔为d个信号帧长dM<sub>i</sub>、相关长度为帧头长度L<sub>i</sub>的自相关,并且把所述固定感知时间内共C<sub>d</sub>个自相关结果累加平均,公式如下:<maths num="0004"><math><![CDATA[<mrow><mi>A</mi><mrow><mo>(</mo><mrow><mi>m</mi><mo>,</mo><mi>d</mi></mrow><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>C</mi><mi>d</mi></msub><msub><mi>L</mi><mi>i</mi></msub></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>c</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>C</mi><mi>d</mi></msub><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><msub><mi>L</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub></munderover><mi>y</mi><mrow><mo>(</mo><mrow><mi>m</mi><mo>+</mo><mi>n</mi><mo>+</mo><msub><mi>cM</mi><mi>i</mi></msub></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><msup><mi>y</mi><mo>*</mo></msup><mrow><mo>(</mo><mrow><mi>m</mi><mo>+</mo><mi>n</mi><mo>+</mo><mrow><mo>(</mo><mrow><mi>c</mi><mo>+</mo><mi>d</mi></mrow><mo>)</mo></mrow><msub><mi>M</mi><mi>i</mi></msub></mrow><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0001109945730000016.GIF" wi="1221" he="154" /></maths>其中,m为自相关起始位置,c为帧序号,n为帧内符号序号;构造第二统计量为:<maths num="0005"><math><![CDATA[<mrow><msub><mi>t</mi><mrow><mi>D</mi><mi>I</mi><mi>P</mi><mi>N</mi><mi>A</mi><mi>C</mi></mrow></msub><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>d</mi><mo>=</mo><mn>1</mn></mrow><mi>D</mi></munderover><msub><mi>C</mi><mi>d</mi></msub><msub><mi>C</mi><mrow><mi>d</mi><mo>+</mo><mn>1</mn></mrow></msub><mi>A</mi><mrow><mo>(</mo><mrow><mi>m</mi><mo>,</mo><mi>d</mi></mrow><mo>)</mo></mrow><msup><mi>A</mi><mo>*</mo></msup><mrow><mo>(</mo><mrow><mi>m</mi><mo>,</mo><mi>d</mi><mo>+</mo><mn>1</mn></mrow><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0001109945730000021.GIF" wi="854" he="143" /></maths>其中,共采用D个不同间隔的累积自相关,第二判决统计量为:<maths num="0006"><math><![CDATA[<mrow><msub><mi>T</mi><mrow><mi>D</mi><mi>I</mi><mi>P</mi><mi>N</mi><mi>A</mi><mi>C</mi></mrow></msub><mo>=</mo><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><mn>0</mn><mo>&le;</mo><mi>m</mi><mo>&le;</mo><msub><mi>M</mi><mi>i</mi></msub></mrow></munder><mo>|</mo><msub><mi>t</mi><mrow><mi>D</mi><mi>I</mi><mi>P</mi><mi>N</mi><mi>A</mi><mi>C</mi></mrow></msub><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>|</mo><mo>,</mo></mrow>]]></math><img file="FDA0001109945730000022.GIF" wi="549" he="102" /></maths>对于指定的虚警概率P<sub>FA</sub>,第二判决门限γ<sub>DIPNAC</sub>为:<maths num="0007"><math><![CDATA[<mrow><msub><mi>&gamma;</mi><mrow><mi>D</mi><mi>I</mi><mi>P</mi><mi>N</mi><mi>A</mi><mi>C</mi></mrow></msub><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mn>0</mn><mn>2</mn></msubsup><mi>l</mi><mi>n</mi><mfrac><mn>1</mn><mrow><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>F</mi><mi>A</mi></mrow></msub><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><msub><mi>M</mi><mi>i</mi></msub></mrow></msup></mrow></mfrac><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>,</mo></mrow>]]></math><img file="FDA0001109945730000023.GIF" wi="686" he="139" /></maths>其中,<img file="FDA0001109945730000024.GIF" wi="414" he="143" />为归一化噪声功率,<img file="FDA0001109945730000025.GIF" wi="59" he="70" />为噪声功率,<img file="FDA0001109945730000026.GIF" wi="446" he="131" />为归一化系数;如果所述DTMB帧头为旋转相位模式,则对于所述固定感知时间内的接收信号y(n)进行相关间隔为d个信号帧长dM<sub>i</sub>、相关长度为帧头长度L<sub>i</sub>的自相关,并且把所述固定感知时间内共C<sub>d</sub>个自相关结果累加平均,公式如下:<maths num="0008"><math><![CDATA[<mrow><msub><mi>A</mi><mi>r</mi></msub><mrow><mo>(</mo><mrow><mi>m</mi><mo>,</mo><mi>d</mi></mrow><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>C</mi><mi>d</mi></msub><msub><mi>L</mi><mi>i</mi></msub></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>c</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>C</mi><mi>d</mi></msub><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><msub><mi>L</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub></munderover><munderover><mi>&Sigma;</mi><mrow><mi>a</mi><mo>=</mo><mn>0</mn></mrow><mn>1</mn></munderover><mi>y</mi><mrow><mo>(</mo><mrow><mi>m</mi><mo>+</mo><mi>n</mi><mo>+</mo><msub><mi>cM</mi><mi>i</mi></msub></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><msup><mi>y</mi><mo>*</mo></msup><mrow><mo>(</mo><mrow><mi>m</mi><mo>+</mo><mi>n</mi><mo>+</mo><msup><mrow><mo>(</mo><mrow><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>a</mi></msup><mrow><mo>(</mo><mrow><mi>c</mi><mo>+</mo><mi>d</mi></mrow><mo>)</mo></mrow><msub><mi>M</mi><mi>i</mi></msub></mrow><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0001109945730000027.GIF" wi="1398" he="150" /></maths>其中,m为自相关起始位置,c为帧序号,n为帧内符号序号;构造第三统计量为:<maths num="0009"><math><![CDATA[<mrow><msub><mi>t</mi><mrow><mi>r</mi><mo>,</mo><mi>D</mi><mi>I</mi><mi>P</mi><mi>N</mi><mi>A</mi><mi>C</mi></mrow></msub><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>d</mi><mo>=</mo><mn>1</mn></mrow><mi>D</mi></munderover><msub><mi>C</mi><mi>d</mi></msub><msub><mi>C</mi><mrow><mi>d</mi><mo>+</mo><mn>1</mn></mrow></msub><msub><mi>A</mi><mi>r</mi></msub><mrow><mo>(</mo><mrow><mi>m</mi><mo>,</mo><mi>d</mi></mrow><mo>)</mo></mrow><msubsup><mi>A</mi><mi>r</mi><mo>*</mo></msubsup><mrow><mo>(</mo><mrow><mi>m</mi><mo>,</mo><mi>d</mi><mo>+</mo><mn>1</mn></mrow><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0001109945730000028.GIF" wi="902" he="142" /></maths>其中,第三判决统计量为:<maths num="0010"><math><![CDATA[<mrow><msub><mi>T</mi><mrow><mi>D</mi><mi>I</mi><mi>P</mi><mi>N</mi><mi>A</mi><mi>C</mi></mrow></msub><mo>=</mo><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><mn>0</mn><mo>&le;</mo><mi>m</mi><mo>&le;</mo><msub><mi>M</mi><mi>i</mi></msub></mrow></munder><mo>|</mo><msub><mi>t</mi><mrow><mi>r</mi><mo>,</mo><mi>D</mi><mi>I</mi><mi>P</mi><mi>N</mi><mi>A</mi><mi>C</mi></mrow></msub><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>|</mo><mo>,</mo></mrow>]]></math><img file="FDA0001109945730000029.GIF" wi="566" he="111" /></maths>对于指定的虚警概率P<sub>FA</sub>,第三判决门限为:<maths num="0011"><math><![CDATA[<mrow><msub><mi>&gamma;</mi><mrow><mi>r</mi><mo>,</mo><mi>D</mi><mi>I</mi><mi>P</mi><mi>N</mi><mi>A</mi><mi>C</mi></mrow></msub><mo>=</mo><msup><mrow><mo>(</mo><mrow><msubsup><mi>&sigma;</mi><mrow><mn>0</mn><mo>,</mo><mi>r</mi></mrow><mn>2</mn></msubsup><mi>ln</mi><mfrac><mn>1</mn><mrow><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>F</mi><mi>A</mi></mrow></msub></mrow><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><msub><mi>M</mi><mi>i</mi></msub></mrow></msup></mrow></mfrac></mrow><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>,</mo></mrow>]]></math><img file="FDA00011099457300000210.GIF" wi="726" he="141" /></maths>其中,<img file="FDA00011099457300000211.GIF" wi="434" he="143" />为归一化噪声功率。
地址 100084 北京市海淀区100084-82信箱