阅读原文还请移步我的知乎专栏:
https://www.zhihu.com/column/c_1287066237843951616
继上篇文章,本篇继续学习,包括内容为:
均匀线列阵常规波束形成原理概述
例3.4 均匀线阵列常规波束图
采用下图所示坐标系统:
图1 均匀线列阵坐标系统
对于该坐标系统,假设由 个均匀分布的阵元组成的线列阵,假设阵元间距为 ,则线列阵总长度为 。这里计算线阵列长度时,将两端阵元向外各延伸了 ,该均匀线列阵相当于对原连续阵列进行了空间采样。
各阵元的位置可以表示为:
运用每个阵元的加权系数表示将其视作连续线阵时的加权函数,即:
其中, 是第个阵元的加权系数; 是 函数。
函数为:
在之前的篇章讲述过,假设一个放置在 轴的连续线阵的频率-波数响应为:
将上述公式替换为 轴的一个均匀分布的线阵列,即为:
改为波束响应的形式为:
还是基于此均匀线列阵坐标系统,阵列的流行项链为:
波束响应可表达为:
对该均匀线列阵进行常规波束形成,假设波束指向角为 ,波束加权向量为:
代入上式可计算出波束方向响应为:
例3.4 均匀线阵列常规波束图
假设波束观测方向为 ,计算常规波束形成获得的波束响应。假设阵元数目为 ,即阵元间隔 。采用上述计算的波束形成公式计算波束响应。
图2 均匀线阵列常规波束图
图2 均匀线阵列常规波束图
另附实现代码如下:
c=340; %声速
theta_d = 0*pi/180; %入射角度
f=1000; %频率
space=c/f/2; %麦克风间距
M=10; %麦克风数量
theta_angle=0:0.1:360;
theta=theta_angle*pi/180;
B=sin((M*pi*f*space*(sin(theta)-sin(theta_d)))/c)...
./(M*sin((pi*f*space*(sin(theta)-sin(theta_d)))/c));
B_db = 20*log10(B);
limit_dB = -50;
index = B_db < limit_dB;
B_db(index) = limit_dB;
plot(theta_angle, B_db, 'linewidth', 1.5);
grid on;
title('均匀线阵列常规波束响应');
xlabel('\theta/(\circ)');ylabel('20lg|B(\theta)|/dB');
figure;
GraphicHandle = polar(theta, B_db);
set( GraphicHandle, 'LineWidth', 1.5);
参考书籍:
《优化阵列信号处理》