matlab 谱峰,MUSIC算法怎么求谱峰对应的角度值

我有个参考的 你研究出来了记得告诉我，我毕设急用

%姓名：吴茜  学号：xxxxxxxxx  08.6.26

%(3)在信噪比变化范围为-20dB－+20dB时，分别采用MUSIC算法,CAPON算法, ESPRIT算法，

%   对上述四个信号源的波达方向进行估计,并画出波达方向估计误差,随着信噪比变化的曲线图

clear all;

clc;

close all;

format short

snr=-20:20;                 %定义信噪比

snap=2000;                  %快拍数

n_x=0:snap-1;

M=7;                        %阵元数

N=4;                        %信源数

K=0:M-1;

fs=100;                     %采样率

c=300000000;

lamda=2;

d=lamda/2;

f0=c/lamda ;                %目标信号的中心频率

t=[0:snap-1]/fs;

angle=[45,60,30,-45];        %信号方位角

fre=[0.015,0.05,0.02,0.035];  %信号归一化频率

for jj=1:length(snr)

%============产生阵元输出数据矩阵=