| ²é¿´: 1771 | »Ø¸´: 6 | |||
¡îð©Ôµ±¿Õ¡îгæ (³õÈëÎÄ̳)
|
[ÇóÖú]
ÇóÒ»¸öÓÃFDTD£¨Ê±ÓòÓÐÏÞ²î·Ö·¨£©Ä£Äâµç´Å²¨ÈýάRCS£¨É¢É䣩µÄÔ´³ÌÐò ÒÑÓÐ2È˲ÎÓë
|
| ÇóÒ»¸öÓÃFDTD£¨Ê±ÓòÓÐÏÞ²î·Ö·¨£©Ä£Äâµç´Å²¨ÈýάRCS£¨É¢É䣩µÄÔ´³ÌÐò¡£±¾ÈËÏÖÔÚ¼±ÓÃÒ»¸öÈýάɢÉäµÄ³ÌÐò£¬ÓÃÓÚ·ÂտģÄâÊÒÄÚ³£ÓòÄÁ϶Եç´Å²¨µÄ·øÉäºÍÉ¢É仹ÓÐ͸Éä¹ý³Ì£¬Ï£ÍûÓÐÄÄλºÃÐÄÈ˸ø¸öÔ´³ÌÐò£¬matlab»òÕßC¶¼ÐУ¨±¾È˲»»ác++£¬fortranµÈÓïÑÔ£©¡£Ò»¶¨ÒªÊÇÈýάµÄ£¬¶øÇÒÊÇ·ÖÎöÉ¢ÉäµÄ¡£ |
» ²ÂÄãϲ»¶
PRBͶ¸å£¬Acknowledgment sent to author״̬ʮÌìÁË
ÒѾÓÐ3È˻ظ´
PVAÈܽâ
ÒѾÓÐ2È˻ظ´
ÎïÀíѧIÂÛÎÄÈóÉ«/·ÒëÔõôÊÕ·Ñ?
ÒѾÓÐ289È˻ظ´
»ù½ðÉêÇë
ÒѾÓÐ44È˻ظ´
CSCÓëÐÂÎ÷À¼Î¬¶àÀûÑÇ´óѧPhD½±Ñ§½ðÏîÄ¿
ÒѾÓÐ0È˻ظ´
ÐÂÎ÷À¼RobinsonÑо¿Ëù ÕÐÆ¸CSC¹«ÅÉ·ÃÎÊÈËÔ±
ÒѾÓÐ0È˻ظ´
°ïÎÒµÄÓ¢Óï¿ÚÓïÀÏʦÕÒѧÉú
ÒѾÓÐ0È˻ظ´
ʲôʱºò¿ª½±£¿
ÒѾÓÐ13È˻ظ´
É¢½ð±ÒÆí¸£
ÒѾÓÐ94È˻ظ´
Çà»ùÒÑÖÐ
ÒѾÓÐ1È˻ظ´
» ±¾Ö÷ÌâÏà¹ØÉ̼ÒÍÆ¼ö: (ÎÒÒ²ÒªÔÚÕâÀïÍÆ¹ã)
» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:
feixiaolin
ÈÙÓþ°æÖ÷ (ÎÄ̳¾«Ó¢)
-

