FDTD参数选择程序

*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;