任意阵列的常规波束形成matlab程序

clear all
clc
close all
% n基元的数量
% b波束的指向，单位为角度
% fs抽样频率
% c波速
% f信号频率
% an为阵元点的位置

%以下为设置可以修改
N=1                         % N波束图的角度分辨力
T=0.01                      % T信号的长度
f=100;c=1500;b=-30;fs=5*f;an=[0 0;1 1; -1 2; 4 5;6  7;-3 -4;-7 -5;-2 -5]   % 信号初始信息



%以下为程序
t=0:1/fs:T;                 % 抽样时间
s0=exp（j*2*pi*f*t）;         % 原始信号
number=length（t）;           % 采样点的长度
[nv]=size（an）;             % n为阵元个数
S=zeros（nnumber）;          % 建立一个N行length（t）列的零矩阵用来表示收到的信号                       
P=zeros（nnumber）;          % 建立一个N行length（t）列的零矩阵用来表示处理后的信号
power=zeros（1（360/N）+1）;   % 用来放功率
for m=1:n
    bn（m:）=an（m:）-an（1:）;
end
cn（:1）=-sin（b*pi/180）*bn（:1）;
cn（:2）=-cos（b*pi/180）*bn（:2）;
dn=cn（:1）+cn（:2）;
tao=dn/c;                     % tao为各阵元点与第一个阵元点的时延
for m=1:n
    kn=tao（m1）;
    S（m:）=exp（j*2*pi*f*（t+kn））;