发明名称 基于循环平稳特性的空时分组码MC-CDMA信号盲识别方法
摘要 本发明请求保护一种基于循环平稳特性的空时分组码MC-CDMA信号盲识别方法,属于号处理技术领域。该方法先根据空时分组码MC-CDMA系统的基带物理模型、信号构造及矩阵表示,建立了信号相关函数的选择标准,估计出来这些相关函数对应的四阶累积量,然后把四阶循环累积量跟信道的参数结合起来去建立一个循环统计量,最后根据循环频率的检测提取检测门限跟这个循环统计量进行比较做判决实现识别。计算机模拟表明,本发明采用的算法能以较小的复杂度有效地识别STBCMC-CDMA信号,且能在较低的输入信噪比条件下良好地工作。同时,跟传统的识别算法比较,本发明采用的算法利用四阶循环累积量的优点减少随机噪声的影响,从而改善系统性能。
申请公布号 CN104393963A 申请公布日期 2015.03.04
申请号 CN201410521619.7 申请日期 2014.09.30
申请人 重庆邮电大学 发明人 张天骐;裴光盅;张刚;高超
分类号 H04L1/06(2006.01)I 主分类号 H04L1/06(2006.01)I
代理机构 重庆华科专利事务所 50123 代理人 康海燕
主权项 基于循环平稳特性的空时分组码MC‑CDMA信号盲识别方法,其特征在于,以空时分组码多载波码分多址信号的循环平稳特性为基础进行分析,先建立信号相关函数的选择标准,估计出来这些相关函数对应的四阶累积量,然后把四阶累积量跟信道的参数结合起来去建立一个循环统计量,最后根据循环频率的检测提取检测门限跟这个循环统计量进行比较做判决实现识别;具体步骤如下:步骤一:将信号Y(t)按T<sub>s</sub>采样形成一个时间序列,提取信号的码率,按照式(29)取得信号的相关函数数目R:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mi>f</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>E</mi><mo>{</mo><mi>g</mi><mo>[</mo><msub><mi>y</mi><mi>i</mi></msub><mrow><mo>(</mo><mn>2</mn><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msub><mi>y</mi><mi>j</mi></msub><mrow><mo>(</mo><mn>2</mn><mi>t</mi><mo>)</mo></mrow><mo>]</mo><mo>}</mo><mi>orE</mi><mo>{</mo><mi>g</mi><msub><mrow><mo>[</mo><mi>y</mi></mrow><mi>i</mi></msub><mrow><mo>(</mo><mn>2</mn><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msubsup><mi>y</mi><mi>j</mi><mo>*</mo></msubsup><mrow><mo>(</mo><mn>2</mn><mi>t</mi><mo>)</mo></mrow><mo>]</mo><mo>}</mo></mtd></mtr><mtr><mtd><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mi>lp</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000581803300000011.GIF" wi="1196" he="184" /></maths>                (29)t是空时码块的编号,t=1,2,...,LT,l是正数,f(t)是周期p的函数,g[·]是一个线性相关函数,定义<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>g</mi><mo>[</mo><mo>&CenterDot;</mo><mo>]</mo><mo>=</mo><msub><mi>y</mi><mi>i</mi></msub><mrow><mo>(</mo><mn>2</mn><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msub><mi>y</mi><mi>i</mi></msub><mrow><mo>(</mo><mn>2</mn><mi>t</mi><mo>)</mo></mrow><mi>or</mi><msub><mi>y</mi><mi>i</mi></msub><mrow><mo>(</mo><mn>2</mn><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msubsup><mi>y</mi><mi>j</mi><mo>*</mo></msubsup><mrow><mo>(</mo><mn>2</mn><mi>t</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000581803300000016.GIF" wi="900" he="97" /></maths>0≤i,j∈I&lt;M,I为整数集;步骤二:在循环频率α=±T/2T<sub>S</sub>,根据式(28)计算出来对应信号的四阶累积量值<img file="FDA0000581803300000012.GIF" wi="235" he="93" /><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mover><mi>C</mi><mo>^</mo></mover><mrow><mn>4</mn><mi>Y</mi></mrow></msub><mrow><mo>(</mo><mi>&alpha;</mi><mo>,</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mi>N</mi></munderover><mi>Y</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><mi>Y</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mi>p</mi><mo>)</mo></mrow><mi>Y</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><mi>Y</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mi>p</mi><mo>)</mo></mrow><mo>.</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>j&alpha;t</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000581803300000013.GIF" wi="1231" he="139" /></maths>(28);步骤三:根据式(33)估计对应每个相关函数的循环稳定第一峰的相位φ<sup>r</sup><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>&phi;</mi><mo>=</mo><mi>angle</mi><mo>{</mo><mfrac><mn>1</mn><mi>T</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>T</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>E</mi><mo>{</mo><mi>g</mi><mo>[</mo><msub><mi>y</mi><mi>i</mi></msub><mrow><mo>(</mo><mn>2</mn><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msub><mi>y</mi><mi>i</mi></msub><mrow><mo>(</mo><mn>2</mn><mi>t</mi><mo>)</mo></mrow><mo>]</mo><mo>}</mo><mo>.</mo><mi>exp</mi><mo>[</mo><mo>-</mo><mfrac><mrow><mi>j</mi><mn>2</mn><mi>&pi;t</mi></mrow><mi>T</mi></mfrac><mo>]</mo><mo>}</mo></mrow>]]></math><img file="FDA0000581803300000014.GIF" wi="1086" he="158" /></maths>    (33);步骤四:根据假设检测与<img file="FDA0000581803300000015.GIF" wi="220" he="102" />收敛于零均值方差渐近为Σ<sub>ac</sub>的多元正态分布这个性质估计一个循环统计量值th<sub>ac</sub>,采用式(31)或式(34)进行估计:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>th</mi><mi>ac</mi></msub><mo>=</mo><mi>N</mi><mover><msub><mi>c</mi><mi>aY</mi></msub><mrow><mo>^</mo><mo>^</mo></mrow></mover><mrow><mo>(</mo><mi>t</mi><mo>;</mo><mi>&tau;</mi><mo>)</mo></mrow><msubsup><mover><mi>&Sigma;</mi><mo>^</mo></mover><mi>ac</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>c</mi><mi>aY</mi></msub><msup><mrow><mo>(</mo><mi>t</mi><mo>;</mo><mi>&tau;</mi><mo>)</mo></mrow><mi>T</mi></msup></mrow>]]></math><img file="FDA0000581803300000021.GIF" wi="578" he="92" /></maths>    (31)式中N为数据长度,<img file="FDA0000581803300000022.GIF" wi="79" he="82" />是<img file="FDA0000581803300000023.GIF" wi="174" he="80" />的协方差矩阵的估计<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>th</mi><mi>ac</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>r</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>R</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mrow><mo>|</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><msup><mi>&phi;</mi><mi>r</mi></msup><mo>)</mo></mrow><msubsup><mover><mi>C</mi><mo>^</mo></mover><mrow><mn>4</mn><mi>Y</mi></mrow><mi>r</mi></msubsup><mrow><mo>(</mo><mi>&alpha;</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0000581803300000024.GIF" wi="662" he="138" /></maths>    (34)式中R是按式(29)分类的相关数目;步骤五:根据循环频率的式(27)),用式(35)取得检测门限值γ<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>H</mi><mn>0</mn></msub><mo>:</mo><mi>&alpha;</mi><mo>&NotElement;</mo><msub><mi>A</mi><mn>4</mn></msub><mo>,</mo><mo>&ForAll;</mo><msubsup><mrow><mo>{</mo><msub><mi>&tau;</mi><mi>n</mi></msub><mo>}</mo></mrow><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><mo>&DoubleRightArrow;</mo><msub><mover><mi>c</mi><mo>^</mo></mover><mrow><mn>4</mn><mi>Y</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>;</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>&epsiv;</mi><mrow><mn>4</mn><mi>Y</mi></mrow><mi>T</mi></msubsup><mrow><mo>(</mo><mi>&alpha;</mi><mo>;</mo><mi>&tau;</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000581803300000025.GIF" wi="921" he="92" /></maths>H<sub>1</sub>:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>&alpha;</mi><mo>&NotElement;</mo><msub><mi>A</mi><mn>4</mn></msub><mo>,</mo></mrow>]]></math><img file="FDA0000581803300000028.GIF" wi="166" he="67" /></maths>存在<img file="FDA0000581803300000026.GIF" wi="140" he="92" />(27)<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mo>&DoubleRightArrow;</mo><msub><mover><mi>c</mi><mo>^</mo></mover><mrow><mn>4</mn><mi>Y</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>;</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>c</mi><mrow><mn>4</mn><mi>Y</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>;</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>&epsiv;</mi><mrow><mn>4</mn><mi>Y</mi></mrow><mi>T</mi></msubsup><mrow><mo>(</mo><mi>&alpha;</mi><mo>;</mo><mi>&tau;</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000581803300000027.GIF" wi="678" he="82" /></maths>其中A<sub>4</sub>是四阶循环累积量的循环频率集合;p<sub>F</sub>=p(χ<sup>2</sup>>γ)   (35)p<sub>F</sub>是虚警概率;步骤六:进行比较判决,即如果th≥γ说明假设H<sub>1</sub>成立,即采用了STBC信号;否则,反之H<sub>0</sub>成立,即采用了BLAST信号。
地址 400065 重庆市南岸区黄桷垭崇文路2号