Соединительная границы 2D FDTD ПМЛ

[Guest]

Привет каждые 1:

Я работаю с FDTD недавно.Я хочу работать в результате, как рис

uplaoded.это добавить в соединительной границе с вытекающим области и разрозненные области.

Как я могу его код?Я уже написал, но это не правильно.Я думаю, это

быть увязаны с добавлением источника.Он отличается, как usuall ~?

кто-нибудь сказать мне ~

<img src="http://www.edaboard.com/images/smiles/icon_cry.gif" alt="Плачущий или очень расстроенный" border="0" />

Добавлено спустя 14 минут:Это на 2 аспекта код ~ (я писал ~~~~)

Код:% 100x100iblcb = 11;% массива индекс соединительной границей, левый%

ibrcb = IB-10;% праве%

jbfcb = 11;% передней%

jbbcb = JB-10;% обратно%для N = 1:7.0 * тау

Источник (N) = 10 * (- (2/tau ^ 2) * N 6 / тау) * EXP (- ((N-задержка) 2/tau ^ ^ 2));% гаусс%

конецEZ (12,50) = источник (N);% Соединительной границе с 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);% передней 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); обратно% CB%Hy (iblcb, jbfcb: jbbcb) = Hy (iblcb, jbfcb: jbbcb) ...

ц / Muz * ((EZ (iblcb, jbfcb: jbbcb)-EZ (iblcb-1, jbfcb: jbbcb)) / DX) -dt/muz/dx * EZ (iblcb, jbfcb: jbbcb);% влево CB%Hy (ibrcb 1, jbfcb: jbbcb) = Hy (ibrcb 1, jbfcb: jbbcb) ...

ц / Muz * ((EZ (ibrcb 1, jbfcb: jbbcb)-EZ (ibrcb, jbfcb: jbbcb)) / DX) DT / Muz / DX * EZ (ibrcb, jbfcb: jbbcb);% CB% права
% Соединительной границе с ЭзEZ (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);% влево 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);% право ЦБ %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);% передней 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); обратно% CB%% 4 углахEZ (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]

Вы можете увидеть FDTD2D кодекса Taflove А. вычислительной электродинамики.Конечных разностей временной области методом (Artech Дом, 1995).

Пока