一维FDTD色散媒质.docxVIP

%%%%%%%% Frequency-dependent media %%%%%%%%%%%%%%%% More on One-dimensional simulation 2_3 %%%%%%%%%%%%%%%% 1D FDTD simulation %%%%%%%%KE=200; %%%%%%%% space long %%%%%%%%ex=zeros(1,KE+1);hy=zeros(1,KE+1);dx=zeros(1,KE+1);ix=zeros(1,KE+1);sx=zeros(1,KE+1);ga=zeros(1,KE+1);gb=zeros(1,KE+1);gc=zeros(1,KE+1);mag=zeros(1,KE+1);arg=zeros(1,3);%%%%%%%% Real and imaginary parts of the Fourier Transform %%%%%%%%real_pt=zeros(3,KE+1);imag_pt=zeros(3,KE+1);%%%%%%%% Amplitude and phase of the Fourier Transform %%%%%%%%ampn=zeros(3,KE+1);phasen=zeros(3,KE+1);%%%%%%%% Amplitude and phase of input pulse %%%%%%%%amp_in=zeros(1,3);phase_in=zeros(1,3);%%%%%%%% Fourier Trans. of input pulse %%%%%%%%real_in=zeros(1,3);imag_in=zeros(1,3);%%%%%%%% Parameters %%%%%%%%ex_low_m1=0;ex_low_m2=0;ex_high_m1=0;ex_high_m2=0; T=0; NSTEP=500; %%%%%%%% time step %%%%%%%%kc=KE/2; %%%%%%%% Medium interface %%%%%%%%t0=50 ; %%%%%%%% Gaussian pulse parameter %%%%%%%%spread=10; %%%%%%%% Gaussian pulse parameter %%%%%%%%ddx=0.01; %%%%%%%% set the cell size to 1cm %%%%%%%%dt=ddx/(2*3e8); %%%%%%%% calculate the time step %%%%%%%%freq_in=7*1e8; %%%%%%%% signal frequency %%%%%%%%epsilon=2; %%%%%%%% relative dielectric constant %%%%%%%%kstart=100; %%%%%%%% position of the dielectric medium %%%%%%%%epsz=8.85419e-12; %%%%%%%% absolute dielectric constant %%%%%%%%sigma=0.01; %%%%%%%% conductivity %%%%%%%%chil=2;tau=1e-9;%%%%%%%% signal frequency %%%%%%%%freq=[50.e6,200.e6,500.e6]; for i=1:3arg(i)=2*pi*freq(i)*dt;end%%%%%%%% initialize to free space %%%%%%%%for n=1:KE+1 ga(n)=1; gb(n)=0; gc(n)=0;enddel_exp=exp(-dt/tau);for n=kstart:KE+1ga(n)=1/(epsilon+(sigma*dt/epsz)+chil*dt/tau);gb(n)=sigma*dt/epsz;gc(n)=chil*dt/tau;end for n=1:NSTEP T=T+1;%%%%%%%% calculate dx field %%%%%%%% for k=2:KE+1 dx(k)=dx(k)+0.5*(hy(k-1)-hy(k));