主权项 |
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>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<M<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>θ</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>θ</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" />内,对空间谱进行谱峰搜索,谱峰位置即为估计出的信号波达方向,谱峰个数即为信号个数。 |