Research Article  Open Access
ChunPing Ren, NengJian Wang, ChunSheng Liu, "Identification of Random Dynamic Force Using an Improved Maximum Entropy Regularization Combined with a Novel Conjugate Gradient", Mathematical Problems in Engineering, vol. 2017, Article ID 9125734, 14 pages, 2017. https://doi.org/10.1155/2017/9125734
Identification of Random Dynamic Force Using an Improved Maximum Entropy Regularization Combined with a Novel Conjugate Gradient
Abstract
We propose a novel mathematical algorithm to offer a solution for the inverse random dynamic force identification in practical engineering. Dealing with the random dynamic force identification problem using the proposed algorithm, an improved maximum entropy (IME) regularization technique is transformed into an unconstrained optimization problem, and a novel conjugate gradient (NCG) method was applied to solve the objective function, which was abbreviated as IMENCG algorithm. The result of IMENCG algorithm is compared with that of ME, MECG, MENCG, and IMECG algorithm; it is found that IMENCG algorithm is available for identifying the random dynamic force due to smaller root meansquareerror (RMSE), lower restoration time, and fewer iterative steps. Example of engineering application shows that Lcurve method is introduced which is better than Generalized Cross Validation (GCV) method and is applied to select regularization parameter; thus the proposed algorithm can be helpful to alleviate the illconditioned problem in identification of dynamic force and to acquire an optimal solution of inverse problem in practical engineering.
1. Introduction
The inverse problem, which is a typical mathematical problem, has been widely applied in many fields of practical engineering, such as vibration isolation, health monitoring, and fault diagnosis [1–3]. The primary concern of mathematical workers and engineering technicians in the inverse problem derives from the development of the theory of solving the problems of universal existence and illconditioned problems in practical engineering applications in the past few years [4–6]. Identification of dynamic force between pick and coalseam has been a major concern in failure analysis of rotary machinery; once accurate force sources are obtained, it is possible to apply some advanced algorithms to analyze the characteristics of the dynamic force signals. However, we can not obtain the accurate dynamic force directly by the traditional methods due to the illposedness of systems [7–9]. Aiming to reduce the above problems, it is necessary to obtain a method for solving illconditioned problem.
Accordingly, regularization techniques, which can treat the illconditioned problem, have been utilized in many fields [10–13]. In the past few years, maximum entropy (ME) regularization techniques have been put forward one after another, which has been applied successfully in the wide areas of image reconstruction, signal processing, force identification problem, and so on [14], and the advantages of using ME regularization technique are as follows: firstly, some important information from the incomplete data can be extracted; then the probability distribution of the constraint is hidden in it; and, finally, entropy term can be changed according to different objects [15]. Hence, ME regularization technique mentioned above can lay a firm foundation for the steady and accurate solutions of illconditioned problem.
However, ME regularization technique is not perfect, which has been made better in many forms, and the means of improvement are generally concentrated in the entropy function combined with some optimization algorithms to solve regularization solution [16]. MohammadDjafari et al. [17] applied the improved entropy function combining the Lagrange method to process mass spectrometry data problems. Haffner et al. [18] discussed an improved maximum entropy modeling to handle with multiclass categorization problems using Newton iteration algorithm. In [19], by combining maximum entropy and Bayesian method, a new maximum entropy regularization is constructed to distribute force acting on boundary structure. Myrheim and Rue [20] presented a novel algorithm combining the ME method with the preconditioned conjugate gradient method to reconstruct image. In [21], the maximum entropy function and the conjugate gradient method were utilized, which can be used to solve absolute value equation problems. However, as for the improved entropy function mentioned, we consider that it is a numerical optimization algorithm, and its disadvantage is that the rate of convergence is not high, which is only suitable for simple objects, but the identified results are not ideal for complex objects. In fact, the entropy function can be redesigned according to the properties of different objects, and identification results of different objects can be achieved when we use different entropy functions combined with some different optimization algorithms [22].
In our study, a fast mathematical algorithm that can identify dynamic force between pick and coalseam is determined. The technique we propose is new to our research field. Significant improvements are made in comparison with our previous work.
First, the force identification procedure is transformed into an unconstrained optimization problem due to being suboptimal for conventional regularization methods. In this work, a novel optimization mathematical algorithm, combining a novel improved maximum entropy (IME) regularization technique and a novel conjugate gradient (NCG) method (for short IMENCG), is proposed to simplify the optimization process and improve the applicability of the algorithm.
Second, as a fast technique of IMENCG, the optimized objective function is redesigned using the 1norm estimation and the elongated entropy function in order to overcome the numerical instability.
Third, the optimal solution (stable solution) of inverse problem in practical engineering is studied. The technique proposed does not require any a priori information on the model for the random dynamic force identification. Only the displacement response measured by experiment is adequate for the identification algorithm.
Fourth, global convergence and stability of the solution for identified force model are proved to improve the precision and detailed information of random dynamic force identification.
Fifth, the Lcurve criteria, combined with the IME regularization technique, are applied to select the suitable regularization parameter for dynamic force identification in complex vibration system.
Finally, the applicability of the proposed technique is effectively demonstrated through an engineering application example. It is found that the performance of random dynamic force identification using the combination of IMENCG algorithm and Lcurve is better than that of other previous works, which can effectively reduce the illposedness of random dynamic force identification.
In the present work, we aim at identifying random dynamic force between pick and coalseam combining a novel improved maximum entropy (IME) regularization technique and a novel conjugate gradient (NCG) method. The structure of the paper is arranged as follows. In Section 2, the identification model of random dynamic force is established. In Section 3, the IMENCG algorithm is proposed. In Section 4, we present the engineering application and results discussion. Some important conclusions are obtained in Section 5.
2. Model of Identified Force
An abridged general view of random dynamic force between pick and coalseam structure is plotted in Figure 1. For deterministic MDOF structures, the kinetic equation can be expressed by the following form [23]:where , , and represent the mass matrix, damping matrix, and stiffness matrix, respectively; denotes the dynamic force column vector in different directions; , and are displacement, velocity, and acceleration response vector, respectively.
The displacement of the structure can be described as where is matrix of mode and ; is the displacement vector and time function in the generalized coordinates; .
By means of (1) and (2), the dynamic equations represented by the physical coordinates are transformed into the modal coordinate system, and the decoupled modal equations are obtained in the following form: in which is the order modal damping ratio; are the generalized mass and the generalized force, respectively, and and .
The initial velocity and initial displacement of the system are assumed to be and , respectively. So in (2) can be expressed as follows:
It is well known that can be obtained by the homogeneous equation of (1), which can be described as follows [24]:where and denote the corresponding values in modal coordinates, which can be described as follows:
can be expressed in the following form due to the fact it has nothing to do with the initial conditions of the system:
The displacement of the system can be expressed as follows:in which and .
The system force is identified according to the known displacement; the displacement is represented as follows:
The following expressions can be obtained through (5) and (7) which are introduced into (9):
In general, can be described as a convolution integral of the modal force and the corresponding Green kernel function:
Hence, Green’s kernel function can be described by comparing (10) and (11): Let , , and then t (11) was discretized, and it is transformed into a matrix form that can be expressed as follows:or simply noted as
Owing to the responses containing large noise, so (14) is expressed by the following form:where denotes the identified dynamic force, is the unknown noise, and , denotes a parameter level, denotes the standard deviation, and denotes a random number which ranges from −1 to 1.
To identify , we need to get and ; moreover, it was found that (15) is a typical illconditioned problem. In some cases, small varieties in the righthand side of (15) may lead to dramatic varieties in the solution [25–28]. In other words, they are very susceptible to error, in the inverse problem, this may directly lead to greater deviations of the identified results [29–31], and thus it is not realistic to solve (15) using traditional mathematical methods [32–34].
3. Algorithm
From the above analysis it can be seen that identification problem of dynamic force is illconditioned, whose solution is usually astaticism [35–38]. Therefore, it is necessary to find a way to achieve a stable solution. The regularization techniques are good method to alleviate the illconditioned problem. In the past few years, maximum entropy (ME) regularization techniques have been put forward one after another, which have been applied successfully in the wide area of practical engineering problem [39]. However, ME regularization technique is not very perfect; it has some disadvantages. Hence, in this work, an improved maximum entropy (IME) regularization technique is studied.
3.1. Maximum Entropy (ME) Regularization Method
It is generally known that maximum entropy (ME) regularization technique has been widely adopted in the past few years. The idea of the ME method is shown that the solving of (15) is transformed into that of an optimization problem [40]:where represents the regularization parameter.
Let , and it is called the entropy item. Thus, (16) can be replaced by the following:
3.2. An Improved Maximum Entropy (IME) Regularization Method
In some circumstances, such as complex identified objects, ME regularization technique is not very perfect, and it has some disadvantages of poor quality of identification. Hence, in this paper, the entropy item of ME regularization technique is redesigned using the 1norm estimation and the elongated entropy function according to the illposedness characteristics of dynamic force identification, which was named as an improved maximum entropy (IME) regularization technique, was proposed to process the dynamic force.
In this section, a function is used to design the entropy item which is expressed as follows according to 1norm estimation method:where represents 1 norm and .
And can be expressed by [41] Hence, (18) is described bywhere denotes the predetermined parameter.
Hence, a novel improved entropy item for dynamic force identification can be obtained, and it can be described as follows:
It is obvious that the proposed entropy item is different with that of ME regularization technique; it mainly considers the 1norm estimation and uses a better entropy item to make regular solutions more stable, thus which may be more suitable for the force identification of complex objects.
3.3. Selection of the Regularization Parameter
In this paper, we discuss Lcurve method and Generalized Cross Validation (GCV) due to the validity and the classical character in choosing the regularization parameter [35, 42].
According to the regularization theory, the regularization parameter λ for an improved maximum entropy (IME) regularization plays an important role in applying the numerical experiment process. First, the Lcurve criterion is introduced to select an optimal regularization parameter.
By reading the relevant literature, we know that the Lcurve criterion is a very powerful way for selecting a suitable value of the regularization parameter for many meaningful engineering and mathematical problems.
It is well known that the Lcurve is represented as a loglog plot of the norm of the residual versus the corresponding norm of a regularized solution.
The curve shows a typical “L” shape, and the optimal value for the regularization parameter is considered to be the one that corresponds to the corner of the curve [43].
Typically, the GCV function can be defined in the following form [44]:where denotes trace of matrix.
3.4. Solving of Objective Function
Equation (21) is considered as a class of unconstrained optimization problems and a series of optimization algorithms, which can treat such unconstrained optimization problem. In this work, a novel conjugate gradient (NCG) method was proposed to minimize (21). It has always been known that iterative form of conjugate gradient (CG) method are described as follows [45]:where denotes the step length, which can be computed using standard Armijo search conditionswhere and .
denotes the search directionwhere denotes the gradient of at the point , , and is conjugate parameter.
A novel search direction based on [46, 47] PRP conjugate gradient method is determined as follows: where
Algorithm 1.
Step 1. Given an initial point , , and , set and .
Step 2. If , then terminate; else go to the next step.
Step 3. is computed based on (24).
Step 4. Set .
Step 5. is computed based on (26).
Step 6. Set and go to Step .
Theorem 2. If the search direction is calculated by (26), and the step length meets the standard Armijo search conditions, for all , and then one has
Proof. When , we can obtain Then Theorem 2 is proved.
When , because the step length is calculated by formula (24), , and then we can obtainHence, the conclusion is established.
3.5. Global Convergence
Under the following assumptions, the global convergence of the algorithm is established.
Assumption A. (i) The objective function is bounded in the following level:(ii) The gradient is Lipschitz continuous; that is, there exists a constant such that, for any , we obtain
Lemma 3. Suppose that Assumption A holds, the step length is calculated by (24) and the search direction is calculated by (26), and there is a constant , then
Proof. Set , because the step length is calculated by (24), then .
When , according to Theorem 2, we can obtain the following:Because , then we obtainThen one hasWe haveHence, (33) is proved.
When , we know that can not meet (24) according to the definition of the step length, and thenAccording to the mean value theorem, we havewhere .
According to the formula, we obtainSo we haveHence, (33) is proved.
We can see that is monotone bounded according to (24) and Assumption A, so we can obtain the following conclusion based on (24) and (33) and Lemma 3:Hence, (34) is proved.
Theorem 4. Suppose that the step length and the search direction are generated by IMENCG algorithm, and then one has that
Proof. Owing to the fact that the step length and the search direction are generated by (24) and (26), respectively, thenSo we can obtainHence, (34) is proved.
4. The Engineering Application
4.1. Experimental SetUp
In this paper, a diagram of experimental setup is organized in Figure 2. The size of the coalseam structure specimen manufactured artificially was mm, the pick is mounted on the rotary cutting arm, it can be driven by the reducer and torque, and its power driven by the motor is rated at 55 kW [48]. The cutting bench achieves free forward and backward movement by the hydraulic pressure drive control system. System output responses can be measured by the corresponding sensors, then converted by a signal amplifier, and finally recorded by V10Dasp data vibration signal acquisition system [49, 50].
4.2. The Engineering Application
All codes were written in MATLAB 7.0 and run on a HP with 2.0 GB RAM and Windows 7 operating system. Stop the iteration if criterion is defined as being satisfied or run time more than 500 seconds. And the step length is calculated by (24), and the parameters are obtained as follows: , , , ,, , and . For performance analysis of dynamic force identification techniques used, the performance measurement metrics for dynamic force are described as follows: restoration time, iterative steps, and root meansquareerror (RMSE), where RMSE is defined bywhere is the actual measured random dynamic force and is the identified random dynamic force.
In this section, the material parameters of experimental conditions are shown in Tables 1 and 2.


