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. php 继承内核中的基类,php – 从基类调用继承类的父方法
  2. 最新版Xshell 6安装教程详解
  3. Git复习(八)之快速理解Git结构
  4. ie下js不执行的几种可能
  5. linux 应用层gpio中断_如何在嵌入式Linux中检测GPIO线路上的中断?
  6. 第五章functions.py中的交叉熵代码解释
  7. cron表达式 每天0点10分和30分_揭开考研阅卷的内幕,注意这些多得20分!
  8. 活动子项父项的复杂CSS选择器[重复]
  9. python字符串格式化符号含义及转义字符含义
  10. CSDN博客专家申请的条件及规则
  11. 只用十行 Python 代码就提取了韦小宝的身份证信息
  12. 微软云中国的服务器在哪,微软云计算平台Windows Azure将落户中国
  13. JavaScript高级04 正则表达式
  14. 亚马逊AVS开发--树莓派3麦克风配置
  15. KDA,token.Kadena链的一些常识
  16. ArcGIS10从入门到精通系列实验图文教程(附配套实验数据持续更新)
  17. Wilson定理推论
  18. HTML5视频标签 video 的 poster 属性
  19. tmux系列:3、tmux窗口和窗格
  20. 实体间的联系存在着( )_医药零售的网络销售与实体渠道如何协同

热门文章

  1. python vba 区别_VBA和Python该学哪个?
  2. 软考软件设计师下午真题-面向对象的程序设计与实现-装饰设计模式(2012年上半年试题六))Java代码讲解
  3. 简易语音助手—python
  4. 软考高级 真题 2016年上半年 信息系统项目管理师 综合知识
  5. 数字信号处理----数字上变频和下变频都FPGA实现
  6. 万能打印之Delphi 2010实现(完结)
  7. 从斐波那契数列讲解算法设计的思路
  8. Ubuntu16.04安装QQ(for linux)
  9. 100个常用的 JS 代码
  10. 老罗Android开发 视频教程