实例介绍
【实例简介】
【实例截图】
【实例截图】
【核心代码】
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 | clear all;close all;clc;source_number=2;%信元数% sensor_number=13;%阵元数 sensor_A=[0 4 10 13 17 23]; %52A天线布局% sensor_A=[0 1 6 13 14 19]; %52A天线布局% sensor_A=[0 2 6 9 11 15]; %博世天线布局% sensor_A=[0 1 6 19 20 25]; %均匀天线布局 % sensor_A=[0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23]; %互质天线布局 sensor_number=length(sensor_A);N_x=1; %信号长度snapshot_number=N_x;%快拍数w=[pi/4 pi/4].';%信号频率l=((2*pi*3e8)/w(1) (2*pi*3e8)/w(2))/2;%信号波长 d=0.5*l;%阵元间距snr=26;%信噪比j=sqrt(-1);for n=1:1:100source_doa=[0 15];%两个信号的入射角度% A=[1*exp(-j*(0:sensor_number-1)*d*2*pi*sin(source_doa(1)*pi/180)/l);1*exp(-j*(0:sensor_number-1)*d*2*pi*sin(source_doa(2)*pi/180)/l)].';%阵列流型;以来波方向为source_doa（k）入射% %M根天线，阵元间距为d的均匀线阵的阵列响应矢量A=[1*exp(-j*sensor_A*d*2*pi*sin(source_doa(1)*pi/180)/l);0.1*exp(-j*sensor_A*d*2*pi*sin(source_doa(2)*pi/180)/l)].';s=sqrt(10.^(snr/10))*exp(j*w*(0:N_x-1));%仿真信号x=awgn(s,snr);xx=A*s;xxx=abs(xx);% plot(xxx);x=A*s (1/sqrt(2))*(randn(sensor_number,N_x) j*randn(sensor_number,N_x));%加了高斯白噪声后的阵列接收信号R=x*x'/snapshot_number;searching_doa=[-20:0.5:20]';for m=1:1:length(searching_doa) for i=1:1:length(searching_doa) th=[searching_doa(m) searching_doa(i)];% th=[35,40]; A1=exp(-j*sensor_A'*d*2*pi*sin(th.*pi/180)/l); p=A1*pinv(A1'*A1)*A1'; xx=real(diag(p*R)); ang(m,i)=sum(xx);% ang=sum(xx); endend[kk1,kk2]=find(ang==max(max(ang)));%寻找最大值所在位置kk1=searching_doa(kk1);%寻找最大值对应的角度坐标kk2=searching_doa(kk2);%寻找最大值对应的角度坐标k1(n)=max(kk1);k2(n)=min(kk2);endfigure(2);mesh(searching_doa,searching_doa,ang);figure(1);plot(k1,'y');hold on;plot(k2,'r'); |
好例子网口号：伸出你的我的手 — 分享！
小贴士
感谢您为本站写下的评论，您的评论对其它用户来说具有重要的参考价值，所以请认真填写。
- 类似“顶”、“沙发”之类没有营养的文字，对勤劳贡献的楼主来说是令人沮丧的反馈信息。
- 相信您也不想看到一排文字/表情墙，所以请不要反馈意义不大的重复字符，也请尽量不要纯表情的回复。
- 提问之前请再仔细看一遍楼主的说明，或许是您遗漏了。
- 请勿到处挖坑绊人、招贴广告。既占空间让人厌烦，又没人会搭理，于人于己都无利。
关于好例子网
本站旨在为广大IT学习爱好者提供一个非营利性互相学习交流分享平台。本站所有资源都可以被免费获取学习研究。本站资源来自网友分享，对搜索内容的合法性不具有预见性、识别性、控制性，仅供学习研究，请务必在下载后24小时内给予删除，不得用于其他任何用途，否则后果自负。基于互联网的特殊性，平台无法对用户传输的作品、信息、内容的权属或合法性、安全性、合规性、真实性、科学性、完整权、有效性等进行实质审查；无论平台是否已进行审查，用户均应自行承担因其传输的作品、信息、内容而可能或已经产生的侵权或权属纠纷等法律责任。本站所有资源不代表本站的观点或立场，基于网友分享，根据中国法律《信息网络传播权保护条例》第二十二与二十三条之规定，若资源存在侵权或相关问题请联系本站客服人员，点此联系我们。关于更多版权及免责申明参见 版权及免责申明
网友评论
我要评论