matlab 程序,采用迎风格式求解波动方程