Kinematic of 4SPRR-SPR parallel robots based on VQTAM

4SPRR-SPR parallel robot is a novel closed-loop mechanism. The research of kinematics is the basis of real-time and robust robot control. This paper aims to proposing a method to address a surrogate model of forward kinematics for PMs (Parallel Mechanisms). Herein, the forward kinematics model is derived by training the VQTAM (Vector-Quantified Temporal Associative Memory) network, which originates from a SOM (Self-Organized Mapping). During the processes of training, testing and estimating this neural network, the priority K-means tree search algorithm is utilized, thus improving the training efficacy. Furthermore, LLR (Local Linear Regression), LWR (Local Weighted Linear Regression) and LLE (Local Linear Embedding) algorithms are respectively combined with VQTAM to get three improvement algorithms, aiming to further optimize the prediction accuracy of the networks. To speed up solving the least squared equation in the three algorithms, SVD (Singular Value Decomposition) is introduced. Finally, Data from inverse kinematics by geometric method is obtained, which is for constructing and validating the VQTAM neural network. Results show that the prediction effect of LLE algorithm is better than others, which could be a potential surrogate model to estimate the output of forward kinematics.


Introduction
Parallel robots offer noticeable advantages in machining applications, especially in the field of heavy equipment manufacture, due to their high rigid-weight ratio and load capacity. Recently, the machine tool industry has discovered the potential advantages of PMs and many parallel machine tools have been developed based on either the 6-DOF (degrees-of-freedom) parallel mechanisms or 5-DOF PMs. 1 Parallel robots or PMs with 5 or 6 DOF are often used as motion simulators and crane devices. Stewart introduced a 6-DOF PM, which is now popularly known as Stewart platform. 2 It is commonly used in flight or driving simulators, vibration isolation platform, space docking mechanisms, and space telescope. 3,4 Hereafter the lower-mobility (2-5 DOFs) PMs have emerged and have become popular owing to their simple structure, low kinematic coupling, low cost and convenient control. [5][6][7] Xie et al., 5 Huang et al. 6 and Joshi et al. 7 studied the type synthesis for different types of lower-mobility PM. Therefore the 2UPR/SPR PM 8 and 3UPS/UP PM 9,10 have been used as position adjusters and space telescope (R, P, U and S stand for revolute, prismatic, universal and spherical joints respectively). As far as 5-DOF PMs are concerned, researchers have mainly worked on the type synthesis. 6,[11][12][13][14][15] Through the type synthesis, a variety of PMs are obtained. 16,17 For instance, M. Garcia-Murillo et al. 18 proposed a 2(3-RRPS) parallel robot and analyzed its workspace. Nicholas Baron, et al. 19 displayed a family of kinematically redundant planar parallel robots, and presented a geometric method of singularity avoidance. In addition, the cable-driven manipulators are a class of PMs, where the end-effector is displaced by means of servo-actuated extendable cables. 20 They have high ratio of payload to robot weight and a much larger workspace which is appropriate for camera positioning in sports stadiums and many others. 21 Arian Bahrami et al. studied the optimal design of various forms of the cable-driven manipulators. 22,23 In summary, the structures of PMs are various and have been applied to many fields. However, there is still no executable method to solve the kinematics, especially the forward kinematics.
The 5SPRR PM is a 6-DOF PM, which has been used as a machining tool. [24][25][26] The structural form of 4SPRR-SPR parallel robot, which has great potential to transform into a new type of mobile machining robot, is derived from the 5SPRR PM. To the best knowledge of the author, up to now, Research, such as its kinematic, workspace analysis et al. has not been conducted yet.
Parallel manipulators do not have formulaic kinematic solutions. Therefore, several studies are focused on methods to simplify the solution of robot kinematics. Intelligent algorithm or machine learning is introduced into kinematics solution of robot. Simpler mapping models of kinematic input and output are established to replace the complex analytical solution. Evolutionary computational methods such as genetic algorithm (GA), 27-31 differential evolution (DE), and particle swarm optimization (PSO) 32 are used to optimize the solutions by formulating the problem in the form of an objective or a fitness function. Furthermore, approximation methods based on artificial neural networks and polynomial regression models [33][34][35][36][37] can also be used to estimate forward kinematic. Robot kinematics can be regarded as a non-linear autoregressive model with exogenous input (NARX). In reference, 38 the forward kinematics of parallel manipulators is solved by combining the wavelet network and the NARX method. However, the validity of this method is still dependent on parameter initialization and back propagation learning process. The simplification of complex kinematics model often leads to decreasing the accuracy. To overcome this problem, efficient and accurate kinematics identification is needed. 39 The VQTAM algorithm does not depend on the adjustment of super parameters, but on SOM to achieve topology preservation. It can be applied to the mapping from input space to output space, as an extension of SOM. The mathematical foundation of VQTAM is non-linear manifold embedding. This method is suitable for time series modeling and prediction of non-linear systems. Depending on the selforganizing mapping structure, the discrete mapping relationship can effectively realize output estimation according to exogenous input. [39][40][41] Since the input space and output space of VQTAM are discrete, a network with more neurons is required to achieve highprecision prediction. In order to improve the computational efficiency, it is necessary to reduce the number of VQTAM neurons, which decreases the efficiency. Improved local linear algorithms of VQTAM are proposed to solve the inverse kinematics of a 6-DOF robot efficiently and accurately.
In this paper, a novel parallel robot that could be used as mobile machining robot is introduced. The architecture and general kinematic properties of the 4SPRR-SPR parallel robots are outlined, and the study on its kinematic follows.

