A Prediction Model of Forest Preliminary Precision Fertilization Based on Improved GRA-PSO-BP Neural Network

+e optimal amount of fertilizer application which was needed by the trees and the factors that influence the fertilization have an intricated nonlinear relationship. According to the problems that the traditional fertilization predictionmodel has, such as lacking of the scalability and practicality, this paper initiates an accurate fertilization prediction model that was based on the GRA-PSOBP neural network which canmake the accurate fertilization come true and improve the economic benefits of forest industry.+is paper uses the GRA method to determine the input of the neural network as the site index and make the forest age, nutrient content of the advantage trees, biomass of the advantage trees, biomass of average trees, and target yield as the output numbers of the Actual amount of fertilizer applied. During the calculation process, the global particle swarm optimization algorithm is used to optimize the initial numbers and threshold numbers of BP neural network which build a phased GRA-PSO-BP accurate fertilization model. Compared with the prediction algorithm of full input variate that is based on the single BP neural network and the prediction algorithm of full input variate that is based on PSO-BP Neural Network, the GRA method can determine the key factors that influence the amount of fertilizer applied in different forest areas and modify the prediction model to improve the scalability and accuracy of the prediction and finally achieve the precision fertilization as the data of different forests updated, so we can see that the prediction result of this paper is more accurate.+e result demonstrates that the GRA-PSO-BP neural network Segment fertilization model is more accurate than the traditional BP neural network and BP Neural Network that was optimized by the PSO algorithm, and specifically, the error of the predicted amount of fertilizer application and the actual amount of fertilizer application is less than 5%, which can effectively guide the fertilization in stages.


Introduction
In the process of forestry production, it is fundamentally required to calculate the amount of fertilizer, such as nitrogen, phosphorus, and potassium fertilizer, for different forest species and different soil conditions [1]. Cunninghamia lanceolata enjoys the characteristics of fast growth, wide application, high economic value, and easy reproduction [2]. Cunninghamia lanceolata is dependent on fertilizer and humidity. In the condition of 25°C-35°C, Cunninghamia lanceolata can grow rapidly. In this regard, fertilization is the key link for high yield of Chinese fir forest. Fertilization generally applies from young forest to middleaged forest, which plays an important role in forestry economy. Forest fertilization functions as a critical technique to enhance soil fertility, improve tree nutrition, and facilitate rapid growth and high yield [3,4]. At present, forestry production is increasingly transformed to precision operations. Cunninghamia lanceolata production is no exception. e premise of precision production of Chinese fir forest is precise fertilization since such an agrotechnique can help balance the nutrients in soil and increase the nutrients required by trees, so as to achieve the target yield. In addition, the implementation of precision fertilization can save the fertilizer and phase down environmental pollution.
Fertilization prediction model, as the key of precision fertilization in Cunninghoya lanceolata forest, has become a hot spot in forestry research. In response to this issue under discussion, Hu Yueli et al. proposed a site-nutrient effect fertilization model [5][6][7][8][9]. However, sometimes this model can hardly help to get satisfactory accuracy and stability. e traditional fertilization model can only describe the static relationship between the growth of forest and the required fertilization nutrients, which can hardly reflect the different growth periods of trees. In view of the differences in nutrient absorption capacity and fertilizer application in different growth stages, this model cannot be applied to guide the stage fertilization of trees. In the past, it was difficult to collect data due to the long growth period of forest. erefore, there were few models of forest fertilization. However, in recent years, sensor network can be used to collect real-time data in forest areas with the burgeoning development of artificial intelligence. It is time to construct a predictive model of accurate fertilization of forest through the collected data. Similarly, as the traditional one, the model can only exert real-time and accurate small-scale prediction in forest areas considering the complexity of different forest areas and site indices, the diversified factors affecting forest fertilization. At the same time, since there is a high degree of nonlinearity between site index and dominant wood-related nutrient content, forest age, and other influencing factors, the traditional fertilization model often has inaccurate results due to manifold parameters or low goodness of fitting. As a result, the large amount of accumulated data can hardly serve to guide production practice. In the later stage, expert experience is still needed, resulting in a waste of manpower, material, and financial resources.
Neural network is considered to be a powerful tool to solve nonlinear problems [10][11][12][13][14][15]. In order to ensure the integrity of information and improve the precision of fertilization prediction, the prediction model of multifactor precision fertilization based on neural network can be used to solve these problems. It can reflect different absorption of nutrients in different growth periods of young and middle ages of trees and apply fertilizer efficiently in the optimal fertilization period of trees. e BP algorithm shows great potential to introduce the prediction model of forest fertilization. With the updating of the data of forest fertilization experiment, the BP algorithm modifies the existing prediction model of fertilization, which makes the forecasts of the model more precise and get closer to the actual situation. However, the accuracy of predicted results and actual values is not fairly reliable, on account of some big defects, such as local optimization, slow convergence speed, and large dependence on training data, when a single BP neural network algorithm is used to predict the amount of fertilizer applied to trees. In this paper, the analyzing method of grey relativity [16][17][18][19] is introduced to extract the original variable data of multiple indexes in advance, and the factors with high correlation degree are taken as the network input layer, and then the global searching ability of PSO [20][21][22][23][24][25] is used to optimize the BP neural network, which greatly avoids the defect of the model falling into the local optimum and further improves the accuracy. Although neural network prediction has been applied in medical engineering, electric power, system, and other fields at home and abroad, and some gratifying results have been achieved, there are only few studies on the application of neural network prediction model to forest fertilizer application. erefore, an improved GRA-PSO-BP algorithm was used to establish the accurate fertilization model of trees in stages.

