技術(shù)特征：

1.一種基于機翼形變擬合的共形陣列幅相誤差校正方法，包括如下步驟：

(1)假設探測區(qū)域內(nèi)信號源數(shù)N為1，共形陣列參考單元的偏移角度φ₁為零或已知，利用接收數(shù)據(jù)計算協(xié)方差矩陣R，并對該協(xié)方差矩陣R進行特征分解，得到白噪聲功率σ²和噪聲子空間U_N，其中協(xié)方差矩陣的表達式為：

$R=\frac{1}{K}\underset{t=1}{\overset{K}{\Σ}}x\left(t\right)x{\left(t\right)}^H,$

其中，K為接收數(shù)據(jù)的快拍數(shù)，x(t)為第t次快拍的接收數(shù)據(jù)，(·)^H為接收數(shù)據(jù)矩陣取共軛轉(zhuǎn)置；

(2)利用協(xié)方差矩陣R和白噪聲功率σ²，計算共形陣列中各陣元的相對接收強度并由計算得到的共形陣列中各陣元的相對接收強度構(gòu)建對角化矩陣，該矩陣即共形陣列的相對幅度誤差矩陣

(3)設定迭代序數(shù)變量j、第j次迭代中波達方向估計角度θ_j和代價函數(shù)J_j，并對序數(shù)變量j、波達方向估計角度θ_j以及代價函數(shù)J_j進行初始化，使j＝1，θ₁＝0，J₁＝0；

(4)由波達方向估計角度θ_j計算共形陣列參考陣元的絕對接收強度F(1)，再利用共形陣列參考陣元的絕對接收強度F(1)計算共形陣列的絕對幅度誤差矩陣F，并利用共形陣列的絕對幅度誤差矩陣F計算共形陣列中各陣元的偏移角度φ₁,φ₂,…,φ_M；

(5)利用共形陣列中各陣元的偏移角度φ₁,φ₂,…,φ_M，計算共形陣列中各陣元的偏移斜率ρ₁ ρ₂ … ρ_M，并根據(jù)計算出的共形陣列中各陣元的偏移斜率ρ₁ ρ₂ … ρ_M，利用多項式曲線擬合算法擬合出共形陣列的斜率曲線f(x)，再對擬合出的共形陣列的斜率曲線f(x)進行積分，得到共形陣列的變形曲線F(x)，最后利用積分得到的共形陣列的變形曲線F(x)，計算共形陣列中陣元的偏移量矢量y；

(6)根據(jù)計算出的共形陣列中陣元的偏移量矢量y，計算共形陣列的相位誤差矩陣W，并根據(jù)共形陣列的相位誤差矩陣W，得到變形后共形陣列的響應矩陣Γ(θ_j)A(θ_j)，將該響應矩陣Γ(θ_j)A(θ_j)代入MUSIC譜函數(shù)，搜索MUSIC譜，得到波達方向估計角度θ_j+1；

(7)利用波達方向估計角度θ_j+1計算代價函數(shù)，將代價函數(shù)的值代入收斂條件，判斷收斂條件是否滿足，若是，則代表共形陣列幅相誤差校正完成，輸出波達方向估計角度θ_j+1；否則，令j＝j+1，執(zhí)行步驟(4)。

2.根據(jù)權(quán)利要求1所述的基于機翼形變擬合的共形陣列幅相誤差校正方法，其特征在于，步驟(2)中所述的共形陣列的相對幅度誤差矩陣其獲取步驟為：

2a)利用協(xié)方差矩陣R與白噪聲功率σ²，計算共形陣列中各陣元的相對接收強度計算公式為：

${\stackrel{\‾}{\mathrm{\α}}}_i=\sqrt{\frac{R_{i,i}-{\mathrm{\σ}}^2}{R_{1,1}-{\mathrm{\σ}}^2}},i=1,2,...,M---\<1\>$

其中，R_i,i為協(xié)方差矩陣R第i個對角元素，M為共形陣列陣元數(shù)目；

2b)將共形陣列中各陣元的相對接收強度通過下式構(gòu)建共形陣列的相對幅度誤差矩陣

$\hat{F}=diag({\stackrel{\‾}{\mathrm{\α}}}_1,{\stackrel{\‾}{\mathrm{\α}}}_2,...{\stackrel{\‾}{\mathrm{\α}}}_M)---\<2\>$

其中，diag(·)表示由序列構(gòu)建對角化矩陣。

3.根據(jù)權(quán)利要求1所述的基于機翼形變擬合的共形陣列幅相誤差校正方法，其特征在于，步驟(4)中所述的計算共形陣列中各陣元的偏移角度φ₁,φ₂,…,φ_M，其實現(xiàn)步驟為：

4a)利用波達方向估計角度θ_j，計算共形陣列參考陣元的絕對接收強度F(1)：

F(1)＝f(θ_j+φ₁)， <3>

其中，φ₁為共形陣列參考陣元的偏移角度；

4b)利用共形陣列參考陣元的絕對接收強度F(1)和共形陣列的相對幅度誤差矩陣根據(jù)計算共形陣列的絕對幅度誤差矩陣F；

4c)利用共形陣列的絕對幅度誤差矩陣F，計算信號源相對于共形陣列中各陣元的波達方向η₁,η₂,…η_M，其計算公式為：

