⛄一、获取代码方式
获取代码方式1:
完整代码已上传我的资源:【方位估计】基于matlab music算法均匀园阵方位估计【含Matlab源码 2668期】
点击上面蓝色字体,直接付费下载,即可。
获取代码方式2:
付费专栏Matlab物理应用(初级版)
备注:
点击上面蓝色字体付费专栏Matlab物理应用(初级版),扫描上面二维码,付费29.9元订阅海神之光博客付费专栏Matlab物理应用(初级版),凭支付凭证,私信博主,可免费获得1份本博客上传CSDN资源代码(有效期为订阅日起,三天内有效);
点击CSDN资源下载链接:1份本博客上传CSDN资源代码
⛄二、部分源代码
%%%%%%%% MUSIC for Uniform Circle Array%%%%%%%%
r=1; % 半径(m)
N=16; % 阵元数目
d=2rsin(pi/N);
M=1; % 信源数量
gamma=2pi/N(0:N-1); % 和参考阵元的角度
fc=16000; % 采样率
c=340; % 声音的速度(m/s)
lambda=c/fc;
a_theta=45; % 仰角
a_phi=90; % 方位角
zeta=2pi/lambdarsin(a_thetapi/180);
A=exp(1izetacos((a_phi-gamma)pi/180)).'; % 导向矢量
K=100; % 快拍数
t=(0:K-1)/1000;
% 构建信号模型
S=sin(2pifct);
X=AS;
snr=10;
X1=awgn(X,snr,‘measured’);
% 计算协方差矩阵
R=X1X1’/K;
% 特征值分解
[Q,D]=eig®;
[D,I]=sort(diag(D),1,‘descend’);
Q=Q(:,I);
Qn=Q(:,M+1:N); % 信号子空间
theta=0:90;
phi=0:1:360;
p_MUSIC=zeros(length(theta),length(phi));
for ii=1:length(theta)
for iii=1:length(phi)
zeta=2pi/lambdarsin(theta(ii)pi/180);
A=exp(1izetacos((phi(iii)-gamma)pi/180)).‘;
p_MUSIC(ii,iii)=(1/(A’(Qn*Qn’)*A));
end
end
⛄三、运行结果
⛄四、matlab版本及参考文献
1 matlab版本
2014a
2 参考文献
[1] 门云阁.MATLAB物理计算与可视化[M].清华大学出版社,2013.
3 备注
简介此部分摘自互联网,仅供参考,若侵权,联系删除