二维时域有限差分PML的结缔组织边界

G

Guest

Guest
您好每1:

我最近与差分。我想订出的结果如峰

uplaoded。它的加入造成的领域和与散射场结缔组织边界。

我怎样才能代码吗?我已经写了一些,但它不是正确的。我猜想这是

与增加的相关源。不同的是它作为usuall〜?

有人告诉我〜

<img src="http://www.edaboard.com/images/smiles/icon_cry.gif" alt="哭或很伤心" border="0" />

14分钟后添加:这2维代码〜(我写~~~~)

代码:100 × 100%iblcb = 11;%的结缔组织边界数组的索引,左%

ibrcb =伊布- 10;%权%

jbfcb = 11;%战线%

jbbcb = jb - 10;%返回%对于n = 1:7.0 *头

源(北)= 10 *( - (2/tau ^ 2)* n的 6 /头)*格式( - ((正延迟)^ 2/tau ^ 2));%高斯%

末端宰(12,50)=源(n)的;%结缔组织地界黄嘌呤和海兰HX(注册iblcb:ibrcb,jbfcb)= HX(注册iblcb:ibrcb,jbfcb)...

-dt/muz *((内斯(iblcb:ibrcb,jbfcb)内斯(iblcb:ibrcb,jbfcb - 1))/ dx的) 胸苷/穆日/ dx的*内斯(iblcb:ibrcb,jbfcb);%前会CB%HX(注册iblcb:ibrcb,jbbcb 1)= HX(注册iblcb:ibrcb,jbbcb 1)...

-dt/muz *((内斯(iblcb:ibrcb,jbbcb 1)-内斯(iblcb:ibrcb,jbbcb))/ dx的)-dt/muz/dx *内斯(iblcb:ibrcb,jbbcb);%回立法会CB%海兰(iblcb,jbfcb:jbbcb)=海兰(iblcb,jbfcb:jbbcb)...

胸苷/穆日*((内斯(iblcb,jbfcb:jbbcb)内斯(iblcb - 1,jbfcb:jbbcb))/ dx的)-dt/muz/dx *内斯(iblcb,jbfcb:jbbcb);%会CB%左海兰(ibrcb 1,jbfcb:jbbcb)=海兰(ibrcb 1,jbfcb:jbbcb)...

胸苷/穆日*((内斯(ibrcb 1,jbfcb:jbbcb)内斯(ibrcb,jbfcb:jbbcb))/ dx的) 胸苷/穆日/ dx的*内斯(ibrcb,jbfcb:jbbcb);%权会CB%
%结缔组织边界与雌二醇宰(iblcb,jbfcb:jbbcb)=内斯(iblcb,jbfcb:jbbcb) 胸苷/ epsz *((海兰(iblcb 1,jbfcb:jbbcb)杂交(iblcb,jbfcb:jbbcb))/ dx的...

- (HX(注册iblcb,jbfcb 1:jbbcb 1)HX(注册iblcb,jbfcb:jbbcb))/ dx的) - (胸苷/ epsz / dx的)*海兰(iblcb,jbfcb:jbbcb);%会CB%左宰(ibrcb,jbfcb:jbbcb)=内斯(ibrcb,jbfcb:jbbcb) 胸苷/ epsz *((海兰(ibrcb 1,jbfcb:jbbcb)杂交(ibrcb,jbfcb:jbbcb))/ dx的...

- (HX(注册ibrcb,jbfcb 1:jbbcb 1)HX(注册ibrcb,jbfcb:jbbcb))/ dx的) (胸苷/ epsz / dx的)*海兰(ibrcb 1,jbfcb:jbbcb);%权文件%宰(iblcb:ibrcb,jbfcb)=内斯(iblcb:ibrcb,jbfcb) 胸苷/ epsz *((海兰(iblcb 1:ibrcb 1,jbfcb)杂交(iblcb:ibrcb,jbfcb))/ dx的.. 。

- (HX(注册iblcb:ibrcb,jbfcb 1)HX(注册iblcb:ibrcb,jbfcb))/ dx的) (胸苷/ epsz / dx的)* HX(注册iblcb:ibrcb,jbfcb);%前会CB%宰(iblcb:ibrcb,jbbcb)=内斯(iblcb:ibrcb,jbbcb) 胸苷/ epsz *((海兰(iblcb 1:ibrcb 1,jbbcb)杂交(iblcb:ibrcb,jbbcb))/ dx的.. 。

- (HX(注册iblcb:ibrcb,jbbcb 1)HX(注册iblcb:ibrcb,jbbcb))/ dx的) - (胸苷/ epsz / dx的)* HX(注册iblcb:ibrcb,jbbcb 1);%回立法会CB%%4角宰(iblcb,jbfcb)=内斯(iblcb,jbfcb) 胸苷/ epsz *((海兰(iblcb 1,jbfcb)杂交(iblcb,jbfcb))/ dx的...

- (HX(注册iblcb,jbfcb 1)HX(注册iblcb,jbfcb))/ dx的) - (胸苷/ epsz / dx的)*海兰(iblcb,jbfcb) (胸苷/ epsz / dx的)* HX(注册iblcb,jbfcb );宰(iblcb,jbbcb)=内斯(iblcb,jbbcb) 胸苷/ epsz *((海兰(iblcb 1,jbbcb)杂交(iblcb,jbbcb))/ dx的...

- (HX(注册iblcb,jbbcb 1)HX(注册iblcb,jbbcb))/ dx的) - (胸苷/ epsz / dx的)*海兰(iblcb,jbbcb) - (胸苷/ epsz / dx的)* HX(注册iblcb,jbbcb 1);宰(ibrcb,jbbcb)=内斯(ibrcb,jbbcb) 胸苷/ epsz *((海兰(ibrcb 1,jbbcb)杂交(ibrcb,jbbcb))/ dx的...

- (HX(注册ibrcb,jbbcb 1)HX(注册ibrcb,jbbcb))/ dx的) (胸苷/ epsz / dx的)*海兰(ibrcb 1,jbbcb) - (胸苷/ epsz / dx的)* HX(注册ibrcb ,jbbcb 1);宰(ibrcb,jbfcb)=内斯(ibrcb,jbfcb) 胸苷/ epsz *((海兰(ibrcb 1,jbfcb)杂交(ibrcb,jbfcb))/ dx的...

- (HX(注册ibrcb,jbfcb 1)HX(注册ibrcb,jbfcb))/ dx的) (胸苷/ epsz / dx的)*海兰(ibrcb 1,jbfcb) (胸苷/ epsz / dx的)* HX(注册ibrcb ,jbfcb);
 
你可以看到Taflove答
: 计算电动力学FDTD2D代码。时域有限差分方法(Artech众议院,1995年)。

附属的

 

Welcome to EDABoard.com

Sponsor

Back
Top