The former two is the main program, the middle two is what went wrong, the last two are the emd subroutine,


(how are you, have problem has been solved, because it is because of the spline interpolation in p take is no, thanks,)
This is the author said, but I don't know how to specific modifications,