MUSIC空间谱估计算法matlab实现

%--------------------------------------------------------------------------
%       Copyright 2012 @ Destiny Corporation.  All rights reserved.
%
%       Description: 均匀圆阵MUSIC仿真
%       By         : Destiny
%       Time       : 2013 05
%       Modify     :
%
%--------------------------------------------------------------------------
clear all;
clc;
close all;

M=10;                                %阵元个数;   
sn=4;                                %信号个数
dw=6;                                %半径波长比   
snr1=[51015 20];
snr=10.^（snr1/20）;                   %信噪比（幅度表示）

If=[2e6 3e6 4e6  6e6 ];              %信号中频率（单位：MHz）
f=[2.0e3 3.0e3 4.0e3 6.0e3];         %信号调制频率（单位：KHz）
fs=20e6;                             %采样频率（满足采样定理）（单位：MHz）

N=1024;                              %采样点数
n=1:N; 
theta = （0:0.5:180）*pi/180;
phi = （0:1:90）*pi/180;

atheta=[70 150 50 120] *pi/180;   %信号方位角
aphi = [15 55 30 60] *pi/180;     %仰角
ss=zeros（snN）;

for k = 1:sn
      A（:k） = exp （-1j*2*pi