Research methodology
4SPRR-SPR parallel robot is a novel closed-loop mechanism. Its inverse kinematics can be achieved using a geometric method. The rod length can be obtained quickly through the known pose of end effector. But the forward kinematics of parallel robot is very complex. It's a non-linear system. Generally, it is difficult to accurately determine the forward kinematics.
Thus, the length of drive limbs corresponding to different poses can be obtained using inverse kinematics conveniently. The data obtained from the inverse kinematics can be used to train the surrogate model of the forward kinematics. In this paper, VQTAM neural network is used to replace the complex calculation of forward kinematics. The research content and structure of this paper are shown in Figure 1.
In Section 2, the inverse kinematics is studied. Using the inverse kinematics algorithm, the length of drive limbs corresponding to poses is calculated as the sample set.
In Section 3, the non-linear system identification method based on VQTAM is studied, and the training method of VQTAM neural network is established. In the process of VQTAM training and testing, it is core to search and activate neuron nodes. Therefore, the priority search K-means tree algorithm is introduced for this step.
In Section 4, the improved method of local linearization of input and output mapping based on VQTAM neural network is studied.

4SPRR-SPR parallel robots
The 4SPRR-SPR parallel robot is one of the new 5-DOF PMs, and its structure is shown in Figure 2. The four limbs of the PM consist of one spherical joint, one moving pair and two rotating pairs, and the other one lack for a rotating pair at one end of the limb. The rotating pairs connected with an end effector are coaxial. The spherical joints connected with the static platform are not arranged in the same plane but are distributed in different positions in space. This provides more possibilities for the optimization of the parallel mechanism.
The arrangements of the kinematic pairs of the 4SPRR-SPR parallel robot are similar, except that the end of the first limb is fixed to the moving platform. The kinematic pair arrangement of the i th chain is shown in Figure 3 (i = 2,3,...5).
A i is the center point of the universal joint connected by the ith limb and the static platform, B i and C i are the two revolving pair centers of the limb, and D is the reference point of the position of the end effector. The screw of the center line S of each pair is used to represent the motion of the pair. 12 Where b i , c i , m i , and n i represent the length of segment The six screws of the i th limb are linearly independent (i = 2,3,...5). The screws of the i th limb constitute a 6-system of screws, which provide no constraint to the end effector.
The screw system of the first limb differs from those of others, because the end of the first limb is fixed to the moving platform. The screws represented in coordinated O 1 -x 1 y 1 z 1 are The reciprocal screw of the 5-system screws can be written as Therefore, five limbs provide one constraint to the end effector. It can be confirmed that 4SPRR-SPR parallel robots have 5 DOFs. The motion screw system which reflects the DOF property of the end-effector can be obtained by finding the secondary inverse screw system of the end-effector.
Formula (4) shows that the end effector can rotate around the X and Y axes and around the line (0,0,1;0,n 1 ,0), which passes through point A 1 and parallels to the Z axis. The end effector can also moving along the X and Z-axes.

