相控阵雷达阵列天线仿真计算

1.基本原理

        阵列天线由多个单元组成，每个单元都是一个天线，它们可以协同工作。当这些元件结合在一起时，可以使得阵列天线的辐射方向更加的集中。如果我们把天线的辐射方向互相叠加，就可以得到阵列天线的发射方向图。

方向图可以用下面这个公式来表示：

       AF(theta) = |sum( w(n)*exp(j(n-1)beta*sin(theta)) )/N|^2

       其中，Theta是天线的方向，w(n)是每个天线元件的权重系数，beta是阵元间的相位差，N是阵列天线的元素数。这个方程刻画了阵列天线在不同方向上的辐射强度。

      说明：^2 表示对整个复数幅值取平方，即对实部和虚部的平方分别相加。这个平方操作可以将幅度转化为单位功率，同时也可以将负数幅度转化为正数。在天线阵列理论中，常常使用平方幅度函数（或称方向图）来描述天线在不同方向的辐射/接收功率情况，而这里的 AF(theta) 就是一个典型的平方幅度函数。

 2.推导过程

      结构示意图如图1所示。

                                            图1 相控阵雷达结构示意图（截取魏老师雷达原理教程的图）

说明：  

       学习魏青老师雷达原理 相位扫描法 章节，把公式全部推导一遍，包括魏老师省略的推导部分 ，有几点说明如下：

1. 为简化推导难度， 以下用E表示电场强度，E中既包含幅度，又有相位信息。
2. 不同单元不考虑接收回波幅度的大小（因为远场对于不同单元可以等效），但相位必须要考虑，因为一个波长就是2pi；
3. 如上图，前端路径（不同阵元接收到的信号）的延迟（相位）1点相当于0点（基准点）相位超前，后端的移相器补偿则滞后。

 

 

 3.matlab程序仿真

% 定义阵列天线参数

theta = -45 : 0.1 :45; % 方位角角度范围

d = 0.05; % 阵列天线元件间距

lambda = 0.05; % 信号波长

N = 10; % 阵列天线元件个数

beta = 2 * pi * d / lambda; % 阵列元件间距所对应的相位差

% 计算阵列权重

w = ones(N, 1); % 全向天线阵列权重为 1

% 计算阵列天线方向图

AF = zeros(size(theta));

for i = 1 : length(theta)

phi = 2 * pi * d / lambda * sin(theta(i) / 180 * pi); % 阵列天线相位差公式

AF(i) = abs(sum(w .* exp(1j * (0 : N-1)' * beta * sin(theta(i) / 180 * pi) + 1j * phi))); % 阵列天线方向图公式

end

% 绘制阵列天线方向图

plot(theta, AF);

title('阵列天线方向图');

xlabel('方位角（度）');

ylabel('增益（dB）');

grid on;

 说明如下：

（1）sin(theta(i) / 180 * pi)：将输入的方向角度值转化为弧度制，并用正弦函数计算该方向相对于阵列中心的角度值。

（2）(0 : N-1)' * beta：生成一个列向量，表示阵列中各个天线对应的相位差值（根据阵列间距和波长计算）。该值会被乘以该方向角度的正弦值，用于计算该方向上每个天线产生的相位差值。

（3）1j * phi：根据之前计算的相位差公式，生成标量相位差值（phi），并使用虚数单位（1j）来表示该相位差值。

（4）exp(1j * (0 : N-1)' * beta * sin(theta(i) / 180 * pi) + 1j * phi)：将上述第2步和第3步计算的值相加，并取指数得到复数值（即各个天线的电场幅度和相位）。

（5）w .* exp(1j * (0 : N-1)' * beta * sin(theta(i) / 180 * pi) + 1j * phi)：将阵列的元素权重和第4步中计算的电场值相乘。

（6）sum(w .* exp(1j * (0 : N-1)' * beta * sin(theta(i) / 180 * pi) + 1j * phi))：求出第5步中的各个元素的和，即为该方向的阵列电场。

（7）AF(i) = abs(sum(w .* exp(1j * (0 : N-1)' * beta * sin(theta(i) / 180 * pi) + 1j * phi)))：将第6步中的结果取绝对值，并将结果存入一个数组，即为该方向的方向图。

（8）上述是接收天线方向图，发射天线方向图只是计算方法类似于接收天线方向图，这句代码有区别：

           AF(i) = sum(w .* exp(1j * (0 : N-1)' * beta * sin(theta(i) / 180 * pi) - 1j * phi));

      a）sum前不需要取绝对值，原因是：在计算过程中使用的是幅度值，即各个天线在指定方向角度上的电场幅度大小的加权和。与接收天线方向图不同的是，发射天线不需要求解接收到的电磁波强度，而只需要计算发射出去的电磁波在各个方向角度上的幅度分布，因此不需要取绝对值。

      b）相位差和正负号的计算公式不同：对于发射天线，由于是向外辐射电磁波，因此一般相位差、正负号等参数要取负号；而对于接收天线，由于是接收入射电磁波，因此这些参数是正号。