mean_var.m

来自「music高阶谱分析算法研究」· M 代码 · 共 39 行

function   [signal_azimuth_12_solved,s_azimuth_12_solved,s_azimuth_12_samplemean,s_azimuth_12_mean]=mean_var(R,antenna_number,samplepoint_number,antenna_space,lambda_bochang,measured_number,signal_azimuth)
           s_azimuth_1_samplemean=0;
           s_azimuth_2_samplemean=0;
  for      measured_index=1:measured_number 
%调用信号方位角signal_azimuth_solved及对应谱值Pmu_solved的子程序fromeigencomposition_lastsolved        
         %?在调用该函数过程中有问题
          [Pmu{measured_index},search_azimuth,AOA_solved{measured_index},Pmu_solved{measured_index},signal_number_juded(measured_index),big_index(measured_index)]=fromeigencomposition_lastsolved(R{measured_index},antenna_number,samplepoint_number,antenna_space,lambda_bochang);
           AOA_solved_size{measured_index}=size(signal_azimuth_solved{measured_index});
           AOA_solved_size_1(measured_index)=signal_azimuth_solved_size{measured_index}(1,2);
           if     signal_azimuth_solved_size_1(measured_index)==2
                  signal_azimuth_1_solved(measured_index)=signal_azimuth_solved{measured_index}(1); 
                  signal_azimuth_2_solved(measured_index)=signal_azimuth_solved{measured_index}(2); 
                  if   abs(signal_azimuth_solved{measured_index}(1)-signal_azimuth(1))<=abs(signal_azimuth_solved{measured_index}(1)-signal_azimuth(2))| abs(signal_azimuth_solved{measured_index}(2)-signal_azimuth(2))<=abs(signal_azimuth_solved{measured_index}(2)-signal_azimuth(1))
                       s_azimuth_1_samplemean=s_azimuth_1_samplemean+signal_azimuth_solved{measured_index}(1,1)/measured_number;
                       s_azimuth_2_samplemean=s_azimuth_2_samplemean+signal_azimuth_solved{measured_index}(1,2)/measured_number;
                       s_azimuth_1(measured_index)=signal_azimuth_solved{measured_index}(1,1);
                       s_azimuth_2(measured_index)=signal_azimuth_solved{measured_index}(1,2);
                   else
                       s_azimuth_1_samplemean=s_azimuth_1_samplemean+signal_azimuth_solved{measured_index}(1,2)/measured_number;
                       s_azimuth_2_samplemean=s_azimuth_2_samplemean+signal_azimuth_solved{measured_index}(1,1)/measured_number;
                       s_azimuth_1(measured_index)=signal_azimuth_solved{measured_index}(1,2);
                       s_azimuth_2(measured_index)=signal_azimuth_solved{measured_index}(1,1);
                   end;     
           else
                  if   abs(signal_azimuth_solved{measured_index}(1)-signal_azimuth(1))<=abs(signal_azimuth_solved{measured_index}(1)-signal_azimuth(2))
                       signal_azimuth_1_solved(measured_index)=signal_azimuth_solved{measured_index}(1);
                       s_azimuth_1(measured_index)=signal_azimuth_solved{measured_index}(1);
                   else
                       signal_azimuth_2_solved(measured_index)=signal_azimuth_solved{measured_index}(1);
                       s_azimuth_2(measured_index)=signal_azimuth_solved{measured_index}(1);
                   end;
           end;
   end;
           signal_azimuth_12_solved=[signal_azimuth_1_solved',signal_azimuth_2_solved']; 
           s_azimuth_12_solved=[s_azimuth_1',s_azimuth_2'];
           s_azimuth_12_samplemean=[s_azimuth_1_samplemean,s_azimuth_2_samplemean];
           s_azimuth_1_mean=mean(s_azimuth_1);
           s_azimuth_2_mean=mean(s_azimuth_2);
           s_azimuth_12_mean=[s_azimuth_1_mean,s_azimuth_2_mean];