Research Article  Open Access
Adaptive PD Control Based on RBF Neural Network for a WireDriven Parallel Robot and Prototype Experiments
Abstract
An adaptive PD control scheme is proposed for the support system of a wiredriven parallel robot (WDPR) used in a wind tunnel test. The control scheme combines a PD control and an adaptive control based on a radial basis function (RBF) neural network. The PD control is used to track the trajectory of the end effector of the WDPR. The experimental environment, the external disturbances, and other factors result in uncertainties of some parameters for the WDPR; therefore, the RBF neural network control method is used to approximate the parameters. An adaptive control algorithm is developed to reduce the approximation error and improve the robustness and control precision of the WDPR. It is demonstrated that the closedloop system is stable based on the Lyapunov stability theory. The simulation results show that the proposed control scheme results in a good performance of the WDPR. The experimental results of the prototype experiments show that the WDPR operates on the desired trajectory; the proposed control method is correct and effective, and the experimental error is small and meets the requirements.
1. Introduction
Neural networks are widely used in engineering, especially in the field of industrial automation [1]. The application is also very popular for robot control. A neural network is based on the biological neural networks of the human brain from the perspective of information processing and represents a simple model that forms different networks according to different connection modes. Why have neural networks attracted broad attention in the field of robot control? The reason is that this approach meets the basic requirements of the control theory and control systems engineering and helps to solve relevant problems [2]. Neural networks fully approximate complex nonlinear mapping relationship, possess the characteristics of learning and adapting to the dynamic characteristics of uncertain systems, and exhibit strong robustness and fault tolerance [3]. As a result, a large number of researchers have applied neural networks for the control of parallel robots in nonlinear systems.
Dinh et al. [4] proposed a dynamic neural network (DNN) based on a robust observer for uncertain nonlinear systems. The observer structure consists of a DNN, a dynamic filter, and a slidingmode feedback term. Simulations and experiments using a twolink robot manipulator were performed to demonstrate the effectiveness of the proposed method in comparison with several other state estimation methods. Wang et al. [5] used a recursive neural network to optimize the structure of complex robots. Melingui et al. [6] proposed two subcontrollers for a class of bionic continuum robots termed the “compact bionic handling arm” (CBHA). One controller is encapsulated in the other and both are implemented in real time to control the CBHA’s endeffector position. The first subcontroller controls the CBHA’s kinematics based on a distal supervised learning scheme. The second subcontroller controls the CBHA’s kinetics based on an adaptive neural control. The experimental results of using a CBHA robot demonstrated the accurate tracking of the CBHA’s endeffector position. Sabahi et al. [7] proposed a new indirect type2 fuzzy neural network predictive (T2FNNP) controller for a class of nonlinear inputdelay systems in the presence of unknown disturbances and uncertainties. The proposed T2FNNP controller was applied to a nonlinear inverted pendulum and singlelink robot manipulator systems with input timevarying delay; the performance was compared with that of a type1 fuzzy sliding predictive (T1FSP) controller. The simulation results showed the efficiency of the proposed T2FNNP controller. Wen et al. [8] proposed an improved radial basis functional (RBF) neural network by means of a hybrid force/position control for robotic manipulation; good stability and transient performance of the system were achieved. Zhao et al. [9, 10] proposed an adaptive outputfeedback control for a class of nonsmooth nonlinear systems and an adaptive fuzzy hierarchical slidingmode control method for a class of multiinput multioutput unknown nonlinear timedelay systems with input saturation. Wang et al. [11, 12] presented an adaptive neural control approach for nonstrictfeedback nonlinear systems in the presence of unmodeled dynamics, unknown control directions, and input deadzone nonlinearity and proposed a robust adaptive tracking control for nonaffine stochastic nonlinear switching systems.
At present, most of the literature on the use of neural networks for robot control is focused on rigid robots (arms are rigid, such as link robots) and papers on multiDOF wiredriven parallel robots (WDPRs) used in wind tunnel tests are rarely seen. In this study, we propose a new adaptive PD control scheme based on an RBF neural network for a WDPR with 6 DOF and driven by 8 wires; the nature and task requirements are based on the WDPR described in reference [13]. Compared with other types of neural networks, RBF neural networks avoid the issue of local minimization and slow training and satisfy the requirements of realtime control. In the RBF neural network compensation control strategy, it is important to ensure the stability of the closedloop control system. At the same time, in the designed control system, the RBF neural network can approximate the unknown function [14]. One advantage of the proposed control method is that the RBF neural network and adaptive PD control are integrated into the WDPR, and other advantages of the proposed control method are verified by the prototype experiment.
In this study, the main purpose is to design a new tracking control scheme for a 6DOF 8WDPR used in a wind tunnel test. The main contributions of this study are the following: (1) the PD control is used to track the trajectory of the end effector; (2) the RBF neural network approximation is used to compensate for the uncertainties of the system, such as external disturbances; (3) the adaptive control is used to achieve a high trajectory tracking performance of the WDPR; (4) the simulation results and the prototype experimental results show that the RBF neural network results in good control performance of the WDPR and meets the requirements of the wind tunnel test.
The rest of the paper is organized as follows. Section 2 presents the problem. Section 3 describes the RBF network design of the WDPR. The control design and the stability analysis of the WDPR are described in Section 4. The simulation results are presented in Section 5. The experimental results of the prototype are provided in Section 6. The conclusions are given in Section 7.
2. Problem Description
2.1. Description of WDPR
A WDPR is a specially designed support system and a parallel device for wind tunnel tests. The robot converts the motion and force of the driver to the motion and force of the endeffector in parallel using the wire as the medium. The system has the advantages of a simple structure, high rigidity, large workspace, low inertia, fast movement, easy assembly and disassembly, and reconfiguration [15, 16]. It is a complex multipleinput and multipleoutput nonlinear system. The structural diagram of the WDPR is shown in Figure 1. Figure 2 shows the WDPR in the wind tunnel. In Figure 2, the end effector is a standard dynamics model (SDM). The end effector (aircraft model) is located on the moving platform and it is driven by eight wires with six DOF.
In this study, the WDPR is composed of a control system, a driving system, an executive system, a measurement system, and other components. The driving system represents the power source to guarantee the movement of the end effector and is used as the input for the motion control. The dynamic analysis and modeling of the end effector are the basis of the motion control and the key to ensuring the system performance. The two dynamic equations are shown as follows:where is equivalent to the inertia matrix of the driver, is equivalent to the viscous friction coefficient matrix of the driver, is the angle of the motor, is the load torque produced by wire tension, , is the transmission coefficient of the ball screw, is the wire tension vector, is the output torque vector of the driver, is the inertial matrix of the end effector, is the pose of the end effector, is the pose velocity of the end effector, the relationship between the angle of the motor and the pose velocity is , is a transition matrix based on the attitude angle, is the Jacobi matrix, is a nonlinear Coriolis centrifugal force matrix, is the gravity vector of the end effector, , and is subject to the external dynamic loads at the end effector.
By combining the dynamic equation of the driving system with the dynamic equation of the model, the dynamic equation of the entire system is as follows [13]:where , . The meanings and derivations of the parameters and are detailed in reference [13], because the derivation of (2) and its parameters are relatively complex.
Equation (2) is multiplied by on both sides and (3) is obtained.Let , , , .
Equation (3) is simplified as
2.2. Problem Formulation of the WDPR
In practical engineering applications, there are many uncertainties, such as the experimental environment, the external disturbances, and other factors. These factors cannot be taken into account when the dynamic model of the system is developed. If those factors are taken into account, the dynamic equation will become too complex to be solved. Thus, the dynamic equations for an established system are generally incomplete.
Because it is impossible to obtain a complete dynamic equation of a complex nonlinear system in practical engineering applications, the nominal model of the WDPR is defined as follows:where , , and are known parameters, and represent matrices with 6 rows and 6 columns, and represents a matrix of 6 rows and 1 column.
Let , , and ; then, combine (4) with (5) to be simplified as
Equation (6) is rearranged and (7) is obtained:
Due to the error between the whole dynamic model and the nominal model, the unknown disturbance term is obtained. Let , then .
In the torque control algorithm, the nominal model is used; if is known, the control can be designed as follows:where and are the proportional gain and differential gain, respectively; they are positive definite matrices.
Equation (8) is substituted into (4) and we obtainwhere and is a matrix with 8 rows and 6 columns.
Equation (9) is rearranged and the error system equation (10) that can be obtained by simplification is as follows:
In (10), is the pose error, is the pose velocity error, and is the pose acceleration error. is the theoretical pose of the end effector.
The purpose of this study is to design a stable robust controller based on a nominal model. In practical engineering applications, is usually unknown; therefore, it is necessary to estimate and compensate for it.
3. RBF Neural Network Design of the WDPR
An inaccurate robot model can reduce the tracking performance of the robot. Therefore, neural networks are widely used in robot control due to their learning linearity and nonlinear mapping [17–20]. In order to compensate for the unknown item , an RBF neural network is used to approximate the unknown item and further improve the control accuracy of the WDPR.
The RBF neural network with compensated control is a network structure with three layers based on the practical application and the input of the WDPR, as shown in where is the input vector comprised of the pose error and the velocity error; is the output of the RBF neural network; is the output of the Gaussian basis function, which is a column vector with 5 rows and 1 column; is the weight of the RBF neural network and is a matrix with 5 rows and 6 columns, as shown below.where is the coordinate vector of the center point of the Gaussian basis function of the jth neuron in the hidden layer; is the matrix with 12 rows and 5 columns; is the width of the Gaussian basis function of the jth neuron in the hidden layer.
The RBF neural network structure for the control of the WDPR is shown in Figure 3 where the input of the RBF neural network structure is 12 and the output is a column vector with 6 rows and 1 column.
Assumption 1. Let , is the ideal output of the RBF neural network, and is the ideal weight. Suppose that a very small positive number is given so that satisfies the following condition:
Assumption 2. Let the approximation error ; it is assumed that the approximation error is bounded; that is,When is bounded, the unknown interference term and the approximation compensation term of the RBF neural network tend to be consistent and when both are equal, the approximation error is zero. At that moment, the control system operates under optimal conditions.
4. Control Design and Stability Analysis
4.1. Control Design
The schematic diagram of the control system of the WDPR is shown in Figure 4.
Considering the experimental environment, the external disturbances, and other factors, the adaptive PD control based on the RBF neural network is as shown in where , is the PD control, is an adaptive control based on the RBF neural network, and is the estimated value of the ideal output of the RBF neural network.
The proposed control defined by (4) is used and the error closedloop system equation (16) is obtained:where , , is the subtraction of the theoretical pose and the actual pose , respectively, and is the subtraction of the theoretical pose velocity and the actual pose velocity , respectively. Let , then (16) can be transformed into where and .
Because of , let ; (17) is transformed into the following form:where and is bounded. is the estimated weight.
4.2. Stability Analysis
(1) The stability of the abovedescribed closedloop system is analyzed. The Lyapunov function iswhere is the symmetric positive definite matrix that satisfies; therefore, the function is a positive definite matrix.
(2) The derivative of the Lyapunov function can be obtained:where is used in (21) and we obtain
Equation (22) is simplified and we obtain
Because , , (24) is simplified as
(3) Adaptive control technology is usually used to deal with dynamic uncertainty [21–23]; therefore, in order to respond to the dynamic characteristics of the WDPR motion process and disturbances and ensure the asymptotic stability of the system, the adaptive control algorithm is designed as follows [24]:
Equation (27) is substituted into (25) and we obtain
To ensure that the closedloop system is stable, first of all, needs to be proven. According to the known conditions, , , , can be obtained:
In (29), is the minimum eigenvalue of the matrix and is the maximum eigenvalue of the matrix .
In order to ensure that , the conditions of (30), (31), and (32) must be satisfied:
That is, it is possible to improve the convergence effect by either increasing the eigenvalue of the matrix , decreasing the eigenvalue of the matrix , or decreasing . When , the closedloop system is stable.
5. Simulation Experiments
In order to verify whether the proposed control method satisfies the requirements for practical applications, it is combined with the dynamic equation of the WDPR and a MATLAB/Simulink simulation is conducted. The 6DOF WDPR is used as the simulation object. Based on the theory from [24], better convergence results can be obtained when and have the following values:
Due to the experimental environment, the external disturbances, and other factors, the external disturbances of the system are set to [24]:
In (31), , , , .
The simulation program of the RBF neural network is established, as shown in Figure 5; the parallel robot has six DOF and the theoretical input value is set as follows:
Due to the nature of the task and the required results, the relative error is used, because it reflects the measurement errors appropriately. The range of the relative error of the pose and pose velocity is from 5% to 5%; the range of the relative error of the pose compensation is from 5% to 5%.
In Figure 1, the position of the connection point P_{i} in the aircraft model (the end effector) is in a moving coordinate system and the position of the hinge point B_{i} in the pulley is in a static coordinate system; the coordinates are shown in Table 1. In the moving coordinate system in Figure 1, P_{1} (P_{6}) and P_{2} (P_{7}), P_{3} and P_{4}, and P_{5} and P_{8} are symmetrical on the xaxis.

