close allclear allclcc=1500;fs=10e3;T = 1;t = 0:1/fs:T;L=length(t);f=2000;w=2*pi*f;k=w/c;M=11; %阵元个数Nmid=1; %参考点d=3;%阵元间距m=0:1:M-1;yi=zeros(M,1);% 返回一个M*1维的零矩阵zi=zeros(M,1);xi=m*d;xi=xi.;%各阵元坐标y1=12;x1=10;z1=10;% 声源位置， y轴指向声源平面Ric1=sqrt(x1-xi).2+(y1-yi).2+(z1-zi).2);%声源至各阵元的距离 M*1维Rn1=Ric1-Ric1(Nmid);%声源至各阵元与参考阵元的声程差矢量 M*1维s1=cos(w*t);%参考阵元接收到的信号 1*L维snr =20;Am= 10(-snr/20);n1=Am*(randn(M,L)+j*randn(M,L);%各阵元噪声矢量p1=zeros(M,L)
