Traditional whole-brain dynamical modeling techniques are typically constrained by static single features, neglecting dynamic fluctuations in brain networks and lacking qualitative analysis of corresponding indicators, which limits modeling accuracy and comprehensibility. In order to address this issue, a multi-objective expectation maximization algorithm based on bifurcation analysis was proposed. This approach integrates a dynamic mean-field model with brain structural-functional features extracted from multi-mode imaging data for modeling purposes. Bifurcation theory was employed to qualitatively analyze multiple constraint indicators of the model, including functional connectivity, dynamic functional connectivity, and metastability for model inversion. Initial parameter values were determined through bifurcation analysis, and parameter combinations were iteratively refined using an expectation maximization algorithm. Quantitative analysis validates the accuracy and stability of this method.