Test Area Overview and Plot Settings.
e experimental data are from the state-run forest farm in Renhua County, Guangdong Province, and Liujiashan Forest Farm in Shixing County, Guangdong Province, which is the main production area of fir trees, with an altitude of less than 500 m, and the soil type is the red and yellow soil developed on slates. Largescale fertilization experiments were carried out for young and mature timber of fir forests, and the experimental forests with different site indexes were classified during the experiment according to the utilization of the last round of forest land and the Nanling Mountains fir trees in the different site conditions; the experimental forest can be divided into four types: (1) triple-tilled woodlands of Cunninghamia lanceolata, with a site index of 5; (2) cultivated Chinese fir double-tilled woodlands of Cunninghamia lanceolata, with a site index of 18; (3) first ploughing woodlands after shrub felling, with a site index of 21; and (4) first ploughing woodlands after logged broad-leaf forests, with a site index of 23. Under the condition of the same site index, the slope direction of the test site is the same, the terrain difference is small, and the trees are arranged at the same height. e young forest and the middle-aged forest test area are, respectively, set in the forest stands with appropriate site index, and the area of each test site is 1 hm 2 .
e initial planting density of young trees in the experimental forests was 3,600 plants/1 hm 2 . In the first 6 to 8 years, the first thinning was carried out, and the forests in the middle-aged forests were retained at 2,700 plants/hm 2 . e size, distribution, and density of the forests in each experimental area were basically the same. ere are fixed stakes around each test area, and the dominant wood and average wood in the test area have specific numbers. e growth and nutrient dynamics of dominant wood and average wood in the test area are measured regularly every year, and nitrogen, phosphorus, and potash fertilizers are applied according to the fertilization method in the test area.
16 plots with age of 5a, 19a, 22a, and 25a were selected in the above four different site types. e selected average and dominant trees in different areas were logged and then divided by leaves, branches, stems, barks, and roots to weigh each organ of them. e age of leaves and branches were weighed, respectively, and the root system was dug to a depth of 60 cm in the range of the canopy. Samples of all organs were evenly sampled and dried at 80°C to determine water and nutrient content. Meanwhile, four soil sampling points of different stand ages were set up in four different site index types, 16 soil sampling points in total. Take soil samples in 0∼60 cm soil layer of each soil sampling point. e physicochemical analysis of plant and soil samples was carried out in accordance with conventional methods, and the determination of the physicochemical properties of the soil in the experimental forest land is shown in Table 1. e discrepancies of site level are due to the different locations of the samples, which results in different climatic zones, slopes, slope positions, and soil types. e differentiated use of forest land at different site levels leads to significant differences in the organic matter content, nutrient content, and fertilizer utilization rate of the humus layer in the woodlands. Consequently, it can be concluded that the availability of soil nutrients is one of the essential factors for different experimental forests.
In the selection of the test plots, the age of the forest stands has been replaced by space instead of time. e accuracy of analysis is, to some extent, affected by the varied environmental conditions between the dominant wood and the average wood in the same site type. However, the effects of random errors on the basic trend of tree growth and nutrient absorption changing with forest age were relatively eliminated due to the significant difference between the growth of dominant and average stands. Table 2 shows the measured values of nitrogen (N), phosphorus (P), and potassium (K) in different organs of dominant Chinese fir trees. e nutrient content of different Chinese firs varies greatly. e order of nutrient content in the various parts of the tree is as leaf, branch, bark, root, and stem. At the same time, it was also shown that the nutrient contents of different organs of Cunninghamia lanceolata with different site indexes were different at different ages. e content of nitrogen in the leaves increased with the increase of the site index and the age of the forest; the content of potassium in the leaves, nitrogen, phosphorus, potassium in the trunk, branches, and bark and the content of the phosphorus and potassium in the roots decreased with the increase of the age. Table 3 shows the biomass range from average tree to dominant tree and the percentage of each organ in Chinese fir forest. Based on the data in Tables 2 and 3, the average nutrient content C of the Chinese fir forest obtained by the weighted average method is calculated as C � ( C i w i )/ w i , where C i and w i represent the nutrient content and biomass of leaves, branches, stems, bark, and roots of a single tree, respectively. Its product sum C i w i is the total amount of nutrients absorbed by a single tree, the total organ biomass w i is the total biomass of a single tree, and the relationship between the total nutrient absorption X per unit area of Chinese fir forest and the average nutrient content C is X � nwC/1000 � CW, where, n is the stand density, w is the average biomass per tree (kg/plant), and W is the biomass per unit area (t/hm 2 ).

