The authors would like to thank the supports of National Natural Science Foundation of China (Grant No. 11371075 and No. 11771453), Hunan Key Laboratory of Mathematical Modeling and Analysis in Engineering, Natural Science Foundation of Hunan Province, China (Grant No. 2020JJ5267) and Scientific Research Funds of Hunan Provincial Education Department (Grant No. 18C0877 and No. 19B301).