技術(shù)特征：

1.一種基于傳播算子的2-L型陣列二維DOA估計(jì)算法，其特征是，步驟如下：

步驟1：構(gòu)造接收信號(hào)矩陣

將位于坐標(biāo)原點(diǎn)的陣元作為參考陣元，線陣X,Y,Z接收到的信號(hào)向量分別為x(t)＝[x₁(t),x₂(t),…,x_N(t)]^T，y(t)＝[y₁(t),y₂(t),…,y_N(t)]^T和z(t)＝[z₁(t),z₂(t),…,z_N(t)]^T，其中x_i(t),y_i(t),z_i(t)分別表示線陣X,Y,Z上的第i個(gè)陣元在t時(shí)刻接收到的信號(hào)，N為子陣陣元數(shù)目，T表示矩陣轉(zhuǎn)置運(yùn)算，構(gòu)造新的接收信號(hào)向量w(t)＝[x^T(t),y^T(t),z^T(t)]^T，且有其中A_x,A_y,A_z分別為線陣X,Y,Z的陣列流型矩陣，A為2-L型陣列的陣列流型矩陣，s(t)為陣列的來波信號(hào)，n(t)為噪聲分量，則對(duì)應(yīng)M快拍的接收數(shù)據(jù)矩陣為W＝[w(1),w(2),…,w(M)]；

步驟2：構(gòu)造傳播矩陣

接收信號(hào)的自相關(guān)矩陣為將其按如下形式分塊，R＝[R₁,R₂]，其中H表示矩陣取共軛轉(zhuǎn)置運(yùn)算，R₁∈C^3N×K，R₂∈C^3N×(3N-K)，C為復(fù)數(shù)，則傳播矩陣為定義一個(gè)新的擴(kuò)展傳播矩陣其中I_K×K為K階的單位矩陣，K為來波信號(hào)的個(gè)數(shù)；

步驟3：估計(jì)旋轉(zhuǎn)矩陣

將P_c按如下形式分塊，其中P_x,P_y,P_z均為N×K的矩陣，定義矩陣其中A₁為A的前K行，P_z,1為P_z的前N-1行，P_z,2為P_z的后N-1行，對(duì)Ψ_z進(jìn)行特征值分解，則其特征值即為的對(duì)角線分量，為Φ_z的估計(jì)值，特征向量矩陣即為A₁的估計(jì)值，其中Φ_z為子陣Z對(duì)應(yīng)的旋轉(zhuǎn)矩陣，其形式為其中diag表示將一個(gè)向量對(duì)角化，d為陣元間距，λ為來波信號(hào)的波長(zhǎng)，θ_i表示第i個(gè)信號(hào)的俯仰角，定義兩個(gè)新的矩陣其中P_x,1為P_x的前N-1行，P_x,2為P_x的后N-1行，則有其中為Φ_x的估計(jì)值，Φ_x的形式為為對(duì)應(yīng)第i個(gè)信號(hào)的方位角，同理定義兩個(gè)新的矩陣其中P_y,1為P_y的前N-1行，P_y,2為P_y的后N-1行，則有其中為Φ_y的估計(jì)值，Φ_y的形式為

步驟4：方位角和俯仰角估計(jì)

設(shè)分別為的第k個(gè)對(duì)角線分量，則方位角和俯仰角的估計(jì)值分別為其中angle表示取幅角運(yùn)算，atan表示取反正切運(yùn)算。