Analysis of Influencing Factors.
ere are a series of important indexes to calculate the amount of fertilizer applied to trees, such as site index, forest age, nutrient concentration of dominant trees, dominant wood biomass, average wood  biomass, fertilizer utilization rate, and target yield. e site index is a natural environmental factor required for forest production and is a collection of conditions such as moisture, temperature, light intensity, and soil fertility. In the forest area with high site index or low site index, the fertilization effect is difficult to show. In the forest area with high site index, the soil fertility, water, temperature, light, and other conditions are very good. After fertilization, the effect on the increase of forest yield is small. On the contrary, the soil fertility, water, temperature, light, and other conditions are very poor in the forest land with low site index. After fertilization, due to natural factors, the fertilization effect is also very poor. In the effective site index interval, the theoretical limit value of forest production is the dominant wood biomass of forest land in practice, so the dominant wood nutrient content is the limit value that can be reached by the average wood nutrient content after fertilization. If the difference between the dominant wood biomass and the average wood biomass is greater, the fertilization effect will be better. In addition, the target yield is also related to the amount of fertilization. e target yield is determined by the input-output ratio of forest land in previous years, and the target yield determines the target increase yield, that is, the amount of fertilization. e relative intensities of the absorbed nutrients of different forest ages are different. e best fertilization period, the more fertilizer is absorbed. At the same time, the utilization rate of the fertilizer is related to the physical and chemical properties of the soil. e physical and chemical properties of the soil in different forest land are different, so the utilization rate of the fertilizer is also different. e higher the utilization rate, the better the fertilizer effect. Generally, three kinds of fixed nitrogen, phosphorus, and potassium fertilizers are applied in forest land. e utilization rate of fertilizer is observed according to the different physical and chemical properties of soil in different forest land. erefore, seven factors including site index, forest age, nutrient concentration of dominant trees, dominant wood biomass, average wood biomass, fertilizer utilization rate, and target yield were selected as the influencing factors of forest fertilizer application.

