阵列非线性产物辐射方向图分析1
问题背景
大概是因为5G通信将大量使用有源相控阵技术，因此关于阵列非线性研究成为这几年的热点1，在几个业内较好期刊比如TMTT、TAP、CL等都能找到相关问题的成果报道2。
非线性问题及相关的线性化技术最早在单通道的收发链路中有广泛的研究，比如主要用于发射链路的数字预失真技术（DPD，Digital Predistortion）和主要用于接收链路的非线性均衡技术（NLEQ，Nonlinear equalization）。
但在对于有源阵列的非线性问题，无论从对系统影响、问题的建模还是从处理方法上都同单通道有所不同，因此很有研究价值和必要（实际上可能是因为好发论文）。
大体上关于阵列非线性的研究可以分为以下三类：
- 对阵列非线性特性的研究；
- 对如何利用阵列非线性的研究；
- 对如何消除阵列非线性的研究。
后续的研究大概都基于对阵列非线性特性的理解上。恰好在浏览文献的过程中，翻到一篇比较基础的文献3，分析了发射阵列谐波及交调产物辐射方向图的特性。这里先学习下单音信号激励的简单情况。
单音非线性辐射方向图推导
用三阶幂级数来模拟发射功率放大器的非线性特性，如式（1），并简单假设所有的放大器具有相同的系数。
v
0
=
k
1
⋅
v
i
+
k
2
⋅
v
i
2
+
k
3
⋅
v
i
3
（
1
）
v_{0}=k_{1} \cdot v_{i}+k_{2} \cdot v_{i}^{2}+k_{3} \cdot v_{i}^{3} （1）
v0=k1⋅vi+k2⋅vi2+k3⋅vi3（1）
考虑均匀线阵列ULA，阵列单元为N，单元间距为d，输入信号可表示为（2），其中An为第n个单元的幅度加权，
n
Δ
φ
n \Delta \varphi
nΔφ为相位梯度。
v
i
n
=
A
n
cos
(
w
t
−
n
Δ
φ
)
（
2
）
v_{i n}=A_{n} \cos (w t-n \Delta \varphi) （2）
vin=Ancos(wt−nΔφ)（2）
对于特定的波束指向 θ 0 \theta_{0} θ0和工作角频率 w 0 w_{0} w0，可以如（3）表示
Δ φ = ( w 0 d sin θ 0 ) / c （ 3 ） \Delta \varphi=\left(w_{0} d \sin \theta_{0}\right) / c （3） Δφ=(w0dsinθ0)/c（3）
将（2）式代入到（1）式中，即可得到第n个单通道的发射输出，包含直流项、基波、二次谐波和三次谐波项，如式（4）所示。
v
0
=
k
2
A
n
2
2
+
(
k
1
A
n
+
3
4
k
3
A
n
3
)
⋅
cos
(
w
t
−
n
Δ
φ
)
+
k
2
A
n
2
2
⋅
cos
(
2
w
t
−
2
n
Δ
φ
)
+
1
4
k
3
A
n
3
⋅
cos
(
3
w
t
−
3
n
Δ
φ
)
（
4
）
\begin{array}{l} v_{0}=\frac{k_{2} A_{n}^{2}}{2}+\left(k_{1} A_{n}+\frac{3}{4} k_{3} A_{n}^{3}\right) \cdot \cos (w t-n \Delta \varphi) \\ +\frac{k_{2} A_{n}^{2}}{2} \cdot \cos (2 w t-2 n \Delta \varphi)\\ +\frac{1}{4} k_{3} A_{n}^{3} \cdot \cos (3 w t-3 n \Delta \varphi) \end{array} （4）
v0=2k2An2+(k1An+43k3An3)⋅cos(wt−nΔφ)+2k2An2⋅cos(2wt−2nΔφ)+41k3An3⋅cos(3wt−3nΔφ)（4）
由（4）可以得到阵因子（5）。
A
F
(
θ
)
=
∑
n
=
1
N
(
k
1
A
n
+
3
4
k
3
A
n
3
)
e
j
(
w
t
+
n
w
d
sin
θ
c
−
n
Δ
φ
)
+
∑
n
=
1
N
k
2
A
n
2
2
e
(
2
w
t
+
2
n
w
d
sin
θ
c
−
2
n
Δ
φ
)
+
∑
n
=
1
N
1
4
k
3
A
n
3
e
j
(
3
w
t
+
3
n
w
d
sin
θ
c
−
3
n
Δ
φ
)
(
5
)
\begin{array}{l} AF(\theta)=\sum_{n=1}^{N}\left(k_{1} A_{n}+\frac{3}{4} k_{3} A_{n}^{3}\right) e^{j\left(w t+n \frac{w d \sin \theta}{c}-n \Delta \varphi\right)} \\ +\sum_{n=1}^{N} \frac{k_{2} A_{n}^{2}}{2} e^{\left(2 w t+2 n \frac{w d \sin \theta}{c}-2 n \Delta \varphi\right)}\\ +\sum_{n=1}^{N} \frac{1}{4} k_{3} A_{n}^{3} e^{j\left(3 w t+3 n \frac{w d \sin \theta}{c}-3 n \Delta \varphi\right)} \end{array} (5)
AF(θ)=∑n=1N(k1An+43k3An3)ej(wt+ncwdsinθ−nΔφ)+∑n=1N2k2An2e(2wt+2ncwdsinθ−2nΔφ)+∑n=1N41k3An3ej(3wt+3ncwdsinθ−3nΔφ)(5)
（5）式的三项分别对用基波方向图、二次谐波方向图和三次谐波方向图。进一步由（6）可知，基波，谐波的方向图主波束都指向，由下式求得的 θ \theta θ方向。
w d sin θ c − Δ φ = 0 （ 6 ） \frac{w d \sin \theta}{c}- \Delta \varphi=0 （6） cwdsinθ−Δφ=0（6）
仿真代码和结果
用被禁的某lab软件编程仿真，非线性参数设置如式（7），代码如下：
v
0
=
0.98
⋅
v
i
+
0.4
⋅
v
i
2
−
0.2
⋅
v
i
3
（
7
）
v_{0}=0.98 \cdot v_{i}+0.4 \cdot v_{i}^{2}-0.2 \cdot v_{i}^{3} （7）
v0=0.98⋅vi+0.4⋅vi2−0.2⋅vi3（7）
clear all; close all; clc;
% 参数初始化
N = 10; % 线阵单元数
Amp = ones(1,N); % 幅度加权向量
Freq0 = 1e6; % 中心频率
Freq = Freq0; % 扫描频率
C = 3*10^8; % 光速
lamda0 = C/Freq0; % 波长
d = lamda0/2; % 单元间距为二分之波长
Steering_Angle = 24; % 波束指向角，单位degree
Phi_Const = 2*pi*Freq0*d*sin(Steering_Angle*2*pi/360)/C;
Elevation_Angle = -90:0.1:90; % 俯仰扫描角，单位degree
K1 = 0.98; K2 = 0.4; K3 = -0.2; % 非线性级数的系数
% 阵因子循环求和
Num = length(Elevation_Angle);
AF_Fund = zeros(1,Num);
AF_Second = zeros(1,Num);
AF_Third = zeros(1,Num);
for k = 1:Num
for n = 1:N
temp = 2*pi*Freq*d*sin(Elevation_Angle(k)*2*pi/360)/C;
AF_Fund(k)= AF_Fund(k)+(K1*Amp(n)+(3/4)*K3*Amp(n)^3)*exp(j*(n*(temp)-n*Phi_Const));
AF_Second(k)= AF_Second(k)+(K2*Amp(n)^2)/2*exp(j*(2*n*(temp)-2*n*Phi_Const));
AF_Third(k)= AF_Third(k)+(K3*Amp(n)^3)/4*exp(j*(3*n*(temp)-3*n*Phi_Const));
end
end
% 分贝转换
AF_Fund_dB = 20*log10(abs(AF_Fund)/N);
AF_Second_dB = 20*log10(abs(AF_Second)/N);
AF_Third_dB = 20*log10(abs(AF_Third)/N);
% 绘制方向图
plot(Elevation_Angle,AF_Fund_dB);grid on;hold on;
plot(Elevation_Angle,AF_Second_dB);
plot(Elevation_Angle,AF_Third_dB);
xlabel('EL(deg.)'),xlim([-90 90]);
ylabel('AF(dB)'),ylim([-50 0]);
legend('Fund','Second','Third');
波束指向24°的运行结果如下，基波、谐波波束指向相同，由于单元间距d按基波波长的一半设定，因此谐波方向图出现栅瓣，如图1所示。
图1 基波、二次谐波、三次谐波扫描方向图（波束指向24°，单元间距lamda/2，ULA单元10）
保持波束指向24°，单元间距按三次谐波波长的一半设置，栅瓣即可消除，如图2所示。
图2 基波、二次谐波、三次谐波扫描方向图（波束指向24°，单元间距3*lamda/2，ULA单元10）
下面进一步看一个稍微复杂点的问题（这部分文献种没有）：波束倾斜时谐波的波束指向。通常相控阵基于窄带假设，因此可用相移代替时延进行波束合成。如果窄带假设不成立，对应下面的代码扫描频率不等于中心频率，就会产生波束倾斜（Beam Squiting）,那么谐波会发生什么变化？
Freq0 = 1e6; % 中心频率
%Freq = Freq0; % 扫描频率
Freq = 1.5*Freq0;
其实通过（5）和（6）可以理论分析出，谐波也会产生相同的波束倾斜，同时波束指向依然与基波相同。通过改变上述代码可以仿真分析如图3。中心频率波束指向24°，而扫描频率波束指向偏离。
图3 中心频率、工作频率基波、二次谐波、三次谐波扫描方向图（中心频率波束指向24°，单元间距3lamda/2，ULA单元20，Freq = 1.5Freq0; ）
简单结论
这里主要学了文章3的第一部分即单频单波束的发射，简单结论如下。后面有空再分享下双频双波束的情况（引用文献）。
- 波束指向 ：对于单音信号，其基波、二次谐波、三次谐波乃至n次谐波的波束指向相同；
- 幅度加权 ：对于单音信号，二次谐波和三次谐波的加权分别是基波加权的2次和3次关系，结论可以类推到n次谐波；
- 谐波栅瓣 ：如果以基波波长的一半来设定天线单元间距，那么二次谐波和三次谐波的辐射方向图同样会出栅瓣。
- 谐波波束倾斜 ：对于宽带信号，谐波也会产生波束倾斜，谐波波束指向与基波波束指向依然保持相同。
N. Peccarelli, B. James, R. Irazoqui, J. Metcalf, C. Fulton and M. Yeary, “Survey: Characterization and Mitigation of Spatial/Spectral Interferers and Transceiver Nonlinearities for 5G MIMO Systems,” in IEEE Transactions on Microwave Theory and Techniques, vol. 67, no. 7, pp. 2829-2846. ↩︎
B. Rupakula and G. M. Rebeiz, “Third-Order Intermodulation Effects and System Sensitivity Degradation in Receive-Mode 5G Phased Arrays in the Presence of Multiple Interferers,” in IEEE Transactions on Microwave Theory and Techniques, vol. 66, no. 12, pp. 5780-5795, Dec. 2018. ↩︎
C. Hemmi, “Pattern characteristics of harmonic and intermodulation products in broadband active transmit arrays,” in IEEE Transactions on Antennas and Propagation, vol. 50, no. 6, pp. 858-865, June 2002. ↩︎ ↩︎
227
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
