参考《声呐技术》第四章，基本的线阵形成
波束扫描技术，（以及电阻相移波束形成）主要适用于窄带信号，即信号带宽与中心频率相比很小的系统。（P100）
clc;close all;clear all;
%%系统参数
c0=1500; %声速
f=1000; %发射信号的中心频率
fs=100000; %采样频率 d=0.2; %系统实孔径0.2
M = 16; % 阵元数
k = 0.5; % k = d/lamada = 1/2
% 形成波束指向doa方向的权值
doa = 30; % doa=0，为静态方向图%doa为线阵法线方向与声源的夹角
m = [0 : M-1]; %阵元序号
w = exp(j*2*pi*k*m'*sin(doa*pi/180)); %每个阵元的时延tow=k*m‘*sin（）不同
% 波束扫描，绘制方向图
theta = linspace(-90,90,360); % 扫描方位角
a = exp(-j*2*pi*k*m'*sin(theta*pi/180)); %theta是扫描角度% 导向适量
Y = abs(w'*a);
Y = 20*log10(Y/max(Y));
subplot(211)
plot(theta, Y); title('未加权方向图'); axis([-90,90,-50,0]); grid on; xlabel('方位角/(\circ)');ylabel( 'P/(dB)');
%%===========================================================================
doa = 30; % doa=0，为静态方向图
m = [0 : M-1];
w = exp(-j*2*pi*k*m'*sin(doa*pi/180));
w = chebwin(M, 30) .* w; % 加-30db切比雪夫权
% 波束扫描，绘制方向图
theta = linspace(-90,90,360); % 方位角
a = exp(j*2*pi*k*m'*sin(theta*pi/180)); % 导向适量
Y = abs(w'*a);
Y = 20*log10(Y/max(Y));
figure; plot(theta, Y); title('加-30dB切比雪夫权方向图'); axis([-90,90,-50,0]); grid on; xlabel('方位角/(\circ)');ylabel( 'P/(dB)');
