1124阅读
3回复

一个关于二维FDTD的问题。

UID ：43237

注册：2009-10-07
登录：2014-01-15
发帖：55
等级：仿真一级

0楼发表于: 2010-04-20 09:20:06

大家好，我是一个新手，我在学习FDTD，下面是一本书上的一个程序运行的结果，附上图，但是我看不出这是哪里出了错。
程序如下：
#include "fdtd2d.h"
int main(void)
{
double ga[IE][JE], dz[IE][JE], ez[IE][JE], hx[IE][JE], hy[IE][JE];
double ddx, dt, T, epsz, pi, epsilon, sigma, eaf;
double t0, spread, pulse;
int i, j, ic, jc, l, n, nsteps;
FILE *fp, *fopen();
ic = IE/2;
jc = JE/2;
ddx = 0.001;
dt = ddx/6e8;
epsz = 8.8e-12;
pi = 3.14159;
for ( j=0; j<JE; j++ )
{
  printf("%2d", j);
  for( i=0; i<IE; i++ )
  {
   dz[j] = 0.; Q+T#J9Y
   ez[j] = 0.; fd&Fn=!
   hx[j] = 0.; q()o|V
   hy[j] = 0.; iY,C0=n5Y
   ga[j] = 1.0; /GIGE##1F
   printf( "%5.2", ga[j] ); |l; Ot=C=
  } B>^6tdz
  printf ("\n"); PN l/}'
} Uz^N6q
t0 = 20.; T>& q8'lD
spread = 6.0; 9j>2C
T = 0; CUI3^;&S
nsteps = 1; ]`&_!T
while(nsteps > 0) K {' atc
{ @-!P1]V|
  printf("nsteps--->"); Fvl\.
  scanf("%d", &nsteps); @|anu&Hm
  printf("%d \n", nsteps); 5qtk#FB
  for ( n=1; n<=nsteps; n++ ) S{+t>en
  { .[O{,r
   T = T + 1; l/'GbuECm
   //start of the main FDTD loop 1'E=R0`pA
    wf\"&xwh?
   //put a gaussian pulse in the middle E(r_mF7:
   pulse = exp( -0.5*(pow((t0-T)/spread, 2.0))); j#5a&Z
   dz[ic][jc] = pulse; x|~8?i$%
   //calculate the Dz field ]!B0= XP
   for ( j=1; j<IE; j++ ) $vegU]-R
   { @*YF!LdU{M
    for ( i=1; i<IE; i++) L@S\ rImw
    { j~-N2b6z
     dz[j] = dz[j] + 0.5 * ( hy[j] - hy[i-1][j] - hx[j] + hx[j-1] ); "/Qz?1>l+
    } 2X@9o4_4q
   } ?s^qWA
   //calculate the Ez field ~v;I>ij
   for ( j=1; j<IE; j++ ) o72r `2
   { ?Ij(B}D
    for ( i=1; i<IE; i++) .`OU\LA
    { JY,$B-l
     ez[j] = ga[j]*dz[j]; dw"Es;^
    } Ek)drt7cy
   } >&PM'k
   //calculate the Hx field OTXZdAv
   for ( j=1; j<IE; j++ ) kXwAw]ogN
   { %} ``:
    for ( i=1; i<IE; i++) *-nO,K>y`
    { 1!v{#w{u7
     hx[j] = hx[j] + 0.5 * ( ez[j] - ez[j+1] ); 3x+lf4"
    } P51M?3&=l
   } @Lnv
   //calculate the Hx field r N$0qo
   for ( j=1; j<IE; j++ ) d?*]/ZiR
   { zr.\7\v
    for ( i=1; i<IE; i++) j&d5tgLB
    { w \b+OW
     hy[j] = hy[j] + 0.5 * ( ez[i+1][j] - ez[j+1] ); v1OVrk>s>
    } PAYw:/(P
   } 0~qc,-)3
  } ,eyh%k*hz
     //----end of the main FDTD loop S0^a)#D &
  for ( j=1; j < jc; j++ ) _Vr- bpAf
  { #DUfEZ
   printf( "%2d", j ); zEI+)|4?r
   for ( i=1; i < ic; i++ ) ^^N|:80
   { {cXr!N^K
    printf ( "%5.2f", ez[2*i][2*j]); "JB4Uaa
   } 5qz,FKx5
    printf("\n"); zC2:c"E I
  } lx:$EJ
  printf( "T = %5.0f \n", T ); )zr*Ecz
  //Write the E field out to a file " .. DhyR
[Yt{h9
未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册