一维线阵干涉仪测角模型及解算

一维线阵干涉仪测角模型

为了方便处理相角变化，我们采用复数的形式进行计算,因为复数的相乘相当于，幅角相加，模长相乘，从而实现相位移动。

设电磁波信号:
S W ( t ) = A c e x p { j ( 2 π f c t + ϕ c ) } S_W(t)=A_c exp\{j(2\pi f_c t + \phi_c)\} SW​(t)=Ac​exp{j(2πfc​t+ϕc​)}
${\varphi }_c$ 代表信号的初始相位，当信号$S_W(t)$的入射角与干涉仪基线之间的夹角为${\theta }_{AOA}$,则各单元天线接收并输出的信号:
$$S_{A,i}(t)=S_W(t)B_{ant}({\theta }_{AOA})exp(j\frac{2\pi d_{i}sin({\theta }_{AOA})}{\lambda }))$$
其中$B_{ant}({\theta }_{AOA})$为单元天线方向图对应角度的增益。

对于具有下边频的超外差接收机，$N_a$个接收通道不仅具有相同的放大增益，并且采用了相同的本振信号：$S_{osc}=A_{osc}exp[-j(2\pi f_{osc}t+{\varphi }_{osc})]$进行下变频。
${\varphi }_{osc}$ ,代表本振源的初相:
$$S_{IF}(t)=Filte{r}_{IF}[S_{A,i}(t)\cdot S_{osc}(t)]$$
$$S_{IF}(t)={\beta }_{IF}{A}_c{B}_{ant}({\theta }_{AOA})\cdot exp[j(2\pi f_{IF}t+{\varphi }_c-{\varphi }_{osc})]\cdot exp(j\frac{2\pi d_{i}sin({\theta }_{AOA})}{\lambda })$$
${\beta }_{IF}$代表变频的各级放大增益，从公式中发现通道间的信号相位关系保持不变。

信号模拟仿真

假设光速3e8m/s，中频信号为400MHz，采样率为20GHz，射频信号为10GHz，设计基线为[0.049,0.106,0.177]，通过Matlab进行仿真如下：

C = 3E8;
RF = 10e9;
IF = 0.4e9
fs = 20e9;
d = [0.049,0.106,0.177];
t = 0:1/fs:0.000000005;

ORGIN_SIG = exp(j*2*pi*IF*t);
lambda = C/RF;

m = exp(j * 2 * pi .*d * sind(12) / lambda);
t_sig = zeros(length(d)+ 1, length(t));
t_sig(1,:) = imag( ORGIN_SIG);

for i = 2:length(d)+1
	t_sig(i,:) = imag(m(i - 1)*ORGIN_SIG);
end

plot(t, (t_sig(1,:)), ...
t, (t_sig(2,:)), ...
t, (t_sig(3,:)), ...
t, (t_sig(4,:)));
legend('主通道信号','0.049','0.106','0.177');

得到如下的相移结果：

![[截图_选择区域_20231127001120.png]]

相位解模糊

将接收通道0的输出的信号$S_{IF},0(t)$作为基准信号,后续的$1-N_a-1$个,接收通道的信号$S_{IF,k}(t)$, k = 1,2,3…$N_a-1$,分别进行相位差测量，即鉴相，从而得到相位差值${\varphi }_{M\Delta ,k}$（归一化到$[-\pi ,\pi ]$之间）：
$${\varphi }_{M\Delta ,k}=mod(2\pi d_{k}sin({\theta }_{AOA})/\lambda +\pi ,2\pi )-\pi$$
由上试可知相位差测量的区间在$[-\pi ,\pi ]$之间，也正因如此，基线越长，导致相位偏差越大，从而形成了相位偏差的$2\pi$模糊，如此意味着下面的式子成立：
$$|2\pi d_{k}sin({\theta }_{AOA})/\lambda |>\pi$$
反过来，如果干涉仪中最短基线长度$d_1$与最大测角范围满足$[-{\theta }_{max},{\theta }_{max}]$时，则不会出现$2\pi$模糊的情况，则$d_1$满足如下关系:
$$d_1<\frac{\lambda }{2sin({\theta }_{max})}$$
θ m a x < a s i n { λ 2 d 1 } \theta_{max}<asin \{ \frac{\lambda}{2d_1} \} θmax​<asin{2d1​λ​}
待续 2023-11-23

