时延波束形成 matlab,有关波束形成程序解释

clc;

clear;

%%  参数

c=1500;

f0=1.5e3;

fs=150e3;

T=50e-3;

SNR=20;

lamda=c./f0;

D=lamda./2;

N=21;

maxdelay=(N-1).*D./c;

Nmaxdelay=round(maxdelay.*fs);

alpha=(-90:90)*pi./180;

px=(0:N-1).*D;

theta0=30*pi./180;

N_delay0=D.*sin(theta0)./c.*fs;

theta0=-30*pi./180;

N_delay1=D.*sin(theta0)./c.*fs;

Wn=2*[1000 2000]/fs;

bb=fir1(128,Wn,'bandpass');                          %带通滤波器系数

%% 信号模型

t=0:1./fs:T;

sig00=sqrt(2).*cos(2.*pi.*f0.*t);

for mm=1:N;

N0=round(N_delay0.*(mm-1));

sig0(mm,:)=[zeros(1,Nmaxdelay.*2-N0) sig00 zeros(1,Nmaxdelay.*1+N0)];

end

for mm=1:N;

N0=round(N_delay1.*(mm-1));

s