Inverse kinematics of 4SPRR-SPR parallel robot
The prismatic joints of the five limbs constitute the drive of a 4SPRR-SPR parallel robot. Therefore, the inverse position solution of the parallel robot can be regarded as the calculation of the rods A i B i about the length l i using the coordinates of the known points D and direction vectors S 1 (i = 1,2,.,5). The directions of end effector are the same for the five limbs in unified coordinates.
The geometric condition of the parallel robot is determined using equations (1) and (2).
As shown in Figure 3, the i th limb is considered for the analysis. E i is the vertical foot from A i to DC i . From the geometric relation, the result is deduced using equations (6) to (8).
The coordinates of points C i and B i can be obtained upon calculating c i . Accordingly the input parameter l i corresponding to the pose can be obtained.
The main structural parameters of 4SPRR-SPR parallel robots are the coordinates of points A i (i = 1,2,...5) and rod lengths c i and d i . The above parameters can determine the corresponding parallel robot for inverse kinematics calculation.
The inverse kinematics solution of the robot is calculated considering the structural parameters shown in Table 1 as an example.
By determining the position and posture of the moving platform, the coordinate of point D and the direction vector S 1 is found to be (0,1119.17,3641.27), and (0,0,-1). The length l i of the rod A i B i , as seen in Table  2, is calculated by substituting these parameters in the aforementioned algorithm. System identification of non-linear dynamical system using VQTAM The forward kinematics of the parallel robot aims at solving the corresponding pose given the length of drive limbs. The pose of the end effector is regarded as output Q and the length of the drive limbs as input p. The forward kinematics problem of the robot can be transformed into the identification problem of a non-linear dynamic system. The time-discrete difference formula is established as shown in Formula (9). 39,41 Where Q(t) =[x(t),y(t),z(t);S 1 (t)], which is the pose of the end effector; p(t) = [l 1 (t),l 2 (t),l 3 (t),l 4 (t),l 5 (t)], which is the length of drive limbs at t time. f(Á) represents a non-linear function reflecting the characteristics of the system. The output Q at t + 1 is determined by the previous n q outputs and n p inputs.
The purpose of forward kinematics is to recognize the mapping relation f(Á) from input to output. For a non-linear problem without formulaic solutions, simplifying the model can guarantee efficiency and ensure favorable real-time performance in engineering applications.
Q(t + 1) = f ½Q(t), :::, Q(t À n q + 1); p(t), :::, p(t À n u + 1) ð9Þ X in = ½Q(t), :::, Q(t À n q + 1); p(t), :::, p(t À n u + 1) VQTAM is an extension of SOM, which can be used as an effective manner to realize non-linear mapping. The concept of non-linear manifold embedding can be used to establish the non-linear mapping. SOM is realized through competition and cooperation among feature neurons. Different neurons are activated for different inputs, and the activated neurons affect the output together with the neurons in their neighborhood. Therefore, time series data can be constructed as the training samples of VQTAM.
VQTAM uses the training samples to build a topological structure embedded in multi-dimensional data space, which can construct the mapping relationship between multi-dimensional input vectors and multidimensional output vectors. The VQTAM network consists of three layers: input space v in , output space v out and lattice space. The establishment of the mapping process is shown in Figure 4. The definitions of X in and X out are shown in equations (10) and (11).
The input space and output space are composed of weight vectors v i in and v i out respectively. i is the index of neuron in lattice space, which reflects the location of neuron in lattice topology structure. The weight vectors v i in and v i out have mapping relationships with neurons in lattice space. v i in and v i out have the same dimensions as X in and X out respectively. In the mapping process of VQTAM, the nearest weight vector v i* in to X in is found in the input space, which corresponds to the activated neuron in the lattice space. In this process, The Euclidean distance is used to measure the distance between two vectors. The definition of i* is shown in Formula (12), where A is the collection of all neuron indexes in lattice space. 41 v i* out is obtained according to the neuron index i* and the output X out is estimated, as shown in Formula (13) Learning strategy of VQTAM The main parameters of VQTAM are neuron weight vector v i in , v i out and lattice topology. The initial lattice topology structure can be determined as a super parameter before training, so the learning strategy mainly aims at updating the weight vectors v i in and v i out in the training process.
There is competition and cooperation among neurons in SOM. Different neuron and their neighborhoods are activated by input to participate in the learning process, and the weight vectors v i in and v i out are updated. The index search of the active neurons in the lattice space is consistent with equation (12).
In order to improve the convergence rate in the learning process, the influence range parameters s(t) of the activation neuron and the learning speed a(t) decrease exponentially with the learning epoch, as shown in equations (14) and (15).
Where t is the current learning epoch, and T is the total learning epoch. a 0 and s 0 are the initial value, and a T and s T are the parameter value of training epoch T.
The Gauss neighborhood function is used to determine the effect of input on the neighbors of activated neuron, as shown in equation (16). [39][40][41] v i in and v i out are updated, as shown in equations (17) and (18).
To sum up, the steps involved in the VQTAM algorithm are as follows: VQTAM algorithm: (training part) Begin 1. Input: X in , X out in training set 2. for j 0 to length(X in ) À 1

