Empirical mode decomposition (EMD) and its variants are adaptive algorithms that decompose a time series into a few oscillation components called intrinsic mode functions (IMFs). They are powerful signal processing tools and have been successfully applied in many applications. Previous research shows that EMD is an efficient algorithm with computational complexity π(π) for a given number of IMFs, where π is the signal length, but its memory is as large as (13 + π πππ )π, where π πππ is the number of IMFs. This huge memory requirement hinders many applications of EMD. A physical or physiological oscillation (PO) mode often consists of a single IMF or the sum of several adjacent IMFs. Let π ππ’π‘ denote the number of PO modes and, by definition, π ππ’π‘ β€ π πππ . In this paper, we will propose a low memory cost implementation of EMD and prove that the memory can be optimized to (2 + π ππ’π‘ )π without aggravating the computational complexity, while gives the same results. Finally, we discuss the optimized memory requirements for different noise-assisted EMD algorithms.INDEX TERMS EMD, EEMD, CEEMD, memory cost.