zyh算法