Collaborative power tracking method of diversified thermal loads for optimal demand response: A MILP-Based decomposition algorithm. (1st December 2022)