signal_noise_complex.m

三种DOA方法：经典的MUSIC,前后向空间平滑和改进的空间平滑算法

   %产生复正弦信号signal与高斯噪声noise的函数
   function   [signal_complex,noise_gauss_complex]=signal_noise_complex(antenna_number,signal_number,samplepoint_number,SNR,signal_frequency,sample_frequency)
            %生成高斯噪声
               sigma_noise=1;
               noise=normrnd(0,sigma_noise,antenna_number,samplepoint_number);
            %产生复正弦连续波信号signal
              signal_amplitude=sigma_noise*10.^(SNR/20);
%                signal_initialphase=[7  50  124   86]*pi/180;
          signal_initialphase=[6   186    8      -174    42.7]*pi/180;
          for row_2=1:signal_number;    
%               signal_initialphase(row_2)=2*pi*rand;
                 for  column_2=1:samplepoint_number;
%                       signal_phase(row_2,column_2)=2*pi*signal_frequency*column_2/sample_frequency+signal_initialphase(row_2);
                      signal_phase(row_2,column_2)=2*pi*signal_frequency*column_2/sample_frequency+2*pi*rand;%signal_initialphase(row_2);
                      signal_complex(row_2,column_2)=signal_amplitude(row_2)*exp(j*signal_phase(row_2,column_2));
                             signal_real(row_2,column_2)=sqrt(2)*signal_amplitude(row_2)*conj(exp(j*signal_phase(row_2)));
                 end;
          end;
    
               %生成复高斯型噪声noise
               for  row_noise=1:antenna_number
                    for column_noise=1:samplepoint_number;
                         p1=rand(1,1);
                         p2=rand(1,1);
                         noise_gauss_real(row_noise,column_noise)=sqrt(-2*sigma_noise^2*log(p1))*cos(2*pi*p2);
                         noise_gauss_imag(row_noise,column_noise)=sqrt(-2*sigma_noise^2*log(p1))*sin(2*pi*p2);
                         noise_gauss_complex(row_noise,column_noise)=(noise_gauss_real(row_noise,column_noise)+j*noise_gauss_imag(row_noise,column_noise))/sqrt(2);
                    end 
               end