【DOA估计】基于Matlab实现L阵型MUSIC算法 |
您所在的位置:网站首页 › 二次雷达代码含义 › 【DOA估计】基于Matlab实现L阵型MUSIC算法 |
✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,matlab项目合作可私信。 🍎个人主页:Matlab科研工作室 🍊个人信条:格物致知。 更多Matlab仿真内容点击👇 智能优化算法 神经网络预测 雷达通信 无线传感器 电力系统 信号处理 图像处理 路径规划 元胞自动机 无人机 ⛄ 内容介绍该文提出一种基于MUSIC算法的L型阵列多输入多输出雷达降维波达方向(DOA)估计算法.该算法首先针对L型阵列导向矢量的结构,构造出一个降维矩阵,将回波信号转换到低维空间.然后利用二次优化方法将2维DOA估计分解为两个1维DOA估计.最后利用MUSIC空间谱估计. ⛄ 代码%阵列信号处理MATLAB编程 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(-1i*twpi*d.'*(sin(theta*rad).*cos(fe*rad))); Ay=exp(-1i*twpi*d1.'*(sin(theta*rad).*cos(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(-1i*twpi*d.'*sin(phim1)*cos(phim2)); a2=exp(-1i*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)') zlabel('magnitude(dB)') ⛄ 运行结果 ⛄ 参考文献[1] 王伟, 王晓萌, 李欣,等. 基于MUSIC算法的L型阵列MIMO雷达降维DOA估计[J]. 电子与信息学报, 2014, 36(8). [2] 张琦. 基于时空二维MUSIC算法的DOA估计及其硬件实现方法研究[D]. 西安电子科技大学. [3] 常文静. 基于MUSIC算法的相干信号DOA估计研究[D]. 河南理工大学. ⛳️ 代码获取关注我❤️部分理论引用网络文献,若有侵权联系博主删除 ❤️ 关注我领取海量matlab电子书和数学建模资料
|
今日新闻 |
点击排行 |
|
推荐新闻 |
图片新闻 |
|
专题文章 |
CopyRight 2018-2019 实验室设备网 版权所有 win10的实时保护怎么永久关闭 |