发明名称 一种高速移动环境下的多符号联合信道估计方法
摘要 一种高速移动环境下的多符号联合信道估计方法,基于OFDM无线通信系统,采用基于导频的信道估计方法,包括步骤1)利用复指数基扩展模型对时间频率双选信道进行建模;2)计算联合估计OFDM符号的个数J,满足连续J个OFDM符号对应信道具有联合稀疏特性;3)对于连续J个OFDM符号,设计稀疏导频模式,推导信道估计模型;4)利用离散随机最优化算法求得最优导频位置分布;5)利用块同步正交匹配追踪(BSOMP)算法重建稀疏系数;6)根据BSOMP算法重建的系数恢复信道抽头系数;7)对估计得到的信道抽头系数进行分段线性光滑处理。本发明能够有效对抗时间频率双重选择性衰落,提高信道估计的准确度和频谱利用率。
申请公布号 CN105227505A 申请公布日期 2016.01.06
申请号 CN201510654644.7 申请日期 2015.10.10
申请人 上海交通大学 发明人 秦启波;宫博;归琳;罗汉文
分类号 H04L25/02(2006.01)I;H04L25/03(2006.01)I 主分类号 H04L25/02(2006.01)I
代理机构 上海新天专利代理有限公司 31213 代理人 张宁展
主权项 一种高速移动环境下的多符号联合信道估计方法,用于正交频分复用系统,利用多符号的联合稀疏特性,基于压缩感知理论对信道进行估计,其特征在于,该方法包括如下步骤:①利用复指数基扩展模型对时间频率双选信道进行建模,模型如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>h</mi><mi>l</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mrow><mo>(</mo><mrow><msub><mi>b</mi><mn>0</mn></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>b</mi><mrow><mi>Q</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow><mo>)</mo></mrow><mfenced open = "(" close = ")"><mtable><mtr><mtd><mrow><msup><mi>c</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msup><mrow><mo>&lsqb;</mo><mrow><mn>0</mn><mo>,</mo><mi>l</mi></mrow><mo>&rsqb;</mo></mrow></mrow></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr></mtable></mtd></mtr><mtr><mtd><mrow><msup><mi>c</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msup><mrow><mo>&lsqb;</mo><mrow><mi>Q</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>l</mi></mrow><mo>&rsqb;</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>+</mo><msubsup><mi>&xi;</mi><mi>l</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msubsup><mo>,</mo><mn>0</mn><mo>&le;</mo><mi>l</mi><mo>&le;</mo><mi>L</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000819089390000011.GIF" wi="1403" he="269" /></maths>其中,j代表OFDM符号的序号,Q代表CE‑BEM阶数,L代表路径数,<img file="FDA0000819089390000012.GIF" wi="84" he="87" />表示第l径对应的信道抽头系数,b<sub>q</sub>(0≤q≤Q‑1)代表CE‑BEM基函数,c<sup>(j)</sup>[q,l]代表CE‑BEM系数,<img file="FDA0000819089390000013.GIF" wi="86" he="87" />代表CE‑BEM建模误差;CE‑BEM基函数表示为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>b</mi><mi>q</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><mrow><mn>1</mn><mo>,</mo><mo>...</mo><mo>,</mo><msup><mi>e</mi><mrow><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>N</mi></mfrac><mi>n</mi><mrow><mo>(</mo><mi>q</mi><mo>-</mo><mfrac><mrow><mi>Q</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac><mo>)</mo></mrow></mrow></msup><mo>,</mo><mo>...</mo><mo>,</mo><msup><mi>e</mi><mrow><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>N</mi></mfrac><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mrow><mo>(</mo><mi>q</mi><mo>-</mo><mfrac><mrow><mi>Q</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac><mo>)</mo></mrow></mrow></msup></mrow><mo>)</mo></mrow><mi>T</mi></msup><mo>,</mo><mn>0</mn><mo>&le;</mo><mi>q</mi><mo>&le;</mo><mi>Q</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000819089390000014.GIF" wi="1438" he="188" /></maths>其中,N表示一个OFDM符号子载波的个数;②计算联合估计OFDM符号的个数J,满足连续J个OFDM符号对应信道具有联合稀疏特性,公式如下:J<0.01c/(N+L<sub>CP</sub>)v         (3)其中c是光速,v是移动终端和基站的相对速度,L<sub>CP</sub>代表循环前缀CP的长度;③对于连续J个OFDM符号,设计稀疏导频模式,推导信道估计模型:导频序列包括G个有效导频和(2Q‑2)G个保护导频,其中K<G<<JL,K表示无线信道延时域的稀疏度;有效导频幅度和保护导频幅度分别取1和0,有效导频序列记为P<sub>eff</sub>,保护导频序列记为P<sub>guard</sub>,保护导频位于有效导频两侧;将导频分为Q个子序列:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>P</mi><mn>0</mn></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>e</mi><mi>f</mi><mi>f</mi></mrow></msub><mo>-</mo><mfrac><mrow><mi>Q</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac></mrow></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr></mtable></mtd></mtr><mtr><mtd><mrow><msub><mi>P</mi><mfrac><mrow><mi>Q</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>e</mi><mi>f</mi><mi>f</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr></mtable></mtd></mtr><mtr><mtd><mrow><msub><mi>P</mi><mrow><mi>Q</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>e</mi><mi>f</mi><mi>f</mi></mrow></msub><mo>+</mo><mfrac><mrow><mi>Q</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000819089390000021.GIF" wi="1082" he="541" /></maths>其中<img file="FDA0000819089390000024.GIF" wi="98" he="108" />表示所有有效导频构成的子序列;得到信道估计模型如下:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mrow><mo>&lsqb;</mo><mi>Y</mi><mo>&rsqb;</mo></mrow><msub><mi>P</mi><mn>0</mn></msub></msub><mo>=</mo><msub><mi>&Phi;s</mi><mn>0</mn></msub><mo>+</mo><msub><mi>W</mi><mn>0</mn></msub></mrow></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr></mtable></mtd></mtr><mtr><mtd><mrow><msub><mrow><mo>&lsqb;</mo><mi>Y</mi><mo>&rsqb;</mo></mrow><msub><mi>P</mi><mfrac><mrow><mi>Q</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac></msub></msub><mo>=</mo><msub><mi>&Phi;s</mi><mfrac><mrow><mi>Q</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac></msub><mo>+</mo><msub><mi>W</mi><mfrac><mrow><mi>Q</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac></msub></mrow></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr></mtable></mtd></mtr><mtr><mtd><mrow><msub><mrow><mo>&lsqb;</mo><mi>Y</mi><mo>&rsqb;</mo></mrow><msub><mi>P</mi><mrow><mi>Q</mi><mo>-</mo><mn>1</mn></mrow></msub></msub><mo>=</mo><msub><mi>&Phi;s</mi><mrow><mi>Q</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>W</mi><mrow><mi>Q</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000819089390000025.GIF" wi="1083" he="464" /></maths>其中,Y代表接收机OFDM解调之后的数据,W<sub>q</sub>包括噪声和建模误差,<img file="FDA0000819089390000026.GIF" wi="153" he="115" />是需要重建的系数向量,测量矩阵Φ=[Φ<sub>0</sub>,...Φ<sub>L‑1</sub>],<img file="FDA0000819089390000027.GIF" wi="595" he="176" />④利用离散随机最优化算法求得最优导频位置分布,即最优导频位置确认原则是使测量矩阵Φ的互相关值μ最小,μ计算如下:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>&mu;</mi><mrow><mo>(</mo><mi>&Phi;</mi><mo>)</mo></mrow><mo>=</mo><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&NotEqual;</mo><mi>j</mi><mo>&le;</mo><mi>L</mi></mrow></munder><mfrac><mrow><mo>|</mo><mo>&lt;</mo><msub><mi>&phi;</mi><mi>i</mi></msub><mo>,</mo><msub><mi>&phi;</mi><mi>j</mi></msub><mo>&gt;</mo><mo>|</mo></mrow><mrow><mo>|</mo><mo>|</mo><msub><mi>&phi;</mi><mi>i</mi></msub><mo>|</mo><msub><mo>|</mo><mn>2</mn></msub><mo>|</mo><mo>|</mo><msub><mi>&phi;</mi><mi>j</mi></msub><mo>|</mo><msub><mo>|</mo><mn>2</mn></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000819089390000028.GIF" wi="1055" he="192" /></maths>其中φ<sub>i</sub>,φ<sub>j</sub>是矩阵Φ的任意两列;⑤利用块同步正交匹配追踪算法重建得到稀疏系数<img file="FDA0000819089390000029.GIF" wi="172" he="115" />⑥由稀疏系数恢复出信道抽头系数<img file="FDA00008190893900000210.GIF" wi="106" he="86" />恢复CE‑BEM系数<img file="FDA00008190893900000211.GIF" wi="106" he="91" />公式如下:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>c</mi><mi>q</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msubsup><mo>=</mo><msup><msub><mi>&Lambda;</mi><mi>q</mi></msub><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mrow><mo>&lsqb;</mo><msub><mi>s</mi><mi>q</mi></msub><mo>&rsqb;</mo></mrow><mrow><mi>j</mi><mo>:</mo><mi>J</mi><mo>:</mo><mi>J</mi><mi>L</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00008190893900000212.GIF" wi="1053" he="106" /></maths>其中对角矩阵<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>&Lambda;</mi><mi>q</mi></msub><mo>=</mo><mi>d</mi><mi>i</mi><mi>a</mi><mi>g</mi><mrow><mo>(</mo><mrow><mn>1</mn><mo>,</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>i</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>N</mi></mfrac><mrow><mo>(</mo><mi>q</mi><mo>-</mo><mfrac><mrow><mi>Q</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac><mo>)</mo></mrow></mrow></msup><mo>,</mo><mo>...</mo><mo>,</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>i</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>N</mi></mfrac><mrow><mo>(</mo><mi>q</mi><mo>-</mo><mfrac><mrow><mi>Q</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mrow><mo>(</mo><mi>L</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></msup></mrow><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000819089390000031.GIF" wi="914" he="175" /></maths>根据式(1)得到信道抽头系数<img file="FDA0000819089390000032.GIF" wi="102" he="85" />⑦对信道抽头系数进行分段线性光滑处理:对第j个符号,第l条径对应的信道抽头系数取平均,公式如下<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msubsup><mi>h</mi><mi>l</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo><mi>a</mi><mi>v</mi><mi>e</mi></mrow></msubsup><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><msubsup><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mi>j</mi><mrow><mo>(</mo><mi>N</mi><mo>+</mo><msub><mi>L</mi><mrow><mi>C</mi><mi>P</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>L</mi><mrow><mi>C</mi><mi>P</mi></mrow></msub></mrow><mrow><mo>(</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo><mo>(</mo><mi>N</mi><mo>+</mo><msub><mi>L</mi><mrow><mi>C</mi><mi>P</mi></mrow></msub><mo>)</mo><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>h</mi><mrow><mi>n</mi><mo>,</mo><mi>l</mi></mrow></msub><mo>,</mo><mn>0</mn><mo>&le;</mo><mi>l</mi><mo>&le;</mo><mi>L</mi><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000819089390000033.GIF" wi="1267" he="133" /></maths>分别利用前后两个相邻OFDM符号的信道抽头系数,对第j个OFDM符号对应的信道抽头系数进行光滑处理,<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msubsup><mi>h</mi><mi>l</mi><mrow><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mi>r</mi><mn>1</mn></mrow></msubsup><mo>=</mo><mfrac><mn>1</mn><mrow><mi>N</mi><mo>+</mo><msub><mi>L</mi><mrow><mi>C</mi><mi>P</mi></mrow></msub></mrow></mfrac><mrow><mo>(</mo><mrow><msubsup><mi>h</mi><mi>l</mi><mrow><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mi>a</mi><mi>v</mi><mi>e</mi></mrow></msubsup><mo>-</mo><msubsup><mi>h</mi><mi>l</mi><mrow><mrow><mo>(</mo><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>a</mi><mi>v</mi><mi>e</mi></mrow></msubsup></mrow><mo>)</mo></mrow><mrow><mo>(</mo><mrow><mi>n</mi><mo>+</mo><msub><mi>L</mi><mrow><mi>C</mi><mi>P</mi></mrow></msub><mo>+</mo><mn>1</mn><mo>+</mo><mfrac><mi>N</mi><mn>2</mn></mfrac></mrow><mo>)</mo></mrow><mo>+</mo><msubsup><mi>h</mi><mi>l</mi><mrow><mrow><mo>(</mo><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>a</mi><mi>v</mi><mi>e</mi></mrow></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>h</mi><mi>l</mi><mrow><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mi>r</mi><mn>2</mn></mrow></msubsup><mo>=</mo><mfrac><mn>1</mn><mrow><mi>N</mi><mo>+</mo><msub><mi>L</mi><mrow><mi>C</mi><mi>P</mi></mrow></msub></mrow></mfrac><mrow><mo>(</mo><mrow><msubsup><mi>h</mi><mi>l</mi><mrow><mrow><mo>(</mo><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow><mo>)</mo></mrow><mi>a</mi><mi>v</mi><mi>e</mi></mrow></msubsup><mo>-</mo><msubsup><mi>h</mi><mi>l</mi><mrow><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mi>a</mi><mi>v</mi><mi>e</mi></mrow></msubsup></mrow><mo>)</mo></mrow><mrow><mo>(</mo><mrow><mi>n</mi><mo>+</mo><mn>1</mn><mo>-</mo><mfrac><mi>N</mi><mn>2</mn></mfrac></mrow><mo>)</mo></mrow><mo>+</mo><msubsup><mi>h</mi><mi>l</mi><mrow><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mi>a</mi><mi>v</mi><mi>e</mi></mrow></msubsup></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000819089390000034.GIF" wi="1427" he="316" /></maths>对上述结果取平均得到第j个符号,第l条径对应的信道抽头系数如下:<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msubsup><mi>h</mi><mi>l</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><msubsup><mi>h</mi><mi>l</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo><mi>r</mi><mn>1</mn></mrow></msubsup><mo>+</mo><msubsup><mi>h</mi><mi>l</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo><mi>r</mi><mn>2</mn></mrow></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000819089390000035.GIF" wi="1110" he="131" /></maths>
地址 200240 上海市闵行区东川路800号