发明名称 一种频域超分辨率多径时延估计方法
摘要 本发明公开了一种频域超分辨率多径时延估计方法,包括如下步骤:(1)频域信道响应采样;(2)利用采样数据构造Hankle矩阵;(3)利用采样数据分块矩阵X1和X2构造矩阵E和矩EI;(4)构造噪声子空间、模式矢量;(5)搜索多径时延。本发明提供的频域超分辨率多径时延估计方法,在相同信噪比的情况下,具有比传统的多重信号分类(MUSIC)算法还高的时间分辨率,并且不再需要奇异值分解或是特征值分解等计算量较大的矩阵运算操作,让多径时延估计更加实用、简单。
申请公布号 CN103200136A 申请公布日期 2013.07.10
申请号 CN201310071995.6 申请日期 2013.03.07
申请人 东南大学 发明人 胡爱群;来;黄毅;姜禹;韩远致;邓小伟
分类号 H04L25/02(2006.01)I 主分类号 H04L25/02(2006.01)I
代理机构 南京瑞弘专利商标事务所(普通合伙) 32249 代理人 杨晓玲
主权项 1.一种频域超分辨率多径时延估计方法,其特征在于:包括如下步骤:第一步:频域信道响应采样在频域利用矢量网络分析仪对信道频域响应进行N点的均匀采样,其采样频率范围为[f<sub>L</sub>,f<sub>H</sub>],采样间隔记为△f=(f<sub>H</sub>-f<sub>L</sub>)/(N-1),接收到的数据采样点分别为R(0),R(1),...,R(N-1);其中f<sub>L</sub>为频率下界,f<sub>H</sub>为频率上界;第二步:利用采样数据构造Hankle矩阵将在频域采样得到的N个频域数据组织成为如下的L×(N-L+1)矩阵形式:<img file="FDA00002892571200011.GIF" wi="992" he="304" />其中,X<sub>1</sub>和X<sub>2</sub>是采样数据分块矩阵,且X<sub>1</sub>是D×(N-L+1)矩阵,X<sub>2</sub>是(L-D)×(N-L+1)矩阵,D是多径的传播路径数目,即:<img file="FDA00002892571200012.GIF" wi="951" he="299" /><img file="FDA00002892571200013.GIF" wi="1025" he="301" />第三步:利用采样数据分块矩阵X<sub>1</sub>和X<sub>2</sub>构造矩阵E和矩E<sub>I</sub><maths num="0001"><![CDATA[<math><mrow><mi>E</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mrow><mo>(</mo><msub><mi>X</mi><mn>1</mn></msub><msubsup><mi>X</mi><mn>1</mn><mi>H</mi></msubsup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>X</mi><mn>1</mn></msub><msubsup><mi>X</mi><mn>2</mn><mi>H</mi></msubsup></mtd></mtr><mtr><mtd><mo>-</mo><mi>I</mi></mtd></mtr></mtable></mfenced></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>E</mi><mi>I</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mrow><mo>(</mo><msup><msub><mi>X</mi><mn>1</mn></msub><mi>H</mi></msup><mo>)</mo></mrow><mo>+</mo></msup></mtd></mtr><mtr><mtd><mo>-</mo><msup><mrow><mo>(</mo><msup><msub><mi>X</mi><mn>2</mn></msub><mo>+</mo></msup><mo>)</mo></mrow><mi>H</mi></msup></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中,矩阵I是(L-D)×(L-D)单位矩阵;第四步:构造噪声子空间、模式矢量构造噪声子空间矩阵Q=E×E<sup>H</sup>或Q<sub>I</sub>=E<sub>I</sub>×E<sub>I</sub><sup>H</sup>以及模式矢量<maths num="0003"><![CDATA[<math><mrow><mi>e</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>L</mi></msub><mi>&tau;</mi></mrow></msup><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>L</mi></msub><mo>+</mo><mi>&Delta;f</mi><mo>)</mo></mrow><mi>&tau;</mi></mrow></msup><mo>.</mo><mo>.</mo><mo>.</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>L</mi></msub><mo>+</mo><mrow><mo>(</mo><mi>L</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>&Delta;f</mi><mo>)</mo></mrow><mi>&tau;</mi></mrow></msup></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>其中τ为传播路径的传播时延;第五步:搜索多径时延采用传统的MUSIC算法中的如下搜索方法来估计多径时延:<maths num="0004"><![CDATA[<math><mrow><mi>MU</mi><mrow><mo>(</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>L</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mrow><mo>|</mo><mo>|</mo><msup><mi>e</mi><mi>H</mi></msup><msub><mi>q</mi><mi>i</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup></mrow></mfrac></mrow></math>]]></maths>其中,q<sub>i</sub>是矩阵Q或者Q<sub>I</sub>的第i列。
地址 211189 江苏省南京市江宁区东南大学路2号