下面是一个使用 **Python** 实现的频控阵天线（FDA）**方向图与能量波束仿真**的完整示例代码。该代码不仅计算了频控阵的方向图响应，还通过引入一个简单的“能量波束”模型，模拟了不同距离下的波束分布。 --- ## ✅ 功能概述： - 模拟一个线性频控阵（Linear FDA） - 每个天线单元使用不同的频率 - 仿真方向图（Angle Response） - 仿真波束能量分布（Angle + Range） --- ## ✅ Python 仿真代码： ```python import numpy as np import matplotlib.pyplot as plt # 频控阵参数设置 num_antennas = 8 # 天线数量 d = 0.5 # 天线间距（单位为参考频率下的波长λ） f0 = 10e9 # 参考频率（Hz） delta_f = 1e6 # 频率偏移步进（Hz） c = 3e8 # 光速（m/s） # 角度和距离范围 theta = np.linspace(-90, 90, 360) # 扫描角度范围 r = np.linspace(100, 1000, 100) # 距离范围（米） # 计算每个天线单元的频率 frequencies = f0 + delta_f * np.arange(num_antennas) # 定义频控阵响应函数（考虑角度和距离） def fda_pattern(theta, r, d, f0, frequencies, c): theta = np.radians(theta) k0 = 2 * np.pi * f0 / c n = np.arange(num_antennas)[:, None, None] # 天线索引 [N, 1, 1] theta = theta[None, :, None] # 角度 [1, Angles, 1] r = r[None, None, :] # 距离 [1, 1, Ranges] # 每个天线单元的波数 k = 2 * np.pi * frequencies[:, None, None] / c # 相位计算：每个单元在不同角度和距离的相位 phase = k * d * n * np.sin(theta) - k * r + k0 * r # 相干叠加所有天线响应 pattern = np.abs(np.sum(np.exp(1j * phase), axis=0)) # [Angles, Ranges] return pattern # 计算方向图 pattern = fda_pattern(theta, r, d, f0, frequencies, c) # 归一化 pattern_norm = pattern / np.max(pattern) # 绘图：角度 vs 距离 响应图 plt.figure(figsize=(12, 6)) plt.imshow(20 * np.log10(pattern_norm), aspect='auto', cmap='jet', extent=[r[0], r[-1], theta[-1], theta[0]], vmin=-40, vmax=0) plt.colorbar(label='幅度 (dB)') plt.title("频控阵天线能量波束仿真（角度-距离响应）") plt.xlabel("距离 (m)") plt.ylabel("角度 (°)") plt.grid(False) plt.show() # 绘图：固定距离下的角度响应 plt.figure(figsize=(10, 5)) plt.plot(theta, 20 * np.log10(pattern_norm[:, 50]), label='Range = 550m') plt.title("频控阵角度响应（固定距离）") plt.xlabel("角度 (°)") plt.ylabel("幅度 (dB)") plt.grid(True) plt.ylim(-40, 0) plt.legend() plt.show() ``` --- ## 🔍 代码解释： ### 1. 参数设置： - `f0`: 参考频率（如 10 GHz） - `delta_f`: 每个天线单元之间的频率偏移 - `d`: 天线间距（通常设为参考频率下波长的一半） - `theta`: 扫描角度范围（-90° 到 90°） - `r`: 距离范围（从 100m 到 1000m） ### 2. 频率分布： - 每个天线单元使用一个频率：`f_i = f0 + i * delta_f` ### 3. 波束响应模型： - 使用公式： $$ \phi_i(\theta, r) = k_i (d i \sin\theta - r) + k_0 r $$ 其中： - $ k_i = 2\pi f_i / c $ - $ k_0 = 2\pi f_0 / c $ - 最终波束响应为所有天线单元响应的相干叠加。 ### 4. 绘图部分： - 第一张图显示角度和距离的联合响应（能量波束图） - 第二张图显示固定距离下的角度响应（类似传统方向图） --- ## ✅ 仿真结果分析： - 主瓣方向会随着距离变化而偏移（角度-距离耦合） - 不同距离下，能量最大值对应的角度不同 - 可用于雷达中实现角度-距离联合估计 --- ##