The simulation results of the proposed control method for the WDPR are shown in Figures 6–18.
(a) Control input of motors 14
(b) Control input of motors 58
(a) Compensation of X, Y, and Z
(b) Compensation of roll, pitch, and yaw
(a) WDPR8 prototype platform
(b) IMAC controller
It can be seen from Figures 6 and 7 that the tracking curves of the theoretical pose and the actual pose are well aligned. In Figure 8, the relative error results show that the designed control is correct and meets the experimental requirements. The reason for the errors in the simulation results may be due to some unknown factors.
Figures 9–11 show the results of the motor control input, the pose compensation, and the relative error, respectively. It can be seen from Figure 9 that the initial control input is slightly larger, which is probably due to the initial start of the motor; however, it does not exceed the rated torque of the motor (640 N.mm) and the control input of the 8 motors tends to stabilize after 2.8 s. In Figure 10, it can be seen that the pose uncertainty () and the approximation compensation () tend to be consistent after 0.02s. The relative error between the pose uncertainty and the approximation compensation is in the range of 0.025 to 0.025 within 10s (Figure 11). These results indicate that the compensation is effective. Figure 12 shows that the tension changes for the 8 wires are relatively small, which minimizes the risk of broken wires. It can be seen from Figure 13 that the RBF neural network has a good approximation performance. The simulation results are correct and represent the actual conditions.
Therefore, the simulation results show that the RBF neural network effectively approximates the error function and the tracking error of the WDPR is greatly improved; in addition, the proposed RBF neural network adaptive PD control is feasible and reliable and meets the basic requirements of the WDPR.
6. Prototype Experiments
In order to verify the correctness and effectiveness of the designed neural network adaptive PD control method, an experimental verification is performed using the WDPR prototype platform (Figure 14(a)).
The experimental platform consists of software and hardware. The software is programmed using C++ in Visual Studio. The hardware includes the PC, servo system, motion controller, mechanical structure, and other components. The PC serves as the main controller; the type of servo system is SGM7J04AFC6S with a power of 400 W. The mechanical structure consists of the frame, the ball screw, the end effector, the wiredriven system, etc. The frame is constructed of aluminum profiles and its specific size is shown in Figure 1. An IMACHX integrated motion controller (Turbo PMAC) is used. The output of the controller is an analog signal and the DC regulated power supply of DRP24024 is used to supply the controller, as shown in Figure 14(b).
In order to carry out the neural network adaptive PD control experiment, the servo motor is set to torque mode and the hardware and software environment are configured and connected. The dynamic control software is designed to complete the control experiment of the WDPR prototype platform. The control software includes the main control program, trajectory planning, inverse kinematics model, dynamic model, and other modules; a screenshot of the control software is shown in Figure 15.
Since the dynamic control software involves many matrix operations for modeling and control, the Eigen library is applied to handle the operation of the matrix vector. The dynamic control software invokes Microsoft Foundation Classes (MFC) to design the visual interface of the program. The software interface of the WDPR prototype experimental platform designed by using the class library MFC is shown in Figure 16.
The Jacobi matrix , the pseudoinverse of the transpose of the Jacobi matrix, and the neural network compensation operation involve many matrix operations and the calculation speed is too slow for realtime online operation. Therefore, the pseudoinverse of the transpose of the Jacobi matrix is calculated offline by MATLAB simulation and the compensation of the neural network is sampled and invoked in real time to improve the control efficiency of the WDPR prototype platform.
Because the WDPR prototype is a complex nonlinear system, the feasibility and control performance of the designed neural network adaptive PD control are verified by using a oneDOF motion and the desired trajectory of the end effector is set as .
The experimental results of the WDPR prototype are shown in Figures 17–23.
The experimental results (Figures 17–23) show that when the proposed control method based on the RBF neural network compensation is used for the dynamic tracking control of the WDPR, the actual pose of the end effector is smooth and continuous, the absolute error of the pose is small (±3mm), and the relative error is also small (±1%). We use the #1 motor as an example; the results show that the trend of the actual angular velocity of the motor is correct, the absolute error is ±0.1 rad/s, and the relative error is ±1.5%; the actual output torque of the motor (±0.08N.m) is also within the range of the rated torque of the #1 motor. Therefore, it can be concluded that the proposed RBF neural network adaptive PD control method ensures that the end effector of the WDPR operates on the basis of the desired trajectory. The errors meet the experimental requirements and the tracking performance is good.
7. Conclusions
In this study, a new tracking control method is proposed for a WDPR and prototype experiments are conducted. The proposed control method is a torque control scheme based on the known dynamic model of the WDPR. The adaptive PD control method is based on RBF neural network compensation and the trajectory tracking performance of the WDPR is improved by an approximation of the unknown parameters. The approach is based on the Lyapunov stability theory and an adaptive control method is used to improve the convergence of the RBF input; the stability of the closedloop system of the WDPR is demonstrated. The simulation results and the experimental results show that it is feasible to use an RBF neural network for the WDPR and fully proves the validity of the proposed control method. The research results provide guidance for the application of RBF neural networks to control a WDPR prototype for wind tunnel tests. In order to ensure the applicability of the WDPR to engineering practice, in future research, we plan to upgrade the WDPR system and enable the WDPR to achieve complex multidegreeoffreedom movement to meet the requirements of wind tunnel test indicators.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (Grants nos. 11072207, 11472234, and 11702232) and the Fundamental Research Funds for the Central Universities (20720180071).
References
 A. K. Salima, G. S. Dwarakish, K. V. Liju, T. Justin, D. Gayathri, and R. Rajeesh, “Weekly prediction of tides using neural networks,” in proceedings of the 8th International Conference on Asian and Pacific Coasts (APAC 2015), vol. 116, pp. 678–682, 2015. View at: Google Scholar
 T. Fukuda and T. Shibata, “Theory and applications of neural networks for industrial control systems,” IEEE Transactions on Industrial Electronics, vol. 39, no. 6, pp. 472–489, 1992. View at: Publisher Site  Google Scholar
 Z. Defeng, Neural Network Application Design of MATLAB, China Machine Press, Beijing, China, 2nd edition, 2011.
 H. T. Dinh, R. Kamalapurkar, S. Bhasin, and W. E. Dixon, “Dynamic neural networkbased robust observers for uncertain nonlinear systems,” Neural Networks, vol. 60, pp. 44–52, 2014. View at: Publisher Site  Google Scholar
 Y. Wang, L. Cheng, Z.G. Hou, J. Yu, and M. Tan, “Optimal formation of multirobot systems based on a recurrent neural network,” IEEE Transactions on Neural Networks and Learning Systems, vol. 27, no. 2, pp. 322–333, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 A. Melingui, O. Lakhal, B. Daachi, J. B. Mbede, and R. Merzouki, “Adaptive neural network control of a compact bionic handling arm,” IEEE/ASME Transactions on Mechatronics, vol. 20, no. 6, pp. 2862–2875, 2015. View at: Publisher Site  Google Scholar
 K. Sabahi, S. Ghaemi, J. Liu, and M. A. Badamchizadeh, “Indirect predictive type2 fuzzy neural network controller for a class of nonlinear input  delay systems,” ISA Transactions®, vol. 71, pp. 185–195, 2017. View at: Publisher Site  Google Scholar
 S. Wen, J. Yuan, and J. Zhu, “Radial basis functional link network and hamilton jacobi issacs for force/position control in robotic manipulation,” Mathematical Problems in Engineering, vol. 2012, Article ID 568796, 10 pages, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 X. Zhao, X. Wang, G. Zong, and H. Li, “Fuzzyapproximationbased adaptive outputfeedback control for uncertain nonsmooth nonlinear systems,” IEEE Transactions on Fuzzy Systems, vol. 26, no. 6, pp. 3847–3859, 2018. View at: Google Scholar
 X. Zhao, H. Yang, W. Xia, and X. Wang, “Adaptive fuzzy hierarchical slidingmode control for a class of MIMO nonlinear timedelay systems with input saturation,” IEEE Transactions on Fuzzy Systems, vol. 25, no. 5, pp. 1062–1077, 2017. View at: Publisher Site  Google Scholar
 H. Wang, H. R. Karimi, P. X. Liu, and H. Yang, “Adaptive neural control of nonlinear systems with unknown control directions and input deadzone,” IEEE Transactions on Systems, Man, and Cybernetics: Systems, vol. 48, no. 11, pp. 1897–1907, 2018. View at: Publisher Site  Google Scholar
 H. Wang, P. X. Liu, and B. Niu, “Robust fuzzy adaptive tracking control for nonaffine stochastic nonlinear switching systems,” IEEE Transactions on Cybernetics, vol. 48, no. 8, pp. 2462–2471, 2018. View at: Publisher Site  Google Scholar
 Y. SuiLu, Key Theoretical Issue Analysis and Applied Research on Wind Tunnel Tests of WDPSS8, [Ph.D. thesis], Xiamen University, China, 2013.
 “Summary of RBF Neural Networks,” Baidu.com, 2012, https://wenku.baidu.com/view/fdfe3e3e0b4c2e3f572763ec.html. View at: Google Scholar
 S. Behzadipour and A. Khajepour, “Design of reduced DOF parallel cablebased robots,” Mechanism and Machine Theory, vol. 39, no. 10, pp. 1051–1065, 2004. View at: Publisher Site  Google Scholar
 C. B. Pham, S. H. Yeo, G. Yang, M. S. Kurbanhusen, and I.M. Chen, “Forceclosure workspace analysis of cabledriven parallel mechanisms,” Mechanism and Machine Theory, vol. 41, no. 1, pp. 53–69, 2006. View at: Publisher Site  Google Scholar
 J. B. Butcher, D. Verstraeten, B. Schrauwen, C. R. Day, and P. W. Haycock, “Reservoir computing and extreme learning machines for nonlinear timeseries data analysis,” Neural Networks, vol. 38, pp. 76–89, 2013. View at: Publisher Site  Google Scholar
 M. Chowdhury, A. Alouani, and F. Hossain, “How much does inclusion of nonlinearity and multipoint pattern recognition improve the spatial mapping of complex patterns of groundwater contamination?” Nonlinear Processes in Geophysics, vol. 16, no. 2, pp. 313–317, 2009. View at: Publisher Site  Google Scholar
 W. Si, X. Dong, and F. Yang, “Adaptive neural control for stochastic purefeedback nonlinear timedelay systems with output constraint and asymmetric input saturation,” IET Control Theory & Applications, vol. 11, no. 14, pp. 2288–2298, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 M. Bhuyan, K. K. Sarma, and N. E. Mastorakis, “Nonlinear mobile link adaptation using modified FLNN and channel sounder arrangement,” IEEE Access, vol. 5, pp. 10390–10402, 2017. View at: Publisher Site  Google Scholar
 W. Huanqing, S. Peng, L. Hongyi, and Q. Zhou, “Adaptive neural tracking control for a class of nonlinear systems with dynamic uncertainties,” IEEE Transactions on Cybernetics, vol. 47, no. 10, pp. 3075–3087, 2017. View at: Google Scholar
 B. S. Park, S. J. Yoo, J. B. Park, and Y. H. Choi, “Adaptive neural sliding mode control of nonholonomic wheeled mobile robots with model uncertainty,” IEEE Transactions on Control Systems Technology, vol. 17, no. 1, pp. 207–214, 2009. View at: Publisher Site  Google Scholar
 M. Ghaemi, S. K. HosseiniSani, and M. H. Khooban, “Direct adaptive general type2 fuzzy control for a class of uncertain nonlinear systems,” IET Science, Measurement & Technology, vol. 8, no. 6, pp. 518–527, 2014. View at: Publisher Site  Google Scholar
 L. Jinkun, RBF Neural Network Control for Mechanical Systems Design, Analysis and MATLAB Simulation, Tsinghua University Press, Beijing, China, 2014. View at: MathSciNet
Copyright
Copyright © 2019 Yuqi Wang 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.