It is generally known that direct numerical treatment fails to solve the present inverse problem in the engineering application. Therefore, the proposed method is to be applied to identify dynamic force. However, the key problem is how to select the suitable regularization parameter to obtain the optimal solution. Moreover, the suitable regularization parameter obtained by Lcurve is really important for numerical experiment process.
First, the graphical illustration of Lcurve was shown in Figure 3(a) according to Lcurve criterion. Let it be noted that it is very possible to obtain the optimal solution by combining the Lcurve criterion with the proposed regularization method.
(a) The graphical illustration of Lcurve
(b) λ = 10−1
(c) λ = 10−2
(d) λ = 10−5
(e) λ = 10−7
It is clearly seen from Figure 3(a) that when the regularization parameter is small, the graph shows very clearly that the norm of regularized solution is large, but the norm of the corresponding residual is small; we can speculate that it is likely to be caused by measurement errors. However, when the regularization parameter is large, it is found that the norm of regularized solution is a poor approximation; on the contrary, the norm of the corresponding residual is large.
The corner of the Lcurve indicates this transition, owing to the fact it shows a compromise between the minimization of the norm of the residual and norm of the solution . The selection of the regularization parameter of IME regularization can be not completely guaranteed to be suitable for all systems with the illposedness. However, many examples show that the Lcurve criterion is a powerful method for determining a suitable value of the regularization parameter for many significant engineering and mathematical problems.
In order to further discuss the influence of Lcurve on numerical experiments, then make sure that the optimal regularization parameter is selected from the Lcurve. According to the Lcurve criterion, the value of different regularization parameters are selected from Figure 3(a) to find the optimal regularization parameter. The measured displacement response was shown in Figure 4. And the identified results were formed in Figures 3(b), 3(c), 3(d), and 3(e).
Figures 3(d) and 3(e) show that the identified result is extremely unsatisfactory due to small regularization parameter. However, we can see from Figure 3(b) that even if when the regularization parameter is large, that is, beyond a certain threshold, the identified result is far away from what is normal. At the same time, the identified result of the value at the corner of the Lcurve is shown in Figure 3(c). Figure 3(c) clearly shows that the most accurate results are obtained under the optimal regularization parameter (λ = 10^{−2}) compared with other cases of regularization parameters. So we conclude that the values at the corner of the Lcurve are selected as the most suitable regularization parameter.
As we can see from the above discussion, the regularization parameter was obtained by Lcurve method, and . Then identified dynamic force can be given based on measured displacement response using the identified modal model and the proposed algorithm, and the identified results were formed in Figure 5 by assembling the analysis results of ME, MECG, MENCG, IMECG, and IMENCG algorithm.
(a) ME algorithm
(b) MECG algorithm
(c) MENCG algorithm
(d) IMECG algorithm
(e) IMENCG algorithm
It was obtained that the random dynamic force is easily identified with application of the above identification algorithms, and the identified random dynamic force is obviously observed using the regularization parameter selected by Lcurve method as seen in Figures 5(a), 5(b), 5(c), 5(d), and 5(e), respectively.
Figure 5(a) shows that the identified result is obtained using ME algorithm, Figure 5(b) shows that the identified result is obtained using MECG algorithm, identified result using MENCG algorithm is obtained from Figure 5(c), Figure 5(d) shows that the identified result is obtained using IMECG algorithm, and the identified result using IMENCG algorithm is shown in Figure 5(e).
Moreover, it can be shown that different identification algorithms can commendably identify the random dynamic force. However, to some extent, the performance of random dynamic force identification has a certain difference.
We analyzed the different identification algorithms at And restoration time, iterative steps, and RMSE were used to evaluate the quality of dynamic force identification.
In Table 3, it was found that the proposed IMENCG algorithm is associated with minimum RMSE, minimum restoration time, and minimum iterative steps in comparison to ME, MECG, MENCG, and IMECG algorithm, which proves that the IMENCG algorithm is superior to the other algorithms in the identification of random dynamic force at , which is closer to the actual dynamic force.

