基于FFT的阵列天线方向图的快速计算———截断

银尘879

已于 2025-03-12 10:59:55 修改

阅读量255

点赞数 3

CC 4.0 BY-SA版权

文章标签：算法机器学习人工智能

于 2025-03-12 10:58:37 首次发布

本文链接：https://blog.csdn.net/qq_61584218/article/details/146199094

一、知识回顾

使用离散傅里叶变换（DFT）及其快速算法（FFT）计算天线阵列方向图具有以下优点：

（1）计算效率高

（2）可进行并行计算

（3）适用于多种阵列天线结构

（4）灵活性高

（5）便于优化和综合

二、理论推导

三、仿真分析

clc
clear
close all 

%% 定义天线的基本参数
N = 16;%阵元数量
d = 0.4;%阵元间距，初始为0.4个波长
lamda = 1;%设定波长等于1
j = sqrt(-1);
A = ones(1,N);%设定均匀激励的权重向量
M = 512;%采样点数
m = -(M/2):(M/2)-1;%采样点索引

%% 相应的计算,并截断多余的点
u = m * lamda/(M*d);
u(u>1) = 1;
u(u<-1) = -1;
theta_0 = acos(u);%计算u值对应的采样角度

%% 初始化方向图
AF = zeros(1,M);
BF = zeros(1,M);

%% 使用IFFT计算均匀线阵的方向图
BF = ifft(A,M);
BF_NEW = ifftshift(BF);
BF_ABS = abs(BF_NEW);
BF_MAX = max(BF_ABS);
BF_MAG = BF_ABS/BF_MAX;
BF_MAG(u<=-1|u>=1) = 0;
BF_DB = 20*log10(BF_MAG);

%% 计算阵列因子
for n = 1:M
    for k = 1:N
        AF(n) = AF(n) + A(k)*exp(j*2*pi*d*(k-1)*u(n)/lamda);
    end
end
omega = (theta_0 *180)/pi;%将频率向量转化为角频率向量

AF_ABS = abs(AF);
AF_MAX = max(AF_ABS);
AF_MAG = AF_ABS/AF_MAX;
AF_MAG(u<=-1|u>=1) = 0;
AF_DB = 20*log10(AF_MAG);

%% 绘制方向图
figure;
plot(omega,AF_DB,'-b','LineWidth',1);%绘制方向图，蓝色实线线宽为1
hold on;
plot(omega,BF_DB,'o','MarkerSize',6,'MarkerEdgeColor','red');%绘制方向图，红色尺寸为6的空心圆
axis([0,180,-50,0]);
xlabel('\theta(deg)');
ylabel('Pattern (dB)');
legend('direct sum','IFFT');
grid on;

四、运行结果与总结

1.运行结果