Searching activated neuron by priority search Kmeans tree algorithm
The weight vector v i* in is the nearest neuron to X in in the input space, and is searched during the training, testing and estimation of VQTAM, as shown in equation (12). The search algorithm applied has a great impact on the efficiency of VQTAM. The global traversal method needs to calculate all the Euclidean distance from the weight vectors v i in to the input X in , and searches the minimum value through traversal, which is very inefficient. The K-Dtree algorithm can reduce the time complexity of searching, but for d-dimensional data, the time complexity of K-Dtree data structure search is O(n d/(d-1) ). For high-dimensional data search, especially when the data dimension d.20, K-Dtree data structure cannot provide high efficiency. v i in is 30 dimension. The searching efficiency is severely restricted. The search time complexity of the priority search k-means tree algorithm is O(Ld logn/logK), where n is the amount of data in the data set, and K is the number of nearest neighbors to be searched, and L is the maximum number of data retrieved in the search process. The priority search k-means tree algorithm is suitable for searching high-dimensional data.
The priority search K-means tree algorithm is to divide the space into B different regions using k-means clustering algorithm. Then, the points in each region are partitioned by the same operation until the number of data points in the region is not more than K. In the process of searching, the idea of ''divide and conquer'' is used to find the nearest point in a smaller area and reduce the amount of data to be searched, so as to improve the efficiency.
The algorithm of establishing a priority search Kmeans tree data structure is as follows: K-means tree data structure building algorithm: 42  Create leaf nodes from data sets 5. else 6.
P uses C alg algorithm to select B points from data set D 7. for j 0 to length(P) 8. do 9.
C clustering data in D centered on P j 10.
P j new finding the mean value of group C data after clustering 11.
while (P = P j new & P j new is the non-leaf node) 12. Output: the entire K-means tree data structure In the establishment of K-means tree, the algorithm C alg used to select B points from data set D can be chosen from random selection algorithm, Gonzale algorithm and K-Means ++ algorithm. The specific algorithm has little effect on data structure establishment. Generally, random algorithm is recommended. 42 In the built K-means tree, leaf nodes are the data points in the original data set, and non-leaf nodes are the central points of the regions after segmentation. The search process can only be executed in a few areas, thus avoiding the global search of data sets. The Kmeans tree data structure is searched from the root node. Sub-nodes are sorted according to the distance between the clustering center and the query data points. The nearest sub-nodes are searched in advance. The priority K-means tree search algorithm is as follows: Priority Search K-means Tree Algorithm 42 : Begin 1. Input: K-means tree, Q, K, L 2. /*Q is the query data X in ; K is the nearest neighbor number; L is the maximum number of search data */ 3. S the root node /*Initialize a stack S and place the root node in the stack*/ 4. do 5. if S(0) is a leaf node 6. P S(0) /*Add S(0) to the retrieved data array P */ 7. else 8.
remove S(0), E sub-nodes /* The top node goes out of the stack*/ 9.

VQTAM local linear improvement algorithms
The VQTAM algorithm is used to quantify the input space v in and the output space v out . The neurons in v in and v out correspond through mapping relations. The input X in is approximated to the nearest neuron v i* in during estimation. The accuracy of the estimation results can be guaranteed when the number of neurons is adequate. However, an increase in the number of neurons increases network size and the decreases the computing efficiency. Thus, local linearization of the activated node is used, which can balance the number of neurons and prediction accuracy. Three improved algorithms for VQTAM are proposed: LLR, LWR and LLE based on local linearization. The inverse kinematics of the robot can be further optimized and ensure the computational efficiency.
Local linearization is performed in VQTAM networks. The K-means tree algorithm is used to search the nearest n data points v i*n in in the input space and the corresponding output data v i*n out are mapped, as shown in Figure 5.

