This research introduces a novel and unified optimization design method for multi-stage non-circular gear transmission (MNCGT) to address the challenges in designing MNCGT for complex motion requirements. The method optimizes non-circular gears comprehensively, reducing design and manufacturing difficulties while ensuring the realization of specified transmission requirements. A unified parameterization method, grounded on periodic B-spline interpolation, is introduced to establish the transmission function of non-circular gears and map it to a finite unified variable space. This innovative approach effectively reduces the difficulty and constraints of MNCGT optimization design. The proposed method takes into account crucial factors such as non-circularity, smoothness, processing conditions, and contact ratio, which significantly impact the transmission performance and manufacturing feasibility of non-circular gears. The effectiveness and superiority of this method are demonstrated through two practical examples and a real-world application in a planetary gear transplant mechanism, highlighting its potential for solving complex engineering problems.