Stolt偏移和Gazdag相移法的编程实现