从理论上讲，在干涉仪的设计中通常最短基线的长度$d_1$控制在上面的公式约束范围以内，以此没有相位模糊的相位差去解下一级较长基线的相位差测量值的模糊，如此进行逐级解模糊获得$N_a-1$个无模糊的相位差测量值${\varphi }_{W\Delta ,k}$ ，($W:无模糊，M：模糊$)
$${\gamma }_{M\Delta ,k}=2\pi \cdot floor[\frac{{\alpha }_k{\varphi }_{W\Delta ,k}}{2\pi }],{\alpha }_k=\frac{d_{k+1}}{d_k},{\varphi }_{W\Delta ,1}={\varphi }_{M\Delta ,1}$$
$${\beta }_k={\varphi }_{M\Delta ,k+1}+{\gamma }_{M\Delta ,k}-{\alpha }_k\cdot {\varphi }_{W\Delta ,k}$$
$${\zeta }_k=\{\begin{array}{c}0,{\beta }_k\in [-\pi ,\pi )\\ \pi ,{\beta }_k<-\pi \\ -\pi ,{\beta }_k\ge \pi \end{array}$$
$${\varphi }_{W\Delta ,k}={\varphi }_{M\Delta ,k+1}+{\gamma }_{M\Delta ,k}+{\zeta }_k$$

侧向误差因素影响

待续 2023-11-24

虚拟基线

根据求${\theta }_{max}$的公式进行计算，当频率为18GHz的情况下，波长$\lambda =0.01667m$, 则不模糊角度为${\theta }_{max}=60$,则两个最短基线长度为:$d_1=0.0096m$，天线的安装间距限制了天线的尺寸，从而导致接收的信号质量下降，一般情况下，天线的尺寸须满足${\lambda }_{max}/2$的量级，而基线长度${\lambda }_{min}/2$的量级，侦察设备的带宽大，也会存在${\lambda }_{max}$比${\lambda }_{min}$大几倍的情况，这就形成了一组矛盾。

为了解决上诉问题，专门引入了虚拟基线的概念，设计由3个单元天线构成的两个长度接近的基线来合成一条虚拟的基线：

记天线$A_{n0}$与$A_{n1}$之间的基线长度为：$d_{fic,1}$, $A_{n1}$与$A_{n2}$之间的基线长度为：$d_{fic,2}$ ，并且满足$d_{fic,2}>d_{fic,2}$ , 则可以$A_{n}1$为对称点，与$A_{n}0$镜像处虚拟一个天线单元，记为$A_n{0}^{\prime }$,则可以合成一个基于$A_n{0}^{\prime }与A_{n}2$ 形成的长度为$\Delta d_{fic}$的虚拟基线。

将这条虚拟基线作为整个干涉仪的最短基线，使之满足干涉仪侧向的最短基线要求，从而解决了干涉仪物理长度过短而无法实现的问题。
$$\Delta d_{fic}=d_{fic,2}-d_{fic,1}<\frac{\lambda }{2sin({\theta }_{max})}$$
以上描述中，当前的最短基线是基于3单元天线进行合成，所以该基线的相位差虚拟测量值${\theta }_{fic,M\Delta }$ ，通过长度为：$d_{fic,1}$和$d_{fic,2}$的实基线相位差测量值${\varphi }_{fic,M\Delta ,1}$和${\varphi }_{fic,M\Delta ,2}$，按照下式进行计算：
$${\varphi }_{fic,M\Delta }={\varphi }_{fic,M\Delta ,2}-{\varphi }_{fic,M\Delta ,1}$$
再进行修正到$[-\pi ,\pi ]$之间：
$${\varphi }_{fic,M\Delta }=\{\begin{array}{c}{\varphi }_{fic,M\Delta },{\varphi }_{fic,M\Delta }\in [-\pi ,\pi )\\ {\varphi }_{fic,M\Delta }+\pi ,{\varphi }_{fic,M\Delta }<-\pi \\ {\varphi }_{fic,M\Delta }-\pi ,{\varphi }_{fic,M\Delta }\ge \pi \end{array}$$
在此基础上能够参照实基线逐级解模糊的方式对剩余的基线进行相位差值模糊解算。

解模糊模拟仿真

根据上一章节的场景，假设在没有测量误差的情况下进行解模糊运算

C = 3e8;
RF = 18e9;
lambda = C/RF;
angle = -70;       %假设角度
d = [0.049,0.106,0.177]; 	%真实基线

real_phase = 2 * pi .*d * sind(angle) / lambda; %真实相差
phis = mod(real_phase, 2*pi) - pi; %归一化后的相差，有模糊

%模拟获得虚拟基线的无模糊相位差
d1 = d(2) - d(1) ;
d0 = d1 - d(1);

%天线1-2之间的相位差
i_phi1 = 2 * pi .*d(1) * sind(angle) / lambda;
%天线2-3之间的相位差
i_phi2 = 2 * pi .*(d1) * sind(angle) / lambda;
%获得虚拟基线无模糊的相位差
phi_m_1 = i_phi2 - i_phi1;
phi_w_1 = phi_m_1;

%计算第一个实基线和虚拟基线之间的基线倍数
a1 = d(1) / d0;
y_m_1=2*pi*floor(a1*phi_w_1/(2*pi));
%测试角度是否正确
theta = asin(phi_w_1*lambda/(2*pi*d0))*180/pi; 

%%%%%%%%%%解第一级实基线相位差模糊%%%%%%%%%
phi_m_2 = phis(1); %有模糊的相位
phi_w_2 = phi_m_2 + y_m_1 ;

if(phi_w_1 < 0) %如果在法线左边
	m_2 = phi_w_2 + a1 * phi_w_1;
else %如果在法线右边
	m_2 = phi_w_2- a1 * phi_w_1;
end

if(m_2 >= pi)
	phi_w_2 = phi_w_2 - pi;
elseif(m_2 <= -pi)
	phi_w_2 = phi_w_2 + pi;
end

theta = asin((phi_w_2)*lambda/(2*pi*d(1)))*180/pi; %%测试计算角度
a2 = d(2) / d(1);
y_m_2=2*pi*floor(a2*phi_w_2/(2*pi));

%%%%%%%%%%解第二级实基线相位差模糊%%%%%%%%%
phi_m_3 = phis(2);%有模糊的相位
phi_w_3 = phi_m_3 + y_m_2 ;
if(phi_w_1 < 0) %如果在法线左边
	m_3 = phi_w_3 + a2 * phi_w_2;
else %如果在法线右边
	m_3 = phi_w_3- a2 * phi_w_2;
end

if(m_3 >= pi)
	phi_w_3 = phi_w_3 - pi;
elseif(m_3 <= -pi)
	phi_w_3 = phi_w_3 + pi;
end

theta = asin((phi_w_3)*lambda/(2*pi*d(2)))*180/pi; %%测试角度
a3 = d(3) / d(2);
y_m_3=2*pi*floor(a3*phi_w_3/(2*pi));

%%%%%%%%%%解第三级实基线相位差模糊%%%%%%%%%
phi_m_4 = phis(3);%有模糊的相位
phi_w_4 = phi_m_4 + y_m_3 ;
if(phi_w_1 < 0) %如果在法线左边
	m_4 = phi_w_4 + a3 * phi_w_3;
else %如果在法线右边
	m_4 = phi_w_4- a3 * phi_w_3;
end

if(m_4 >= pi)
	phi_w_4 = phi_w_4 - pi;
elseif(m_4 <= -pi)
	phi_w_4 = phi_w_4 + pi;
end
theta = asin((phi_w_3)*lambda/(2*pi*d(2)))*180/pi; %%测试角度

%根据虚拟基线求得最大测量角度
max_theta = asin(lambda/(2*d0))*180/pi;

经过以上模拟程序，解归一化后的相位差模糊，第一、二、三级的计算结果相位差与理论相位差一致，测试角度也与输入角度-70度一致。

误差因素

但此处有一个需要注意的地方，虚拟基线的相位差值是通过两条长度接近的实基线相位测量后通过上面的公式进行计算得到。

待续 2023-11-27