MATLAB|基于3D FDTD的微带线馈矩形天线分析[用于模拟超宽带脉冲通过线馈矩形天线的传播，以计算微带结构的回波损耗参数]

 👨‍🎓个人主页：研学社的博客   

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

---

💥1 概述

📚2 运行结果

部分代码：

close all; clear; clc;
%% Physical constants
   epsilon0 = 8.85418782e-12; mu0 = 1.25663706e-6;
   c = 1.0/sqrt(mu0*epsilon0);

%% Gaussian half-width
   t_half = 15.0e-12;

%% Microstrip fed parameters
   lineW = 1.6604e-3; 
   lineH = 0.4000e-3;
   % Roger's 5880 Duroid parameters
   lineEr = 2.2;      % eps_r 
   lineTan = 0.00041; % loss tangent 
   
%% End time
   t_end = 2.0e-9;

%% Total mesh dimensions and grid cells sizes (without PML)
   nx = 70; ny = 100; nz = 6;
   dx = 0.2372e-3; dy = 0.2265e-3; dz = 0.1000e-3;

%% Number of PML layers (ten or more!)
   PML = 10;

% Add PML layers
   nx = nx + 2*PML; ny = ny + 2*PML; nz = nz + 2*PML;
% Calculate dt, number of FDTD iterations and Z0 microstrip line matched load    
   dt = (1.0/c/sqrt( 1.0/(dx^2) + 1.0/(dy^2) + 1.0/(dz^2)))*0.999;
   number_of_iterations = ceil(t_end/dt);
   eps_eff = 0.5*(lineEr+1) + 0.5*(lineEr-1)/sqrt(1+12*lineH/lineW);
   Z0 = 377/sqrt(eps_eff)*1/(lineW/lineH + 1.393 + 0.667*log(lineW/lineH) + 1.444);

%% Matrix of material's constants
   number_of_materials = 4;
   % For material of number x = 1,2,3... :
   % Material(x,1) - relative permittivity, Material(x,2) - relative permeability,
   % Material(x,3) - specific conductivity
   % Vacuum
   Material(1,1) = 1.0;   Material(1,2) = 1.0;   Material(1,3) = 0.0;
   % Metal (Copper)
   Material(2,1) = 1.0;   Material(2,2) = 1.0;   Material(2,3) = 5.88e+7;
   % Substrate material (RT/Duroid 5880)
   Material(3,1) = lineEr;   Material(3,2) = 1.0;
   % Calculate conductivity of Duroid at 20 GHz from loss tangent and eps_r
   Material(3,3) = 2*pi*10e9*lineTan*lineEr*epsilon0;
   % Matched load material is calculated from transmission line parameters
   Material(4,1) = 1.0;   Material(4,2) = 1.0;   Material(4,3) = lineH/(Z0*lineW*dy);

%% 3D array for geometry
   Index = ones(nx, ny, nz);

🎉3 参考文献

部分理论来源于网络，如有侵权请联系删除。

[1]井甜甜. 基于FDTD算法的微带滤波天线的设计[D].曲阜师范大学,2020.DOI:10.27267/d.cnki.gqfsu.2020.000943.

[2]井甜甜,赵建平,张月,杨君,徐娟.基于FDTD算法的微带滤波天线的设计[J].通信技术,2019,52(04):991-995.

[3]张敏. 基于并行FDTD的大型天线及阵列模拟与优化研究[D].西南交通大学,2017.

[4]Vasily Kozhevnikov (2023). Microstrip line-fed rectangular antenna analysis using 3D FDTD.

🌈4 Matlab代码实现