程序

doas=[-30 -5 40]*pi/180; %DOA’s of signals in rad.

P=[1 1 1]; %Power of incoming signals

N=10; %Number of array elements

K=1024; %Number of data snapshots

d=0.5; %Distance between elements in wavelengths

noise_var=1; %Variance of noise

r=length(doas); %Total number of signals

% Steering vector matrix. Columns will contain the steering vectors

% of the r signals

A=exp(-i*2*pi*d*(0:N-1)'*sin([doas(:).']));

% Signal and noise generation

sig=round(rand(r,K))*2-1; % Generate random BPSK symbols for each of the

% r signals

noise=sqrt(noise_var/2)*(randn(N,K)+i*randn(N,K)); %Uncorrelated noise

X=A*diag(sqrt(P))*sig+noise; %Generate data matrix

R=X*X'/K; %Spatial covariance matrix

[Q ,D]=eig(R); %Compute eigendecomposition of covariance matrix

[D,I]=sort(diag(D),1,'descend'); %Find r largest eigenvalues

Q=Q (:,I); %Sort the eigenvectors to put signal eigenvectors first

Qs=Q(:,1:r); %Get the signal eigenvectors

Qn=Q(:,r+1:N); %Get the noise eigenvectors

% MUSIC algorithm

% Define angles at which MUSIC “spectrum” will be computed

angles=(-90:0.1:90);

%Compute steering vectors corresponding values in angles

a1=exp(-i*2*pi*d*(0:N-1)'*sin([angles(:).']*pi/180));

for k=1:length(angles)

%Compute MUSIC “spectrum”

music_spectrum(k)=(a1(:,k)'*a1(:,k))/(a1(:,k)'*Qn*Qn'*a1(:,k));

end

figure(1)

plot(angles,abs(music_spectrum))

title('MUSIC Spectrum')

xlabel('Angle in degrees')

matlab显示峰值点横坐标,怎么求出峰值对应横坐标问题,急相关推荐

  1. matlab转换为部分分式,用MATLAB的residuez函数,求出下列各式的部分分式展开式和反变换f[n]。 说明:residuez函数...

    用MATLAB的residuez函数,求出下列各式的部分分式展开式和反变换f[n]. 说明:residuez函数 更多相关问题 需求价格弹性反映需求量与价格之间的敏感程度.( ) (2.0分) 实验( ...

  2. matlab根据末端位姿p求出xyz,Matlab在智能采摘机械手中的应用

    Matlab在智能采摘机械手中的应用 Matlab在智能采摘机械手中的应用 --RB03机器人的运动空间分析及其避障规划仿真 在Matlab中用数值法对广州数控RB03机器人进行工作空间的分析,然后根 ...

  3. matlab 固有频率 振型,请教如何求出有阻尼系统的固有频率及相应振型矩阵

    用矩阵形式表示方程:M*X''+C*X'+K*X=F(t) 方程实际上是六自由度系统的振动响应方程,就是说X是一个6×1的矩阵,M,C,K为质量,阻尼,刚度矩阵.而且C,K不是对角阵,各个自由度之间有 ...

  4. matlab怎么求两个数的和,matlab怎么求出两个函数的交点

    matlab中,两个自变量的函数怎么求最大值(急!1) 需求:利用matlab求解二元函数y=f(x1,x2)=(339-0.01*x1-0.003*x2)*x1+(399-0.004*x1-0.01 ...

  5. c语言编程序按下列公式计算e的值,编写程序,利用公式 ,求出e的近似值

    已知一个三角形中三条边的长度分别为a,b和c,编写程序利用公式求出三角形的面积,其中S=(a+b+c)/2.要求边长a. usingSystem;usingSystem.Collections.Gen ...

  6. matlab中多项式拟合如何给出r方,matlap拟合函数后r^2怎么求

    matlab拟合函数 求助 尝试用二次多项式拟合: clearall; x=2:2:20; y=[0.31.232.416.267.958.529.049.179.299.37]; scatter(x ...

  7. Java黑皮书课后题第7章:7.9(找出最小元素)使用下面的方法头编写一个方法,求出一个整数数组中的最小元素。编写测试程序,提示用户输入10个数字,调用这个方法返回最小值,并显示这个最小值

    7.9(找出最小元素)使用下面的方法头编写一个方法,求出一个整数数组中的最小元素.编写测试程序,提示用户输入10个数字,调用这个方法返回最小值,并显示这个最小值 题目 题目描述与运行示例 破题 代码 ...

  8. Java黑皮书课后题第2章:2.12(物理:求出跑道长度)编写程序,提示用户输入以米/秒为单位的速度v和加速度a,然后显示最短跑道长度

    2.12(物理:求出跑道长度)编写程序,提示用户输入以米/秒为单位的速度v和加速度a,然后显示最短跑道长度 题目 题目描述 运行示例 补充 代码块 题目 题目描述 2.12(物理:求出跑道长度) 假设 ...

  9. Matlab 散点 拟合 曲率,有数据点,希望得到一条拟合曲线,再求出这条曲线的曲率,求助!...

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 得到X,Y的数据点后,我尝试用origin7.5的analysis-fit sigmoidal功能拟合出曲线,由于方程太复杂,用公式计算很困难,我想用or ...

最新文章

  1. 《从Excel到R 数据分析进阶指南》一3.4 更改数据格式
  2. c语言字母表关键字,读书笔记-C语言关键字
  3. Python 技术篇-操作oracle数据库执行SQL语句报错,提示ORA-00911: 无效字符解决方法
  4. 西安交通大学17年9月课程考试计算机应用基础,西安交通大学17年9月课程考试《计算机应用基础》作业考核试题 (含主观题)...
  5. DuplicateHandle复制句柄保护文件
  6. 光纤交换机光纤通道协议介绍
  7. Yii1.1 CGridView 简单使用
  8. OpenCV的工具函数
  9. 获取到的数字证书如何配置在自己的Apache中?
  10. resnet结构_来聊聊ResNet及其变种
  11. Spring MVC 基础及相关概念(基础一)
  12. 学好java再学c 可以吗_再论学好C的重要性!!!
  13. win10 中的系统变量 windir、userProfile 等
  14. 吃了核辐射食物怎么办_尿酸过高怎么办?这几种食物,平时可多吃
  15. TF-tf.keras.layers.Dense
  16. css背景颜色设置后没生效
  17. 祝新年快乐html小程序,小程序网页开发真诚的祝福大家2018...
  18. 你不知道的原生js广播频道接口
  19. html美食网站源码
  20. osd 显示 png格式图片水印

热门文章

  1. 在photoshop中如何设置一张a4纸上打印多个照片,如小二寸图片排版
  2. [Swift]判断手机设备是否插入SIM卡
  3. 网络安全应急响应具体操作流程
  4. webgl坐标转换_WebGL笔记_绘制流程以及三维坐标变换(一)
  5. Word设置页码格式为共 页 第 页,附带页码不从第一页开始的教程
  6. Clover使用过程中遇到的问题
  7. qPCR技术发展史及最新进展
  8. 计算机变慢的十大原因
  9. PAT(乙级) 1032 挖掘机技术哪家强 (20 分)(C语言)
  10. 那些害惨我们的伪科学(下)