二维DOA算法

clc;
clear all;

Pd=2000;
Fd=1;
Fs=4*Fd;
R=0.5;
Delay=5;
M=4;
x1=randint（Pd1M）; 
x2=randint（Pd1M）;
x3=randint（Pd1M）;
y1=modmap（x1FdFs‘qask‘M）;
y2=modmap（x2FdFs‘qask‘M）;
y3=modmap（x3FdFs‘qask‘M）;
[rcv_a1ti]=rcosflt（y1FdFs‘fir/sqrt/Fs‘RDelay）;
[rcv_a2ti]=rcosflt（y2FdFs‘fir/sqrt/Fs‘RDelay）;
[rcv_a3ti]=rcosflt（y3FdFs‘fir/sqrt/Fs‘RDelay）;
s1=amodce（rcv_a110‘qam‘）; 
s2=amodce（rcv_a210‘qam‘）;
s3=amodce（rcv_a310‘qam‘）;
save sig3 s1 s2 s3;

m=8; % 阵元数
p=3; % 信号数
%俯仰角为入射方向与圆阵轴线之间的夹角方位角为x轴到入射线在圆阵上投影的夹角
st1=10; % 俯仰角
dt1=10; % 方位角
st2=20;
dt2=20;
st3=30;
dt3=30;

st=[st1;st2;st3];
dt=[dt1;dt2;dt3];

nn=4096; 
SN1=20; 
SN2=20;
SN3=20;
sn=[SN1;SN2;SN3];
degrad=pi/180;

load sig3