说明：  波动方程有限差分数值模拟，采用二阶差分，AEA自由表面边界条件
(Finite difference numerical simulation of wave equation using second order difference and AEA free surface boundary conditions)