[发明专利]一种基于极化修正的短波测向方法

权利要求书

1.一种基于极化修正的短波测向方法，用于短波天线阵，天线阵由M个水平极化天线阵元按照圆形阵排列组成，M≥3，天线阵元间隔为θ，其特征在于：

设有一个具有中心频率ω₀的空间窄带平面波以角Θ₁入射到该天线阵，其中Θ₁＝(θ₁,φ₁)，θ₁和φ₁分别是该来波信号的仰角和相位角，0≤θ₁＜90°，0≤φ₁＜360°，则天线阵第m个阵元的输出x_m(t)表示为：

1≤m≤M，s(t)为入射天线阵的源信号、为来波信号幅度因子，jω₀为来波信号的相位因子，n_m(t)为第m个阵元的加性噪声，τ_m(Θ₁)为来波信号投射到第m个阵元时的相对时延；

步骤1，记各天线阵元的实测值cv_M，cv_M＝[cv₁,…,cv_N,…,cv_M]，该cv_M为M维复数数组，并求得各天线阵元实测接收信号的幅度值，比较选择幅度最大值及其对应的天线方向θ_N；

步骤2，记最大幅度值对应的天线阵元为N，选定N-i，…，N-1，N，N+1，…，N+i,共2*i+1个天线阵元作为参与运算的测向天线阵元；

步骤3，以天线阵元N为中心，设定扫描角fai_scan，扫描范围θ_N-i*θ≤fai_scan≤θ_N+i*θ，扫描精度Δθ取0.1°、0.5°、1°或2°；

步骤4，采用数值方法计算上述2*i+1个天线阵元的理论主极化响应，记为cv，cv＝[x_n-i，…，x_n-1，x_n，x_n+1，…，x_n+i]，该cv为2*i+1维复数数组；

步骤5，从步骤1所得天线阵元的实测值数组cv_m中，取接收幅度最大的天线阵元N-i对应的天线实际接收数据，记为cv_N-i；

步骤6，通过插值计算第N-i天线阵元关于扫描角fai_scan对称点的响应值，记为cv'_N-i，并计算(cv_N-i+cv'_N-i)/2，作为第N-i天线阵元的修正响应数据；

步骤7，按照步骤6的方法依次计算步骤2中选定的其余2*i个天线阵元的修正响应数据，将这2*i+1个数据按照顺序存储，并记为cv1，该cv1为2*i+1维复数数组；

步骤8，计算cv和cv1的相关系数R：

上式中，x是cv数据，是指cv数据均值，y是cv1数据，是指cv1数据均值；

步骤9，在扫描范围内，每间隔Δθ取一次fai_scan，重复步骤5-8，依次计算cv和cv1的相关系数，得到扫描角-相关系数曲线，相关系数最大点对应的扫描角即为来波方向，该频点测向结束。