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