In this thesis, we propose a method for penalized estimation of multiple changepoints in functional data sequences. We introduces a multiple changepoint detection method enhanced with a penalty function. This method eliminates the need to pre-specify the number of changepoints, enabling simultaneous detection of changepoint locations and quantities. Additionally, to mitigate the occurrence of overly close or consecutive changepoints, we introduce an additional constraint. We utilize two methods, the signal-to-noise ratio and sample splitting, to choose the optimal penalty parameter. Furthermore, we compare our method with others, including Multiple Changepoint Isolation (MCI) and the binary segmentation, demonstrating superior performance under some common scenarios.