However, as is well known, different regularization parameters are obtained by different regularization parameter selection methods, which will lead to different identified results. In this paper, the GCV is applied to compute the regularization parameter compared with Lcurve method.
Figure 6 shows that the regularization parameter was obtained by GCV method, when , reaches the minimum, and the identified results were formed in Figure 7 by assembling the analysis results of different algorithms at .
(a) ME algorithm
(b) MECG algorithm
(c) MENCG algorithm
(d) IMECG algorithm
(e) IMENCG algorithm
In Figures 7(a), 7(b), 7(c), 7(d), and 7(e), it can be shown that five different identification algorithms can identify the random dynamic force, and it was found that the proposed IMENCG algorithm is associated with minimum RMSE, minimum restoration time, and minimum iterative steps in comparison to the other algorithm from Table 4; hence, we can draw a conclusion that the IMENCG algorithm is closer to the actual force at , which shows that the IMENCG algorithm is superior to the other algorithms whether the regularization parameter is selected by using Lcurve or GCV method.

From Tables 3 and 4, we can see that RMSE, restoration time, and iteration steps using Lcurve method are lower than those of using GCV method whether using IMENCG algorithm or the other algorithms.
Hence, we can conclude that the IMENCG algorithm combined with Lcurve method is available to identify the dynamic force.
5. Conclusions
In this paper, the identification model for random dynamic force between pick and coalseam is established through the relationship between the identified force vector and the displacement responses measured by machinery vibration system. Then the entropy item of maximum entropy (ME) regularization is redesigned using the 1norm estimation method and the elongated entropy function according to the illconditioned characteristics of force identification, which was named as an improved maximum entropy (IME) regularization technique, was proposed to process the random dynamic force. Then the force identification problem is transformed into an unconstrained optimization problem and a novel conjugate gradient (NCG) method was applied to solve the objective function, which was called IMENCG algorithm.
The results show that the proposed IMENCG algorithm is associated with minimum RMSE, minimum restoration time, and minimum iterative steps in comparison to ME, MECG, MENCG, and IMECG algorithm, and Lcurve method is introduced which is better than GCV method and is applied to select the regularization parameter. Thus we can safely conclude that the combination of IMENCG algorithm and Lcurve method can reduce the illposedness of random dynamic force identification and achieve the optimal solution of inverse problem in practical engineering.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the Chinese National Natural Science Foundation (Contract nos. 51674106 and 51274091).
References
 T. Poggio, V. Torre, and C. Koch, “Computational vision and regularization theory,” Nature, vol. 317, no. 6035, pp. 314–319, 1985. View at: Publisher Site  Google Scholar
 T. S. Jang, H. G. Sung, S. L. Han, and S. H. Kwon, “Inverse determination of the loading source of the infinite beam on elastic foundation,” Journal of Mechanical Science & Technology, vol. 22, no. 12, pp. 2350–2356, 2008. View at: Publisher Site  Google Scholar
 J. M. Varah, “A practical examination of some numerical methods for linear discrete illposed problems,” SIAM Review, vol. 21, no. 1, pp. 100–111, 1979. View at: Publisher Site  Google Scholar
 Y. Dai, L. Liu, and S. Feng, “On the identification of coupled pitch and heave motions using oppositionbased particle swarm optimization,” Mathematical Problems in Engineering, vol. 2014, Article ID 784049, 10 pages, 2014. View at: Publisher Site  Google Scholar
 W. Luo, “Parameter identifiability of ship manoeuvring modeling using system identification,” Mathematical Problems in Engineering, vol. 2016, Article ID 8909170, 10 pages, 2016. View at: Publisher Site  Google Scholar
 L. Shi, H. He, and W. Yan, “Prestress force identification for externally prestressed concrete beam based on frequency equation and measured frequencies,” Mathematical Problems in Engineering, vol. 2014, Article ID 840937, 13 pages, 2014. View at: Publisher Site  Google Scholar
 J. Gao, D. Wang, and J. Peng, “A Tikhonovtype regularization method for identifying the unknown source in the modified Helmholtz equation,” Mathematical Problems in Engineering, vol. 2012, Article ID 878109, 13 pages, 2012. View at: Publisher Site  Google Scholar
 X. Liu, S. Liu, and H. Ji, “Mechanism of rock breaking by pick assisted with water jet of different modes,” Journal of Mechanical Science and Technology, vol. 29, no. 12, pp. 5359–5368, 2015. View at: Publisher Site  Google Scholar
 J. M. Varah, “On the numerical solution of illconditioned linear systems with applications to illposed problems,” SIAM Journal on Numerical Analysis, vol. 10, pp. 257–267, 1973. View at: Publisher Site  Google Scholar  MathSciNet
 K. M. Hanke and L. P. C. Hansen, “Regularization methods for largescale problems,” Surveys on Mathematics for Industry, vol. 3, no. 4, pp. 253–315, 1993. View at: Google Scholar  MathSciNet
 E. Talic, A. Schirrer, M. Kozek, and S. Jakubek, “Multiobjective parameter identification of EulerBernoulli beams under axial load,” Journal of Sound and Vibration, vol. 341, pp. 86–99, 2015. View at: Publisher Site  Google Scholar
 T. S. Jang, H. Baek, S. L. Hana, and T. Kinoshita, “Indirect measurement of the impulsive load to a nonlinear system from dynamic responses: inverse problem formulation,” Mechanical Systems and Signal Processing, vol. 24, no. 6, pp. 1665–1681, 2010. View at: Publisher Site  Google Scholar
 Y. Jia, Z. Yang, N. Guo, and L. Wang, “Random dynamic load identification based on error analysis and weighted total least squares method,” Journal of Sound and Vibration, vol. 358, pp. 111–123, 2015. View at: Publisher Site  Google Scholar
 S. J. Phillips and R. E. Schapire, “Maximum entropy density estimation with generalized regularization and an application to species distribution modeling,” Journal of Machine Learning Research (JMLR), vol. 8, no. 2, pp. 1217–1260, 2007. View at: Google Scholar  MathSciNet
 U. Amato and W. Hughes, “Maximum entropy regularization of Fredholm integral equations of the first kind,” Inverse Problems, vol. 7, no. 6, 793 pages, 1991. View at: Publisher Site  Google Scholar
 Y.W. Chiang, P. P. Borbat, and J. H. Freed, “Maximum entropy: A complement to Tikhonov regularization for determination of pair distance distributions by pulsed ESR,” Journal of Magnetic Resonance, vol. 177, no. 2, pp. 184–196, 2005. View at: Publisher Site  Google Scholar
 A. MohammadDjafari, J.F. Giovannelli, G. Demoment, and J. Idier, “Regularization, maximum entropy and probabilistic methods in mass spectrometry data processing problems,” International Journal of Mass Spectrometry, vol. 215, no. 13, pp. 175–193, 2002. View at: Publisher Site  Google Scholar
 P. Haffner, S. Phillips, and R. Schapire, “Efficient multiclass implementations of L1regularized maximum entropy,” Computer Science, abs/cs/0506101, 2005. View at: Google Scholar
 T. W. Grandy, “Maximumentropy and Bayesian methods in inverse problems,” Acta Applicandae Mathematica, vol. 20, no. 12, pp. 189–191, 1990. View at: Google Scholar
 J. Myrheim and H. Rue, “New algorithms for maximum entropy image restoration,” Graphical Models & Image Processing, vol. 54, no. 3, pp. 223–238, 1992. View at: Publisher Site  Google Scholar
 W. Zhu and T. Yan, “The conjugate gradient algorithm for solving the absolute value equations based on maximum entropy function,” Journal of Hefei Normal University, 2016. View at: Google Scholar
 G. Nicolas, J. Andrew, and C. C. Finlay, “Maximum entropy regularization of timedependent geomagnetic field models,” Geophysical Journal of the Royal Astronomical Society, vol. 171, no. 3, pp. 1005–1016, 2010. View at: Google Scholar
 J. Liu, X. Sun, X. Han, C. Jiang, and D. Yu, “Dynamic load identification for stochastic structures based on Gegenbauer polynomial approximation and regularization method,” Mechanical Systems and Signal Processing, vol. 56–57, pp. 35–54, 2015. View at: Publisher Site  Google Scholar
 J. Liu, X. Sun, X. Meng et al., “A novel shape function approach of dynamic load identification for the structures with interval uncertainty,” International Journal of Mechanics & Materials in Design, no. 3, pp. 12, 2015. View at: Google Scholar
 L. Wang, X. Han, J. Liu, X. He, and F. Huang, “A new regularization method and application to dynamic load identification problems,” Inverse Problems in Science and Engineering, vol. 19, no. 6, pp. 765–776, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 W. Chen, “A new extrapolation algorithm for bandlimited signals using the regularization method,” IEEE Transactions on Signal Processing, vol. 41, no. 3, pp. 1048–1060, 1993. View at: Publisher Site  Google Scholar
 V. Melicher and V. Vrabel, “On a continuation approach in Tikhonov regularization and its application in piecewiseconstant parameter identification,” Inverse Problems, vol. 29, no. 11, pp. 193–212, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 Z. Zhao and Z. Meng, “A modified Tikhonov regularization method for a backward heat equation,” Inverse Problems in Science and Engineering, vol. 19, no. 8, pp. 1175–1182, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 H. Cheng, X.L. Feng, and C.L. Fu, “A mollification regularization method for the Cauchy problem of an elliptic equation in a multidimensional case,” Inverse Problems in Science and Engineering, vol. 18, no. 7, pp. 971–982, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 A. Qian and Y. Gui, “Determining an unknown source in the heat equation by a mollification regularization method,” in Proceedings of the 2010 International Conference on Computational Intelligence and Software Engineering (CiSE '10), pp. 1–4, IEEE, Wuhan, China, December 2010. View at: Publisher Site  Google Scholar
 F. Yang, C.L. Fu, and X.X. Li, “A mollification regularization method for identifying the timedependent heat source problem,” Journal of Engineering Mathematics, vol. 522, no. 18, pp. 1–14, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 S. Marchesiello and L. Garibaldi, “A time domain approach for identifying nonlinear vibrating structures by subspace methods,” Mechanical Systems and Signal Processing, vol. 22, no. 1, pp. 81–101, 2008. View at: Publisher Site  Google Scholar
 E. Lourens, E. Reynders, G. De Roeck, G. Degrande, and G. Lombaert, “An augmented Kalman filter for force identification in structural dynamics,” Mechanical Systems and Signal Processing, vol. 27, no. 1, pp. 446–460, 2012. View at: Publisher Site  Google Scholar
 M. S. Allen and T. G. Carne, “Delayed, multistep inverse structural filter for robust force identification,” Mechanical Systems and Signal Processing, vol. 22, no. 5, pp. 1036–1054, 2008. View at: Publisher Site  Google Scholar
 T. Reginska, “A regularization parameter in discrete illposed problems,” Society for Industrial and Applied Mathematics, 1996. View at: Google Scholar  MathSciNet
 B. Movahedian and B. Boroomand, “Inverse identification of timeharmonic loads acting on thin plates using approximated Green's functions,” Inverse Problems in Science and Engineering, vol. 24, no. 8, pp. 1475–1493, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 H. Du and M. Cui, “Approximate solution of the Fredholm integral equation of the first kind in a reproducing kernel Hilbert space,” in Proceedings of the International Conference on Mathematics and Engineering Techniques in Medicine and Biological Scienes (METMBS '04), DBLP, pp. 322–328, Las Vegas, Nevada, USA, June 2124, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 C. Xiang, Y. Tong, Z. Wu et al., “A new method for tsvd regularization truncated parameter selection,” Mathematical Problems in Engineering, vol. 2013, Article ID 161834, 9 pages, 2013. View at: Publisher Site  Google Scholar
 H. Gzyl, “Illposed linear inverse problems and maximum entropy in the mean,” Acta Científica Venezolana, vol. 53, no. 2, pp. 74–93, 2002. View at: Google Scholar  MathSciNet
 J.M. Loubes and B. Pelletier, “Maximum entropy solution to illposed inverse problems with approximately known operator,” Journal of Mathematical Analysis and Applications, vol. 344, no. 1, pp. 260–273, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 L. Jing, S. Liu, L. Zhihong, and S. Meng, “An image reconstruction algorithm based on the extended Tikhonov regularization method for electrical capacitance tomography,” Measurement, vol. 42, no. 3, pp. 368–376, 2009. View at: Publisher Site  Google Scholar
 M. E. Hochstenbach, L. Reichel, and G. Rodriguez, “Regularization parameter determination for discrete illposed problems,” Journal of Computational and Applied Mathematics, vol. 273, pp. 132–149, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 P. C. Hansen and D. P. O'Leary, “The use of the Lcurve in the regularization of discrete illposed problems,” SIAM Journal on Scientific Computing, vol. 14, no. 6, pp. 1487–1503, 1993. View at: Publisher Site  Google Scholar  MathSciNet
 G. H. Golub, M. Heath, and G. Wahba, “Generalized crossvalidation as a method for choosing a good ridge parameter,” Technometrics. A Journal of Statistics for the Physical, Chemical and Engineering Sciences, vol. 21, no. 2, pp. 215–223, 1979. View at: Publisher Site  Google Scholar  MathSciNet
 X.P. Wu, L.Y. Liu, F.J. Xie, and Y.F. Li, “A new conjugate gradient algorithm with sufficient descent property for unconstrained optimization,” Mathematical Problems in Engineering, vol. 2015, Article ID 352524, 8 pages, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 L. Zhang, W. Zhou, and D.H. Li, “A descent modified PolakRibièrePolyak conjugate gradient method and its global convergence,” Ima Journal of Numerical Analysis, vol. 26, no. 4, pp. 629–640, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 C. S. Liu, “Novel algorithms based on the conjugate gradient method for inverting illconditioned matrices, and a new regularization method to solve illposed linear systems,” Computer Modeling in Engineering and Sciences, vol. 60, no. 3, pp. 279–308, 2010. View at: Google Scholar  MathSciNet
 C. P. Ren and C. S. Liu, “Experimental research on properties of analog materials of coal and rock,” Journal of Heilongjiang University of Science and Technology, vol. 24, no. 6, pp. 581–584, 2014. View at: Publisher Site  Google Scholar
 C.S. Liu, C.P. Ren, and D.G. Li, “Reconstruction and deduction of cutting coal and rock load spectrum on modified discrete regularization algorithm,” Journal of China Coal Society, vol. 39, no. 5, pp. 981–986, 2014. View at: Publisher Site  Google Scholar
 C. S. Liu, C. P. Ren, and F. Han, “Study on timefrequency spectrum characteristic of dynamic cutting load based on wavelet regularization,” Applied Mechanics and Materials, vol. 577, pp. 196–200, 2014. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 ChunPing Ren et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.