ר¼Ò¾Ñé: +518 - Ó¦Öú: 942 (²©ºó)
- ¹ó±ö: 1.275
- ½ð±Ò: 3891
- É¢½ð: 58785
- ºì»¨: 532
- ɳ·¢: 11
- Ìû×Ó: 24217
- ÔÚÏß: 2601.8Сʱ
- ³æºÅ: 2139575
- ×¢²á: 2012-11-21
- רҵ: ¹âѧÐÅÏ¢»ñÈ¡Óë´¦Àí
- ¹ÜϽ: Êýѧ
2Â¥2013-11-02 18:04:26
¡îð©Ôµ±¿Õ¡î
гæ (³õÈëÎÄ̳)
|
3Â¥2013-11-03 10:58:58
¡¾´ð°¸¡¿Ó¦Öú»ØÌû
|
3-D FDTD code with PEC boundaries clear %*********************************************************************** % Fundamental constants %*********************************************************************** cc=2.99792458e8; %speed of light in free space muz=4.0*pi*1.0e-7; %permeability of free space epsz=1.0/(cc*cc*muz); %permittivity of free space %*********************************************************************** % Grid parameters %*********************************************************************** ie=50; %number of grid cells and Ex samples along x-direction je=24; %number of grid cells and Ey samples along y-direction ke=10; %number of grid cells and Ez samples along z-direction ib=ie+1; %grid boundary coordinates jb=je+1; kb=ke+1; is=26; %location of z-directed current source js=13; %location of z-directed current source kobs=5; dx=0.002; %space increment of cubic lattice dt=dx/(2.0*cc); %time step nmax=500; %total number of time steps %*********************************************************************** % Differentiated Gaussian pulse excitation %*********************************************************************** rtau=50.0e-12; tau=rtau/dt; ndelay=3*tau; J0=-1.0; %*********************************************************************** % Material parameters %*********************************************************************** eps=1.0; sig=0.0; %*********************************************************************** % Updating coefficients %*********************************************************************** ca=(1.0-(dt*sig)/(2.0*epsz*eps))/(1.0+(dt*sig)/(2.0*epsz*eps)); cb=(dt/epsz/eps/dx)/(1.0+(dt*sig)/(2.0*epsz*eps)); da=1.0; db=dt/muz/dx; %*********************************************************************** % Field arrays %*********************************************************************** ex=zeros(ie,jb,kb); ey=zeros(ib,je,kb); ez=zeros(ib,jb,ke); hx=zeros(ib,je,ke); hy=zeros(ie,jb,ke); hz=zeros(ie,je,kb); %figure %set(gcf,'DoubleBuffer','on') %*********************************************************************** % BEGIN TIME-STEPPING LOOP %*********************************************************************** for n=1:nmax %*********************************************************************** % Update electric fields %*********************************************************************** ex(1:ie,2:je,2:ke)=ca*ex(1:ie,2:je,2:ke)+... cb*(hz(1:ie,2:je,2:ke)-hz(1:ie,1:je-1,2:ke)+... hy(1:ie,2:je,1:ke-1)-hy(1:ie,2:je,2:ke)); ey(2:ie,1:je,2:ke)=ca*ey(2:ie,1:je,2:ke)+... cb*(hx(2:ie,1:je,2:ke)-hx(2:ie,1:je,1:ke-1)+... hz(1:ie-1,1:je,2:ke)-hz(2:ie,1:je,2:ke)); ez(2:ie,2:je,1:ke)=ca*ez(2:ie,2:je,1:ke)+... cb*(hx(2:ie,1:je-1,1:ke)-hx(2:ie,2:je,1:ke)+... hy(2:ie,2:je,1:ke)-hy(1:ie-1,2:je,1:ke)); ez(is,js,1:ke)=ez(is,js,1:ke)+... J0*(n-ndelay)*exp(-((n-ndelay)^2/tau^2)); %*********************************************************************** % Update magnetic fields %*********************************************************************** hx(2:ie,1:je,1:ke)=hx(2:ie,1:je,1:ke)+... db*(ey(2:ie,1:je,2:kb)-ey(2:ie,1:je,1:ke)+... ez(2:ie,1:je,1:ke)-ez(2:ie,2:jb,1:ke)); hy(1:ie,2:je,1:ke)=hy(1:ie,2:je,1:ke)+... db*(ex(1:ie,2:je,1:ke)-ex(1:ie,2:je,2:kb)+... ez(2:ib,2:je,1:ke)-ez(1:ie,2:je,1:ke)); hz(1:ie,1:je,2:ke)=hz(1:ie,1:je,2:ke)+... db*(ex(1:ie,2:jb,2:ke)-ex(1:ie,1:je,2:ke)+... ey(1:ie,1:je,2:ke)-ey(2:ib,1:je,2:ke)); %*********************************************************************** % Visualize fields %*********************************************************************** timestep=int2str(n); tview(:, =ez(:,:,kobs);sview(:, =ez(:,js, ;subplot('position',[0.15 0.45 0.7 0.45]),imagesc(tview'); shading flat; caxis([-1.0 1.0]); colorbar; axis image; axis xy; title(['Ez(i,j,k=5), time step = ',timestep]); xlabel('i coordinate'); ylabel('j coordinate'); subplot('position',[0.15 0.10 0.7 0.25]),imagesc(sview'); shading flat; caxis([-1.0 1.0]); colorbar; axis image; axis xy; title(['Ez(i,j=13,k), time step = ',timestep]); xlabel('i coordinate'); ylabel('k coordinate'); pause(0.05) %*********************************************************************** % END TIME-STEPPING LOOP %*********************************************************************** end |
4Â¥2014-10-07 16:11:09
ѩӰÎ趯
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 229.5
- Ìû×Ó: 6
- ÔÚÏß: 1.4Сʱ
- ³æºÅ: 4139955
- ×¢²á: 2015-10-13
- ÐÔ±ð: MM
- רҵ: ¹âѧ

5Â¥2015-10-21 15:27:14
1099797132
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 16.5
- Ìû×Ó: 10
- ÔÚÏß: 3.6Сʱ
- ³æºÅ: 3490643
- ×¢²á: 2014-10-22
6Â¥2015-11-02 21:33:42
ѩӰÎ趯
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 229.5
- Ìû×Ó: 6
- ÔÚÏß: 1.4Сʱ
- ³æºÅ: 4139955
- ×¢²á: 2015-10-13
- ÐÔ±ð: MM
- רҵ: ¹âѧ

7Â¥2015-11-03 09:36:37









»Ø¸´´ËÂ¥
=ez(:,:,kobs);
ͶƱ:
10