Data Normalization.
ree hundred tree growth and nutrient measurement records were obtained. e data of each group were X 1 site index, X 2 forest age, X 3 dominant wood-related nutrient concentration, X 4 dominant wood biomass, X 5 average wood biomass, X 6 target production, and X 7 fixed fertilizer utilization rate. Since the input parameters have different dimensions and the orders of magnitude of units are also greatly different, these will have a great impact on the training of neural network. erefore, the input parameters need to be normalized. In order to improve the training efficiency, an accurate fertilization model is established. e minimum maximum transformation method is used to normalize different data and deal with the dimensional influence between data. e parameters are in the range of [−1, 1], and then the parameters of neural network are trained. e following formula is adopted for normalization: Map the data to (−1, +1) and replace the formula with where y is the normalized value, x is the original value of a parameter, x max is the maximum value in the set of values, and x min is the minimum value in the set of values.

Basic BPNN Model.
e structure of the basic BPNN consists of the input layer, the hidden layer, and the output layer, as shown in Figure 1. e BPNN algorithm consists of two parts: the forward transfer of the input data and the back propagation of the error between the output data and the expected data [26,27]. One part is the input data through the input layer to the hidden layer, the hidden layer to the output layer. In the process, the calculation will be based on the given initial weight and threshold, and finally the output data will be obtained. e other part is to calculate the error change value δj (k) between the actual output and the expected value and then turn to the back propagation, and the error signal δj (k) is propagated back through the original connection path through the network to modify the weight of each layer of neurons. e value ωij (jk) and the threshold θj (k) until the target accuracy are reached. e hidden layer of BPNN can be one or more layers. It is proved theoretically that the neural network of a single hidden layer can approximate the nonlinear function with arbitrary precision, so that the model can realize the nonlinear mapping from input to output. With the increase of the number of hidden layers, the output error of the network will decrease [28]. erefore, the increase of hidden layers will improve the accuracy of the network but will make the network structure become complex, the running time become longer, and even lead to the overfitting phenomenon. erefore, after the number of neurons in the hidden layer is determined by empirical formula, the structure of BPNN tree fertilization model is 3-layer structure of 7-12-1. e BPNN has strong nonlinear mapping ability, selfadaptive and good self-learning ability, and strong fault tolerance. In the absence of a mathematical model, the nonlinear mapping between the amount of fertilizer applied to trees and the factors affecting the effect of fertilization can be realized. By learning and training the relationship between input and output data, when the actual output value and the expected value have a large error, the weights and thresholds of each layer of the neural network are constantly updated until the accuracy of network error is reached, so as to better improve the model.
In the figure, i, j, and k are the number of neurons in each layer; x i is any input signal of each group of data in the input p in the input layer; ω ij , ω jk and θ j , θ k , respectively, represent the weight and threshold of each layer; φ (·) and τ (·) represent the hidden layer and output layer activation functions, respectively. In the BP neural network algorithm, the hidden layer neurons input the signal net j when the sample p acts.
Mathematical Problems in Engineering 5 e hidden layer neuron node outputs a signal o j when the sample p acts: e output layer neuron node inputs the signal net k when the sample p acts: e output layer neurons output the signal y k when the sample p acts: In the formula, η is the error back propagation learning rate (0 < η < 1), δ j (k) is the error signal of each layer, and the weight and threshold are updated and iterated according to e error variation value of each layer gradually corrects the fertilization model during the back propagation process.

