波束形成matlab程序

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

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



%以下为程序
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（180/N）+1）;   % 用来放功率
coef=sqrt（10^（SNR/10））; 
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