corr_coefficient.m

来自「三种DOA方法：经典的MUSIC,前后向空间平滑和改进的空间平滑算法」· M 代码 · 共 25 行

   %该函数求出了信号（噪声）的相关系数矩阵  
   function  [signal_corr,noise_corr]=corr_coefficient(signal,noise,signal_number,antenna_number)             
              for  m=1:signal_number-1
                   for  n=m+1:signal_number  
                        signal_corr_coefficient=corrcoef(signal(m,:),signal(n,:));
                        signal_corr(m,n)=signal_corr_coefficient(1,2);
                        signal_corr(n,m)=signal_corr_coefficient(2,1);
                    end;
                    signal_corr(m,m)=signal_corr_coefficient(1,1);
                    signal_corr(n,n)=signal_corr_coefficient(2,2);
              end;
              for  m2=1:antenna_number-1
                   for  n2=m2+1:antenna_number
                        noise_corr_coefficient=corrcoef(noise(m2,:),noise(n2,:));
                        noise_corr(m2,n2)=noise_corr_coefficient(1,2);
                        noise_corr(n2,m2)=noise_corr_coefficient(2,1);
                    end;
                    noise_corr(m2,m2)=noise_corr_coefficient(1,1);
                    noise_corr(n2,n2)=noise_corr_coefficient(2,2);
              end;
                        
%                    correlation_coefficient_2=corrcoef(signal(2,:),noise(1,:));
%                    correlation_coefficient_3=corrcoef(signal(3,:),noise(1,:));,sig_noi_corr
           co_12=(mean(signal(1,:).*signal(2,:))-mean(signal(1,:))*mean(signal(2,:)))/sqrt(var(signal(1,:))*var(signal(2,:)))