FDTD參數選擇估計程序

　　delay=(time-rickper);

　　excitation(i)=exp((-rickth)*delay*delay);

　　time=time+dt;

　　i=i+1;

　　end;

　　end;

　　if(strcmp(type,'cont_sine')==1)

　　i=1;

　　while(time=total_time)

　　ramp=time*RAMPD*freq;

　　if(ramp>1.0)

　　ramp=1.0;

　　end;

　　excitation(i)=ramp*sin(2.0*pi*freq*time);

　　time=time+dt;

　　i=i+1;

　　end;

　　end;

　　if(strcmp(type,'sine')==1)

　　i=1;

　　while(time=total_time)

　　excitation(i)=sin(2.0*pi*freq*time);

　　if(time*freq>1.0)

　　excitation(i)=0.0;

　　end;

　　time=time+dt;

　　i=i+1;

　　end;

　　end;

　　excitation=excitation.*amp;