coefficient_h.m

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

% coefficient of H field for FDTD formula
function [Hxhe,Hxhh,Hyhe,Hyhh,Hzhe,Hzhh]=Coefficient_H(u,segema,dt,Nx,Ny,Nz)
u0=pi*4e-7;
% u_Hx(1:Nx,1:Ny,1:Nz)=u0;
% u_Hx(2:Nx-1,:,:)=0.5*(u_Hx(1:Nx-2,:,:)+u_Hx(3:Nx,:,:));
% temp=segema*dt./(2*u_Hx);
temp=segema*dt./(2*u);
Hxhh=(1-temp)./(1+temp);
% Hxhe=dt./u_Hx./(1.+temp);
Hxhe=dt./u./(1+temp);

% u_Hy(1:Nx,1:Ny,1:Nz)=u0;
% u_Hy(:,2:Ny-1,:)=0.5*(u_Hy(:,1:Ny-2,:)+u_Hy(:,3:Ny,:));
% temp=segema*dt./(2*u_Hy);
Hyhh=(1-temp)./(1+temp);
% Hyhe=dt./u_Hy./(1.+temp);
Hyhe=dt./u./(1+temp);

% u_Hz(1:Nx,1:Ny,1:Nz)=u0;
% u_Hz(:,:,2:Nz-1)=0.5*(u_Hz(:,:,1:Nz-2)+u_Hz(:,:,3:Nz));
% temp=segema*dt./(2*u_Hz);
Hzhh=(1-temp)./(1+temp);
% Hzhe=dt./u_Hz./(1.+temp);
Hzhe=dt./u./(1+temp);