利用中心差分格式求解一阶波动方程(附Matlab代码)