bindevæv grænseforholdene af 2D FDTD PML

[Maggie]

Hi hver 1:

Jeg arbejder med FDTD nylig.Jeg ønsker at finde ud af det resultat, ligesom pic

uplaoded.det tilføje i bindevæv grænse med resultanten område og spredt felt.

Hvordan kan jeg kode det?Jeg har allerede skrevet nogle, men det
er ikke korrekt.Jeg tror det
er
være korreleret med tilføjer kilden.Er det anderledes som normalt ~?

nogen fortælle mig ~

<img src="http://www.edaboard.com/images/smiles/icon_cry.gif" alt="Gråd eller meget trist" border="0" />

Lagt efter 14 minutterette er for 2 dimension kode ~ (jeg skrev ~~~~)

Code:% 100x100iblcb = 11% array indekset for bindevæv grænseforholdene, venstre%

ibrcb = ib-10;% ret%

jbfcb = 11% foran%

jbbcb = jb-10;% tilbage%for n = 1:7.0 * tau

kilde = 10 * (- (2/tau ^ 2) * n 6 / tau) * exp (- ((n-forsinkelse) ^ 2/tau ^ 2));% gauss%

endeEZ (12,50) = kilde ;% Bindevæv grænse med Hx & HyHX (iblcb: ibrcb, jbfcb) = hx (iblcb: ibrcb, jbfcb) ...

-dt/muz * ((EZ (iblcb: ibrcb, jbfcb)-EZ (iblcb: ibrcb, jbfcb-1)) / dx) dt / muz / dx * EZ (iblcb: ibrcb, jbfcb);% foran CB%HX (iblcb: ibrcb, jbbcb 1) = hx (iblcb: ibrcb, jbbcb 1) ...

-dt/muz * ((EZ (iblcb: ibrcb, jbbcb 1)-EZ (iblcb: ibrcb, jbbcb)) / dx) -dt/muz/dx * EZ (iblcb: ibrcb, jbbcb);% tilbage CB%hy (iblcb, jbfcb: jbbcb) = hy (iblcb, jbfcb: jbbcb) ...

dt / muz * ((EZ (iblcb, jbfcb: jbbcb)-EZ (iblcb-1, jbfcb: jbbcb)) / dx) -dt/muz/dx * EZ (iblcb, jbfcb: jbbcb);% venstre CB%hy (ibrcb 1, jbfcb: jbbcb) = hy (ibrcb 1, jbfcb: jbbcb) ...

dt / muz * ((EZ (ibrcb 1, jbfcb: jbbcb)-EZ (ibrcb, jbfcb: jbbcb)) / dx) dt / muz / dx * EZ (ibrcb, jbfcb: jbbcb);% ret CB%
% Bindevæv grænse med EzEZ (iblcb, jbfcb: jbbcb) = EZ (iblcb, jbfcb: jbbcb) dt / epsz * ((hy (iblcb 1, jbfcb: jbbcb)-hy (iblcb, jbfcb: jbbcb)) / dx ...

- (HX (iblcb, jbfcb 1: jbbcb 1)-HX (iblcb, jbfcb: jbbcb)) / dx) - (dt / epsz / dx) * hy (iblcb, jbfcb: jbbcb);% venstre CB%EZ (ibrcb, jbfcb: jbbcb) = EZ (ibrcb, jbfcb: jbbcb) dt / epsz * ((hy (ibrcb 1, jbfcb: jbbcb)-hy (ibrcb, jbfcb: jbbcb)) / dx ...

- (HX (ibrcb, jbfcb 1: jbbcb 1)-HX (ibrcb, jbfcb: jbbcb)) / dx) (dt / epsz / dx) * hy (ibrcb 1, jbfcb: jbbcb);% ret CB %EZ (iblcb: ibrcb, jbfcb) = EZ (iblcb: ibrcb, jbfcb) dt / epsz * ((hy (iblcb 1: ibrcb 1, jbfcb)-hy (iblcb: ibrcb, jbfcb)) / dx .. .

- (HX (iblcb: ibrcb, jbfcb 1)-HX (iblcb: ibrcb, jbfcb)) / dx) (dt / epsz / dx) * HX (iblcb: ibrcb, jbfcb);% foran CB%EZ (iblcb: ibrcb, jbbcb) = EZ (iblcb: ibrcb, jbbcb) dt / epsz * ((hy (iblcb 1: ibrcb 1, jbbcb)-hy (iblcb: ibrcb, jbbcb)) / dx .. .

- (HX (iblcb: ibrcb, jbbcb 1)-HX (iblcb: ibrcb, jbbcb)) / dx) - (dt / epsz / dx) * HX (iblcb: ibrcb, jbbcb 1);% tilbage CB%% 4 hjørnerEZ (iblcb, jbfcb) = EZ (iblcb, jbfcb) dt / epsz * ((hy (iblcb 1, jbfcb)-hy (iblcb, jbfcb)) / dx ...

- (HX (iblcb, jbfcb 1)-HX (iblcb, jbfcb)) / dx) - (dt / epsz / dx) * hy (iblcb, jbfcb) (dt / epsz / dx) * HX (iblcb, jbfcb );EZ (iblcb, jbbcb) = EZ (iblcb, jbbcb) dt / epsz * ((hy (iblcb 1, jbbcb)-hy (iblcb, jbbcb)) / dx ...

- (HX (iblcb, jbbcb 1)-HX (iblcb, jbbcb)) / dx) - (dt / epsz / dx) * hy (iblcb, jbbcb) - (dt / epsz / dx) * HX (iblcb, jbbcb 1);EZ (ibrcb, jbbcb) = EZ (ibrcb, jbbcb) dt / epsz * ((hy (ibrcb 1, jbbcb)-hy (ibrcb, jbbcb)) / dx ...

- (HX (ibrcb, jbbcb 1)-HX (ibrcb, jbbcb)) / dx) (dt / epsz / dx) * hy (ibrcb 1, jbbcb) - (dt / epsz / dx) * HX (ibrcb , jbbcb 1);EZ (ibrcb, jbfcb) = EZ (ibrcb, jbfcb) dt / epsz * ((hy (ibrcb 1, jbfcb)-hy (ibrcb, jbfcb)) / dx ...

- (HX (ibrcb, jbfcb 1)-HX (ibrcb, jbfcb)) / dx) (dt / epsz / dx) * hy (ibrcb 1, jbfcb) (dt / epsz / dx) * HX (ibrcb , jbfcb);

 

[pippopamp]

Du kan se FDTD2D kode af Taflove A. Computational electrodynamics.Finite Difference Time Domain Metode (Artech House, 1995).

Farvel