Abstract: Taking Bregman D-Function as a regularization function, we present a general numerical model to identify the multi-variables, including thermal parameters and boundary conditions. A time stepping scheme is used for the transient analysis and a homotopy method is adopted to solve the process. Both the data with measurement errors and the initial value are used to indentify these parameters. The numerical validation exhibits good results.