【使用 3D FDTD 代码和 UPML 进行微带分支线耦合器分析】三维有限差分时域方法在平面微带电路分析中的应用（Matlab代码实现）

创新优化代码学习

于 2024-09-13 11:37:50 发布

阅读量362

点赞数 5

文章标签： 3d 平面 matlab

本文链接：https://blog.csdn.net/m0_58086403/article/details/142206705

版权

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

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

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

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

目录

⛳️赠与读者

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码、文章下载

⛳️赠与读者

👨‍💻做科研，涉及到一个深在的思想系统，需要科研者逻辑缜密，踏实认真，但是不能只是努力，很多时候借力比努力更重要，然后还要有仰望星空的创新点和启发点。当哲学课上老师问你什么是科学，什么是电的时候，不要觉得这些问题搞笑。哲学是科学之母，哲学就是追究终极问题，寻找那些不言自明只有小孩子会问的但是你却回答不出来的问题。建议读者按目录次序逐一浏览，免得骤然跌入幽暗的迷宫找不到来时的路，它不足为你揭示全部问题的答案，但若能让人胸中升起一朵朵疑云，也未尝不会酿成晚霞斑斓的别一番景致，万一它居然给你带来了一场精神世界的苦雨，那就借机洗刷一下原来存放在那儿的“躺平”上的尘埃吧。

或许，雨过云收，神驰的天地更清朗.......🔎🔎🔎

💥1 概述

FDTD 3D 与 UPML 用于计算平面分支线耦合器的散射系数 S_{11}、S_{21}、S_{31} 和 S_{41}，参考了 D. Sheen、S. Ali、M. Abouzahra 和 J. Kong 的原始论文《三维有限差分时域方法在平面微带电路分析中的应用》，发表在 IEEE 微波理论与技术期刊（IEEE Trans. on Microwave Theory and Techniques），DOI: [10.1109/22.55775](http://dx.doi.org/10.1109/22.55775)。当前的代码在原始计算的基础上做了一些改进：
1) 使用 UPML 替代 Mur ABCs；
2) 使用真实金属（铜）作为补丁导体材料，替代 PEC；
3) 在滤波器传输微带线的端部应用匹配负载，以防止物理反射；
4) 在 Ez 源面不应用“磁壁”或“电壁”条件；
5) 使用具有损耗的介电材料的实际属性（Duroid）。

参考文献：

摘要：
直接三维时域有限差分法（FDTD）被应用于各种微带结构的全波分析。该方法被证明是模拟复杂微带电路元件和微带天线的高效工具。根据时域结果，计算了线馈矩形贴片天线的输入阻抗以及低通滤波器和分支线耦合器的频率相关散射参数。这些电路被制造出来，对它们的测量结果与FDTD结果进行了比较，结果表明它们之间具有良好的一致性

📚2 运行结果

部分代码：

%% 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 transmission lines parameters
lineW = 2.413e-3;
lineH = 1.0e-3;
% Roger's 5880 Duroid parameters
lineEr = 2.2; % eps_r
lineTan = 0.0009; % loss tangent
Z0 = 49.2526; % Matched load transmission line

%% End time
t_end = 1.5e-9;

%% Total mesh dimensions and grid cells sizes (without PML)
nx = 60; ny = 100; nz = 16;
dx = 0.406e-3; dy = 0.406e-3; dz = 0.265e-3;

%% Number of PML layers
PML = 5;

%% 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*20e9*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);

% Add PML layers
nx = nx + 2*PML; ny = ny + 2*PML; nz = nz + 2*PML;
% Calculate dt
dt = (1.0/c/sqrt( 1.0/(dx^2) + 1.0/(dy^2) + 1.0/(dz^2)))*0.9999;
number_of_iterations = ceil(t_end/dt);

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

%% Define of low-pass filter geometry
% Ground plane
Index((1+PML):(nx-PML-1), (1+PML):(ny-PML-1), PML+1) = 2;
% Rectangular patch (thickness is equal to one cell)
Index((nx/2-14):(nx/2+15), (ny/2-14):(ny/2-9), PML+5) = 2;
Index((nx/2-14):(nx/2+15), (ny/2+9):(ny/2+14), PML+5) = 2;
Index((nx/2-14):(nx/2-9), (1+PML):(ny-PML-1), PML+5) = 2;
Index((nx/2+10):(nx/2+15), (1+PML):(ny-PML-1), PML+5) = 2;
Index((nx/2-16):(nx/2-7), (ny/2-12):(ny/2+12), PML+5) = 2;
Index((nx/2+8):(nx/2+17), (ny/2-12):(ny/2+12), PML+5) = 2;
% Dielectric substrate between ground plane and filter patch
Index((1+PML):(nx-PML-1), (1+PML):(ny-PML-1), (PML+2):(PML+4)) = 3;
% Matched load before port 1 and 2
Index((nx/2-14):(nx/2-9), PML+1, (PML+2):(PML+4)) = 4;
Index((nx/2+10):(nx/2+15), PML+1, (PML+2):(PML+4)) = 4;
% Matched load after port 2
Index((nx/2-14):(nx/2-9), ny-PML-1, (PML+2):(PML+4)) = 4;
Index((nx/2+10):(nx/2+15), ny-PML-1, (PML+2):(PML+4)) = 4;

%% 3D FDTD physical (fields) and additional arrays are defined as 'single'
%% to increase performance
Ex = zeros(nx, ny+1, nz+1, 'single');