resolution_mmusic_s_1.m

music高阶谱分析算法研究

                             % function    resolution_mmusic_s_1       
% antenna_number                                天线阵元个数
% A0A                                           信号的方位角/度
% A0A_radian                                    设置信号的角度为弧度
% signal_number                                 信号源个数
% signal_frequency                              信号的中心频率
% sample_frequency                              抽样频率
% lambda_bochang                                信号源波长
% antenna_space                                 天线阵元的间距
% samplepoint_number                            对天线阵输出X的采样个数
% SNR                                           信号与噪声的对数信噪比
clc;                                      % MATLAB语言实现MUSIC算法在天线测向上的应用
clear;
close all;
antenna_number=5;                               %天线阵元个数
AOA=[-2   5];                                  %信号的方位角/度
AOA_radian=AOA*pi/180;                          %设置信号的角度为弧度
signal_number=length(AOA);                      %信号源个数
signal_frequency=6*1e9;
sample_frequency=2.5*signal_frequency;
lambda_bochang=3*1e11/signal_frequency;    
antenna_space=0.5*lambda_bochang;               %天线阵元的间距
samplepoint_number=1000;                        %采样个数
SNR=[7   7   5];                           %信号与噪声的信噪比
%#####设置信号的方向矩阵azimuth_matrice
for row_1=1:antenna_number;                 
    for column_1=1:signal_number;
        t=(row_1-1)*2*pi*antenna_space/lambda_bochang*sin(AOA_radian(column_1));
        azimuth_matrice(row_1,column_1)=exp(-j*t);
    end;
end;
%调用产生信号源和噪声的子程序signal_noise
        [signal,noise]=signal_noise_real(antenna_number,signal_number,samplepoint_number,SNR,signal_frequency,sample_frequency);
%调用求信号源（噪声）相关系数矩阵的子程序corr_coefficient    
        [signal_corr,noise_corr]=corr_coefficient(signal,noise,signal_number,antenna_number);      
         Jm=zeros(antenna_number,antenna_number);
                 for  row_3=1:antenna_number
                      Jm(row_3,antenna_number-(row_3-1))=1;
                 end; 
%#####求天线阵输出信号X、及其协方差矩阵R
measured_number=200;
  for   measured_index=1:measured_number 
         X{measured_index}=azimuth_matrice*signal+normrnd(0,1,antenna_number,samplepoint_number);
         %MUSIC法相关矩阵R
         R{measured_index}=1/samplepoint_number*X{measured_index}*X{measured_index}';
         Rb_m{measured_index}=Jm*conj(R{measured_index})*Jm;
         %MMUISC法相关矩阵Rmmusic
         Rmmusic{measured_index}=R{measured_index}+Rb_m{measured_index};
  end;
%   resoluted_propability_mmusic=reso_propa_1_mmusic(Rmmusic,antenna_number,samplepoint_number,antenna_space,lambda_bochang,measured_number);
  [AOA_solved_mmusic,AOA_12_solved_mmusic,AOA_samplemean_mmusic,AOA_mean_mmusic]=mean_var(Rmmusic,antenna_number,samplepoint_number,antenna_space,lambda_bochang,measured_number,AOA);
 
   resoluted_propability_music=reso_propa_1_mmusic(R,antenna_number,samplepoint_number,antenna_space,lambda_bochang,measured_number);
  [AOA_solved_music, AOA_12_solved_music, AOA_samplemean_music, AOA_mean_music]=mean_var(R,antenna_number,samplepoint_number,antenna_space,lambda_bochang,measured_number,AOA);