FDTD參數選擇程序
加入技術交流群
delay=(time-rickper);本文引用地址：http://www.j9360.com/article/149079.htm
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;
評論