⛄一、获取代码方式
获取代码方式1:
完整代码已上传我的资源:【DOA估计】基于matlab常规波束和CAPON算法DOA估计【含Matlab译码 2661期】
点击上面蓝色字体,直接付费下载,即可。
获取代码方式2:
付费专栏Matlab信号处理(初级版)
备注:
点击上面蓝色字体付费专栏Matlab信号处理(初级版),扫描上面二维码,付费29.9元订阅海神之光博客付费专栏Matlab信号处理(初级版),凭支付凭证,私信博主,可免费获得1份本博客上传CSDN资源代码(有效期为订阅日起,三天内有效);
点击CSDN资源下载链接:1份本博客上传CSDN资源代码
⛄二、部分源代码
clear all;
%close all;
clc;
source_number=3;%信元数
sensor_number=16;%阵元数
N_x=1024; %信号长度
snapshot_number=N_x;%快拍数
w=[pi/4 pi/6 pi/3].';%信号频率
l=sum(2pi3e8./w)/3;%信号波长
d=0.5*l;%阵元间距
snr=15;%信噪比
source_doa=[-45 0 60];%两个信号的入射角度
A=[exp(-1j*(0:sensor_number-1)d2pisin(source_doa(1)pi/180)/l);exp(-1j(0:sensor_number-1)d2pisin(source_doa(2)pi/180)/l);exp(-1j(0:sensor_number-1)d2pisin(source_doa(3)*pi/180)/l)].';%阵列流型
s=sqrt(10.^(snr/10))exp(1jw*[0:N_x-1]);%仿真信号
%x=awgn(s,snr);
x=As+(1/sqrt(2))(randn(sensor_number,N_x)+1j*randn(sensor_number,N_x));%加了高斯白噪声后的阵列接收信号
R=xx’/snapshot_number;
iR=inv®;
%[V,D]=eig®;
%Un=V(:,1:sensor_number-source_number);
%Gn=UnUn’;
[U,S,V]=svd®;
Un=U(:,source_number+1:sensor_number);
Gn=Un*Un’;
⛄三、运行结果
⛄四、matlab版本及参考文献
1 matlab版本
2014a
2 参考文献
[1] 沈再阳.精通MATLAB信号处理[M].清华大学出版社,2015.
[2]高宝建,彭进业,王琳,潘建寿.信号与系统——使用MATLAB分析与实现[M].清华大学出版社,2020.
[3]王文光,魏少明,任欣.信号处理与系统分析的MATLAB实现[M].电子工业出版社,2018.
3 备注
简介此部分摘自互联网,仅供参考,若侵权,联系删除