Files
josh-blog-source/source/_posts/2023-04-20-josh-opa-part-1-2.md
T
2023-04-20 10:12:07 +08:00

96 lines
3.4 KiB
Markdown
Raw Blame History

This file contains ambiguous Unicode characters
This file contains Unicode characters that might be confused with other characters. If you think that this is intentional, you can safely ignore this warning. Use the Escape button to reveal them.
---
title: 'Josh''s Notes: 最优阵列处理(Part 1.2 —— 均匀线阵模型)'
mathjax: true
comments: true
copyright: true
toc:
enable: false
number: true
wrap: true
expand_all: true
max_depth: 4
tags:
- 阵列信号处理
- 均匀线阵
- ULA
categories:
- - Josh 的学习笔记
- 最优阵列处理
- 阵列和空域滤波器
abbrlink: de20fd09
date: 2023-04-20 00:14:31
---
我们感兴趣的线阵如[图 1-2-1](#fig.1-2-1) 所示,其中有 $N$ 个位于 $z$ 轴上的阵元,阵元间距均匀,记为 $d$。我们把这种阵列称为**{% label primary @均匀线性阵列(uniform linear array, ULA) %}**。我们已经把阵列的中心放在了坐标系的原点上,这种阵列中心放置方法将在计算量上具有优点。
<a id="fig.1-2-1"></a>
![图 1-2-1 放置在 $z$ 轴上的线阵](https://josh-blog-1257563604.cos.ap-beijing.myqcloud.com/img/2023-04-20-josh-opa-part-1-2/2023-04-20-josh-oap-part-1-2-010-LinearArrayAlongZAaxis.png!sign){width=600px}
<!-- more -->
&emsp;&emsp;阵元的位置为
$$\begin{equation} \label{ElementLocation}
\begin{cases}
p_{x_n} = 0 \\
p_{y_n} = 0 \\
p_{z_n} = \left( n - \frac{N-1}{2} \right) d
\end{cases}
, \, n = 0,1,\cdots,N-1
\end{equation}$$
为了确定阵列流形矢量 $\vec{v}_{\vec{k}}\left(\vec{k}\right)$,把式 $\eqref{ElementLocation}$ 代入到[阵列流形矢量的定义][],得到
$$\begin{equation} \label{ManifoldVectorOfULA}
\vec{v}_{\vec{k}}\left(k_z\right) = \left[\begin{array}{c}
e^{j \left( \frac{N-1}{2} \right) k_z d} \\
e^{j \left( \frac{N-1}{2} - 1\right) k_z d} \\
\vdots \\
e^{-j \left( \frac{N-1}{2} \right) k_z d} \\
\end{array}\right]
\end{equation}$$
这里
$$\begin{equation}
k_z = -\frac{2\pi}{\lambda}\cos\theta = -k_0 \cos\theta
\end{equation}$$
其中波数的幅度
$$\begin{equation}
k_0 \triangleq \left| \vec{k} \right| = \frac{2\pi}{\lambda}
\end{equation}$$
意线阵在 0 方向是没有分辨能力的。将[复权矢量][] 和式 $\eqref{ManifoldVectorOfULA}$ 代入[用复权矢量表示的波数-频率响应函数][],得到 ULA 在角度域的频率-波数响应函数
$$\begin{equation} \label{WFRFOfULAInAngleDomain}
\begin{aligned}
\vec{\varUpsilon}(\omega,k_z) &= \vec{w}^\mathrm{H}\vec{v}_{\vec{k}}\left(k_z\right) \\
&= \sum_{n=0}^{N-1} w_n^\ast e^{-j\left( n - \frac{N-1}{2} k_z d \right)}
\end{aligned}
\end{equation}$$
定义下面的变量也是很有用的(波数在 $z$ 轴上的投影,乘以阵元间距,也即扫过一个阵元间距时,波数的变化量在 $z$ 轴上的体现):
$$\begin{equation} \eqref{DefinitionOfWavenumberDiffOfTwoElements}
\psi = -k_z d = \frac{2\pi}{\lambda}\cos\theta\cdot d = \frac{2\pi}{\lambda}u_z d
\end{equation}$$
其中 $u_z$ 是针对 $z$ 轴的方向余弦,
$$\begin{equation}
u_z = \cos\theta
\end{equation}$$
将式 $\eqref{DefinitionOfWavenumberDiffOfTwoElements}$ 代入式 $\eqref{WFRFOfULAInAngleDomain}$,得到
$$\begin{equation}
\vec{\varUpsilon}_\psi(\psi) = e^{-j \frac{N-1}{2} \psi } \sum_{n=0}^{N-1} w_n^\ast e^{j n\psi}
\end{equation}$$
[阵列流形矢量的定义]: https://josh-gao.top/posts/8b61f5a7.html#DefinitionOfArrayManifoldVector
[复权矢量]: https://josh-gao.top/posts/8b61f5a7.html#ComplexWeightVector
[用复权矢量表示的波数-频率响应函数]: https://josh-gao.top/posts/8b61f5a7.html#FrequencyWavenumberResponseWithComplexWeigh