一维线阵干涉仪测角模型
为了方便处理相角变化,我们采用复数的形式进行计算,因为复数的相乘相当于,幅角相加,模长相乘,从而实现相位移动。
设电磁波信号:
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)=Acexp{j(2πfct+ϕc)}
ϕ
c
\phi_c
ϕc 代表信号的初始相位,当信号
S
W
(
t
)
S_W(t)
SW(t)的入射角与干涉仪基线之间的夹角为
θ
A
O
A
\theta_{AOA}
θAOA,则各单元天线接收并输出的信号:
S
A
,
i
(
t
)
=
S
W
(
t
)
B
a
n
t
(
θ
A
O
A
)
e
x
p
(
j
2
π
d
i
s
i
n
(
θ
A
O
A
)
λ
)
)
S_{A,i}(t)=S_W(t)B_{ant}(\theta_{AOA})exp(j \frac{ 2\pi d_isin(\theta_{AOA})} {\lambda}))
SA,i(t)=SW(t)Bant(θAOA)exp(jλ2πdisin(θAOA)))
其中
B
a
n
t
(
θ
A
O
A
)
B_{ant}(\theta_{AOA})
Bant(θAOA)为单元天线方向图对应角度的增益。
对于具有下边频的超外差接收机,
N
a
N_a
Na个接收通道不仅具有相同的放大增益,并且采用了相同的本振信号:
S
o
s
c
=
A
o
s
c
e
x
p
[
−
j
(
2
π
f
o
s
c
t
+
ϕ
o
s
c
)
]
S_{osc}=A_{osc} exp[-j(2\pi f_{osc}t+\phi_{osc})]
Sosc=Aoscexp[−j(2πfosct+ϕosc)]进行下变频。
ϕ
o
s
c
\phi_{osc}
ϕosc ,代表本振源的初相:
S
I
F
(
t
)
=
F
i
l
t
e
r
I
F
[
S
A
,
i
(
t
)
⋅
S
o
s
c
(
t
)
]
S_{IF}(t) = Filter_{IF}[S_{A,i}(t) \cdot S_{osc}(t)]
SIF(t)=FilterIF[SA,i(t)⋅Sosc(t)]
S
I
F
(
t
)
=
β
I
F
A
c
B
a
n
t
(
θ
A
O
A
)
⋅
e
x
p
[
j
(
2
π
f
I
F
t
+
ϕ
c
−
ϕ
o
s
c
)
]
⋅
e
x
p
(
j
2
π
d
i
s
i
n
(
θ
A
O
A
)
λ
)
S_{IF}(t) = \beta_{IF} A_c B_{ant}(\theta_{AOA}) \cdot exp[j(2\pi f_{IF}t+ \phi_c - \phi_{osc})] \cdot exp(j\frac{2\pi d_i sin(\theta_{AOA}) }{ \lambda})
SIF(t)=βIFAcBant(θAOA)⋅exp[j(2πfIFt+ϕc−ϕosc)]⋅exp(jλ2πdisin(θAOA))
β
I
F
\beta_{IF}
β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
I
F
,
0
(
t
)
S_{IF},0(t)
SIF,0(t)作为基准信号,后续的
1
−
N
a
−
1
1-N_a-1
1−Na−1个,接收通道的信号
S
I
F
,
k
(
t
)
S_{IF,k}(t)
SIF,k(t), k = 1,2,3…
N
a
−
1
N_a - 1
Na−1,分别进行相位差测量,即鉴相,从而得到相位差值
ϕ
M
Δ
,
k
\phi_{M \Delta,k}
ϕMΔ,k(归一化到
[
−
π
,
π
]
[-\pi,\pi]
[−π,π]之间):
ϕ
M
Δ
,
k
=
m
o
d
(
2
π
d
k
s
i
n
(
θ
A
O
A
)
/
λ
+
π
,
2
π
)
−
π
\phi_{M \Delta,k} = mod(2 \pi d_k sin(\theta_{AOA})/\lambda+\pi,2\pi) - \pi
ϕMΔ,k=mod(2πdksin(θAOA)/λ+π,2π)−π
由上试可知相位差测量的区间在
[
−
π
,
π
]
[-\pi,\pi]
[−π,π]之间,也正因如此,基线越长,导致相位偏差越大,从而形成了相位偏差的
2
π
2\pi
2π模糊,如此意味着下面的式子成立:
∣
2
π
d
k
s
i
n
(
θ
A
O
A
)
/
λ
∣
>
π
|2 \pi d_k sin(\theta_{AOA})/\lambda| > \pi
∣2πdksin(θAOA)/λ∣>π
反过来,如果干涉仪中最短基线长度
d
1
d_1
d1与最大测角范围满足
[
−
θ
m
a
x
,
θ
m
a
x
]
[-\theta_{max},\theta_{max}]
[−θmax,θmax]时,则不会出现
2
π
2\pi
2π模糊的情况,则
d
1
d_1
d1满足如下关系:
d
1
<
λ
2
s
i
n
(
θ
m
a
x
)
d_1<\frac{\lambda}{2sin(\theta_{max})}
d1<2sin(θ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
d_1
d1控制在上面的公式约束范围以内,以此没有相位模糊的相位差去解下一级较长基线的相位差测量值的模糊,如此进行逐级解模糊获得
N
a
−
1
N_{a}-1
Na−1个无模糊的相位差测量值
ϕ
W
Δ
,
k
\phi_{W\Delta,k}
ϕWΔ,k ,(
W
:
无模糊,
M
:模糊
W:无模糊,M:模糊
W:无模糊,M:模糊)
γ
M
Δ
,
k
=
2
π
⋅
f
l
o
o
r
[
α
k
ϕ
W
Δ
,
k
2
π
]
,
α
k
=
d
k
+
1
d
k
,
ϕ
W
Δ
,
1
=
ϕ
M
Δ
,
1
\gamma_{M\Delta,k} = 2\pi \cdot floor[\frac{\alpha_k \phi_{W \Delta,k}}{2 \pi}],\alpha_k=\frac{d_{k+1}}{d_k},\phi_{W\Delta,1}= \phi_{M\Delta,1}
γMΔ,k=2π⋅floor[2παkϕWΔ,k],αk=dkdk+1,ϕWΔ,1=ϕMΔ,1
β
k
=
ϕ
M
Δ
,
k
+
1
+
γ
M
Δ
,
k
−
α
k
⋅
ϕ
W
Δ
,
k
\beta_k = \phi_{M\Delta,k + 1} + \gamma_{M\Delta,k} - \alpha_k \cdot \phi_{W\Delta,k}
βk=ϕMΔ,k+1+γMΔ,k−αk⋅ϕWΔ,k
ζ
k
=
{
0
,
β
k
∈
[
−
π
,
π
)
π
,
β
k
<
−
π
−
π
,
β
k
≥
π
\zeta_k = \left\{\begin{matrix} 0, \beta_k \in[-\pi,\pi)\\ \pi,\beta_k < -\pi \\ -\pi,\beta_k \geq \pi \\ \end{matrix}\right.
ζk=⎩
⎨
⎧0,βk∈[−π,π)π,βk<−π−π,βk≥π
ϕ
W
Δ
,
k
=
ϕ
M
Δ
,
k
+
1
+
γ
M
Δ
,
k
+
ζ
k
\phi_{W\Delta,k} = \phi_{M\Delta,k + 1} + \gamma_{M\Delta,k} + \zeta_k
ϕWΔ,k=ϕMΔ,k+1+γMΔ,k+ζk
侧向误差因素影响
待续 2023-11-24
虚拟基线
根据求 θ m a x \theta_{max} θmax的公式进行计算,当频率为18GHz的情况下,波长 λ = 0.01667 m \lambda=0.01667m λ=0.01667m, 则不模糊角度为 θ m a x = 60 \theta_{max}=60 θmax=60,则两个最短基线长度为: d 1 = 0.0096 m d_1= 0.0096m d1=0.0096m,天线的安装间距限制了天线的尺寸,从而导致接收的信号质量下降,一般情况下,天线的尺寸须满足 λ m a x / 2 \lambda_{max}/2 λmax/2的量级,而基线长度 λ m i n / 2 \lambda_{min}/2 λmin/2的量级,侦察设备的带宽大,也会存在 λ m a x \lambda_{max} λmax比 λ m i n \lambda_{min} λmin大几倍的情况,这就形成了一组矛盾。
为了解决上诉问题,专门引入了虚拟基线的概念,设计由3个单元天线构成的两个长度接近的基线来合成一条虚拟的基线:
记天线 A n 0 A_{n0} An0与 A n 1 A_{n1} An1之间的基线长度为: d f i c , 1 d_{fic,1} dfic,1, A n 1 A_{n1} An1与 A n 2 A_{n2} An2之间的基线长度为: d f i c , 2 d_{fic,2} dfic,2 ,并且满足 d f i c , 2 > d f i c , 2 d_{fic,2}>d_{fic,2} dfic,2>dfic,2 , 则可以 A n 1 A_n1 An1为对称点,与 A n 0 A_n0 An0镜像处虚拟一个天线单元,记为 A n 0 ′ A_n0' An0′,则可以合成一个基于 A n 0 ′ 与 A n 2 A_n0' 与 A_n2 An0′与An2 形成的长度为 Δ d f i c \Delta d_{fic} Δdfic的虚拟基线。
将这条虚拟基线作为整个干涉仪的最短基线,使之满足干涉仪侧向的最短基线要求,从而解决了干涉仪物理长度过短而无法实现的问题。
Δ
d
f
i
c
=
d
f
i
c
,
2
−
d
f
i
c
,
1
<
λ
2
s
i
n
(
θ
m
a
x
)
\Delta d_{fic} = d_{fic,2} - d_{fic,1}<\frac{\lambda}{2sin(\theta_{max})}
Δdfic=dfic,2−dfic,1<2sin(θmax)λ
以上描述中,当前的最短基线是基于3单元天线进行合成,所以该基线的相位差虚拟测量值
θ
f
i
c
,
M
Δ
\theta_{fic,M \Delta}
θfic,MΔ ,通过长度为:
d
f
i
c
,
1
d_{fic,1}
dfic,1和
d
f
i
c
,
2
d_{fic,2}
dfic,2的实基线相位差测量值
ϕ
f
i
c
,
M
Δ
,
1
\phi_{fic,M\Delta,1}
ϕfic,MΔ,1和
ϕ
f
i
c
,
M
Δ
,
2
\phi_{fic,M\Delta,2}
ϕfic,MΔ,2,按照下式进行计算:
ϕ
f
i
c
,
M
Δ
=
ϕ
f
i
c
,
M
Δ
,
2
−
ϕ
f
i
c
,
M
Δ
,
1
\phi_{fic,M\Delta}=\phi_{fic,M\Delta,2}-\phi_{fic,M\Delta,1}
ϕfic,MΔ=ϕfic,MΔ,2−ϕfic,MΔ,1
再进行修正到
[
−
π
,
π
]
[-\pi,\pi]
[−π,π]之间:
ϕ
f
i
c
,
M
Δ
=
{
ϕ
f
i
c
,
M
Δ
,
ϕ
f
i
c
,
M
Δ
∈
[
−
π
,
π
)
ϕ
f
i
c
,
M
Δ
+
π
,
ϕ
f
i
c
,
M
Δ
<
−
π
ϕ
f
i
c
,
M
Δ
−
π
,
ϕ
f
i
c
,
M
Δ
≥
π
\phi_{fic,M\Delta}= \left\{\begin{matrix} \phi_{fic,M\Delta},\phi_{fic,M\Delta} \in[-\pi,\pi) \\ \phi_{fic,M\Delta}+\pi,\phi_{fic,M\Delta} < -\pi \\ \phi_{fic,M\Delta}-\pi,\phi_{fic,M\Delta} \geq \pi \\ \end{matrix}\right.
ϕfic,MΔ=⎩
⎨
⎧ϕfic,MΔ,ϕfic,MΔ∈[−π,π)ϕfic,MΔ+π,ϕfic,MΔ<−πϕfic,MΔ−π,ϕfic,MΔ≥π
在此基础上能够参照实基线逐级解模糊的方式对剩余的基线进行相位差值模糊解算。
解模糊模拟仿真
根据上一章节的场景,假设在没有测量误差的情况下进行解模糊运算
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