coefficient_e.m

FDTD法计算平行耦合双微带线电容电感参数的程序

% coefficient of E field for FDTD formula
function [Exee,Exeh,Eyee,Eyeh,Ezee,Ezeh]=Coefficient_E(eps,segema,dt,Nx,Ny,Nz)
eps0=8.854e-12;  %真空介电常数
% eps_Ex(1:Nx,1:Ny,1:Nz)=eps0;
% eps_Ex(:,2:Ny-1,2:Nz-1)=0.25*(eps(:,3:Ny,3:Nz)+eps(:,3:Ny,1:Nz-2)+eps(:,1:Ny-2,3:Nz)+eps(:,1:Ny-2,1:Nz-2));
% temp=segema*dt./(2*eps_Ex);
temp=segema*dt./(2*eps);
Exee=(1-temp)./(1+temp);
Exeh=(dt./eps)./(1+temp);

% eps_Ey(1:Nx,1:Ny,1:Nz)=eps0;
% eps_Ey(2:Nx-1,:,2:Nz-1)=0.25*(eps(3:Nx,:,3:Nz)+eps(3:Nx,:,1:Nz-2)+eps(1:Nx-2,:,3:Nz)+eps(1:Nx-2,:,1:Nz-2));
% temp=segema*dt./(2*eps_Ey);
Eyee=(1-temp)./(1+temp);
Eyeh=(dt./eps)./(1+temp);

% eps_Ez(1:Nx,1:Ny,1:Nz)=eps0;
% eps_Ez(2:Nx-1,2:Ny-1,:)=0.25*(eps(3:Nx,3:Ny,:)+eps(3:Nx,1:Ny-2,:)+eps(1:Nx-2,3:Ny,:)+eps(1:Nx-2,1:Ny-2,:));
% temp=segema*dt./(2*eps_Ez);
Ezee=(1-temp)./(1+temp);
Ezeh=(dt./eps)./(1+temp);