Particle Swarm Optimization BPNN Model.
When a single BPNN model is used to predict the precise fertilization amount of forest trees, the error between the actual output and the expected value is still large, because in the actual use, the convergence speed of the network is relatively slow, and it is easy to fall into the local optimization, even when there are few training samples, and there may be a fitting problem. Particle swarm optimization algorithm has better global optimization ability. erefore, the particle swarm optimization algorithm was proposed to optimize the BPNN model. e particle swarm optimization (PSO) algorithm is a swarm intelligence optimization algorithm [29,30]. e velocity characteristics, fitness value, and position were used to represent the motion characteristics of the particle. After initializing the potential optimal solution particles, the speed of the particles determines the search direction and distance of the particles, the fitness value determines the quality of the particles. When the particles move in the preset space, the position is changed according to the individual optimal solution and the global optimal solution. By updating the fitness value of the particle, the fitness value is the smallest.
e position corresponding to the particle is the optimal solution, as shown in Figure 2, the motion of particles in the particle swarm leaving its current position in the search area, the velocity v i+1 that this particle is about to update is a combination of various factors, It includes v i , v p , and v G , which are the current particle speed, the past best experience of individual particles, and the best experience of population particles.
e weights and thresholds of each layer in the BPNN can be represented by individual particle swarms. e individual swarms are coded, represented by the position vector of the particle swarm, and the optimal population particles are output by the iterative algorithm. After decoding, the BPNN global optimal is obtained. Weight and threshold establish a PSO-BPNN algorithm model. After decoding, the global optimal weight and threshold of BPNN are obtained, and the pso-bpnn algorithm model with 3layer structure of 7-12-1 is established.

GRA-PSO-BPNN Model Construction Process.
e degree of correlation can delicate the correlation in mutual matters and factors, and the Grey correlation degree analysis method can express the changeable situation in mutual matters and factors. Taking the timeliness and regionality of  the fertilization model into consideration, different forests have different factors that affect the accuracy of fertilization, resulting that there is no certain model to implement the fertilization model to achieve the accurate prediction so that the accurate rate of the prediction is not so well. is research adopts the analysis method of Grey Relational degree which can analyze the key factors that influence the fertilizer application amount in different areas and on the basis of the annual forest data updated by artificial forest, which can make the accurate fertilization come true. e grey relation analysis method is used to calculate the grey relation between input variables and output results. e specific steps are as follows: Step 1: Using the formula (9), establish a data matrix x i where X i (1), X i (2), X i (3),. . .,X i (300) represents the number of sequences. Among them, i represents site index, forest age, nutrient concentration related to dominant wood, biomass of dominant wood, average biomass of wood, fixed fertilizer utilization rate, and target yield, which are different factors affecting the fertilization effect.
Step 2: Establish an initialization change matrix x i ′ using Step 3: Use formula (11) to calculate the difference sequence, where k represents the sequence number: Step 4: Calculate the correlation coefficient and the grey correlation degree by using equations (12) and (13).
Step 4: Use formulas (12) and (13) to calculate correlation coefficient and grey correlation degree. φ is the resolution coefficient, 0 < φ < 1. In this study, φ � 0.5: c oi � 1 300 According to the calculation formula of grey correlation coefficient, the minimum difference min i min k Δ oi (k) and the maximum difference max i max k Δ oi (k) of N, P, and K fertilizers were 8.822e − 06, 3.483; 7.1719e − 05, 7.307; 2.388e − 04, 1.201, respectively; Table 4 shows the grey correlation coefficient and grey correlation degree between the input variables and the output of the prediction model of forest precise fertilization of nitrogen, phosphorus, and potassium.
According to the above analysis, the grey correlation degree of X 1 site index, X 2 forest age, X 3 dominant woodrelated nutrient concentration, X 4 dominant wood biomass, X 5 average wood biomass, and X 6 target yield is higher, and X 7 fixed fertilizer utilization rate is smaller; this study uses data reduction, grey correlation analysis, to reduce the influencing factors of different forest areas for predicting and calculating forest fertilization, identifying key factors affecting forest fertilization and providing more effective input for the prediction model.
Step 1 defines the input and output of the GRA-PSO-BPNN model. Compared with the PSO-BPNN model, using the grey correlation analysis, the site index, the forest age, the dominant wood-related nutrient content, the dominant wood biomass, the average wood biomass, and the target yield are used as the network input, and the actual fertilization amount is used as the output. erefore, there are 300 records, 200 for training and 100 for validation. In the formula, the number of neurons in the input layer x is 6, and the number of neurons in the output layer y is 1. e empirical formula is used to determine the number m of hidden neurons: e reduction of network error can not only adjust the number of hidden layer neurons, but also increase the number of hidden layer. However, increasing the number of hidden layers will complicate the network, increase the training time, and even overfitting. It is biased towards selecting the single hidden layer BPNN [31]. As mentioned above, the structure of GRA-PSO-BPNN model is determined to be a 3-layer structure of 6-12-1.
Step 2: Particle Swarm Initialization. e weights and thresholds in the BPNN model are the parameters to be optimized by the PSO algorithm. At the beginning, a certain number of particle swarm individuals can be randomly generated to represent, Mathematical Problems in Engineering and the global optimal position gbest and the individual optimal position pbest of the particle are initialized. e sum of ownership value and threshold value in network structure is the dimension d of particle swarm individual search space, that is, represents the position vector of the i th particle in the t th generation population.
Step 3: Particle Velocity Location Update. According to the following formula, the particle velocity is updated according to the individual optimal solution and the global optimal solution.
Among them, pbest (t) and gbest (t) are the individual optimal solution of the i th particle in the t th generation of the d-dimensional space, and the global optimal solution in the t th generation, r 1 and r 2 are in [0, 1], c 1 , c 2 are learning factors, w is the inertia weight, and the particle position is updated according to the following formula: e particle group can be evaluated by calculating the fitness value. e smaller the fitness value of the particle group is, the higher the fitness is. e speed and position of each particle are adjusted based on the fitness value. e particle fitness is calculated according to the following formula. Value F is where N is the number of samples in the test area; abs is the absolute value function; y i is the actual value of the sample i; and t i is the predicted value of the sample i.
Step 4: Optimal Population Particles: When the number of population evolution reaches the upper limit T or the iteration error reaches the set precision e, when the algorithm stops, the global optimal solution can be obtained and mapped to the weight and threshold of the BP network.
Step 5: GRA-PSO-BPNN Model Training. After determining the GRA-PSO-BPNN weights and thresholds, the training data are input to train the GRA-PSO-BPNN model. e entire model flow is shown in Figure 3.

