signal_noise_real.m

music高阶谱分析算法研究

   %产生实余弦信号signal与高斯噪声noise的函数
   function   [signal_real,noise]=signal_noise_real(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=[6  90]*pi/180;
          for row_2=1:signal_number;           
                  for column_2=1:samplepoint_number;
                      signal_phase(row_2)=2*pi*signal_frequency*column_2/sample_frequency+signal_initialphase(row_2);
                                   signal_complex(row_2,column_2)=signal_amplitude(row_2)*exp(j*signal_phase(row_2));
                      signal_real(row_2,column_2)=signal_amplitude(row_2)*real(exp(j*signal_phase(row_2)));
                 end;
          end;
    
   %产生复高斯信号源signal与高斯噪声noise的函数
               %生成复高斯型信号源
               sigma_signal=sigma_noise*10.^(SNR/20);
               for  row_signal=1:signal_number
                    for column_signal=1:samplepoint_number;
                         p1=rand(1,1);
                         p2=rand(1,1);
                         signal_gauss_real(row_signal,column_signal)=sqrt(-2*sigma_signal(row_signal)^2*log(p1))*exp(j*(2*pi*p2));
                    end 
               end