clear all

close all

derad = pi/180; % deg -> rad

radeg = 180/pi;

twpi = 2*pi;

kelm = 8; % 阵列数量

dd = 0.5; % space

d=0:dd:(kelm-1)*dd; %

iwave = 4; % number of DOA

theta = [-60 -30 30 60]; % 角度

snr = 10; % input SNR (dB)

n = 500; %

A=exp(-j*twpi*d.'*sin(theta*derad));%%%% direction matrix

S=randn(iwave,n);

X=A*S;

X1=awgn(X,snr,'measured');%在信号X中加入高斯白噪声,信噪比SNR,'measured'函数在加入噪声前测定信号强度

Rxx=X1*X1'/n;

InvS=inv(Rxx); %%%%

[EV,D]=eig(Rxx);%%%% [V,D]=eig(A):求矩阵A的全部特征值,构成对角阵D,并求A的特征向量构成V的列向量。

EVA=diag(D)';%返回矩阵D的主对角线上的元素

[EVA,I]=sort(EVA);%其中I是一个大小等于size(EVA)的数组,其每一列是EVAA中列向量的元素相对应的置换位置记号。

EVA=fliplr(EVA);%将矩阵A的列绕垂直轴进行左右翻转,如果A是一个行向量,fliplr(A)将A 中元素的顺序进行翻转。如果A是一个列向量,fliplr(A)还等于A。

EV=fliplr(EV(:,I));

% MUSIC

foriang = 1:361

angle(iang)=(iang-181)/2;

phim=derad*angle(iang);

a=exp(-j*twpi*d*sin(phim)).';

L=iwave;

En=EV(:,L+1:kelm);

SP(iang)=(a'*a)/(a'*En*En'*a);

end

%画图

SP=abs(SP);

SPmax=max(SP);

SP=10*log10(SP/SPmax);

h=plot(angle,SP);

set(h,'Linewidth',2)

music的matlab程序,MUSIC算法matlab程序相关推荐

  1. 视频教程-三十八课时零基础matlab精通优化算法-Matlab

    三十八课时零基础matlab精通优化算法 图像和算法等领域有多年研究和项目经验:指导发表科技核心期刊经验丰富:多次指导数学建模爱好者参赛. 宋星星 ¥100.00 立即订阅 扫码下载「CSDN程序员学 ...

  2. pid算法matlab仿真程序和c程序,pid算法matlab仿真程序和c程序.doc

    pid算法matlab仿真程序和c程序.doc 1 增量式 PID 控制算 法 Matlab 仿真程序 设一被控对象 G(s) =50/(0.125s^2+7s),用增量式 PID 控制 算法编写仿真 ...

  3. 极值滤波matlab,极值滤波算法MATLAB程序及处理结果对比

    <极值滤波算法MATLAB程序及处理结果对比>由会员分享,可在线阅读,更多相关<极值滤波算法MATLAB程序及处理结果对比(5页珍藏版)>请在人人文库网上搜索. 1.极值滤波算 ...

  4. bfgs算法matlab程序,bfgs算法matlab代码

    (对 Large -scale 问题) 对应文件 \\toolbox\\matlab\\funfun\\fminbnd.m \\toolbox\\optim\\sfminbx.m \\toolbox\ ...

  5. idw matlab 程序_IDW 算法MATLAB 实现 -

    中国Unix/Linux软件开发联盟 http://www.lisdn.com IDW 算法MATLAB 实现 linux软件开发 %IDW(反距离加权插值法) %其中x,y,z为已知坐标及其函数值, ...

  6. idw matlab 程序_IDW 算法MATLAB 实现

    IDW 算法MATLAB 实现 linux软件开发 %IDW(反距离加权插值法) %其中x,y,z为已知坐标及其函数值,X,Y为要插值的坐标 %x,y,z,X,Y最高为二维的,不可为三维 %不考虑x, ...

  7. matlab 图像白平衡算法,Matlab常用白平衡算法

    <Matlab常用白平衡算法>由会员分享,可在线阅读,更多相关<Matlab常用白平衡算法(21页珍藏版)>请在人人文库网上搜索. 1.1 灰色世界法灰色世界法(grey wo ...

  8. 用matlab实现理查森外推算法,Matlab数值积分(2)

    实验目的: 掌握理查森外推法 实验要求: 1. 给出理查森外推算法 2. 用Matlab实现理查森外推算法 3. 用Matlab实现自适应积分算法 实验内容: 1. 理查森外推算法,数学知识:利用Ri ...

  9. 压缩感知算法matlab,压缩感知算法matlab

    MATLAB代码如下: 压缩感知仿真实例 压缩感知仿真实例在MATLAB R2001b中的计算结果如下: 原图像 采样率0.7 采样率0.5 采样率0.3 压缩感知仿真实例采用均方误差...... M ...

  10. 关于精英蚁群算法matlab,蚁群算法MATLAB解VRP问题

    Excel  exp12_3_2.xls内容: ANT_VRP函数: function [R_best,L_best,L_ave,Shortest_Route,Shortest_Length]=ANT ...

最新文章

  1. 折腾Java设计模式之责任链模式
  2. CCNA实验之---单臂路由实现VLAN间路由
  3. php代码优化 -- array_walk 和 foreach, for 的效率的比较
  4. C++CLR类库封装Native类库并用C#调用
  5. 命令行安装Pillow
  6. 六爪机器人_六爪机器人
  7. SQL Server 索引中的碎片和填充因子
  8. UVM-sequence
  9. Arduino使用ESP8266安装问题(包括附加开发板管理网址)
  10. 在Markdown文件中快速插入本地图片
  11. 计算机核心期刊(北大核心)投稿意见
  12. 矩形中分出正方形,递归优化
  13. Connection id “0HMGNTPDIIA4M“.Request id “0HMGNTPDIIA4M:00000001“:An unhandled exception was thrown
  14. C++入门编程之成语接龙
  15. 一个菜鸟实习生的月总结
  16. 解决Cleartext HTTP traffic to 01.minipic.eastday.com not permitted
  17. win7无法保存打印机设置(错误0x000006d9)解决方法
  18. week03-面向对象
  19. 【研究方法】好的研究想法从哪里来--刘知远
  20. Python:随机生成手机号码

热门文章

  1. Eclipse安装Lombok
  2. 一个app管理平台-app-host的搭建
  3. 股票交易数据下载 | 下载股票历史交易数据到本地Excel
  4. java如何开发安卓软件_Java 程序员如何快速入门 Android 开发?
  5. 多智能体MAPPO代码环境配置以及代码讲解
  6. 细数实现流水灯的三种代码形式
  7. IndexedDB封装
  8. 圈圈USB开发板 IDE40
  9. delphi 组件安装教程详解
  10. Java Web学习视频