⛄一、获取代码方式
获取代码方式1:
完整代码已上传我的资源:【DOA估计】基于matlab MVDR高分辨方位估计【含Matlab源码 4189期】
点击上面蓝色字体,直接付费下载,即可。
获取代码方式2:
付费专栏Matlab信号处理(初级版)
备注:
点击上面蓝色字体付费专栏Matlab信号处理(初级版),扫描上面二维码,付费29.9元订阅海神之光博客付费专栏Matlab信号处理(初级版),凭支付凭证,私信博主,可免费获得1份本博客上传CSDN资源代码(有效期为订阅日起,三天内有效);
点击CSDN资源下载链接:1份本博客上传CSDN资源代码
⛄二、部分源代码
%方位估计实验 MVDR
clear all;
close all;
%导入数据
load(‘data_831.mat’);
x0(1,:) = Data1_AI_0___U’;
x0(2,:) = Data1_AI_1___U’;
x0(3,:) = Data1_AI_2___U’;
x0(4,:) = Data1_AI_3___U’;
x0(5,:) = Data1_AI_4___U’;
x0(6,:) = Data1_AI_5___U’;
x0(7,:) = Data1_AI_6___U’;
x0(8,:) = Data1_AI_7___U’;
x0(9,:) = Data1_AI_8___U’;
x0(10,:) = Data1_AI_9___U’;
x0(11,:) = Data1_AI_10___U’;
x0(12,:) = Data1_AI_11___U’;
%圆弧基阵信息
N = 12;
R = 0.24;
BETA = 30pi/180;
%上述基阵尺寸近似满足lambda=d/2
%声信号信息
F = 5700; %5.7khz信号
FS = 200000; %200kHz采样
T = 100/F; %100个CYCLE
T0 = 1; %1s周期
NS = TFS; %脉宽内的采样点
C = 1450; %声速
THETA_D = 0.1; %多波束的波束指向角增量,必须是360的约数
THETA_T = (0:THETA_D:360-THETA_D)*pi/180; %假定信号来得方向
%% 定义空矩阵
w = ones(N,length(THETA_T));
R_M=zeros(N,N);%协方差矩阵
j=sqrt(-1);
a=ones(N,length(THETA_T));
P=ones(1,length(THETA_T)); %功率响应即方位谱图
%% 数据降采样,归一化,数据归一化和有指向性阵元不能同时用
for i = 1:N
x(i,:)=resample(double(x0(i,:)),50,200);
x(i,:)=x(i,:)/max(x(i,:));
end
M = length(x(1,:));
figure(2);
for i=1:N
subplot(3,4,i)
plot(x(i,:));
axis([0 M-1 -5 5])
title(strcat(num2str(i),‘号阵元’));
end
⛄三、运行结果
⛄四、matlab版本及参考文献
1 matlab版本
2014a
2 参考文献
[1] 沈再阳.精通MATLAB信号处理[M].清华大学出版社,2015.
[2]高宝建,彭进业,王琳,潘建寿.信号与系统——使用MATLAB分析与实现[M].清华大学出版社,2020.
[3]王文光,魏少明,任欣.信号处理与系统分析的MATLAB实现[M].电子工业出版社,2018.
3 备注
简介此部分摘自互联网,仅供参考,若侵权,联系删除