Beamforming

最新推荐文章于 2024-09-23 17:33:56 发布

373955482

最新推荐文章于 2024-09-23 17:33:56 发布

阅读量1.8w

点赞数 25

分类专栏： Microphone Array 音频处理数字信号处理文章标签： DOA SRP-PHAT microphone array circular array 宽带波束形成

本文链接：https://blog.csdn.net/u010592995/article/details/81289461

版权

数字信号处理同时被 3 个专栏收录

29 篇文章 45 订阅

订阅专栏

音频处理

19 篇文章 19 订阅

订阅专栏

Microphone Array

12 篇文章 36 订阅

订阅专栏

1. 阵列简介

1.1 均匀线阵

假设0、1、…M-1阵元的均匀线阵中，参考阵元接收的信号为 $s (t)$ ，则第 $m$ 个阵元接收到的信号为 $s(t-\tau_m)$ ，傅里叶变换如下

$S_m(\omega)=\int_{-\infty }^{+\infty} s_m(t)e^{-j\omega t}dt$

　　 $=\int_{-\infty }^{+\infty} s(t-\tau_m)e^{-j\omega t}dt$
　　　
$=S(\omega)e^{-j\omega\tau_m}$

则阵列接收的信号可以表示为

$X(\omega)=S(\omega)[1 e^{-j\omega\tau_１} e^{-j\omega\tau_２} … e^{-j\omega\tau_{M-1}} ]$

其中，定义

$\vec{a}=[1 e^{-j\omega\tau_１} e^{-j\omega\tau_２} … e^{-j\omega\tau_{M-1}} ]$

为阵元的响应向量（array response vector），也称作阵列流形(array manifold vector)。
同时，也定义波数（wave number）如下
$\kappa=\frac{\omega}{c} \vec{a}$
其中 $\omega=2\pi kf_s/N_{fft}$ 为数字角频率
　　在delauy-and-sum beamformer中，每个阵元的权值即为该阵元相对参考阵元需要补偿的相位延迟，因此权值即为阵列流形的共轭转置
　　根据上面的分析，写出窄带线性阵列的delay-and-sum波束图代码如下

delay-and-sum beamformer

M=8;                 % 阵元数
c=340;                         % 声速
f = 2125;                      % 信号频率
d = 0.08;                      % 阵元间距
theta=linspace(-pi/2,pi/2,200);% 入射信号角度范围
theta0=30*pi/180;                      % 注视方向
w=exp(1j*2*pi*f*sin(theta0)*[0:M-1]'*d/c);       % 导向向量
p = zeros(length(1:length(theta)),1);
for  j=1:length(theta)                                     % 角度扫描
    a=exp(-1j*2*pi*f*sin(theta(j))*[0:M-1]'*d/c);% 入射信号方向向量
    p(j)=sum(w.*a)/M;                                        % 延时-求和
end
% 画图
figure;
plot(theta/pi*180,abs(p)),grid on
xlabel('degree')
ylabel('amplitude')
title('8阵元均匀线阵方向图')

方向图如下
在这里插入图片描述

1.2 均匀圆阵

跟上面的分析相同，关键地方就是写出阵列流形向量和权值向量
　　假设 $M$ 阵元均匀圆阵，半径为 $R$ ，阵元在 $x - y$ 平面上，则各个阵元相对 $x$ 轴的夹角 $\gamma_m=2\pi m/M$ ，均匀圆阵中，一般以圆心为参考，要写出阵列方向向量，就是要计算当单位平面波以夹角 $(\phi，\theta)$ 入射时，各阵元相对圆心接收到的信号时间差。

计算上图¹

中的距离 $d$ ，就要计算这个余弦角 $\beta$ ，由三余弦定理（查百度才知道叫这个名字，逃~）有

$cos(\beta)=cos(\phi-\gamma_m)sin(\theta)$

然后，就可以计算得到距离差,得到时间差

$\tau_m=Ｒcos(\phi-\gamma_m)sin(\theta)／ｃ$

最后，根据上面线阵的相同分析，可以直接写出均匀圆阵的的阵列流形向量如下

$\vec{a}=\begin{bmatrix}e^{-j\omega Ｒcos(\phi-\gamma_0)sin(\theta)/c} \\ e^{-j\omega Ｒcos(\phi-\gamma_1)sin(\theta)/c} \\ ... \\ e^{-j\omega Ｒcos(\phi-\gamma_{M-1})sin(\theta)/c} \end{bmatrix}$

然后，写出均匀圆阵的delay-and-sum波束图就简单了

M=8;                 % 阵元数
c=340;                         % 声速
f = 2125;                      % 信号频率
r = 0.08;                      % 阵元间距
phi=linspace(0,2*pi,360);% 入射信号角度范围

theta = 80*pi/180; %固定一个俯仰角
gamma = (0:360/M:360-360/M)*pi/180;%麦克风位置

tao = r*sin(theta)*cos(angle(1)-gamma)/c;     %方位角 0 < angle <360

phi0=60*pi/180;                      % 注视方向
w=exp(1j*2*pi*f*r*cos(phi0-gamma)*sin(theta)/c);       % 导向向量
p = zeros(length(1:length(phi)),1);
for  j=1:length(phi)                                     % 角度扫描
    a=exp(-1j*2*pi*f*r*cos(phi(j)-gamma)*sin(theta)/c);% 入射信号方向向量
    p(j)=sum(w.*a)/M;                                        % 延时-求和
end
% 画图
figure;
plot(phi/pi*180,abs(p)),grid on
xlabel('degree')
ylabel('amplitude')
title('均匀圆阵方向图')

%画极坐标图
figure,polarplot(phi,abs(p))

在这里插入图片描述

另外,圆阵是面阵,可以进行二维空间搜索,简单改下程序,画出三维图

M=8;                 % 阵元数
c=340;                         % 声速
f = 2125;                      % 信号频率
r = 0.08;                      % 阵元间距
phi=linspace(0,2*pi,360);% 入射信号方位角范围

theta = linspace(0,pi/2,90);% 入射信号俯仰角范围

gamma = (0:360/M:360-360/M)*pi/180;%麦克风位置


phi0=130*pi/180;                      % 注视方向
theta0 = 80*pi/180;

w=exp(1j*2*pi*f*r*cos(phi0-gamma)*sin(theta0)/c);       % 导向向量
p = zeros(length(1:length(phi)),length(theta));
for i = 1:length(theta)
    for  j=1:length(phi)                                     % 角度扫描
        a=exp(-1j*2*pi*f*r*cos(phi(j)-gamma)*sin(theta(i))/c);% 入射信号方向向量
        p(j,i)=sum(a.*w)/M;                                        % 延时-求和
    end
end
% 画图

figure,mesh(abs(p))

注视方向在(130,80)时的阵列方向图如下所示

在这里插入图片描述