Improvement algorithm of VQTAM with LLR
It's assumed that the relationship between local input and output can be expressed linearly, so the LLR can be performed as shown in equation (19).
Y is the nearest neighborhood v i*n in of X in , and Z is the corresponding data point v i*n out , As shown in equations (20) and (21).
The linear relationship between the local input and output can be obtained by solving W. Equation (21) is an over-determined equation, which is generally treated as a least squares problem and solved by a generalized inverse matrix, as shown in equation (22).
Y T Y is an n-dimensional square matrix. It needs great amount of calculation when the dimension increases. Meanwhile, Y T Y may be a singular matrix or ill-conditioned matrix, so it is difficult to find the inverse matrix directly. The SVD (singular value decomposition) can be used to solve the least squares problem. Y can be decomposed as shown in equation (23).
The column vectors of U are the left singular vectors of Y, that is, the eigenvectors of YY T ; the column vectors of V T are the right singular vectors of Y, that is, the eigenvectors of Y T Y.
Where S is a diagonal matrix whose value is the singular value y of matrix Y, that is, the eigenvalue of YY T is l = y 2 , and U and V T are orthogonal matrices.
U can be disassembled as [U n ,U m-n ]. For the least squares problem as shown in equation (21), the solution obtained by SVD is shown in Formula (25).
The input X in is used to estimate the output, as shown in Formula (26).

Improvement algorithm of VQTAM with LWR
The impact of all neighbor points on the prediction results is consistent in LLR algorithm. Considering the different influence of neighbor points, a distance weight is added in the regression process. A local weighted linear regression (LWR) is performed. The cost function is defined as shown in equation (27). The weight w k in equation (27) is determined using the Gauss neighborhood function, as shown in equation (28), where t is the bandwidth parameter.
The analytical solution of the cost function minimization is equivalent to the solution of the overdetermined equation shown in equation (29).
W is an n-dimensional diagonal matrix and the diagonal element is w k .
The SVD for matrix WY is also used to solve the over-determined equation.

Improvement algorithm of VQTAM with LLE
It's assumed that input data can be represented by a linear combination of several samples in its neighborhood. That is to say, the input X in to be predicted is represented by a linear combination of n data points in its neighborhood in the input space v i*n in .
The output has the same linear combination.
c k is the coefficients of the linear combination. The output can be estimated using equation (33) after solving c k . Define the cost function as shown in equation (34) and rewrite it as a matrix.
For the least squares problem as shown in equation (33), it can also be solved by SVD of Y.
Simulation results and discussion

Standard VQTAM network test results
The poses for the sample set are generated according to the robot parameters listed in Table 1, using equations (25) to (30). The generated data can effectively cover the entire workspace of the manipulator.
x(t) = 1500(1 À e Àpt ) cos 1:88pt ð36Þ y(t) = 1500(1 À e Àpt ) sin 1:88pt + 1500 ð37Þ u 2 (t) = 2p(1 À e Àpt ) sin 0:86pt ð40Þ S 1 = cos u 1 sin u 2 sin u 1 sin u 2 cos u 2 ½ ð 41Þ where t2[0,1000] and step t = 0.01. Thus Q(t) = [x(t),y(t),z(t);S 1 (t)] is calculated as the output of the sample set. The three components of S 1 are s x , s y and s z , respectively. p(t) = [l 1 (t),l 2 (t),l 3 (t),l 4 (t),l 5 (t)] is calculated using the inverse kinematic algorithm. p(t) and Q(t) are used as the sample set. Thus the VQTAM network is trained. The data is selected randomly from the sample set in each epoch. At the same time, in order to test the robustness of the algorithm, random errors are added to the test set data according to N(0,1) distribution. Table 3 shows the hyper-parameter setting of the VQTAM network training. A VQTAM network of dimension 80 3 80 (M x 3 M y ) is trained based on the hyper-parameters, and the effect of training is evaluated using the test set.
Root mean squared error (RMSE), R squared (R 2 ), and root mean absolute error (RMAE) are used to evaluate the prediction accuracy of VQTAM for forward kinematics. Table 4 shows the prediction accuracy of the six parameters in output calculated using RMSE, R 2 , and RMAE.
The value of R 2 is between 0 and 1. The closer it is to 1, the better is the effect of the regression fitting, which means that the learning VQTAM network is more accurate as an approximate model. Generally, R 2 . 0.95 can be used in engineering applications. For VQTAM networks of dimension 80 3 80, the R 2 values after learning are all above 0.95. This shows that this network is suitable for engineering applications. The RMSE and RMAE reflect the prediction errors in the VQTAM network, and their small values also indicate the prediction accuracy of the network. The convergence result in neural network training process is shown in Figure 6.
It can be seen from Figure 6 that the R 2 of the VQTAM neural network tested using the test set in the training process gradually increases and converges near 1. This confirms the effectiveness of the VQTAM neural network training algorithm.
The consistency between the actual data and the estimated data can be displayed more intuitively through box diagrams, as shown in Figure 7. The box diagrams of the actual values and predicted values are notably similar, which verifies the excellent prediction ability of the VQTAM network of dimension 80 3 80.

