%該函數求出了信號（噪聲）的相關系數矩陣  
   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,:)))