This work was supported by the National Natural Science Foundation of China (Grant No. 11571052, 11731012), the Hunan Provincial Natural Science Foundation of China (Grant No. 2018JJ2417), and the Open Fund of Hunan Provincial Key Laboratory of Mathematical Modeling and Analysis in Engineering (Grant No. 2018MMAEZD02). The authors are grateful to anonymous referees and Professor Quansheng Liu for their very valuable comments and remarks, which significantly contributed to improving the quality of the paper.