VQTAM local linear improvement algorithms
The output of the inverse kinematics problem is estimated using the VQTAM local linear improvement algorithms. The priority search K-means tree search algorithm is used to search the neighbor data of the input. The influence of parameter k on the prediction accuracy is analyzed, considering a VQTAM network with dimension 70 3 70 as an example, as shown in Figure 8. Figure 8 shows that R 2 of each parameter increases continuously with the increase in k, at first; and it gradually increases to more than 0.95. This indicates the influence on the prediction accuracy. When k = 7, the overall prediction accuracy is the best.   Figure 9 shows that R 2 for all variables is more than 0.95 when k = 3. Comparing Figure 8 with Figure 9, the overall prediction accuracy of LWR algorithm is higher than that of LLR algorithm. Comparing among Figures 8 to 10, the overall prediction accuracy of improvement algorithm of VQTAM with LLE is the highest. When k = 7, the overall prediction accuracy is the best. In conclusion, when the number of neurons is

Overall examples and test results
According to the VQTAM network parameters shown in Table 3, the prediction accuracy of four algorithms under different network dimensions (30 3 30, 40 3 40, 50 3 50, 60 3 60, 70 3 70, 80 3 80) is studied, as shown in Figure 11. The number of the nearest nodes of the local linear algorithm is set to k = 7. Figure 11 shows that an increase of the number of neurons results in a significant increase in the prediction accuracy of the four algorithms. The overall prediction accuracy of improvement algorithm of VQTAM with LLE is the highest. Local linearization of VQTAM yields remarkable results. Using the improved VQTAM algorithm, the prediction accuracy of a network of dimension 50 3 50 is found to be close to that of the standard VQTAM network of dimension 70 3 70. Considering improvement algorithm of VQTAM with LLE as an example, the output of the forward kinematics is estimated, as shown in Figure 12. The estimated poses are compared with the actual poses. This shows that the VQTAM neural network can effectively estimate the output The error of VQTAM can be obtained by comparing the result with the actual pose of the robot. The RSME for position estimation is less than 0.1mm, as shown in   Table 5. This precision has met the control requirements of general robots in industrial applications.

Conclusion
The kinematic of 4SPRR-SPR parallel robot is studied, and the inverse kinematics solution by a geometric method is derived. Given the input pose of the robot, this method can be used to calculate the lengths of the drive limbs conveniently and quickly. It can also be used as the computational basis of VQTAM network to establish sample sets. A fast forward kinematic mapping method for parallel robots using VQTAM for non-linear dynamic systems identification is proposed. The VQTAM algorithm is constructed to train and test the neuron network. Based on the priority K-means tree search algorithm, the training, testing and estimating process of VQTAM network are improved to enhance the search efficiency of nearest neighbors.
According to local neurons activation, the VQTAM network is improved based on local linearization. To further optimize the prediction accuracy of the network, reduce the dimension of the network in application, the improved local linear algorithms for VQTAM networks are proposed.
The sample data set for VQTAM training is constructed. The VQTAM algorithm and its improved algorithm are tested. The test results show that the increase of network dimension can improve the prediction accuracy of VQTAM, but the computational efficiency is affected. By using VQTAM local linearization improvement algorithms, the estimation accuracy can be optimized for low-dimensional network. The above algorithm can show good prediction accuracy, and the prediction effect of LLE algorithm is the best among the three algorithms.