%%%%%%%%%%%%%%%基于线性调频信号的MUSIC算法
clear;
close all;
j=sqrt(-1);
c=3e8;
tao=10e-6;%采样间隔
fs=100e6;%采样频率
M=8;
f1_l=29.9e6;
f1_h=30.1e6;%信号的低频和高频
f2_l=30.1e6;
f2_h=29.9e6;%此种情况下,中频相同且不相关,且为窄带信号
f1_c=(f1_l+f1_h)/2;%中频,且相同
w0=2*pi*f1_c;
B1=f1_h-f1_l;%带宽
kc1=B1/tao;%调制斜率
kc2=-B1/tao;
N_fs=fix(tao*fs);% 采样点数
t0=(0:N_fs-1)/fs;
snr=10;
theta=[-10 20];%入射角度
N_source=length(theta);
lambda=c/f1_c;%波长及阵元间距
d=lambda/2;
dd=(0:M-1)*d;
mecnt=100;
F = fs/N_fs;
f_xlabe = ((0:N_fs-1)-N_fs/2)*F;
%%%%%%%%%%%流行矩阵
delay=zeros(M,N_source);
for i=1:N_source
delay(:,i)=dd.'*sin(theta(i)*pi/180)/c;
end
% A=zeros(M,N_source);
% for i2=1:N_source
% A(:,i2)=exp(j*2*pi*dd.'*sin(theta(i)*pi/180)/c);
% delay(:,i2)=dd.'*sin(theta(i2)*pi/180)/c
% end
%%%%%%%%%%%产生线性调频信号
X1=zeros(M,N_fs);
X2=zeros(M,N_fs);
for h1=1:M;
t1=t0-delay(h1,1);
X1(h1,:)=exp(j*2*pi*(f1_l*t1+1/2*kc1*t1.^2));
end
for h2=1:M;
t2=t0-delay(h2,2);
X2(h2,:)=exp(j*2*pi*(f2_l*t2+1/2*kc2*t2.^2));
end
X=X1+X2;
X_out=awgn(X,snr,'measured');
mesh(1:1000,1:8,abs(X_out));
%%%%%%%music算法————————————————————
R=(X_out*(X_out)')/N_fs;
[edv,d]=eig(R);
d=diag(d);%取出对角矩阵
[d_sort pre]=sort(abs(d),'ascend');%对特征值d进行从小到大排列并取出其索引
edv_sort = edv(:,pre); %对应特征值从小到大的特征向量
Un = edv_sort(:,1:M-N_source); %取出噪声的特征向量
%对角度进行搜索
ang = (-90:1:90);
p = zeros(1,length(ang));
for loop = 1:length(ang) %-90~90度搜索
a = exp(-j*w0*dd.'*sin(ang(loop)*pi/180)/c);
tem = a'*Un;
p(1,loop) = 1/(tem*tem');
end
p_max = max(p);
p=10*log10(p/p_max);
figure;
plot(ang,p);
xlabel('angle(degree)');
axis([-90 90 -60 5]);
ylabel('magnitude(db)');
title('MUSIC算法')
grid on;
没有合适的资源?快使用搜索试试~ 我知道了~
DOA估计基于matlab music算法宽带信号(线性调频信号)DOA估计【含Matlab源码 2677期】.zip
共3个文件
jpg:2个
m:1个
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 192 浏览量
2023-06-07
10:28:51
上传
评论 1
收藏 77KB ZIP 举报
温馨提示
CSDN海神之光上传的代码均可运行,亲测可用,直接替换数据即可,适合小白; 1、代码压缩包内容 主函数:MUSIC.m; 调用函数:其他m文件;无需运行 运行结果效果图; 2、代码运行版本 Matlab 2019b;若运行有误,根据提示修改;若不会,可私信博主; 3、运行操作步骤 步骤一:将所有文件放到Matlab的当前文件夹中; 步骤二:双击打开MUSIC.m文件; 步骤三:点击运行,等程序运行完得到结果; 4、仿真咨询 如需其他服务,可私信博主或扫描博主博客文章底部QQ名片; 4.1 CSDN博客或资源的完整代码提供 4.2 期刊或参考文献复现 4.3 Matlab程序定制 4.4 科研合作 功率谱估计: 故障诊断分析: 雷达通信:雷达LFM、MIMO、通信、成像、定位、干扰、检测、信号分析、脉冲压缩 滤波估计:SOC估计 目标定位:WSN定位、滤波跟踪、目标定位 生物电信号:肌电信号EMG、脑电信号EEG、心电信号ECG 通信系统:DOA估计、编码译码、变分模态分解、管道泄漏、滤波器、数字信号处理+传输+分析+去噪、数字信号调制、误码率、信号估计、信号检测信号识别融合
资源推荐
资源详情
资源评论
收起资源包目录
【DOA估计】基于matlab music算法宽带信号(线性调频信号)DOA估计【含Matlab源码 2677期】.zip (3个子文件)
【DOA估计】基于matlab music算法宽带信号(线性调频信号)DOA估计【含Matlab源码 2677期】
运行结果1.jpg 60KB
MUSIC.m 2KB
运行结果2.jpg 28KB
共 3 条
- 1
资源评论
海神之光
- 粉丝: 3w+
- 资源: 2091
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- 基于matlab实现文档+程序边缘计算任务卸载与资源调度的算法,是论文的源代码,具有价值.rar
- 什么是学生成绩管理系统c++以及学习学生成绩管理系统的意义
- 什么是词向量-以及学习关于了解词向量的意义
- 什么是mybatis动态sql以及学习mybatis动态sql的意义
- 华为数据治理方法论,包括:数据治理框架、数据治理组织架构、数据治理度量评估体系以及华为数据治理案例分享
- 基于matlab实现对表面肌电信号进行归一化处理,并对归一化后的图形显示 .rar
- 基于matlab实现单级倒立摆的 T-S 模型 包括 LMI 程序源码
- 图书管理系统(struts+hibernate+spring+ext).rar
- 基于matlab实现此压缩包包含语音信号处理中的语音变声代码加音频.rar
- STM32使用PWM驱动舵机并通过OLED显示
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功