matlab实现圆阵MVDR波束形成

%%%%%%非相干信号，矢量阵，
clc;
close all; 
clear all; 
M=16;           %阵元数

m=[0:1:M-1];
N=1000;      %采样点数
n=[1:1:N];
f=1000;         %频率
fs=5000;        %采样频率
c=1500;         %声速
lamd=c/f;       %波长
r=2*lamd;

gamam=m*2*pi/M;
Pxn = r*cos（gamam）;
Pyn = r*sin（gamam）;
APos = [Pxn‘ Pyn‘ ];  %%各阵元位置矢量


sita1=90*pi/180;    %目标1方位
gama1=[cos（sita1）sin（sita1）]‘;
tao1=APos*gama1/c;

sita2=95*pi/180;    %目标2方位
gama2=[cos（sita2）sin（sita2）]‘;
tao2=APos*gama2/c;

SN=10;  %信噪比

s1=exp（j*（2*pi*f*n/fs+2*pi*rand（1N）））; 
s2=exp（j*（2*pi*f*n/fs+2*pi*rand（1N）））; 

a1=exp（-j*2*pi*f*tao1）; 
a2=exp（-j*2*pi*f*tao2）; 

u1=[1 cos（sita1） sin（sita1）].‘;
u2=[1 cos（sita2） sin（sita2）].‘;

Av1=kron（a1u1）  ;
Av2=kron（a2u2）  ;

 %s=[s1‘ ]‘;  %%单目标
%A=[a1 ];
%Av=[Av1 ];

s=[s1‘