经典chang算法,maglab实现