麦克风阵列方向图模型

本节讨论宽带声源信号模型下的信号模型。
考虑一个白噪声源$s$（它涵盖了整个频谱）,其方差为${\sigma }_s^2=1$。
在这一信号模型中，考虑第n个传感器为
$$y_n(k)=s\left[k-(n-1){\tau }_s)\right],n=1,2,\cdots ,N式（1.1）$$
对于远场信号${\tau }_s=\frac{d\sin \theta }{c}$，为信号到达相邻麦克风的时延。式（1.1）可以写为向量形式

$$y_n(k)=\vec{g}\left[(n-1){\tau }_s\right]\vec{s}(k),n=1,2,\cdots ,N式（1.2）$$
其中
$$\vec{g}\left[(n-1)\cdot {\tau }_s\right]={\left[\begin{array}{lllllll}0& \cdots & 0& 1_{（n-1）\cdot {\tau }_s}& 0& \cdots & 0\end{array}\right]}_{L_g}^{\mathrm{T}}$$
为长度为$L_g⩾(N-1)\cdot {\tau }_s+1$的滤波器，其中第$(n-1)\cdot {\tau }_s+1$个元素等于1.
其中
$$\vec{s}(k)={\left[s(k)s(k-1)\cdots s\left[k-(n-1){\tau }_s\right]\cdots s\left(k-L_g+1\right)\right]}_{L_g}^{\mathrm{T}}$$

现在考虑长度为$L_h$的宽带波束形成器，其中第$(n-1)\cdot \tau +1$个元素等于1.
$$\vec{h}\left[(n-1)\cdot \tau \right]=\frac{1}{N}{\left[\begin{array}{lllllll}0& \cdots & 0& 1_{（n-1）\cdot \tau }& 0& \cdots & 0\end{array}\right]}_{L_h}^{\mathrm{T}}$$
则波束形成后的输出$$\begin{array}{rl}z(k)& =\sum _{n=1}^N{\vec{h}}^{\mathrm{T}}\left[(N-n)\tau \right]{\vec{y}}_n(k)\\ & =\frac{1}{N}\sum _{n=1}^N{y}_n\left[k-(N-n)\tau \right]\\ & =\frac{1}{N}\sum _{n=1}^{N}s\left[k-(N-n)\tau -(n-1){\tau }_s\right]\end{array}$$
其中
$${\vec{y}}_n(k)={\left[y_n(k)y_n(k-1)\cdots y_n\left(k-L_h+1\right)\right]}^{\mathrm{T}},n=1,2,\cdots ,N式(1.3)$$
为对信号的进行$L_h$采样后的向量。
对于 $\tau ={\tau }_s$,$z(k)=s\left[k-(N-1){\tau }_s\right]$。

1.如何理解上面得波束形成器？
信号从到达第一个麦克风，到到达最后一个麦克风结束。此过程中第一个麦克风接收到的信号持续了了$N\cdot {\tau }_s$，第二个麦克风接收到的信号持续了$(N-1){\tau }_s$…第N个麦克风接收到的信号持续了了$1\cdot {\tau }_s$。，所以将第一个麦克风信号右移$（N-1）\cdot {\tau }_s$对应$y_1[k-(N-1){\tau }_s]$，将第二个麦克风信号右移$(N-2){\tau }_s$，对应$y_2[k-(N-2){\tau }_s]$…最后一个麦克风不动。这样N个麦克风的接收信号经过右移后就会在时间上对齐，这时候再相加除以N就可以得到滤波后的结果。
2.与传统Delay-Sum的区别
传统DS滤波器是在让接收信号与第一个麦克风接收信号对齐，而上面的波束形成方法是与最后一个麦克风接受信号对齐。

上式也可以表示成如下形式：$$z(k)={\mathbf{H}}^{\mathrm{T}}(\tau )\mathbf{Y}(k)$$
其中
$$\mathbf{H}(\tau )={\left[{\vec{h}}^{\mathrm{T}}\left[(N-1)\tau \right]{\vec{h}}^{\mathrm{T}}\left[(N-2)\tau \right]\cdots {\vec{h}}^{\mathrm{T}}\left[0\cdot \tau \right]\right]}_{N\cdot L_h}^{\mathrm{T}}$$
$$\mathbf{Y}(k)={\left[{\vec{y}}_1^{\mathrm{T}}(k){\vec{y}}_2^{\mathrm{T}}(k)\cdots {\vec{y}}_N^{\mathrm{T}}(k)\right]}_{N\cdot L_h}^{\mathrm{T}}$$
将式（1.3）写成矩阵乘以向量的形式。
$${\vec{y}}_n(k)=\mathbf{G}\left[(n-1){\tau }_s\right]{\vec{s}}_L(k),n=1,2,\cdots ,N$$
式中$$\mathbf{G}\left[(n-1){\tau }_s\right]=\left[\begin{array}{ccccc}{g}^{\top }\left[(n-1){\tau }_s\right]& 0& 0& \cdots & 0\\ 0& g^{\mathrm{T}}\left[(n-1){\tau }_s\right]& 0& \cdots & 0\\ ⋮& ⋮& ⋮& \ddots & ⋮\\ 0& 0& 0& \cdots & g^{\mathrm{T}}\left[(n-1){\tau }_s\right]\end{array}\right],n=1,2,\cdots ,N$$
是大小为$L_h\times L$的Sylvester矩阵，$L=L_g+L_h-1$
Sylvester矩阵第二行是第一行往右移一列，第二行第一列的元素是零。下面的$(L_h-1)$行也是用这种方法得出，每次都往右移一列。
$${\vec{s}}_L(k)=[s(k-t)s(k-1)\cdots s(k-L+1){]}^{\mathrm{T}}$$
则$z(k)$可以写为：
$$z(k)=H^{\mathrm{T}}(\tau )G\left(N{\tau }_s\right)s_L(k)$$
其中
$$\mathbf{G}\left(N{\tau }_s\right)={\left[\begin{array}{c}\mathbf{G}\left[0\cdot {\tau }_s\right]\\ \mathbf{G}\left[1\cdot {\tau }_s\right]\\ ⋮\\ \mathbf{G}\left[(N-1)\cdot {\tau }_s\right]\end{array}\right]}_{N{L}_h\times L}$$
可以看作方向矩阵，包含了有关目标信号位置的所有信息。因此波束形成输出的方差为$$E\left[z^2(k)\right]={\Vert G^{\mathrm{T}}\left(N\cdot {\tau }_s\right)H(\tau )\Vert }_2^2$$定义麦克风阵列的模式函数为：$$\mathcal{A}(\tau )=1-{\Vert {\mathbf{G}}^{\mathrm{T}}\left(N\cdot {\tau }_s\right)H(\tau )-\mathbf{u}\left[(N-1)\cdot \tau \right]\Vert }_2$$
其中
$$u\left[(N-1)\tau \right]=[0\cdots 010\cdots 0{]}^{\top }$$
是长度为L的矢量，其中第$(N-1)\cdot \tau$个元素等于1。表示滤波的结果是逼近第N个麦克风接收的信号。