发明名称 一种基于空域自适应滤波的波达方向估计方法
摘要 一种基于空域自适应滤波的波达方向估计方法是针对均匀线阵,实现步骤为:1)接收数据进行K次采样,则参考阵元接收信号为d,辅助阵元接收信号为x(1),......,x(K);2)所有数据进行正交变换得解析信号<img file="DDA00002289303400011.GIF" wi="46" he="58" /><img file="DDA00002289303400012.GIF" wi="95" he="57" />......,<img file="DDA00002289303400013.GIF" wi="120" he="57" />3)对解析信号采用基于功率倒置的空域自适应滤波算法进行M次迭代计算,最后一次迭代计算出最优权值w(M);4)把所有N+1个阵元视为一个新的阵列,构造一个新的权值矢量为w<sub>e</sub>及方向矢量v<sub>e</sub>(θ);5)定义空间谱p(θ),在可见区间<img file="DDA00002289303400014.GIF" wi="259" he="123" />内,对空间谱进行谱峰搜索,完成波达方向估计。该方法实现了窄带信号的DOA估计,无需进行信源数估计和矩阵特征分解,计算量小,分辨率高,实现简单。
申请公布号 CN102866383B 申请公布日期 2014.01.22
申请号 CN201210404468.8 申请日期 2012.10.22
申请人 重庆大学 发明人 曾浩;王娅;赵静;凤林锋;刘陆军
分类号 G01S3/00(2006.01)I 主分类号 G01S3/00(2006.01)I
代理机构 代理人
主权项 1.一种基于空域自适应滤波的波达方向估计方法,其实现步骤为:第一步,由N+1个阵元构成的均匀线阵中,阵元间距l为半个波长,第一个阵元为参考阵元,其余N个阵元为辅助阵元,对阵元接收的窄带信号进行K次采样,采样频率为f<sub>s</sub>,则参考阵元接收信号为d=[d(1) ... d(K)],辅助阵元接收信号为<maths num="0001"><![CDATA[<math><mrow><mi>x</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mn>1</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>x</mi><mi>N</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>......,<maths num="0002"><![CDATA[<math><mrow><mi>x</mi><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>x</mi><mi>N</mi></msub><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>其中N和K都为自然数,K&gt;2N;第二步,对所有接收数据矢量进行正交变换,得到相应的解析信号<maths num="0003"><![CDATA[<math><mrow><mover><mi>d</mi><mo>^</mo></mover><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mover><mi>d</mi><mo>^</mo></mover><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mover><mi>d</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mover><mi>x</mi><mo>^</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><msub><mover><mi>x</mi><mo>^</mo></mover><mi>N</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>......,<maths num="0005"><![CDATA[<math><mrow><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mover><mi>x</mi><mo>^</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><msub><mover><mi>x</mi><mo>^</mo></mover><mi>N</mi></msub><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>第三步,对解析信号采用基于功率倒置的空域自适应滤波算法,进行M次迭代计算求权值,第k次迭代权值为w(k),收敛后的最后一次迭代计算出最优权值w(M),该权矢量为N行列矢量,其中,迭代次数满足k&lt;M&lt;K,k=1,2,…,M-1;第四步,根据得到的最优权值矢量w(M),把所有N+1个阵元视为一个新的阵列,构造一个新的权值矢量为<maths num="0006"><![CDATA[<math><mrow><msub><mi>w</mi><mi>e</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mo>-</mo><mi>w</mi><mrow><mo>(</mo><mi>M</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>及方向矢量<maths num="0007"><![CDATA[<math><mrow><msub><mi>v</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mi>v</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>v(θ)为入射信号进入N个辅助阵元的方向矢量;第五步,定义空间谱<img file="FDA0000372397430000018.GIF" wi="392" he="148" />在可见区间<img file="FDA0000372397430000019.GIF" wi="276" he="145" />内,对空间谱进行谱峰搜索,谱峰位置即为估计出的信号波达方向,谱峰个数即为信号个数。
地址 400044 重庆市沙坪坝区沙正街174号
您可能感兴趣的专利