Case Study
300 groups of tree growth and nutrient dynamic test data were measured in the test area, and the test data were divided into two parts. Among them, 200 groups of experimental data were used as training data to train the stage accurate fertilization model of forest based on neural network, and 100 groups of experimental data were used as the verification set of test fertilization model. BPNN models optimized by different algorithms are compared and verified in the experiment. Tansig function and purelin function are, respectively, assigned as the activation functions of hidden layer and output layer. In the process of adjusting parameters of BPNN model, the accuracy is found higher when the maximum training times of BPNN model reaches about 200 times, while big error and inconsistencies occur when the training times are more than or less than 200 times. erefore, the max of training times of the model is set to 200, the target error is 0.00001, and the learning rate is 0.1. Even so, the target accuracy is still unachievable, and the training time is about 5∼8 s. Based on the BPNN model, the PSO algorithm was used for optimization. When adjusting the parameters of the PSO-BPNN model, it is found that accuracy was higher when the population size was about 200 and the training time was about 8∼14 s. When the population size is significantly lower than 200, the accuracy becomes lower; when the population size is significantly higher than 200, the accuracy remains the same, but the training time becomes longer. At the same time, the accuracy is higher when the max of training times of the PSO-BPNN model is in the range of about 200 times; the accuracy however becomes lower when it is significantly lower than 200 times. When the training times greatly ex-   It can be seen from Figures 4-6 that the basic BPNN model can only predict the change level of forest fertilization amount and not accurately predict the actual fertilization amount, the accuracy is poor, the error range is about 20%, and the PSO-BPNN model predicts the amount of forest fertilization. e accuracy is obviously improved, the error between predicted fertilization amount and actual fertilization amount is reduced, and the error range is about 10%, which indicates that the optimization of PSO algorithm has a great influence on accurately predicting the amount of forest fertilization; the main influencing factors of determining the amount of forest fertilization by using grey correlation analysis method. After that, the GRA-PSO-BP prediction model of this paper further enhances the ability to predict the amount of forest fertilization. e predicted value of the GRA-PSO-BP neural network model has small fluctuations near the actual value. e model has a good prediction effect. e error between the predicted fertilization amount and the actual fertilization amount is within 5%, which can reflect the different forests. e change in nutrient demand during the growth phase can well guide the staged precision fertilization. e three models were used to predict the nitrogen, phosphorus, and potassium fertilizer application rates of 8 experimental sites and compared with the actual fertilization use. e results are shown in Table 5.
It can be seen from Table 5 that the error percentage between the predicted fertilization amount and the actual fertilization use amount of the GRA-PSO-BP prediction model is the smallest and the error is within 5%. e prediction accuracy of the GRA-PSO-BP prediction model is better than other fertilization models. e prediction accuracy is high, the error between the predicted fertilization amount and the actual fertilization usage is small, and the e model can fully describe the mapping relationship between input and output. According to the target production requirements of Chinese fir forest biomass, the practical amount of nitrogen, phosphorus, and potassium fertilizers can be predicted, which can guide the staged precision fertilization.
On the basis of the above data of actual fertilizer application of nitrogen, phosphorus, and potassium and predicted fertilizer application of nitrogen, phosphorus, and potassium, the paper calculated the mean square error between fertilization prediction results and expected values of different neural network algorithms, different absolute error values, and correlations one by one. e formula for calculating the mean square error is shown in (19). Besides the paper also compared the error and correlation to analyze the traditional BP Neural Network, BP Neural Network that was optimized by PSO algorithm, and the BP Neural Network Prediction Model that was optimized by GRA-PSO algorithm: In the formulation, n represents the number of calculated mean square error data, ε 1 , ε 2 , . . ., ε n indicates the difference between the actual fertilizer application amount of nitrogen, phosphorus, and potassium and the predicted fertilizer application amount of nitrogen, phosphorus, and potassium and σ means the square error. Tables 6-8 show the prediction models which are the traditional BP Neural Network, BP Neural Network that was optimized by PSO algorithm, and the BP Neural Network that was optimized by GRA-PSO algorithm, which predict the correlation and error value between the predicted fertilizer application of nitrogen, phosphorus, and potassium and the actual fertilizer application rate of nitrogen, phosphorus, and potassium.
According to the data in above rough the error value and correlation between the predicted results of N, P, and K fertilizer amount and the actual fertilization amount, it can be seen that the prediction accuracy of predication mode of BP neural network Finally, the mean and standard deviation of D ab and D ab are calculated, and DM statistics are calculated. e formulas are as follows: e DM test theory holds that the distribution of DM is in accordance with the requirements of the standard normal distribution. erefore, by querying the confidence values corresponding to DM in the standard normal distribution table, the confidence values of the significant differences between BP network model and PSO-BP network model and between PSO-BP network model and BP network prediction model optimized by GRA-PSO algorithm can be obtained while the predication of the fertilization amount of nitrogen, phosphorus, and potassium is advancing. It is shown in the results that when predicting the fertilization amount of   e prediction data of the network training are limited, so the different trained fertilization models are timeliness and regionality, while as for this research, because it is on the basis of the annual forest data that was updated by artificial forest and adopted the Analysis method of Grey Relational degree to analyze the key factors of fertilizer application amount in different areas, this can make the precision fertilization come true. After comparing these three prediction algorithms, the result came out that the square error, maximum absolute error, minimum absolute error, and average absolute error numbers of the BP Neural Network that were optimized by GRA-PSO algorithm are small, while the correlation number of it is large, so that the accuracy rate of the BP Neural Network that was optimized by GRA-PSO algorithm is high, and the square error, maximum absolute error, minimum absolute error, and average absolute error numbers of the traditional BP Neural Network and BP Neural Network that were optimized by PSO algorithm are relatively large, but the correlation number is relatively small, which draws the conclusion that the accuracy rate of these two prediction methods is relatively low.