η_i＝f^-1(F_i,i),i＝1,2,…,M， <4>

其中，F(xiàn)_i,i表示共形陣列的絕對幅度誤差矩陣F的第i個對角元素，f(θ)為共形陣列單元天線的方向圖函數(shù)，(·)^-1表示該方向圖函數(shù)的反函數(shù)，如果F_i,i≥1，則強制設置f^-1(F_i,i)為0；

4d)取共形陣列中各陣元沿X軸正方向依次排列，根據(jù)波達方向η₁,η₂,…η_M的增減情況，計算共形陣列中各陣元的偏移角度φ₁,φ₂,…φ_i…,φ_M，其中，φ_i為共形陣列中第i個陣元偏離X軸的角度，且偏移角度φ_i偏向Y軸正方向時角度為正：

(1)如果η₁,η₂,…η_M是遞減排列，則共形陣列中各陣元的偏移角度：

φ_i＝η_i-θ_j，i＝1,2,…M； <6>

(2)如果η₁,η₂,…η_M是遞增排列，則共形陣列中各陣元的偏移角度：

φ_i＝θ_j-η_i，i＝1,2,…M； <7>

(3)如果η₁,η₂,…η_M是先增后減排列，取η₁,η₂,…η_p為遞增排列，η_p+1,η_p+2,…η_M為遞減排列，則前p項共形陣列中各陣元的偏移角度：

φ_i＝θ_j-η_i，i＝1,2,…p； <8>

后(M-p)項共形陣列中各陣元的偏移角度：

φ_i＝η_i-θ_j，i＝p+1,p+2,…M， <9>

4.根據(jù)權(quán)利要求1所述的基于機翼形變擬合的共形陣列幅相誤差校正方法，其特征在于，步驟(5)中所述的計算共形陣列中陣元的偏移量矢量y，其實現(xiàn)步驟為：

5a)利用共形陣列中各陣元的偏移角度φ₁,φ₂,…,φ_M，計算共形陣列中各陣元的偏移斜率ρ₁,ρ₂,…,ρ_M：

ρ_i＝tan(φ_i)i＝1,2…,M； <9>

5b)取共形陣列中各陣元的偏移斜率ρ₁,ρ₂,…,ρ_M滿足以X軸坐標變量為自變量的(M-1)階多項式曲線，即共形陣列的斜率曲線：

f(x)＝a₁x^M-1+a₂x^M-2+…+a_M-1x+a_M， <10>

其中，a₁,a₂,…a_M為常數(shù)，x為陣元的X軸坐標變量；將共形陣列陣元的X軸坐標點x₁,x₂,…,x_M代入式<10>中，得到關(guān)系式：

其中，[x₁,x₂,…,x_M]^T為共形陣列中陣元的X軸坐標點矢量，(·)^T表示轉(zhuǎn)置運算，ρ＝[ρ₁ρ₂ … ρ_M]^T為共形陣列中陣元的偏移斜率矢量；將關(guān)系式<11>簡化為

aX₁＝ρ <12>

其中，a＝[a₁ a₂ … a_M]^T為共形陣列的斜率曲線f(x)的系數(shù)矢量，表示以共形陣列中陣元的X軸坐標點為元素的范德蒙矩陣；將式<12>兩邊右乘矩陣得到計算共形陣列的斜率曲線f(x)的系數(shù)矢量a的表達式：

$a={\mathrm{\&rho;X}}_1^{-1},---\&lt;13\&gt;$

其中，為矩陣X₁的逆矩陣，并利用該式與共形陣列中陣元的偏移斜率矢量ρ、矩陣計算共形陣列的斜率曲線f(x)的系數(shù)矢量a；

5c)共形陣列變形曲線F(x)滿足M階多項式：

F(x)＝b₁x^M+b₂x^M-1+…+b_Mx+b_M+1 <14>

其中，b₁,…b_M,b_M+1為共形陣列的變形曲線F(x)的系數(shù)，取常量b_M+1＝0，由共形陣列的變形曲線F(x)與共形陣列的斜率曲線f(x)的積分關(guān)系得到：

$\begin{array}{c}{F}^{\&prime;}\left(x\right)=f\left(x\right)\\ =b_1{\mathrm{Mx}}^{M-1}+b_2(M-1)x^{M-2}+\mathrm{...}+b_M\end{array}---\&lt;15\&gt;$

其中，(·)'表示對共形陣列的變形曲線取導數(shù)，對比式<10>與式<15>中的系數(shù)，得到

a＝b⊕[M,M-1,…,1]^T <16>

其中，b＝[b₁,b₂,…,b_M]^T為共形陣列的變形曲線F(x)的系數(shù)矢量，⊕為矩陣的Hadamard積，將式<16>簡化，得到計算共形陣列的變形曲線F(x)的系數(shù)矢量b的表達式：

b＝a⊕n， <17>

其中，并利用該式與矢量n、共形陣列的斜率曲線f(x)的系數(shù)矢量a計算共形陣列變形曲線F(x)的系數(shù)矢量b；

5d)利用共形陣列變形曲線F(x)的系數(shù)矢量b與式<14>，計算共形陣列中陣元的偏移量矢量y：

y＝bX₂， <18>

其中，y＝[y₁ y₂ … y_M]^T，