reso_propa_mmusic.m

music高阶谱分析算法研究

   function    resoluted_propability=reso_propa_mmusic(R,initial_length,antenna_number,samplepoint_number,antenna_space,lambda_bochang,measured_number)
               resoluted_number=zeros(1,initial_length);
               for   initialphase_index=1:initial_length
                     for  measured_index=1:measured_number
      %调用求信号方位角signal_azimuth_solved及对应谱值Pmu_solved的子程序fromeigencomposition_lastsolved        
                     %？在调用该函数过程中有问题
                         [Pmu,search_azimuth,AOA_solved{initialphase_index,measured_index},Pmu_solved{initialphase_index,measured_index},signal_number_juded(initialphase_index,measured_index),big_index(initialphase_index,measured_index)]=fromeigencomposition_lastsolved(R{initialphase_index,measured_index},antenna_number,samplepoint_number,antenna_space,lambda_bochang)
                          AOA_solved_length(initialphase_index,measured_index)=length(AOA_solved{initialphase_index,measured_index});
                         if    AOA_solved_length(initialphase_index,measured_index)==2
                               AOA_1_solved(initialphase_index,measured_index)=AOA_solved{initialphase_index,measured_index}(1);
                               AOA_2_solved(initialphase_index,measured_index)=AOA_solved{initialphase_index,measured_index}(2);
                               if    abs(AOA_1_solved(initialphase_index,measured_index)-AOA_2_solved(initialphase_index,measured_index))>=2
                                     resoluted_number(initialphase_index)=resoluted_number(initialphase_index)+1;
                               end;
                         end;      
                      end;
                 end;
                 resoluted_propability=resoluted_number/measured_number;