• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Research on Forecasting Flowering Phase of Pear Tree Based on Neural Network

    2021-12-14 06:05:14ZhenzhouWangYinuoMaPingpingYuNingCaoandHeinerDintera
    Computers Materials&Continua 2021年9期

    Zhenzhou Wang,Yinuo Ma,Pingping Yu,*,Ning Cao and Heiner Dintera

    1School of Information Science and Engineering,Hebei University of Science and Technology,Shijiazhuang,050000,China

    2School of Internet of Things and Software Technology,Wuxi Vocational College of Science and Technology,Wuxi,214028,China

    3German-Russian Institute of Advanced Technologies,Karan,420126,Russia

    Abstract:Predicting the blooming season of ornamental plants is significant for guiding adjustments in production decisions and providing viewing periods and routes.The current strategies for observation of ornamental plant booming periods are mainly based on manpower and experience,which have problems such as inaccurate recognition time,time-consuming and energy sapping.Therefore,this paper proposes a neural network-based method for predicting the flowering phase of pear tree.Firstly,based on the meteorological observation data of Shijiazhuang Meteorological Station from 2000 to 2019,three principal components (the temperature factor,weather factor,and humidity factor)with high correlation coefficient with the flowering phase of pear tree are obtained by using the principal component analysis method.Then,the three components are used as input factors for the BP neural network.A BP neural network prediction model is constructed based on genetic algorithm optimization.The crossover operator and mutation operator in the adaptive genetic algorithm are improved.Finally,the meteorological sample data from 2013 to 2019 are used to test and verify the algorithm in this paper.The results demonstrate that,the model can solve the local optimization problem of the BP neural network model.The prediction results of the flowering phase of pear tree are evaluated in terms of relevance and prediction accuracy.Both are superior to the traditional effective accumulated temperature and the prediction results of the stepwise regression method.This method can provide more reliable forecast information for the blooming period,which can provide decision-making reference for improving the development of tourism industry.

    Keywords:Pear flower;flowering phase;principal component analysis;BP neural network;prediction model

    1 Introduction

    Flowering period data of ornamental plants are important for the organization and development of flower-viewing tourism activities.These data can provide guidance for tourists to formulate flower viewing tourism plans,and provide reference for the planning of flower-viewing activities in scenic locations.Ornamental plants have a strong seasonal flowering period,short duration,and are easily affected by the external environment.The earlier the flowering period is monitored,the earlier best views and routs can be got.Traditional manual observation is timeconsuming and laborious.It cannot effectively predict and release information in a timely manner.

    There are many studies on the relationships among the plant flowering period and climate,climate change and prediction technology of flowering period.Gonsamo et al.[1]simulated the anomalous changes in the initial flowering of 19 Canadian plants since 1948.Morin et al.[2]predicted the changes in the initial leaf development of 22 woody plants in North America in the next 100 years.Zhong et al.[3]established a temperature-based temporal and spatial prediction model for the flowering period of Chinese ornamental plants.Liu et al.[4]used wavelet analysis,correlation analysis and other methods to assess the impacts of climate change on flowering phase of peach blossom.Aono [5]investigated the phenological data of cherry blossoms in Edo since the 17th century and their application in March temperature estimation.Ahn et al.[6]conducted a seasonal forecast of surface temperatures and the first flowering period in Korea.At present,the effective accumulated temperature rule is mainly used to predict the flowering phase,that is,the effective accumulated temperature for a certain period of time after the plant’s flowering threshold temperature is met.Dong et al.[7]analyzed the relationships between the flowering phase of peach trees and temperature elements by using the phenological observation data of peach trees in the blooming period and the daily average temperature observation data of the climate station in the Hekou area of Dandong City.The active accumulated temperature,effective accumulated temperature,and sliding accumulated temperature were used to predict the blooming periods of peach trees in the Dandong area.Through analysis,Yao et al.[8]found that the temporal accumulated temperature phase had a small correlation with the variation coefficient of daily accumulated temperature,but had a large correlation with the absolute value of the correlation coefficient of daily ordinal number at flowering stage.These results indicate that temporal accumulated temperature can better reveal the internal relationship between accumulated temperature and crop growth on a smaller time scale.Judging the flowering phase of plants by accumulated temperature is simple to perform and has a certain biological basis.However,because this method disregards the influences of meteorological factors such as sunlight and precipitation on the flowering period,the accuracy of the predicted flowering period is low.

    With the development of computational science,neural networks have been widely investigated and applied in various fields of science and engineering,and relatively satisfactory prediction results have been achieved.Long et al.[9]employed a 3-layer BP neural network to establish a fog forecast model.Based on the air pollutant data of Delhi,India,from 2014 to 2016.Ganesh et al.[10]established an artificial neural network prediction model that is based on the conjugate gradient descent method to predict the air quality index in a specific area.Zhang et al.[11]proposed a BP neural network earthquake damage prediction model that is based on the LM algorithm,which can realize that only part of the key data of a single building needs to be extracted.He et al.[12]established a BP neural network soil temperature prediction model that is optimized by a genetic algorithm to accurately predict the soil temperature in winter orchards.Aiming at the problem of single sources of forecast data and low accuracy in current atmospheric visibility forecasting algorithms.Wang et al.[13]constructed a BP neural network model optimized based on a genetic algorithm.Pooya et al.[14]employed artificial neural networks (ANNs),genetic algorithms and artificial neural network integration (GAINN) to predict the efficiency of fiber filters.Sun et al.[15]established a rape blossom forecast model based on the BP neural network and explored the application of the BP neural network in the field of flowering forecast.Neural network technology provides a new way to simulate and predict flowering.

    Aimed at forecasting the flowering phase of ornamental plants,existing research focuses on cherry blossoms,while there are few simulations of the flowering phase of other ornamental plants.This paper propose a prediction method based on neural network for flowering phase of pear trees.Through principal component analysis,three principal components with higher correlation coefficient with the flowering phase of pear trees were obtained as input factors of the neural network.The problem of low accuracy of the forecasted flowering period is improved.This paper is organized as follows:In the second section,the sources of pear phenological data were introduced,and three principal components were obtained by using principal component analysis.In the third section,a BP neural network model is established based on genetic algorithm optimization,and the model parameters are instantiated.In the fourth section,the model is trained and tested,and compared with the traditional effective accumulated temperature method.The fifth section uses the trained model to forecast the flowering stage of pear trees in 2020.In the last,a conclusion is provided.

    2 Data Sources and Research Methods

    2.1 Flowering Data of Pear Trees

    The phenological data of pear flowers were obtained from the observation data of pear trees planting area provided by Shijiazhuang Meteorological Bureau from 2009 to 2020.The observation basis and standards are the “China Phenological Observation Network” observation standards and the China Meteorological Administration “Agricultural Meteorological Observation Specifications.” The initial flowering period of a plant is defined as the data when the first fully open flower begins to appear on the observation plant.The full bloom period is defined as the data when more than half of the flower buds on the observation plant unfold and show the petals [16].

    The flowering day ordinal number is the conversion of the observation date into a diurnal sequence.January 1 is taken as the beginning data of the flowering diurnal sequence,and the diurnal number is 1,and so on,such as February 5,the diurnal number is 36 [17].

    By analyzing the location observation phenological data and parallel observation meteorological data of the pear trees planting area from 2009 to 2020,it is calculated that the average initial flowering time of pear trees is April 3;the average blooming period is April 10;and the average end of flowering is April 20.The earliest flowering date is March 20,and the latest data is April 16.The flowering phase begins as early as March 28 and as late as April 24.The earliest date of the end of flowering is April 8,and the latest date is May 4.The flowering phase of pear trees lasts for approximately 20 days on average from the beginning to the end of flowering.According to years of observation and analysis,pear blossoms will enter a suitable viewing period after 2-3 days of initial flowering [18],as shown in Tab.1.

    2.2 Meteorological Factors Determining by Principal Component Analysis Method

    2.2.1 Selection of Meteorological Factors

    Meteorological factors,including temperature,precipitation,sunshine,etc.,impose obvious constraints on the growth of pear trees.However,the restrictive effects of different meteorological factors are divided into strong and weak factors that are relatively independent.

    This paper selects the daily maximum temperature,daily minimum temperature,daily average temperature,daily precipitation,daily sunshine duration,daily average ground temperature and daily average relative humidity from 2000 to 2019 at Shijiazhuang City Meteorological Station.These data have been strictly controlled by artificial quality.The average values of each meteorological element in winter from 2001 to 2019 are obtained by calculation.

    Table 1:Design index of pear blossom index

    2.2.2 Principal Component Analysis Method

    Principal Component Analysis (PCA) aims to analyze the characteristics of the covariance matrix to reduce the dimensionality of the data while maintaining the largest contribution to the variance of the data set [19].By dimensionality reduction,the original index is transformed into one or several comprehensive indexes:namely,principal components.The principal components are not correlated.The steps of the principal component analysis to analyze the meteorological factors related to pear flower blooming are described as follows:

    Step 1:Select the meteorological factors that are related to the flowering period of pear blossoms,and perform z-score normalization (zero-mean normalization) on the original data.The formula is expressed as follows:

    In the formula,is the standardized meteorological factor,xijis the original data of each meteorological factor,andandsjare the mean and standard deviation,respectively,of each meteorological factor.

    Step 2:Calculate the correlation coefficient matrix R of the standardized meteorological factor data:

    In the formula,Rijis the correlation coefficient between the i-th standardized index and the j-th standardized index.

    Step 3:According to the meteorological factor correlation coefficient matrix R,find the eigenvalue,the contribution rate of the principal component and the contribution rate of the cumulative variance.Determine the number of principal components,and solve the characteristic equation:

    Find the eigenvalueλi (i=1,2,...,p),wherepis the number of main components.The eigenvalue is the variance between the principal components,which reflects the influence of each principal component.The sharing rate of the principal componentWiis:

    Cumulative contribution rate:

    According to the principle of selecting the number of principal components,eigenvalues greater than 1 and cumulative contribution rates greater than 90% are selected.

    2.2.3 Results of Principal Component Analysis Method

    The meteorological observation data of meteorological stations from 2000 to 2019 were used,and the eigenvalues of the correlation matrix and the contribution rate of each principal component were obtained.The results are shown in Tab.2.M principal components with characteristic roots greater than 1 and cumulative contribution rates greater than 90% are selected.From Tab.2,the characteristic roots of the first three are greater than 1,and the cumulative contribution rate reaches 91.94%.The three principal components can reflect most of the information of the original meteorological indicators.This process can reduce the complexity of the original data and achieve the purpose of dimensionality reduction.

    Table 2:2000-2019 Shijiazhuang weather station sample correlation coefficient matrix eigenvalues

    The eigenvectors of the three principal component factors are listed in Tab.3.Tab.3 shows that for the eigenvector of the first principal component,the factors with larger and positive eigenvalues are the average daily minimum temperature in winter,average daily maximum temperature in winter,average daily average temperature in winter,and average daily average ground temperature in winter.The flowering period is closely related to temperature,so the first principal component can be referred to as the temperature factor.For the eigenvector of the second principal component,the factor with a larger and positive eigenvalue is the average sunshine duration in winter,the larger and negative factor is the average daily precipitation in winter and the average daily average low temperature in winter.Therefore,the flowering period is obviously positively correlated with sunshine duration and negatively correlated with precipitation and relative humidity.The second principal component is referred to as the weather factor.Among the eigenvectors of the third principal component,the factor with a large and positive eigenvalue is the average daily average relative humidity in winter,which can be referred to as the humidity factor.

    Table 3:Principal component feature vector of meteorological data from Shijiazhuang Meteorological Station 2000-2019

    According to the characteristic vector of the principal components (Tab.3),the linear equations between the three principal components and the meteorological factors can be obtained:

    In addition,studies have pointed out that environmental parameters such as soil temperature,cold demand and hourly accumulated temperature,as well as management measures such as tree age,fertilization and irrigation,are also related to the flowering phase of trees.This article forecasts the overall flowering period in the region from the perspective of weather forecasting services,instead of a single fruit tree or orchard.Taking into account the availability of data,these parameters are not involved.

    3 BP Neural Network Prediction Model Based on An Improved Genetic Algorithm

    The research and application of neural networks in various fields of science and engineering are extensive,and relatively satisfactory prediction results have been achieved.The BP neural network,which is the most commonly employed from of artificial neural network [19],has a nonlinear mapping ability and excellent fault tolerance and can be applied to the flowering period prediction model of this article.The learning process of the BP neural network is actually a process of training the network with training samples [20].Aimed at problems such as difficulty in forecasting the flowering period of ornamental plants and inaccurate prediction of long-term data,this paper uses an improved genetic algorithm to iteratively identify the structure and designs a BP neural network prediction model that is based on genetic algorithm optimization for flowering prediction.

    3.1 Neural Network Structure

    The BP neural network is generally composed of an input layer,a hidden layer and an output layer.And a large number of neurons are connected to each other as network nodes [21].Considering a 3-layer neural network as an example,set the input layer toxi,the hidden layer toyh,the output layer tool,and the expected value of output totl.The structure diagram is shown in Fig.1.

    Figure 1:BP neural network structure diagram

    The error function between the output and the expected value is:

    whereθlis the threshold andf′is the inverse function of the transfer function.

    whereθhis the threshold,and the weight in thek+1 iteration is

    Because Δvlh=-η·?E/?vlh=ηδlyh,Δwhi=-η′·?E/?whi=η′δ′hxi,the threshold of the output layer and the change in the threshold of the hidden layer can be expressed as

    where η and η′are the learning rates.

    Meanwhile the King showed the Princess inside all his gold wares, every single bit of it--dishes, goblets, bowls, the birds and game, and all the wonderful beasts

    3.2 Improved Adaptive Crossover Operator and Mutation Operator

    Traditional genetic algorithms use fixed crossover probability and mutation probability [22].The crossover probability is generally selected between 0.3 and 0.7,and the mutation probability is generally selected between 0.1 and 0.3 [23].However,it is difficult to optimize the crossover and mutation probability.If a larger crossover probability is selected,new individuals will be generated too quickly,and individuals with high fitness values are easily destroyed,which may cause the algorithm to evolve into a random search algorithm.If a smaller crossover probability is selected,new individuals will be generated.If the speed is too slow,the algorithm can easily fall into the local optimum.If a larger mutation probability is selected,although the diversity of the population can be maintained,the advanced mode of population inheritance will be destroyed.If a small mutation probability is selected,the diversity of the population will gradually decrease,which leads to the rapid loss of outstanding individuals and failure to achieve the desired result [24].Due to the various shortcomings of traditional genetic algorithms,many scholars have continuously improved and proposed many improved algorithms,such as the adaptive genetic algorithm (AGA).The formula of the AGA is presented as

    wheref′is the larger fitness value of the two individuals to be crossed;fmaxis the maximum fitness in the population;favgis the average fitness of the population;frepresents the fitness of the variant individual;andk1-k4are adaptive control parameters.However,AGA evolves slowly in the initial stage and is prone to stagnation.Excellent individuals are in a static state:that is,no crossover and mutation operations occur.At this time,individuals with high fitness values in the population are likely to converge locally and fail to reach the global optimum.The genetic algorithm can easily evolve in the direction of local convergence.Therefore,this article improves adaptive crossover probability and mutation probability as

    wherepcminandpcmaxrepresent the lower limit and upper limit,respectively,of the crossover rate;pmminandpmmaxrepresent the lower limit and upper limit,respectively,of the mutation rate.The individual’s crossover rate and mutation rate are still linearly transformed between the average fitness and the maximum fitness according to the fitness of the individual.

    3.3 Establishment of the BP Neural Network Model Based on Genetic Algorithm Optimization

    To improve the stability of the BP algorithm and solve the problem of local optimal solutions,this paper constructs a BP neural network model that is based on genetic algorithm optimization [25].The genetic algorithm does not rely on gradient information and uses its global search advantages to optimize the BP neural network to obtain the best solution to the problem [26].The genetic algorithm uses binary coding to assign a real number string to each individual of the population,and the population individual represents all the weights and thresholds of the BP neural network layer.According to the designed fitness value,the coded individual is employed to train the BP neural network,and then by the processes of selection,crossover,mutation,etc.,the structure is iteratively identified by the improved genetic algorithm.Assume that there are a total of N fuzzy rules and the improved genetic algorithm is utilized for identification.In the identification process,the fitness of all individuals in the population P(t)isF1=1/J.The specific process is described as follows:

    (1) Initialize the population P(0).Refer to the range of the data vector to perform coding operations to generateNCindividuals as the initial population P(t)and set the maximum number of iterations T.The initial chromosomes are composed of regular codes that are generated by the BP algorithm.

    (2) Calculate the fitness function.The fitness function F,which is also known as the evaluation function,can realize the measurement of the pros and cons of individuals in the group.The fitness function can be expressed by the sum of the absolute value of the error between the predicted output of the BP neural network and the actual output.The calculation formula is

    (3) The selection operation is performed after completing the evaluation of the fitness function.The commonly employed selection methods of genetic algorithms include the fitness ratio method,partial selection method and roulette selection method.This article chooses the roulette selection method.The probability of an individual that is selected in this method is proportional to the fitness value that is calculated in Step (2).Assume that the fitness of individualiisfi;then

    wherenis the size of the group andpiis the probability that individualiwill be selected.

    (4) Crossover operation.Crossover operation refers to the operation of partially reorganizing the structure of the two selected first-generation individuals according to the principle of crossexchange of biological chromosomes to form a new individual,and the crossover operation acts on the selected group.The crossover probability isPc.

    (5) In the mutation operation,the mutation operator judges the mutation probability of the individual in the group and then randomly selects the mutation position of the mutated individual to obtain the result of the mutation.The mutation operation is applied to the selected population.The mutation probability isPm.

    (6) After the populationP(t)completes Steps (3)-(5),the next generation populationP(t+1)is obtained.

    (7) If the termination condition is met,proceed to the decoding operation;otherwise,return to Step (3).

    The overall design process is shown in Fig.2:

    Figure 2:Improved genetic algorithm BP neural network process

    3.4 Model Parameter Setting

    This study selects a 3-layer BP neural network and uses principal component analysis to obtain three principal component factors that affect the flowering phase of pear trees as the input layer of the BP neural network.Coupled with the input layer threshold,the model has 4 input vectors.The daily ordinal number of the full-blossom period of pear trees (best viewing period)is utilized as the output vector of the output layer,and thus the number of nodes in the output layer is 1,and the output vector is 1.The determination of the number of neurons in the hidden layer is related to the actual problem that needs to be solved,and the optimal number of neurons is generally determined by a formula.

    wherenis the number of neurons in the input layer,qis the number of neurons in the output layer,αis an integer in [1,10],andmis the number of neurons in the hidden layer.The determination of the number of hidden layers has an important role in the accuracy of the prediction model [26].In the case of the same sample test,this paper compares the training errors using different hidden layer numbers and obtains the results in Tab.4 via 6 trials.It can be seen from the test results that when the number of hidden layers is 6,the training error is 0.1868,and the best training result is obtained.

    Table 4:Network training error of different hidden layer nodes

    The relevant parameters of the genetic algorithm are listed as follows:the number of iterations is 10;the randomly generated population size is 100;the crossover ratePc=0.3,and the mutation ratePm=0.1.

    4 Result Analysis

    4.1 Model Training and Test Results

    This article uses the MATLAB R2014a neural network development toolbox and Sheffield genetic algorithm toolbox for training.In the process of training the determined network structure,to prevent “transition training,” the network adopts the following convergence rules:if the samples with an absolute error of less than 2 reach 85% of the total samples,the training is stopped.Otherwise,the maximum number of training iterations (10,000 iterations) is specified,and the samples are used to verify the network when the training meets the requirements.

    In this paper,the flowering phenology data of pear trees from 2000 to 2012 in the pear planting area of Shijiazhuang are employed as the training input samples of the BP neural network.A trial report is conducted using the 2013-2019 sample.When using the three main component factors as the network input,the comprehensive situation of the meteorological factors is considered,the main factor is enlarged and the secondary factor is weakened.Network training and testing were carried out according to the forecast model,and the training effect map (Fig.3)and prediction effect map (Fig.4) were drawn via training.

    Comparing Figs.3 and 4,it can be seen that the average error of the BP neural network training,which is based on genetic algorithm optimization with 3 principal component factors as input,is 1.5 days,and the training effect is reasonable.This finding shows that the forecast model has practical application value.

    Figure 3:Training effect diagram

    Figure 4:Forecast effect diagram

    4.2 Analysis of Results of Different Forecasting Models

    Using the data samples of the flowering phase of pear trees from 2016 to 2019,the improved BP neural network method proposed in this article,the traditional effective accumulated temperature law method and stepwise regression method are verified.The verification results of the three methods are shown in Tab.5.From the table,we can see that the error of the flowering period forecast using the BP neural network model,which is optimized based on the genetic algorithm,is 1 day.The result for the effective accumulated temperature prediction method is 3.25 days,and that of the stepwise regression prediction method is 2.75 days.The flowering date predicted by this paper is closer to the actual flowering date.The meteorological factors filtered by the principal component analysis are used as model inputs,and the modeling method based on the BP neural network optimized by the genetic algorithm improves the accuracy of flowering forecasting.

    Table 5:Forecast results of different forecast models

    4.3 Accuracy of Flowering Fitting

    To intuitively compare the applicability and accuracy of the method in this paper for the prediction of the flowering phase of pear trees in Shijiazhuang,a combination of internal inspection and cross-checking is employed.First,internal inspection is carried out.The parameters fitted from 2000 to 2019 are used to simulate the phenological sequence of the flowering phase of pear trees.The observed sequence of the flowering phase of pear trees is compared with the simulated sequence,and the variance (R2),significance and root mean squared error (RMSE)are calculated.Second,cross-check (also referred to as external check) is conducted.Cross-check adopts the method of elimination one by one,that is,after eliminating the observed value of flowering phase of pear trees in a certain year,the data of other years are used to fit parameters to simulate the phenological period of the eliminated year (full blooming period).After each year’s observations are separately eliminated,the series for cross-checking can be obtained.Additionally,we compare and analyze this sequence with respect to the phenological observation sequence.

    wheredi(x)is the original date sequence of the flowering phase of pear trees;di(x)2is the square of the original date sequence of the flowering phase of pear trees;diobsis the date sequence of the flowering phase of pear trees simulated by the model;andnis the number of test samples.

    We compare and analyze the actual observation sequence of the flowering phase of pear trees with the simulated sequence of a certain year.If the difference between the two sequences is less than 3 days,the prediction is considered accurate and recorded as “1”;otherwise,the prediction is considered wrong and recorded as “0”.The accuracy of the forecast is the percentage of the number of accurate forecasts in the total number of years.It can be seen from Tab.6 that the prediction accuracy of the algorithm in this paper is higher.

    Table 6:The prediction accuracy of the three algorithms

    5 Trial Report of the Flowering Phase of Pear Trees

    After the training is completed,the topology of the neural network is readjusted.The number of nodes in the input layer includes the number of weather factors,the number of days in the initial flowering period and the accumulated temperature since the initial flowering period.The number of nodes in the output layer is the number of forecasted objects (1,use 1/0 to indicate whether flowering occurs the next day).

    Accumulated temperature refers to the daily average temperature that is accumulated in the corresponding period of time when a plant completes a certain period or all-growth period.Accumulated temperature is an important indicator to measure the requirements of crop growth for thermal conditions and to evaluate thermal resources.According to the statistics of the accumulated temperature from January 1 to the full flowering period of pear blossom in Shijiazhuang City,the daily average temperature is greater than or equal to the threshold temperature of pear trees (any value in the range of 0.1-20.0°C).The relationship between the accumulated temperature and the flowering period is analyzed to determine the most stable accumulated temperature threshold of 6.9°C in this experiment.

    A sample of a certain year is selected to form a calculation sample of the model in this paper within 10 days before and after flowering.The forecast value is 1 from the day before the flowering period,which indicates that the next day has entered the flowering period,and the number of days before that is set to 0.For example,the initial of flowering period of a certain year is March 20,the starting data of the sample is March 10-29,the forecast result before March 19 is 0,and the forecast result from March 19 is 1.The flowchart is shown in Fig.5.

    Figure 5:Forecast flow chart

    Using the flowering period forecast model in this paper,the flowering phase of pear trees in 2020 in the Shijiazhuang area is forecasted.The forecast will be carried out every day beginning on March 1st.The maximum forecast time limit is 10 days,that is,the cycle in Fig.5 is iterated 10 times.If the forecast result is 0,it is considered that no flowering will occur within 10 days.If 1 appears on the first time in the cycle,the next day is considered to be the flowering period.After the initial flowering data is forecasted,the blooming day forecast will continue.The first cycle means that March 25 predicts whether blooming will occur on March 26.If the result is 0,the second cycle begins to predict whether blooming will occur on March 27.If the result is 1 for three cycles,then the forecast result indicates that blooming will occur on March 28 (error of 0 day).If the results of 1 appears after 8 cycles,then the forecast result indicates that blooming will occur on April 4 (with an error of 1 day).In terms of service,the model provides a reasonable reference for the public to schedule their visits and has achieved excellent social benefits.

    6 Conclusion

    This paper selects the pear trees planting area in Shijiazhuang as the target area for prediction.Three principal component factors that affect the flowering period of pear are obtained through principal component analysis,namely,the temperature factor,weather factor and humidity factor.The BP neural network optimized based on a genetic algorithm is applied to the forecasting of the pear flowering period,and the crossover operator and mutation operator in the adaptive genetic algorithm are improved.The experimental results show that the error of the flowering period forecast using the forecast model in this paper is 1 day,the value of the effective accumulated temperature forecast method is 3.25 days,and that of the stepwise regression forecast method is 2.75 days.It can be seen that the algorithm proposed in this paper achieves better results than the traditional methods in predicting the flowering period,and also verifies the correlation between the flowering data and meteorological factors.The algorithm presented in this paper can provide more reliable forecast information of the flowering phase and provide reasonable reference for the public’s travel arrangement.

    Funding Statement:This research was funded by the Science and Technology Support Plan Project of Hebei Province (Grant Number 19273703D) and the Science and Technology Research Project of Hebei Province (Grant Number ZD2020318).

    Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

    国产欧美日韩一区二区三| 亚洲精品国产区一区二| 国产高清国产精品国产三级| 久久精品成人免费网站| 男女床上黄色一级片免费看| 热re99久久国产66热| 国产单亲对白刺激| 中文字幕另类日韩欧美亚洲嫩草| 女性生殖器流出的白浆| 天堂动漫精品| 一区二区三区激情视频| 国产成人啪精品午夜网站| tocl精华| 亚洲av电影在线进入| 亚洲精品久久成人aⅴ小说| 一级a爱视频在线免费观看| 亚洲av成人一区二区三| 亚洲欧洲精品一区二区精品久久久| 欧美人与性动交α欧美精品济南到| 亚洲aⅴ乱码一区二区在线播放 | av有码第一页| 99re在线观看精品视频| 久久精品91蜜桃| 欧美日韩瑟瑟在线播放| 亚洲色图综合在线观看| 高清av免费在线| 丝袜人妻中文字幕| 一个人免费在线观看的高清视频| 精品国产一区二区三区四区第35| 又紧又爽又黄一区二区| 国产精品电影一区二区三区| 欧美一级毛片孕妇| 国产av又大| 亚洲av第一区精品v没综合| 十八禁网站免费在线| 久热这里只有精品99| 一进一出抽搐gif免费好疼 | 午夜福利在线免费观看网站| 久久国产乱子伦精品免费另类| 国产精品1区2区在线观看.| 国产亚洲精品一区二区www| 久久久久精品国产欧美久久久| 热99re8久久精品国产| 嫁个100分男人电影在线观看| 视频区欧美日本亚洲| av视频免费观看在线观看| 少妇粗大呻吟视频| 水蜜桃什么品种好| netflix在线观看网站| xxxhd国产人妻xxx| 成人影院久久| 热99re8久久精品国产| 国产精品偷伦视频观看了| 精品国产乱码久久久久久男人| www.999成人在线观看| 久久狼人影院| 亚洲伊人色综图| 久久香蕉精品热| 日韩国内少妇激情av| 看片在线看免费视频| 免费av中文字幕在线| 亚洲欧美激情综合另类| 一区二区三区精品91| 大型av网站在线播放| 久久中文字幕一级| 女性生殖器流出的白浆| 成年人免费黄色播放视频| 国产精品香港三级国产av潘金莲| 91字幕亚洲| 国产亚洲av高清不卡| 热99国产精品久久久久久7| 亚洲午夜理论影院| 级片在线观看| 久久久国产一区二区| 亚洲一区高清亚洲精品| 国产精品成人在线| 最近最新免费中文字幕在线| 交换朋友夫妻互换小说| 精品国产乱码久久久久久男人| 国产精品二区激情视频| 国产av一区二区精品久久| 日本五十路高清| 好男人电影高清在线观看| 搡老熟女国产l中国老女人| 人人澡人人妻人| 欧美黄色淫秽网站| 亚洲色图 男人天堂 中文字幕| 中国美女看黄片| 国产精品久久久久久人妻精品电影| 黑丝袜美女国产一区| 国产精品久久久人人做人人爽| 欧美av亚洲av综合av国产av| 久久久久亚洲av毛片大全| 一区二区三区激情视频| 欧美成人午夜精品| 一级a爱视频在线免费观看| 亚洲 国产 在线| 成年版毛片免费区| 国产亚洲精品一区二区www| 成人永久免费在线观看视频| 男男h啪啪无遮挡| 我的亚洲天堂| www.精华液| 成人亚洲精品av一区二区 | 一区二区三区激情视频| 亚洲人成网站在线播放欧美日韩| 美女国产高潮福利片在线看| 免费在线观看视频国产中文字幕亚洲| 亚洲欧美一区二区三区久久| 高清欧美精品videossex| av天堂在线播放| 麻豆成人av在线观看| 国产精品久久久久成人av| 色播在线永久视频| 一级毛片高清免费大全| 亚洲欧美一区二区三区久久| 亚洲一区二区三区不卡视频| 黑人操中国人逼视频| 香蕉久久夜色| 国产成+人综合+亚洲专区| 超碰成人久久| 高清欧美精品videossex| 黑人操中国人逼视频| 看黄色毛片网站| 热re99久久国产66热| 看黄色毛片网站| netflix在线观看网站| 亚洲av片天天在线观看| 午夜日韩欧美国产| 黑人操中国人逼视频| 一本综合久久免费| 一本大道久久a久久精品| 色婷婷久久久亚洲欧美| 在线十欧美十亚洲十日本专区| 免费在线观看日本一区| 成人精品一区二区免费| 国产亚洲精品久久久久久毛片| 不卡一级毛片| 午夜精品国产一区二区电影| 久久久久亚洲av毛片大全| 99re在线观看精品视频| 久久精品国产清高在天天线| 国产精品一区二区三区四区久久 | 亚洲av电影在线进入| 欧美激情高清一区二区三区| 亚洲av第一区精品v没综合| 狠狠狠狠99中文字幕| 中国美女看黄片| 国产成人影院久久av| 成人黄色视频免费在线看| 午夜老司机福利片| 视频区图区小说| 最近最新免费中文字幕在线| 青草久久国产| 超色免费av| 女生性感内裤真人,穿戴方法视频| 午夜福利一区二区在线看| 国产欧美日韩一区二区三区在线| www.自偷自拍.com| 精品电影一区二区在线| 国产av在哪里看| 亚洲国产精品999在线| 国产不卡一卡二| 亚洲成人精品中文字幕电影 | 欧美日韩亚洲综合一区二区三区_| 国产一区二区在线av高清观看| 色在线成人网| 性欧美人与动物交配| 久久影院123| 久久精品91蜜桃| 亚洲在线自拍视频| 麻豆一二三区av精品| 国产精品99久久99久久久不卡| 亚洲精品在线观看二区| 亚洲av日韩精品久久久久久密| 一级作爱视频免费观看| 国产精品一区二区免费欧美| 人人妻人人添人人爽欧美一区卜| 99精品在免费线老司机午夜| 女人爽到高潮嗷嗷叫在线视频| 岛国在线观看网站| 亚洲欧美日韩另类电影网站| 999久久久精品免费观看国产| 99国产极品粉嫩在线观看| 国产不卡一卡二| 美女高潮喷水抽搐中文字幕| 夜夜爽天天搞| 成人三级做爰电影| 午夜久久久在线观看| 丰满人妻熟妇乱又伦精品不卡| 91大片在线观看| 国产99久久九九免费精品| 99久久99久久久精品蜜桃| 热re99久久精品国产66热6| 天天躁狠狠躁夜夜躁狠狠躁| av视频免费观看在线观看| 亚洲成av片中文字幕在线观看| 19禁男女啪啪无遮挡网站| 黄色视频,在线免费观看| 国产精品久久久久久人妻精品电影| 麻豆av在线久日| 69精品国产乱码久久久| 免费不卡黄色视频| 91麻豆av在线| 校园春色视频在线观看| 搡老乐熟女国产| xxxhd国产人妻xxx| 高清欧美精品videossex| 9191精品国产免费久久| √禁漫天堂资源中文www| 欧美老熟妇乱子伦牲交| 日韩精品青青久久久久久| 亚洲 欧美 日韩 在线 免费| 在线国产一区二区在线| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品 国内视频| 国产一卡二卡三卡精品| 亚洲av成人不卡在线观看播放网| 欧美日韩精品网址| 黑人操中国人逼视频| 欧洲精品卡2卡3卡4卡5卡区| 日本wwww免费看| cao死你这个sao货| 国产一区二区激情短视频| 久久狼人影院| 亚洲一区二区三区色噜噜 | 久久人妻av系列| 亚洲成a人片在线一区二区| 热99re8久久精品国产| 欧美一级毛片孕妇| 欧美最黄视频在线播放免费 | 一进一出抽搐gif免费好疼 | 别揉我奶头~嗯~啊~动态视频| 两人在一起打扑克的视频| 午夜久久久在线观看| 一级毛片高清免费大全| 国产乱人伦免费视频| 纯流量卡能插随身wifi吗| 免费少妇av软件| 亚洲精品粉嫩美女一区| 亚洲精品中文字幕一二三四区| 99riav亚洲国产免费| 狂野欧美激情性xxxx| 亚洲在线自拍视频| 亚洲午夜精品一区,二区,三区| 一进一出抽搐gif免费好疼 | av片东京热男人的天堂| 黄色丝袜av网址大全| 日韩精品中文字幕看吧| 97碰自拍视频| 欧美激情 高清一区二区三区| 久久久精品欧美日韩精品| 免费看十八禁软件| 国产单亲对白刺激| 啦啦啦在线免费观看视频4| 国产精品国产高清国产av| 男人舔女人的私密视频| 亚洲精品久久成人aⅴ小说| 琪琪午夜伦伦电影理论片6080| 母亲3免费完整高清在线观看| 在线观看日韩欧美| 亚洲精品国产一区二区精华液| 日日爽夜夜爽网站| 国产精品乱码一区二三区的特点 | 欧美亚洲日本最大视频资源| 国产xxxxx性猛交| 啦啦啦 在线观看视频| 午夜福利在线观看吧| 欧美中文日本在线观看视频| 最近最新中文字幕大全电影3 | 天天影视国产精品| 国产一区二区三区综合在线观看| 中国美女看黄片| 大码成人一级视频| 国产精品野战在线观看 | 久久精品国产亚洲av高清一级| 一级作爱视频免费观看| 欧美中文综合在线视频| 啦啦啦在线免费观看视频4| 欧美日韩精品网址| 国产精品久久久av美女十八| 国产黄a三级三级三级人| 人成视频在线观看免费观看| 欧美人与性动交α欧美精品济南到| 午夜成年电影在线免费观看| 国产伦一二天堂av在线观看| 亚洲三区欧美一区| 在线观看一区二区三区| videosex国产| 午夜精品久久久久久毛片777| 伦理电影免费视频| 欧美日韩瑟瑟在线播放| 97碰自拍视频| 国产精品国产av在线观看| 亚洲视频免费观看视频| 久久香蕉激情| 亚洲在线自拍视频| 欧美激情久久久久久爽电影 | 97人妻天天添夜夜摸| 亚洲精品一区av在线观看| 成人亚洲精品一区在线观看| 69精品国产乱码久久久| 久久久久精品国产欧美久久久| 亚洲国产欧美网| 三上悠亚av全集在线观看| 99精品欧美一区二区三区四区| 亚洲五月天丁香| 老司机福利观看| 欧美 亚洲 国产 日韩一| 久久狼人影院| 日韩av在线大香蕉| 免费久久久久久久精品成人欧美视频| 亚洲欧美一区二区三区久久| 90打野战视频偷拍视频| 欧美日韩视频精品一区| 亚洲欧洲精品一区二区精品久久久| 精品福利观看| 巨乳人妻的诱惑在线观看| 久久人妻福利社区极品人妻图片| 操出白浆在线播放| av国产精品久久久久影院| 丰满饥渴人妻一区二区三| 1024香蕉在线观看| 亚洲五月天丁香| 欧美另类亚洲清纯唯美| 欧美性长视频在线观看| 人人妻,人人澡人人爽秒播| 日韩三级视频一区二区三区| 一边摸一边抽搐一进一小说| 久久伊人香网站| 岛国在线观看网站| 啦啦啦免费观看视频1| 成人午夜高清在线视频| 国产69精品久久久久777片| 日日摸夜夜添夜夜添小说| 老熟妇仑乱视频hdxx| 九色成人免费人妻av| 成人性生交大片免费视频hd| 亚洲中文日韩欧美视频| 不卡一级毛片| 精品久久久久久久久av| 精品一区二区三区视频在线| 高潮久久久久久久久久久不卡| 国产视频一区二区在线看| 一进一出好大好爽视频| 国产精华一区二区三区| 国产大屁股一区二区在线视频| 亚洲经典国产精华液单 | 日韩有码中文字幕| 亚洲欧美日韩高清在线视频| 小说图片视频综合网站| 久久精品久久久久久噜噜老黄 | 内地一区二区视频在线| 亚洲欧美日韩东京热| 日本五十路高清| 国产精品自产拍在线观看55亚洲| 一进一出抽搐动态| 免费在线观看日本一区| 国产午夜福利久久久久久| 欧美黑人欧美精品刺激| 男插女下体视频免费在线播放| 尤物成人国产欧美一区二区三区| 亚洲专区国产一区二区| 中文字幕人成人乱码亚洲影| 欧美成人性av电影在线观看| 国产69精品久久久久777片| 赤兔流量卡办理| 网址你懂的国产日韩在线| 久久婷婷人人爽人人干人人爱| 深爱激情五月婷婷| 欧美黄色淫秽网站| 免费看光身美女| 一级a爱片免费观看的视频| 淫秽高清视频在线观看| 一进一出抽搐动态| 午夜免费男女啪啪视频观看 | 97热精品久久久久久| 尤物成人国产欧美一区二区三区| 国产探花极品一区二区| 精品久久久久久久久亚洲 | 免费搜索国产男女视频| 美女高潮的动态| 国内精品一区二区在线观看| www日本黄色视频网| 天美传媒精品一区二区| 欧美成人免费av一区二区三区| 欧美日韩国产亚洲二区| 色噜噜av男人的天堂激情| 欧美激情在线99| 直男gayav资源| 网址你懂的国产日韩在线| 高清日韩中文字幕在线| 国产成人影院久久av| 看免费av毛片| 成人毛片a级毛片在线播放| 欧美一区二区国产精品久久精品| 十八禁人妻一区二区| 国产色婷婷99| 中文字幕av成人在线电影| 琪琪午夜伦伦电影理论片6080| 可以在线观看的亚洲视频| 小蜜桃在线观看免费完整版高清| 久久欧美精品欧美久久欧美| 成人特级av手机在线观看| 极品教师在线视频| h日本视频在线播放| 国产蜜桃级精品一区二区三区| 久久久久国内视频| 国内精品久久久久精免费| 成人毛片a级毛片在线播放| 国产伦人伦偷精品视频| 亚洲欧美精品综合久久99| 亚洲精品一卡2卡三卡4卡5卡| 国产成年人精品一区二区| 亚洲欧美日韩高清专用| 国产在线男女| 3wmmmm亚洲av在线观看| 午夜老司机福利剧场| 久久久久免费精品人妻一区二区| 男女做爰动态图高潮gif福利片| 亚洲五月婷婷丁香| 看免费av毛片| 中文字幕av成人在线电影| 嫩草影院新地址| 国产激情偷乱视频一区二区| 91av网一区二区| 久久久色成人| 免费看光身美女| 国产精品,欧美在线| 久久久国产成人免费| 熟妇人妻久久中文字幕3abv| 一进一出好大好爽视频| 能在线免费观看的黄片| 少妇熟女aⅴ在线视频| 亚洲,欧美精品.| 国产69精品久久久久777片| 亚洲专区国产一区二区| 国产精品三级大全| 国产v大片淫在线免费观看| 欧美日韩综合久久久久久 | 99在线人妻在线中文字幕| 无人区码免费观看不卡| 亚洲欧美日韩东京热| 一进一出好大好爽视频| 51国产日韩欧美| 看十八女毛片水多多多| 国产精品久久久久久精品电影| 午夜激情福利司机影院| 久久久久国产精品人妻aⅴ院| 他把我摸到了高潮在线观看| 在线免费观看不下载黄p国产 | a级一级毛片免费在线观看| 少妇被粗大猛烈的视频| 国产成人a区在线观看| 国产探花在线观看一区二区| 一区二区三区激情视频| 久久99热这里只有精品18| 少妇的逼好多水| 天美传媒精品一区二区| 久9热在线精品视频| 亚洲天堂国产精品一区在线| 国产午夜精品论理片| av国产免费在线观看| 精品一区二区三区视频在线| 精品久久久久久久人妻蜜臀av| 中文字幕人成人乱码亚洲影| 亚洲欧美日韩高清专用| 特级一级黄色大片| 美女高潮喷水抽搐中文字幕| 成人鲁丝片一二三区免费| 动漫黄色视频在线观看| 一夜夜www| 小蜜桃在线观看免费完整版高清| 国产真实伦视频高清在线观看 | 色精品久久人妻99蜜桃| netflix在线观看网站| 国产91精品成人一区二区三区| 中文字幕精品亚洲无线码一区| 三级男女做爰猛烈吃奶摸视频| 欧美性猛交╳xxx乱大交人| 最近最新中文字幕大全电影3| 亚洲国产精品999在线| 午夜精品在线福利| 91九色精品人成在线观看| 日韩欧美一区二区三区在线观看| 国产人妻一区二区三区在| 一个人看的www免费观看视频| 精品久久久久久久久av| 在线播放国产精品三级| 久久国产乱子免费精品| bbb黄色大片| 国产精品久久久久久久电影| 欧美zozozo另类| 免费av观看视频| 97超视频在线观看视频| 老司机福利观看| 午夜福利在线在线| 久久久久国内视频| 久久久久免费精品人妻一区二区| 国产高清视频在线观看网站| 亚洲美女搞黄在线观看 | 男人和女人高潮做爰伦理| 少妇熟女aⅴ在线视频| 九九热线精品视视频播放| 亚洲av成人精品一区久久| 国产高清激情床上av| 两人在一起打扑克的视频| 一区二区三区高清视频在线| 国产成人影院久久av| 少妇的逼水好多| 亚洲av成人av| 久久伊人香网站| 亚洲欧美日韩无卡精品| 一进一出抽搐动态| 亚洲av中文字字幕乱码综合| 我的女老师完整版在线观看| 亚洲精品456在线播放app | 婷婷精品国产亚洲av| 欧美性猛交╳xxx乱大交人| 亚洲人成网站高清观看| 久久精品久久久久久噜噜老黄 | 亚洲第一区二区三区不卡| 久久6这里有精品| 成人精品一区二区免费| 99在线人妻在线中文字幕| 国产成人啪精品午夜网站| av在线观看视频网站免费| 男插女下体视频免费在线播放| 99国产精品一区二区三区| 国产黄a三级三级三级人| 乱码一卡2卡4卡精品| 婷婷色综合大香蕉| 狂野欧美白嫩少妇大欣赏| 国产午夜福利久久久久久| 18禁裸乳无遮挡免费网站照片| 欧美精品啪啪一区二区三区| 噜噜噜噜噜久久久久久91| 成人欧美大片| 亚州av有码| 国产aⅴ精品一区二区三区波| 日本 av在线| 国产精品伦人一区二区| av天堂在线播放| 最近在线观看免费完整版| 成人毛片a级毛片在线播放| bbb黄色大片| 女同久久另类99精品国产91| 日韩高清综合在线| 亚洲中文字幕日韩| 夜夜爽天天搞| 精品福利观看| 99久久无色码亚洲精品果冻| 日韩国内少妇激情av| 观看美女的网站| 9191精品国产免费久久| 极品教师在线免费播放| 亚洲中文日韩欧美视频| 精品久久国产蜜桃| 日日摸夜夜添夜夜添小说| 亚洲内射少妇av| 成人欧美大片| 国产一区二区亚洲精品在线观看| 天堂网av新在线| 精品午夜福利在线看| 国产色爽女视频免费观看| 国产蜜桃级精品一区二区三区| 欧美成人免费av一区二区三区| 搡老岳熟女国产| 91麻豆av在线| www日本黄色视频网| 婷婷丁香在线五月| 91av网一区二区| 亚洲欧美精品综合久久99| 美女大奶头视频| 国产v大片淫在线免费观看| 国产精品影院久久| 午夜精品久久久久久毛片777| 免费在线观看影片大全网站| 亚洲国产色片| 91字幕亚洲| 国产日本99.免费观看| 人人妻人人澡欧美一区二区| 日本黄色片子视频| 美女大奶头视频| 99热这里只有精品一区| 久久精品国产亚洲av香蕉五月| 性欧美人与动物交配| 精华霜和精华液先用哪个| 国产v大片淫在线免费观看| 少妇丰满av| 欧美国产日韩亚洲一区| 国产精品一及| 精品国内亚洲2022精品成人| 国产一区二区在线av高清观看| 亚洲三级黄色毛片| 有码 亚洲区| 在线国产一区二区在线| 国内毛片毛片毛片毛片毛片| 久久久久免费精品人妻一区二区| 日本a在线网址| 三级男女做爰猛烈吃奶摸视频| 午夜免费激情av| 精品国产亚洲在线| 熟女人妻精品中文字幕| 久久久久亚洲av毛片大全| 免费看日本二区| 一区福利在线观看| 久久人妻av系列| 国产高清视频在线观看网站| 成人一区二区视频在线观看| 色综合站精品国产| 国产精品人妻久久久久久| 搡老妇女老女人老熟妇| 激情在线观看视频在线高清| 日韩高清综合在线|