1维严格耦合波方法matlab代码