mur2ndformula.m

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

function [W]=Mur2ndformula(W00x1,Wx1,W00x0,W0x0,W0x1,W0x0addy,W0x0suby,W0x1addy,W0x1suby,W0x0addz,W0x0subz,W0x1addz,W0x1subz,dx,dy,dz,dt)  % W00表示n-1时刻，W0表示n时刻，当前时刻为n+1
c=3e8;
c1=(c*dt-dx)/(c*dt+dx);
c2=2*dx/(c*dt+dx);
c3=(c*dt)^2*dx/(2*dy^2*(c*dt+dx));
c4=(c*dt)^2*dx/(2*dz^2*(c*dt+dx));
W=-W00x1+c1*(Wx1+W00x0)+c2*(W0x0+W0x1)+c3*(W0x0addy-2*W0x0+W0x0suby+W0x1addy-2*W0x1+W0x1suby)+c4*(W0x0addz-2*W0x0+W0x0subz+W0x1addz-2*W0x1+W0x1subz);