pec.m

来自「FDTD法计算平行耦合双微带线电容电感参数的程序」· M 代码 · 共 23 行

% 该函数确定导体边界条件
function [Exc,Ezc]=PEC(Nx,Ny,Nz,dx,dy,dz,W,S,TT,NH)
Exc(1:Nx,1:Ny,1:Nz)=1;  Ezc(1:Nx,1:Ny,1:Nz)=1; 
Nxsta=round((Nx-(2*W+S)/dx)/2); Nxend=round(Nx/2+(2*W+S)/(2*dx));
Exc(Nxsta:Nxsta+W/dx,NH:NH+TT/dy,1:Nz)=0;
Exc(Nxend-W/dx:Nxend,NH:NH+TT/dy,1:Nz)=0;
Ezc(Nxsta:Nxsta+W/dx,NH:NH+TT/dy,1:Nz)=0;
Ezc(Nxend-W/dx:Nxend,NH:NH+TT/dy,1:Nz)=0;
Exc(:,1,:)=0;  %接地板 
Ezc(:,1,:)=0;