matlab，fdtd情况下的分析衍射等

The source is located at (20,200), and the total simulation time is 10 ns. The excitation waveform is plotted in Fig. 2, and the pulse response at 10 ns is shown in Fig. 3. The propagating wavefront is clearly seen in Fig. 3, together with evidence of a small reflection which has arisen due to the imperfect nature of the ABC on the left-hand side of the computa- tional domain. It is not meaningful to extract the time-harmonic response from this result, as the solution had not reached steady state. 

fdtd_1.m