The authors would like to thank the supports of the National Natural Science Foundation of China under Grant no. 11371075 , the Hunan Key Laboratory of mathematical modeling and analysis in engineering, China, and the Portuguese Funds through FCT–Fundacão para a Ciência e a Tecnologia, within the Project UID/MAT/00013/2013.