demo0.m

来自「这是一个 York大学开发的FDTD代码。很好用的工具」· M 代码 · 共 20 行

disp('\nSetup for demos')

disp('\nSet Pulse width and Period')
Tw=0.25
Tp=1.0

disp('\nSet number of points and determine time step')
N=100
NX=400

disp('\n set number of harmonics for Fourier series')
K=15

disp('\n Time constants for decaying exponentials')
Tor=[0.2,1,1,0.2,0.2,0.5]
f0=[5,6,19,20.5,38.5,40.2]
%f0=[10,10,10,10,10,10]
Phi=[0,pi/4,pi/2,pi*3/4,pi/3,pi/5]
zeta=[0.01,0.02,0.05,0.1,0.2,0.5]