Conclusions
It is inefficient to determine the amount of fertilization by traditional experience or basic theoretical model, which may also cause environmental pollution, so it is difficult to adapt to the new situation of precision forestry advocated by the contemporary world. As the current application of sensor network is getting more diversified, data acquisition in forest area is easier than before. rough real-time data update, real-time prediction can be carried out using high-tech to build fertilizer model. But different fertilization models will have timeliness and regionality. In different forest ecological areas, it seems to be impossible to accurately make a prediction at present using the definite fertilization model with an eye to different factors that affect the precision fertilization of trees. e traditional methods to determine the amount of fertilizer application are empirical method, fertilizer effect function method, and nutrient balance method. e empirical method means that forest farmers fertilize trees on their willing which is according to their experience and without basis. is method is simple and sometimes feasible, but in most cases, it is not accurate enough, which often results in insufficient or excessive fertilizer application and make it not fit the needs of high yield and efficiency of forestry construction. e fertilizer effect function method needs a lot of experimental data, but the successful rate of the fitting is not high, and a lot of data resources are just wasted, so that a lot of financial wealth, material, and human labors are wasted too. e nutrient balance method often needs to come up with expert experience, because it requires too many parameters and also there is a relative error in the calculation of fertilization. It is important to find an effective and accurate fertilization model. Artificial neural network has a strong ability to solve nonlinear problems. In order to solve the problem of nonlinear precision fertilization, this paper introduces the neural network modeling method and uses BPNN, PSO-BPNN, and GRA-PSO-BPNN to construct three kinds of forest precise fertilization models for comparative verification and result analysis. rough comparison, the accuracy of the BPNN model is low and the error range is about 20%, and it is unstable. Although the PSO-BPNN model has high accuracy, the error range is about 10%, considering that the fertilization model mentioned above is time-sensitive and regional, the factors that may affect the fertilization are different in different time and different regions. In some geographically complex areas, more factors need to be considered to ensure accurate fertilization. Since it is difficult for the PSO-BPNN model to identify the key influencing factors of forest fertilizer application in different regions, the PSO-BPNN model has poor scalability. However, the key factors influencing forest fertilizer application in different regions can be determined by grey relational analysis. e constructed GRA-PSO-BPNN forest tree precision fertilization model not only has high accuracy, the error range is about 5%, and it has a faster calculation speed, but also applicable to the determination of forest fertilizer application in different regions and can constantly update the forest fertilizer model with good expansibility.

