信号基础知识---线阵

参考《声呐技术》第四章，基本的线阵形成

波束扫描技术，（以及电阻相移波束形成）主要适用于窄带信号，即信号带宽与中心频率相比很小的系统。（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)');