%實現(xiàn)算法：信號特征矢量重排法（算法作者：葉中付）
%程序作者：章希睿
%編寫日期：2010年10月12日
%文件名稱：eigenvectorrearrange.m
%算法簡要說明：該算法是針對相干信號源提出的一種解相干方法，其實現(xiàn)步驟如下：
%1.采取N陣元均勻直線標(biāo)量陣列獲取M個相干信號源，假設(shè)信號源全部相干（M%2.求陣列接收數(shù)據(jù)的最大似然協(xié)方差矩陣Rx，并進行特征值分解，確定特征矢量的個數(shù)，進而得到重排矩陣的維數(shù)L；
%3.根據(jù)特征矢量重排的法則確定重排矩陣Rr；
%4.采取MUSIC算法實現(xiàn)信號源數(shù)和DOA的估計（進行100次獨立實驗）。

%------clear&clear------%
clc;
clear?all;

%------確定均勻直線標(biāo)量陣列和信號的相關(guān)參數(shù)------%
K=100;
N=12;
lamda=1;
d=0.5*lamda;
SNR=10;
Ps=1;
Pn=Ps/10^（0.1*SNR）;
theta=[520255070]*（pi/180）;
M=length（theta）;
for?trial=1:100
????
%------產(chǎn)生全部相干的信號源------%
Sk=sqrt（Ps/2）*（randn（1K）+1i*randn（1K））;
Nk=sqrt（Pn/2）*（randn（NK）+1i*randn（NK））;

%------明確陣列的導(dǎo)向矢量矩陣------%
a1=exp（-1i*2*pi*d*[0:N-1]‘*sin（theta（1））/lamda）;
a2=exp（-1i*2*pi*d*[0:N-1]‘*sin（theta（2））/lamda）;
a3=exp（-1i*2*pi*d*[0:N-1]‘*sin（theta（3））/lamda）;
a4=exp（-1i*2*pi*d*[0:N-1]‘*sin（theta（4））/lamda）;
a5=exp（-1i*2*pi*d*[0:N-1]‘*sin（theta（5））/lamda）;
A=[a1a2a3a4a5];

%------明確陣列接收數(shù)據(jù)和數(shù)據(jù)協(xié)方差矩陣，并進行特征分解----