clear
close all
clc
twpi = 2*pi;
rad = pi/180;
deg = 180/pi;
kelm = 8;
snr = 10;
iwave = 3;
theta = [10 30 50];
fe = [15 25 35];
n = 100;
dd = 0.5;
d = 0:dd:(kelm-1)*dd;
d1 = dd:dd:(kelm-1)*dd;
Ax = exp(-j*twpi*d.'*(sin(theta*rad).*cos(fe*rad)));
Ay = exp(-j*twpi*d1.'*(sin(theta*rad).*sin(fe*rad)));
A = [Ax;Ay];
S = randn(iwave,n);
X = A*S;
X1 = awgn(X,snr,'measured');
Rxx = X1*X1'/n;
[EV,D] = eig(Rxx);
[EVA,I] = sort(diag(D).');
EV = fliplr(EV(:,I));
Un = EV(:,iwave+1:end);
for ang1 = 1:90
for ang2 = 1:90
thet(ang1) = ang1-1;
phim1 = thet(ang1)*rad;
f(ang2) = ang2-1;
phim2 = f(ang2)*rad;
a1 = exp(-j*twpi*d.'*sin(phim1)*cos(phim2));
a2 = exp(-j*twpi*d1.'*sin(phim1)*sin(phim2));
a = [a1;a2];
SP(ang1,ang2) = 1/(a'*Un*Un'*a);
end
end
SP=abs(SP);
SPmax=max(max(SP));
SP=SP/SPmax;
h = mesh(thet,f,SP);
%set(h,'Linewidth',2)
xlabel('elevation(degree)')
ylabel('azimuth(degree)')
axis([0 90 0 90 0 1])
zlabel('magnitude(dB)')
Matlab实现L形阵列的二维DOA估计
版权申诉
166 浏览量
2022-06-17
11:45:29
上传
评论 2
收藏 73KB ZIP 举报
天天Matlab科研工作室
- 粉丝: 2w+
- 资源: 7257
最新资源
- 电子设计竞赛的单相不间断电源设计
- cutcamera1715961370938.png
- 基于MATLAB的图像处理课程设计报告.doc
- tensorflow-gpu-2.6.0-cp38-cp38-manylinux2010-x86-64.whl
- mmexport1715960553858.png
- tensorflow-gpu-2.6.0-cp37-cp37m-manylinux2010-x86-64.whl
- 通过 .NET 应用程序中的源代码查找 SQL 注入
- 电子设计竞赛2007年B题 无线识别装置.doc
- Wox全局搜索工具,一款win下的全局搜索软件
- 使用高级集群管理 (ACM) 模板来管理用户、组和命名空间
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