发明名称 一种基于滤波器组的时延估计方法及其系统
摘要 本发明涉及一种基于滤波器组的时延估计的方法及其系统,利用滤波器组分别估计两路目标信号的互谱和相应的加权谱,以实现时延估计,1)采集两路目标信号的采样序列x1(n),x2(n);2)预先生成L×N维傅里叶矩阵F=[f0,f1,…,fN-1];3)分别计算协方差矩阵R11,R22和相应的各次逆矩阵R11-1,R22-1,以及互协方差矩阵R12;4)计算得到N个长度为L相应数据子段的滤波器组;5)计算序列x1(n),x2(n)的互谱 <mrow> <msubsup> <mi>P</mi> <mn>12</mn> <mi>p</mi> </msubsup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mi>g</mi> <mrow> <mn>1</mn> <mo>,</mo> <mi>i</mi> </mrow> <mi>H</mi> </msubsup> <msub> <mi>R</mi> <mn>12</mn> </msub> <msub> <mi>g</mi> <mrow> <mn>2</mn> <mo>,</mo> <mi>i</mi> </mrow> </msub> </mrow>和加权函数谱6)计算采用广义互谱相位法的时延估计。本方案中还可以进行数据分段和加窗,目的是保证数据序列的平稳性,提高估计的稳健度。本发明原理明确,易于工程实现;不仅适合窄带和功率谱平坦的情况,也适合功率谱有较大起伏的情况,在信道起伏导致功率谱不平坦的情况下,仍可实现稳健的时延估计。
申请公布号 CN101645701B 申请公布日期 2011.07.27
申请号 CN200910091453.9 申请日期 2009.08.25
申请人 中国科学院声学研究所;北京中科海跃科技有限责任公司 发明人 田彪;李宇;黄海宁
分类号 H03H17/00(2006.01)I 主分类号 H03H17/00(2006.01)I
代理机构 北京法思腾知识产权代理有限公司 11318 代理人 杨小蓉
主权项 1.一种基于滤波器组的时延估计方法,该方法利用滤波器组分别估计两路目标信号的互谱和相应的加权谱,以实现时延估计,所述的基于滤波器组的时延估计方法包括以下步骤:1)采集两路目标信号的采样序列x<sub>1</sub>(n),x<sub>2</sub>(n);所述的两路目标信号的到达接受传感器时存在时延差,即x<sub>1</sub>(n),x<sub>2</sub>(n)之间存在时延τ,采样序列x<sub>1</sub>(n),x<sub>2</sub>(n)的采样频率为F<sub>s</sub>;2)预先生成L×N维傅里叶矩阵F=[f<sub>0</sub>,f<sub>1</sub>,L,f<sub>N-1</sub>],傅里叶矩阵的列向量<img file="FSB00000458660500011.GIF" wi="669" he="143" />ω<sub>n</sub>=2πn/N,其中,n=0,1,L,N-1;Fs/N为频率分辨率;3)分别计算两路目标信号的采样序列x<sub>1</sub>(n),x<sub>2</sub>(n)的各子段相应的协方差矩阵R<sub>11</sub>,R<sub>22</sub>和相应的各次逆矩阵<img file="FSB00000458660500012.GIF" wi="214" he="59" />以及两路目标信号的采样序列x<sub>1</sub>(n)和x<sub>2</sub>(n)各子段相应的互协方差矩阵R<sub>12</sub>;所述协方差矩阵R<sub>11</sub>,R<sub>22</sub>和互协方差矩阵R<sub>12</sub>的计算公式为:<img file="FSB00000458660500013.GIF" wi="440" he="58" /><img file="FSB00000458660500014.GIF" wi="449" he="75" /><img file="FSB00000458660500015.GIF" wi="443" he="59" />其中上标H表示共轭转置操作,E{·}表示求数学期望;4)设定参数α,由设定参数α、协方差矩阵估计及其逆矩阵,得到相应的N个长度为L相应数据子段的滤波器组:<img file="FSB00000458660500016.GIF" wi="465" he="63" /><img file="FSB00000458660500017.GIF" wi="470" he="65" />其中i=0,1,L,N-1;f<sub>i</sub>为频率点;所述参数α的取值范围为0≤α≤1;所述滤波器<img file="FSB00000458660500018.GIF" wi="54" he="67" />和<img file="FSB00000458660500019.GIF" wi="60" he="57" />分别表示第p个子段数据对应的滤波器系数;5)计算序列x<sub>1</sub>(n),x<sub>2</sub>(n)的互谱<img file="FSB000004586605000110.GIF" wi="358" he="67" />和加权函数谱<img file="FSB000004586605000111.GIF" wi="140" he="82" />6)计算采用广义互谱相位法的时延估计。
地址 100190 北京市海淀区北四环西路21号中国科学院声学研究所