Discussion
With the in-depth development of modern forestry, precision forestry is receiving unprecedented attention. e emergence of forest precise fertilization model will greatly promote the forestry production. Forest economy plays an important role in the society improvement, and to improve the forest economy, only the development of artificial fastgrowing and high-yielding forest can be depended, and the effective way to improve the productivity of artificial forest is just precision fertilization. e prediction of precision fertilization is a scientific and innovative technology which can achieve five goals: ensuring high quality, increasing the production, improving soil quality and ecology, and high efficiency; taking the factors which affect the amount of fertilizer application as the basis and making the fertilizer application rate and fertilization target clear before implementing it; linking the fertilization with economic benefit; achieving the quantitative fertilization of forest trees under the realization of target benefit, making a prediction of the benefit that earned by the fertilization, and according to the investigation of basic data and soil chemical analysis in forest area, so that the ratio of tree fertilization under the target benefit can be calculated; making forest fertilization become an important safeguard measure to improve the quality and yield of commercial forest. At the same time, as a direction of precision forestry technology, precision fertilization model will attract more and more attention of forestry science and technology personnel and forestry production managers. In this paper, the GRA-PSO-BPNN prediction model of forest periodic precise fertilization is of far-reaching significance to solve the problem of forest fertilization that has been troubling people for a long time, which can be applied by forest-related workers in forestry ecological zones where conditions are basically the same. If the generalization is carried out in a larger and more complex area to improve the generalization ability of the model, it is necessary to consider other factors affecting the amount of forest fertilization, such as soil type, organic matter, pH, and other factors. In some cases, the GRA-PSO-BPNN precise fertilization model for trees may need to be updated based on the experience of experts, which is what needs to be further studied.
Data Availability e raw/processed data required to reproduce these findings cannot be shared at this time as the data also forms part of an ongoing study.

Conflicts of Interest
e authors declare no conflicts of interest.

Authors' Contributions
Chen Zuxing and Wang Dian contributed analytical tools and ideas, collected and processed the data, conducted experiments and wrote papers, and finally finalized the final version of the paper.