Research article

Nonlocal delay gives rise to vegetation patterns in a vegetation-sand model

  • The vegetation pattern generated by aeolian sand movements is a typical type of vegetation patterns in arid and semi-arid areas. This paper presents a vegetation-sand model with nonlocal interaction characterized by an integral term with a kernel function. The instability of the Turing pattern was analyzed and the conditions of stable pattern occurrence were obtained. At the same time, the multiple scales method was applied to obtain the amplitude equations at the critical value of Turing bifurcation. The spatial distributions of vegetation under different delays were obtained by numerical simulation. The results revealed that the vegetation biomass increased as the interaction intensity decreased or as the nonlocal interaction distance increased. We demonstrated that the nonlocal interaction between vegetation and sand is a crucial mechanism for forming vegetation patterns, which provides a theoretical basis for preserving and restoring vegetation.

    Citation: Jichun Li, Gaihui Guo, Hailong Yuan. Nonlocal delay gives rise to vegetation patterns in a vegetation-sand model[J]. Mathematical Biosciences and Engineering, 2024, 21(3): 4521-4553. doi: 10.3934/mbe.2024200

    Related Papers:

    [1] Cuicui Cai, Maosheng Fu, Xianmeng Meng, Chaochuan Jia, Mingjing Pei . Indoor high-precision visible light positioning system using Jaya algorithm. Mathematical Biosciences and Engineering, 2023, 20(6): 10358-10375. doi: 10.3934/mbe.2023454
    [2] Cheng Zhu, Yong Zhang, Xuhua Pan, Qi Chen, Qingyu Fu . Improved Harris Hawks Optimization algorithm based on quantum correction and Nelder-Mead simplex method. Mathematical Biosciences and Engineering, 2022, 19(8): 7606-7648. doi: 10.3934/mbe.2022358
    [3] Shuang Wang, Heming Jia, Qingxin Liu, Rong Zheng . An improved hybrid Aquila Optimizer and Harris Hawks Optimization for global optimization. Mathematical Biosciences and Engineering, 2021, 18(6): 7076-7109. doi: 10.3934/mbe.2021352
    [4] Yuheng Wang, Yongquan Zhou, Qifang Luo . Parameter optimization of shared electric vehicle dispatching model using discrete Harris hawks optimization. Mathematical Biosciences and Engineering, 2022, 19(7): 7284-7313. doi: 10.3934/mbe.2022344
    [5] Shenghan Li, Linlin Ye . Multi-level thresholding image segmentation for rubber tree secant using improved Otsu's method and snake optimizer. Mathematical Biosciences and Engineering, 2023, 20(6): 9645-9669. doi: 10.3934/mbe.2023423
    [6] Wang Cai, Jianzhuang Wang, Longchao Cao, Gaoyang Mi, Leshi Shu, Qi Zhou, Ping Jiang . Predicting the weld width from high-speed successive images of the weld zone using different machine learning algorithms during laser welding. Mathematical Biosciences and Engineering, 2019, 16(5): 5595-5612. doi: 10.3934/mbe.2019278
    [7] Hong Gao, Cuiyun Wu, Dunnian Huang, Dahui Zha, Cuiping Zhou . Prediction of fetal weight based on back propagation neural network optimized by genetic algorithm. Mathematical Biosciences and Engineering, 2021, 18(4): 4402-4410. doi: 10.3934/mbe.2021222
    [8] Feng Li, Zhongao Bi, Hongzeng Xu, Yunqi Shi, Na Duan, Zhaoyu Li . Design and implementation of a smart Internet of Things chest pain center based on deep learning. Mathematical Biosciences and Engineering, 2023, 20(10): 18987-19011. doi: 10.3934/mbe.2023840
    [9] Chun-Chao Yeh, Ke-Jia Jhang, Chin-Chun Chang . An intelligent indoor positioning system based on pedestrian directional signage object detection: a case study of Taipei Main Station. Mathematical Biosciences and Engineering, 2020, 17(1): 266-285. doi: 10.3934/mbe.2020015
    [10] Xinli Hu . Threshold dynamics for a Tuberculosis model with seasonality. Mathematical Biosciences and Engineering, 2012, 9(1): 111-122. doi: 10.3934/mbe.2012.9.111
  • The vegetation pattern generated by aeolian sand movements is a typical type of vegetation patterns in arid and semi-arid areas. This paper presents a vegetation-sand model with nonlocal interaction characterized by an integral term with a kernel function. The instability of the Turing pattern was analyzed and the conditions of stable pattern occurrence were obtained. At the same time, the multiple scales method was applied to obtain the amplitude equations at the critical value of Turing bifurcation. The spatial distributions of vegetation under different delays were obtained by numerical simulation. The results revealed that the vegetation biomass increased as the interaction intensity decreased or as the nonlocal interaction distance increased. We demonstrated that the nonlocal interaction between vegetation and sand is a crucial mechanism for forming vegetation patterns, which provides a theoretical basis for preserving and restoring vegetation.



    With the rapid development of innovations in communication technology and the notoriety of remote systems, the desire for enhanced indoor positioning continues to expand [1]. However, when facing complex indoor environments, such as with the real-time positioning of personnel, equipment and materials in office buildings, factories and construction sites [2], real-time care and management of the elderly [3], precise electronic positioning guidance inside hospitals [4] and precise location trajectories of supervised persons in prisons [5], there is an expanding requirement for high-precision situating innovation. Thus, strategies for advancing the exactness of indoor situating tasks have become of increasing interest to a range of researchers.

    Indoor positioning technology is developing rapidly today and is classified as building-dependent or building-independent, depending on the sensors used for positioning [6]. The most common building-independent indoor positioning techniques are pedestrian dead reckoning (PDR) systems [7] and image-based techniques that use cameras [8]. Although the noise is small and stable in the short term when using PDR systems, such approaches cannot be used for long times due to the accumulation of errors, and toward solving this problem, the mainstream research direction is currently the fusion algorithm. Various fusion algorithms are proposed in [9,10,11], such as fusion algorithms with Wi-Fi, Bluetooth, RSS fingerprinting, etc., which effectively reduce positioning error and improve robustness. The four aspects of a typical inertial measurement unit(IMU)-based PDR system framework are direction updating, step detection, step estimation and map matching [12]. Navigation systems include techniques such as orientation updating and map matching, whereas step detection and estimation are specialized to human motion since they exploit the constraint of pedestrian steps. The camera is used for positioning in image-based indoor positioning tasks. The camera used in such applications can be monocular [13], stereo [14] or RGB-D [15]. Because it is difficult to obtain accurate positioning results with a single method, a hybrid positioning method in such a positioning system is introduced as the reliable option. The existing techniques for IMU-based positioning, camera-based positioning and various existing hybrid methods for indoor positioning are discussed in [16], and a hybrid indoor positioning system that uses IMU sensors and smartphone cameras is proposed, which effectively improves indoor positioning accuracy.

    Building-independent positioning technologies, of which there are many types, fall into two main categories: namely, those that utilize the existing building infrastructure and those that require dedicated infrastructure. The former category includes Wi-Fi, cellular and Bluetooth. Most Wi-Fi-based positioning systems can be used as RSS-based measurement systems. Various Wi-Fi restrictions must also be considered: for example, signal attenuation in static contexts such as walls and floors, as well as changes in interior environments [17]. Although the cellular network has a wide coverage area, it still suffers from low positioning accuracy. One of the disadvantages of employing Bluetooth technology for localization is that it performs the device discovery operation for each location finding event; as a result, the localization latency (10-30 s) and consumption are greatly increased. As a result, the latency of a Bluetooth device is inadequate for real-time positioning applications. However, for the previously described Wi-Fi localization, a hybrid deep learning model (HDLM)-based indoor localization system was developed in [18]. Instead of raw RSSI signals from APs, the proposed HDLM-based localization system uses RSSI heatmaps. As a result, the Wi-Fi-RSSI signal-based localization system's localization performance is improved. For the latter, the main positioning technologies are RFID, UWB, infrared, ultrasonic, ZigBee, VLC, acoustic, etc. RFID technology requires a large number of tags, which do not meet the requirements of real-time positioning, and its positioning accuracy is related to the deployment density of tags [19]. Indoor positioning systems based on infrared possess several drawbacks, such as security and privacy concerns. The position determination of infrared signals is limited due to interference from fluorescent light and sunlight. In addition, the hardware and maintenance expenses of IR-based indoor systems are high [20]. Reflected waves propagated by other sources can cause severe interference for ultrasound devices (e.g., metal collisions) [21]. The ZigBee positioning accuracy is low, as is its data transfer rate [22].

    UWB is a remote correspondence innovation that transfers a bandwidth of more than 1 GHz. Due to its strong anti-interference performance, resistance to multipath interference, high transmission rate, low transmission power and increased penetration capability, UWB positioning technology is widely used for positioning in indoor environments. UWB can provide high positioning accuracy [23,24,25]. Currently, the various UWB ranging methods can be classified into time of arrival (TOA), received signal strength (RSS) and time difference of arrival (TDOA) [26]. However, in complex indoor environments, where walls or other signals block objects, the signal transmission distance increases, and the transmission time lengthens, significantly reducing UWB positioning accuracy. A single indoor positioning method can no longer meet the growing requirements for localization accuracy. Therefore, mixed localization algorithms are becoming increasingly attractive [27].

    The localization equation is usually a nonlinear equation based on TOA or TDOA; as a result, the equation must be converted into a linear equation for solving [28]. Currently, nonlinear problems are mainly solved by the least-squares estimation (LSE) algorithm, Fang algorithm, Chan algorithm and Taylor algorithm. The weighting matrix of the LSE algorithm should be decided according to the precision of the measurement data; in practical applications, the selection of the weighting matrix is difficult. The best estimation performance of the algorithms developed by Chan and Fang is achieved when the time difference measurement error is small, although the performance of both algorithms is seriously affected by the large time difference measurement error in the real channel environment [29]. In recent years, deep learning approaches have been commonly employed for indoor positioning because of their excellent nonlinear mapping capabilities and self-adaptive and self-learning abilities [30]. For example, Zhuo et al. suggested a UWB locating technique based on BPNN that is applied to automatic parking, showing good real-time performance and accuracy [31]. In [32], a long short-term memory (LSTM) network is proposed to reduce the position error in UWB systems; this network solves the LOS problem associated with UWB systems and improves the positioning accuracy. Shuai et al. applied the 60 GHz pulsed-IR algorithm to a BPNN, which was used to improve the ability of conventional algorithms to accurately estimate signal fading due to the complexity of indoor channels, thus improving localization accuracy [33]. Lian et al. proposed an optimized K-means clustering BPNN UWB positioning approach, which enhances the original data in the UWB positioning system and preprocesses the K-means clustering algorithm to filter the data with a significant positional deviation, and the experimental measurement accuracy can reach 26 cm with good stability [34]. Zhang and Jin proposed a BP neural network indoor positioning algorithm based on VLC. Through simulation experiments, the authors showed that the algorithm has a lower relative positioning error and better indoor positioning accuracy than the MDS-MAP and MDS-MAP(P) algorithms [35].

    However, the BPNN is sensitive to random weights and thresholds, which enhances the likelihood of being trapped in nearby minima, thus generating significant localization errors [36]. Therefore, numerous optimization algorithms have been proposed to reduce the downsides of the BPNN. Among them, population intelligence optimization algorithms are the most widely applied algorithms toward solving the inherent defects of traditional training algorithms [37]. For example, a location fusion method of genetic algorithm (GA) optimization of the BP neural network that effectively reduces the localization error was proposed by Yang et al. [38]. Li et al. [39] proposed a mixed model that combines the cuckoo search optimization algorithm (CS) with the BPNN, where the random weights and thresholds of the BPNN are optimally trained, and the measured time differences are corrected after training. Finally, the Chan algorithm is used for localization calculation according to the corrected time difference values. In addition to these studies, the current swarm intelligence algorithm is developing rapidly, and better results have been demonstrated [40,41,42,43,44]. Previous literature has investigated indoor localization methods using the particle swarm optimization algorithm and BPNN hybrid algorithm (PSO-BP) combined with RFID using artificial intelligence techniques [45] and UWB indoor localization methods using the immune algorithm and BPNN hybrid algorithm (IA-BP) [46]. The above research shows that artificial intelligence techniques show strong capability in UWB positioning. It is essential to study new artificial intelligence techniques and prediction models to improve UWB positioning accuracy [47].

    As indicated by the no free lunch (NFL) theorem [48], theoretically, a certain algorithm cannot be treated as a generic global optimal optimization algorithm. Therefore, although a range of algorithms is used to optimize UWB indoor positioning, the presented algorithms are far from being able to solve all problems or be applied to all indoor environments. To the authors' best knowledge, UWB indoor localization has not yet been attempted using a combination of the HHO algorithm and BPNN.

    This study proposes a new hybrid algorithm, called HHO-BP, to improve UWB indoor positioning accuracy and resistance to non-line-of-sight ranges. The algorithm uses the HHO algorithm to optimize the BPNN to find a set of optimal weights and thresholds, thereby improving the positioning accuracy. Finally, the optimization performance of the HHO-BP algorithm and the conventional BPNN are compared and analysed through experimental validation.

    The major contribution of this paper is to propose a new hybrid algorithm, HHO-BP, for optimizing ultrawideband indoor positioning, outlined as follows:

    (1) A new algorithm for hybrid localization, namely, HHO-BP, is proposed.

    (2) Combining the BP neural network with an intelligent population algorithm, the structure of the BP neural network is optimized, thus solving the local optimal solution problem.

    (3) The algorithm greatly improves the accuracy and anti-jamming capability of ultrawideband indoor positioning, with centimetre-level accuracy for both line-of-sight and non-line-of-sight positioning.

    The structure of the paper is as follows: First, the principles of UWB localization and BPNN are explained in Section 2. Section 3 introduces the principle of the proposed swarm intelligence algorithm (HHO), as well as the constructed new model. Section 4 offers the findings as well as a discussion and analysis. Finally, in Section 5, the study's findings are presented.

    The principle of the trilateral positioning algorithm is as follows: The three base stations are located at A0(x1, y1), A1(x2, y2) and A2(x3, y3), and it is known that the distance values of the tag to be measured, Tn(x0, y0), to the 3 base stations are d1, d2 and d3, as displayed in Figure 1.

    Figure 1.  Trilateral positioning algorithm.

    In the case of two-dimensional positioning, according to the Pythagorean theorem, the formula for calculating the position of the tag, Tn, and the focus are attained using Eq (1).

    d1=(x1x0)2+(y1y0)2d2=(x2x0)2+(y2y0)2d3=(x3x0)2+(y3y0)2 (1)

    Under ideal circumstances, the distance values of tag Tn to base Stations A0, A1 and A2 are accurate; tag Tn has only one correct solution, as shown by the red dot in Figure 2.

    Figure 2.  Ideal circumstances.

    Under actual circumstances, real distance measurement is subject to error; the reasons for the error include the error of the range itself and the error of the actual scene, such as obstruction or metal, which affects the penetration and transmission of electromagnetic waves. In a line-of-sight environment, the measured value fluctuates above and below the actual value; however, in a non-line-of-sight environment, the distance measured is theoretically higher than the actual value due to the presence of occlusion between the base station and the tag, resulting in the effect shown below, where the solved Tn coordinate is no longer a point. Nevertheless, it is a region, so we need to find an optimal solution in the red area, as shown in Figure 3.

    Figure 3.  Actual circumstances.

    There are many models of artificial neural networks, but the foremost broadly utilized, natural and simple-to-understand is the backpropagation algorithm, called BPNN [49]. Training is performed through signal forward propagation and error backward propagation.

    The network generally has a multilayer structure, including an input layer for data input, an output layer for result output and several hidden layers, and each hidden layer contains one or more neurons. The topological diagram of the BPNN is shown in Figure 4.

    Figure 4.  The topology of BPNN.

    In the structure, x1, x2, …, xm are inputs, and y1, y2, …, yn are outputs. A hidden layer contains s neurons, the output function of the hidden layer is bj, and the threshold values for the hidden and output layers are θj and θk, respectively. The most frequently used activation function is the sigmoid function, which is set to S(x) here. f1 and f2 are the transfer functions of the hidden and output layers, respectively. The output of the model is set to yk, and the expected output is tk [50]. The above can be modelled as follows.

    The output of the Jth neuron of the hidden layer is [50]

    bj=f1(mi=1wijxiθj) (i=1,2,...,m;j=1,2,...,s). (2)

    a. Calculate the output yk of the output layer as in Eq (3) [50]:

    yk=f2(si=1wjkbjθk) (j=1,2,...,s;k=1,2,...,n). (3)

    b. Define the error function by the actual output of the network as in Eq (4) [50]:

    e=nk=1(tkyk)2. (4)

    c. The sigmoid activation function is given as Eq (5) [50]:

    S(x)=11+ex. (5)

    A new bionic intelligent optimization approach, the Harris hawks optimization (HHO) algorithm, was developed in 2019. The method is a metaheuristic approach that combines Lévy flights to solve large multidimensional problems by simulating the collaborative group behaviour of Harris hawks while hunting. Compared with previous algorithms, this algorithm has a higher search capability, acceptable development and exploration and thus has been successfully applied to many practical engineering problems [51,52].

    As a result, since BPNN is more sensitive to initial weights and thresholds, this research proposes the Harris hawks optimization method for BPNN random weights and threshold optimization.

    The whole optimization process includes a global exploration phase, a transition phase and a local exploitation phase [52].

    1) Exploration phase:

    In HHO, the Harris hawks perch randomly on some locations in a waiting form, monitoring the area [lb, ub] to search for prey. The prey is simultaneously searched at random locations according to two different strategies, and the location update is performed with probability q during the iteration. The mathematical formula for the Harris hawks's location, which is closest to the prey, is the following:

    X(t+1)={Xrand(t)r1|Xrand(t)2r2X(t)|,q0.5(Xrabbit(t)Xb(t))r3(lb+r4(ublb)),q<0.5. (6)

    In Eq (6), X(t+1) is the hawk's location in the following iteration, t. The location of the rabbit is Xrabbit(t), and the current position of the hawk is X(t). r1, r2, r3, r4 and q are all random numbers ranging from 0 to 1. The upper and lower bounds of the variables are represented by ub and lb, respectively, while Xrand(t) is the hawk chosen at random from the population, which is updated in each iteration. Xb(t) is the current hawk average location in Eq (7).

    Xb(t)=1NNi=1Xi(t)

    where N is the number of hawk, and in iteration t, Xi(t) shows the location of each hawk.

    2) Transition phase:

    The HHO algorithm moves from exploration to exploitation and then between multiple exploitation behaviours depending on the energy, E, of prey escape, as shown below in Eq (8):

    E=2E0(1tT)

    where E0 is the energy starting stage, and T is the highest number of iterations [52].

    3) Exploitation phase:

    Four distinct tactics are employed to recreate the assault in this phase, depending on how the prey flees and how the Harris hawk pursues it, combined with Lévy flights, representing the prey's chance of escaping [52].

    a)

    r0.5&|E|0.5

    In this case, the rabbit still has enough energy, and the hawk aim to exhaust its energy before launching the assault. The modelling rules given in Eqs (9) and (10) [52]:

    X(t+1)=ΔX(t)E|JXrabbit(t)X(t)|,
    ΔX(t)=Xrabbit(t)X(t).

    To replicate the nature of the rabbit's movements, the J value varies randomly in each iteration, and the rabbit's random hop intensity throughout the escape phase is denoted by J=2(1r5). r5 is a number in the range from 0 to 1. The disparities between the rabbit's prior location vector and the current position vector are ΔX(t).

    b)

    r0.5&|E|<0.5

    Under this condition, the rabbit is fatigued, and the hawk can launch a surprise assault. Equation (11) expresses how the model updates the rabbit's position [52]:

    X(t+1)=Xrabbit(t)E|ΔX(t)|. (11)

    c)

    r<0.5&|E|0.5

    The rabbit has enough energy in this situation to effectively flee while still constructing a soft envelope before the attack. The following rule, given in Eq (12), can be used to update the Harris hawks' positions [52]:

    Y=Xrabbit(t)E|JXrabbit(t)X(t)|. (12)

    The hawk will dive utilizing the following rule, Eq (13), using LF-based patterns:

    Z=Y+S×LF(D)

    where D is the dimension, S is a one-dimensional random vector, and LF is the Levy flight function and is based on the equation below (Eq (14)).

    LF(x)=0.01×μ×σ|ν|1β,σ=(Γ(1+β)×sin(πβ2)Γ(1+β2)×β×2(β12))1β (14)

    Here, β is 1.5, and μ and ν are random numbers between 0 and 1.

    By calculation, Eq (15) can execute the ultimate plan to update the locations of the hawks:

    X(t+1)={Y F(Y) < F(H(t))Z F(Z) < F(H(t)). (15)

    d)

    r<0.5&|E|<0.5

    In this situation, the rabbit is unable to flee due to a lack of energy, so the hawks set up a hard siege. The modelling rules are given in Eqs (16) and (17) [52]:

    =Xrabbit(t)E|JXrabbit(t)Xb(t)|, (16)
    Z =Y+S×LF(D). (17)

    Here, Y references Eq (15).

    Figure 5.  A method for creating a Harris hawk optimization-BP neural network (HHO-BP) model is described.

    Figure 5 depicts the development path of the HHO-BP model. First, the obtained database is divided into training and test datasets. The BP neural network is first set to random weights and thresholds, and then the training set data are input for training, where the Harris hawk optimization algorithm develops the final BP model by searching for the best combination of weights and thresholds in the BP model in a preset number of iterations. Finally, the test set data are fed into the final constructed BP neural network for predicting the UWB localization results. Performance metrics such as the mean error and position comparison can be used to evaluate the optimized BP model, and the actual and measured coordinates of the test set and the prediction results are used to test its performance and calculate the prediction error.

    The computational complexity of HHO in this approach is determined using three key processes: initial optimization, fitness evaluation and candidate update. The computing complexity of the initialization process is O(N) for all N. The update mechanism has a computational complexity of O(T×N)+O(T×N×D), which comprises updating the search agent location and finding the optimal solution, where T is the maximum number of iterations, set to 200 to assure convergence, and D is the problem's dimensionality. As a result, HHO's computational complexity is O(N×(T+TD+1)) [52].

    Through the following site environment and hardware equipment, UWB measurement data are acquired. Then, they are compared with the actual position, and visualized positioning points and error curves are plotted to see if the algorithm pair is effective and accurate.

    To test the accuracy of UWB localization and the algorithm's efficiency, a combination of offline data acquisition and MATLAB simulation is adopted. For the acquisition, an LOS environment and an NLOS environment were created, and the experimental sites are shown in Figures 6 and 7.

    LOS environment: To create an indoor positioning scenario, four UWB anchors with heights of 2.2 m and positions of (0, 0), (0, 6.45), (5.58, 0) and (5.58, 6.45) were placed in the region illustrated in Figure 6. UWB anchor 0 is connected via Wi-Fi to a computer with the developed interactive interface installed, and the real-time position of the UWB tag is obtained through the host computer software to obtain indoor positioning data for the LOS environment.

    NLOS environment: Because of the indoor environment's complexity, the height of the four UWB anchors was reduced to 1.4 m, and the positioning test points were built around the table using the shelter of the table to build a non-line-of-sight testing environment, as shown in Figure 7.

    The DWM1000 UWB module used in the range circuit is made by DECAWAVE. The Wi-Fi is used to communicate with the PC to obtain UWB measurement data. Figure 8 depicts the system structure and hardware circuit diagram [53].

    Figure 6.  The UWB test field has a size of 6.45 × 5.58 m2. The height of the UWB anchors is 2.2 m. The coordinate's origin is set to A0. (a) LOS distance actual view. (b) LOS distance plan and measurement trajectory.
    Figure 7.  NLOS experimental scene diagram. The height of the UWB anchors is 1.4 m. (a) NLOS distance actual view. (b) NLOS distance plan and measurement trajectory.
    Figure 8.  The UWB ranging system. (a) System structure diagram. (b) Hardware circuit diagram.

    The positioning data were collected separately at the same height at the site. Four label heights were tested (1.3, 1.55, 1.7 and 1.95 m) to ensure the feasibility of experimental rigor. The actual locations and the measured locations were used to form the sample set, which was divided into the training and test sets.

    In Figure 9, the four panels a-d represent the total dataset of the four altitude UWB measurements and the actual coordinate positions. The black asterisks represent the actual coordinate positions measured with the metric ruler; the red diamonds represent the upper computer positioning software data. The data in the graph are all shifted by different magnitudes. To verify the effectiveness and accuracy of the algorithm, the simulation results were subsequently compared with the redraw data. Additionally, to better analyse the error of the original UWB positioning, the error curve is plotted using MATLAB 2020 software, as shown in Figure 10.

    As seen from Figure 10, for different altitudes, there are positioning anomalies, and some points do not achieve centimetre-level positioning accuracy. In contrast, some positions are more volatile in comparison. In other words, conventional UWB positioning either has large fluctuations or large deviations. For this purpose, from the aforementioned sample set, a training dataset and a testing dataset are constructed, which are fed into the BPNN and the HHO-BP neural network, which predicts the localization sites. Finally, the error analysis is performed with the actual coordinate points. To analyse the effect of the algorithm, we choose nonidentical locus points as the test set.

    The black points in Figure 11 are the actual locations (i.e., the metric measurement data). The red points are the original UWB positioning measurement data. As the figure shows, the initial UWB positioning is not ideal, and there are significant errors.

    In Figure 12, the black points are the actual coordinates (measured in metres), the blue triangles are the unoptimized BP neural network (random weight threshold) localization points, and the red circles are the HHO-BP neural network localization points.

    When comparing Figures 11 and 12, the findings reveal that both the BPNN and the HHO-BP have greatly increased accuracy. To more accurately evaluate the difference in effectiveness between the BPNN and the optimized HHO-BP, each test point was compared with the actual location, and three sets of error data were obtained. The initial UWB measurement error, the BPNN prediction result error, and the HHO-BP prediction result error were all identified, as demonstrated in Figures 13 and 14.

    In Figure 13, the black line segment indicates the original UWB error, the blue line segment shows the BP neural network optimization error, and the red line segment indicates the HHO-BP neural network optimization error.

    The experimental findings in Figures 13 and 14 reveal that the HHO-BP algorithm has a superior optimization impact relative to the BP neural algorithm at all four heights of the LOS environment, which effectively reduces the error of UWB indoor positioning and achieves centimetre-level accuracy. At the same time, HHO-BP can effectively decrease outliers and reduce accuracy variation to a certain amount.

    Figure 9.  The LOS environment has four total altitude test positioning points and actual coordinate points. (a) 1.3 m. (b) 1.55 m. (c) 1.7 m. (d) 1.95 m.
    Figure 10.  The original error between the actual position and the test position for each point at four different heights. (a) 1.3 m. (b) 1.55 m. (c) 1.7 m. (d) 1.95 m.
    Figure 11.  The test set actual position and measured position. (a) 1.3 m. (b) 1.55 m. (c) 1.7 m. (d) 1.95 m.
    Figure 12.  The test set actual position and optimized points. (a) 1.3 m. (b) 1.55 m. (c) 1.7 m. (d) 1.95 m.
    Figure 13.  The test set original error and optimized error. (a) 1.3 m. (b) 1.55 m. (c) 1.7 m. (d) 1.95 m.
    Figure 14.  Test set average error and optimized error. (a) 1.3 m. (b) 1.55 m. (c) 1.7 m. (d) 1.95 m.

    In this environment, to achieve the NLOS environmental conditions, four UWB anchors were adjusted to a height of 1.4 m, and the test points were located around the table with the table positioned as an occluder. The localization coordinates of the UWB tag heights of 0.24, 0.4, 0.55 and 0.67 m were measured separately as the sample set for this environment. From the aforementioned sample set, a training dataset and a testing dataset were constructed.

    From Figure 15, from sample sets of different heights, it can be found that in the NLOS environment, the outliers of sample points increase significantly. Because of the shielding effect of the table, the signal's arrival time increases, such that the measured coordinates deviate from the shielding object. As the height of the UWB tag decreases gradually, the complexity of the NLOS environment increases, and the location deviation becomes more chaotic, causing the optimization task of this environment to be more difficult.

    Figure 16 shows the error curves of the sample set for each altitude measurement versus the actual coordinates, indicating a significant increase in error in the NLOS environment, along with greater volatility. Figure 16 shows that UWB positioning has a large error in the NLOS environment, with some test points even showing an error of more than 1.5 m. Therefore, optimization for this situation is essential.

    Figure 15.  NLOS environment with four total altitude test positioning points and actual coordinate points. (a). 0.24 m. (b). 0.4 m. (c). 0.55 m. (d). 0.67 m.
    Figure 16.  The original error between the actual position and the test position for each point at four different heights. (a). 0.24 m. (b). 0.4 m. (c). 0.55 m. (d). 0.67 m.
    Figure 17.  The training set actual position and measured position. (a). 0.24 m. (b). 0.4 m. (c). 0.55 m. (d). 0.67 m.
    Figure 18.  The test set actual position and measured position. (a). 0.24 m. (b). 0.4 m. (c). 0.55 m. (d). 0.67 m.
    Figure 19.  The test set original error and optimized error. (a). 0.24 m. (b). 0.4 m. (c). 0.55 m. (d). 0.67 m.
    Figure 20.  The test set average error and optimized error. (a). 0.24 m. (b). 0.4 m. (c). 0.55 m. (d). 0.67 m.

    Comparing Figures 17 and 18, the results reveal that under the NLOS setting, the accuracy of the BPNN and the HHO-BP has improved significantly. To more accurately evaluate the difference in effectiveness between the BPNN and the optimized HHO-BP, the UWB measurement error is compared to the error of each point in the testing dataset, as shown in Figures 19 and 20.

    Figures 19 and 20 show that HHO-BP effectively reduces errors and improves positioning accuracy at all altitudes in the complex NLOS environment.

    Table 1.  Results of the algorithm.
    BP HHO-BP
    Environment H Max E Min E ME Max E Min E ME
    LOS 130 17.63 0.97 8.63 11.78 0.21 5.31
    155 14.77 2.99 7.39 10.38 0.50 4.94
    170 9.22 2.68 5.15 7.72 1.74 4.18
    195 11.14 4.77 4.89 6.78 2.05 3.56
    NLOS 24 48.60 4.25 18.19 28.82 3.46 13.75
    40 30.12 3.45 11.65 19.57 1.47 6.65
    55 71.08 3.90 20.41 61.10 2.24 12.45
    67 19.35 1.53 9.04 13.25 0.87 5.41
    Note: H: Height, Max E: Maximum Error, Min E: Minimum Error, ME: Mean Error.

     | Show Table
    DownLoad: CSV

    Table 1 indicates the maximum errors, minimum errors and mean errors of all test data. These details are obtained from the graph above. Here, the mean error is counted separately for each height result. According to experimental data, the optimization of typical BP neural networks is not ideal, especially in nonvisual range environments, and rarely has centimetre-level accuracy.

    Table 2.  Results of promoting.
    UWB BP HHO-BP
    Environment H ME ME O P ME O P
    LOS 130 9.15 8.63 5.70 5.31 41.99
    155 8.13 7.39 9.08 4.94 39.23
    170 8.22 5.15 37.29 4.18 49.07
    195 11.11 4.89 55.91 3.56 67.95
    NLOS 24 27.68 18.19 34.28 13.75 50.33
    40 34.80 11.65 66.52 6.65 80.90
    55 38.05 20.41 46.36 12.45 67.28
    67 19.30 9.04 53.15 5.41 71.97
    Note: H: Height (cm), ME: Mean Error (cm); O P: Optimization percentage (%).

     | Show Table
    DownLoad: CSV

    Table 2 uses the statistical mean error for each height and further calculates the optimization percentages of the conventional BPNN and the HHO-BP algorithm. The mathematical formula is as follows in Eq (18):

    YHP=(YUYH)÷YU×100%YBP=(YUYB)÷YU×100% (18)

    where YHP is the optimization percentage of the HHO-BP algorithm, YBP is the optimization percentage of the conventional BPNN, YU is the mean error of the UWB measurement, YH is the mean error of the HHO-BP algorithm, and YB is the mean error of the conventional BPNN.

    According to Table 2, the HHO-BP algorithm outperforms the optimization results of the conventional BPNN at every height, both for the LOS and NLOS distances. In addition, for the HHO-BP algorithm in non-line-of-sight distance conditions, compared with UWB measurement data, the degree of optimization is more than 50%, and the effect is better than the line-of-sight distance environment. For 130 cm data, the conventional BP neural network is only optimized by 5.70%, while HHO-BP performs 36.29% better than the conventional BP neural network.

    Table 3.  Total data comparison.
    UWB BP HHO-BP HHO-BP greater than BP
    Environment MTE MTE O P MTE O P O P
    LOS 9.15 6.52 26.99 4.50 49.56 22.57
    NLOS 29.96 14.82 50.08 9.56 67.62 17.54
    Note: MTE: Mean total error (cm), O P: Optimization percentage (%).

     | Show Table
    DownLoad: CSV

    Table 3 shows the total mean errors for LOS and NLOS and the optimization percentages. Here, the mean errors in Tables 1 and 2 are averaged to obtain the total mean error of the corresponding environment. According to the data in Table 3, the HHO-BP algorithm achieves centimetre-level positioning accuracy in both LOS and NLOS environments. Under line-of-sight conditions, the HHO-BP algorithm improves the localization accuracy by 22.57% compared to the conventional BP algorithm; under non-line-of-sight conditions, the HHO-BP algorithm improves the localization accuracy by 17.54% in comparison to the conventional BPNN.

    In this study, an intelligent group algorithm (HHO) combined with a BP neural network is proposed to calibrate UWB indoor locations. First, the original and UWB coordinates are collected to obtain the actual errors. Next, BPNN training is used to generate the relevant weights and thresholds. Then, the HHO method optimizes the weights and thresholds of the BPNN until the minimum error is achieved. Finally, a new BPNN is constructed to calibrate the UWB measurement coordinates. The original BPNN and HHO-BP errors are compared. The experimental results show that, to some extent, BPNN calibration of the UWB indoor location will fall into the optimal local solution and that the HHO-BP method has higher accuracy. Finally, the proposed HHO-BP method is applied to indoor NLOS environments for calibrating indoor positioning environments with larger errors. The findings demonstrate that the HHO-BP algorithm's overall positioning error in an LOS environment is approximately 4.5 cm, which is 22.57% higher than the conventional BPNN, and for 130 cm data, the maximum improvement is 36.29% compared to the conventional BPNN. In an NLOS scenario, the HHO-BP algorithm has a positioning error of approximately 9.6 cm, which is 17.54% better than the conventional BPNN method. Experiments have confirmed the accuracy and efficacy of the suggested indoor positioning technique, indicating that it has greater calibration precision and stability, as well as strong future application possibilities.

    CXH conceived the experimental study methodology, participated in the design and drafted the manuscript. FMS and LZY conducted the research and helped refine the experimental ideas. JCC provided technical support and specific directions for the algorithms and helped revise the manuscript. LY collected the data and analysed the results. All authors read and approved the final manuscript.

    The datasets used and analysed during the current study are available from the corresponding author upon reasonable request.

    We are very grateful to the Robotics Research Laboratory of Shandong University of Science and Technology for providing an experimental site for our research. This work was supported by the Natural Science Research Project of Universities in Anhui Province (NO. KJ2021A0953), Natural Science Key Scientific Research Project of West Anhui University (NO. 0041021003, WXZR201903), Key Research Programs of Shandong Province (NO. 2018GGX103011) and Key Science and Technology Innovation Programs in Shandong Province (NO. 2017CXGC0919).

    There are no conflicts of interest declared by the authors.



    [1] UNCED, The United Nations conference on environment and development, 1992.
    [2] A. J. Bach, Assessing conditions leading to severe wind erosion in the Antelope Valley, California, 1990–1991, Prof. Geogr., 50 (2010), 87–97. https://doi.org/10.1111/00330124.00106 doi: 10.1111/00330124.00106
    [3] J. Leys, G. Mctainsh, Soil loss and nutrient decline by wind erosion-cause for concern, Aust. J. Soil Water Conserv., 7 (1994), p30–35.
    [4] D. P. C. Peters, K. M. Havstad, Nonlinear dynamics in arid and semi-arid systems: Interactions among drivers and processes across scales, J. Arid Environ., 65 (2006), 196–206. https://doi.org/10.1016/j.jaridenv.2005.05.010 doi: 10.1016/j.jaridenv.2005.05.010
    [5] D. D. Breshears, J. J. Whicker, M. P. Johansen, J. E. Pinder, Wind and water erosion and transport in semi-arid shrubland, grassland and forest ecosystems: quantifying dominance of horizontal wind-driven transport, Earth Surf. Processes Landf., 28 (2003), 1189–1209. https://doi.org/10.1002/esp.1034 doi: 10.1002/esp.1034
    [6] J. F. Weltzin, M. E. Loik, S. Schwinning, D. G. William, P. A. Fay, B. M. Haddad, et al., Assessing the response of terrestrial ecosystems to potential changes in precipitation, BioScience, 53 (2003), 941–952. https://doi.org/10.1641/0006-3568(2003)053[0941:ATROTE]2.0.CO;2 doi: 10.1641/0006-3568(2003)053[0941:ATROTE]2.0.CO;2
    [7] N. English, J. Weltzin, A. Fravolini, L. Thomas, D. G. Williams, The influence of soil texture and vegetation on soil moisture under rainout shelters in a semi-desert grassland, J. Arid Environ., 63 (2005), 324–343. https://doi.org/10.1016/j.jaridenv.2005.03.013 doi: 10.1016/j.jaridenv.2005.03.013
    [8] G. H. Guo, J. J. Wang, Pattern formation and qualitative analysis for a vegetation-water model with diffusion, Nonlinear Anal. Real World Appl., 76 (2024), 104008. https://doi.org/10.1016/j.nonrwa.2023.104008 doi: 10.1016/j.nonrwa.2023.104008
    [9] G. H. Guo, S. H. Zhao, J. J. Wang, Y. X. Gao, Positive steady-state solutions for a water-vegetation model with the infiltration feedback effect, Discrete Contin. Dyn. Syst. B, 29 (2023), 426–458. https://doi.org/10.3934/dcdsb.2023101 doi: 10.3934/dcdsb.2023101
    [10] G. H. Guo, Q. J. Qin, D. F. Pang, Y. H. Su, Positive steady-state solutions for a vegetation-water model with saturated water absorption, Commun. Nonlinear Sci. Numer. Simul., 131 (2024), 107802. https://doi.org/10.1016/j.cnsns.2023.107802 doi: 10.1016/j.cnsns.2023.107802
    [11] R. Bastiaansen, M. Chirilus-Bruckner, A. Doelman, Pulse solutions for an extended Klausmeier model with spatially varying coefficients, SIAM J. Appl. Dyn. Syst., 19 (2020), 1–57. https://doi.org/10.1137/19M1255665 doi: 10.1137/19M1255665
    [12] C. A. Klausmeier, Regular and irregular patterns in semiarid vegetation, Science, 284 (1999), 1826-1828. https://doi.org/10.1126/science.284.5421.1826 doi: 10.1126/science.284.5421.1826
    [13] E. Meron, E. Gilad, J. Von Hardenberg, M. Shachak, Y. Zarmi, Vegetation patterns along a rainfall gradient, Chaos Soliton. Fract., 19 (2004), 367–376. https://doi.org/10.1016/S0960-0779(03)00049-3 doi: 10.1016/S0960-0779(03)00049-3
    [14] G. McTainsh, C. Strong, The role of aeolian dust in ecosystems, Geomorphology, 89 (2006), 39–54. https://doi.org/10.1016/j.geomorph.2006.07.028 doi: 10.1016/j.geomorph.2006.07.028
    [15] P. P. Hesse, L. R. Simpson, Variable vegetation cover and episodic sand movement on longitudinal desert sand dunes, Geomorphology, 81 (2006), 276–291. https://doi.org/10.1016/j.geomorph.2006.04.012 doi: 10.1016/j.geomorph.2006.04.012
    [16] K. Burri, C. Gromke, M. Lehning, F. Graf, Aeolian sediment transport over vegetation canopies: A wind tunnel study with live plants, Aeolian Res., 3 (2011), 205–213. https://doi.org/10.1016/j.aeolia.2011.01.003 doi: 10.1016/j.aeolia.2011.01.003
    [17] Y. C. Yan, X. L. Xu, X. P. Xin, G. X. Yang, X. Wang, R. Yan, et al., Effect of vegetation coverage on aeolian dust accumulation in a semiarid steppe of northern China, Catena, 87 (2011), 351–356. https://doi.org/10.1016/j.catena.2011.07.002 doi: 10.1016/j.catena.2011.07.002
    [18] F. F. Zhang, L. Yao, W. J. Zhou, Q. J. You, H. Y. Zhang, Using shannon entropy and contagion index to interpret pattern self-organization in a dynamic vegetation-sand model, IEEE Access, 8 (2020), 17221–17230. https://doi.org/10.1109/ACCESS.2020.2968242 doi: 10.1109/ACCESS.2020.2968242
    [19] F. F. Zhang, Y. X. Li, Y. L. Zhao, Z. Liu, Vegetation pattern formation and transition caused by cross-diffusion in a modified vegetation-sand model, Int. J. Bifurcat. Chaos, 32 (2022), 2250069. https://doi.org/10.1142/S0218127422500699 doi: 10.1142/S0218127422500699
    [20] F. F. Zhang, H. Y. Zhang, M. R. Evans, T. Huang, Vegetation patterns generated by a wind driven sand-vegetation system in arid and semi-arid areas, Ecol. Complex., 31 (2017), 21–33. https://doi.org/10.1016/j.ecocom.2017.02.005 doi: 10.1016/j.ecocom.2017.02.005
    [21] M. Alfaro, H. Izuhara, M. Mimura, On a nonlocal system for vegetation in drylands, J. Math. Biol., 77 (2018), 1761–1793. https://doi.org/10.1007/s00285-018-1215-0 doi: 10.1007/s00285-018-1215-0
    [22] L. Eigentler, J. A. Sherratt, Analysis of a model for banded vegetation patterns in semi-arid environments with nonlocal dispersal, J. Math. Biol., 7 (2018), 739–763. https://doi.org/10.1007/s00285-018-1233-y doi: 10.1007/s00285-018-1233-y
    [23] Y. Maimaiti, W. B. Yang, J. H. Wu, Turing instability and coexistence in an extended Klausmeier model with nonlocal grazing, Nonlinear Anal. Real World Appl., 64 (2022), 103443. https://doi.org/10.1016/j.nonrwa.2021.103443 doi: 10.1016/j.nonrwa.2021.103443
    [24] E. Siero, Nonlocal grazing in patterned ecosystems, J. Theor. Biol., 436 (2018), 64–71. https://doi.org/10.1016/j.jtbi.2017.10.001 doi: 10.1016/j.jtbi.2017.10.001
    [25] S. Zaytseva, L. B. Shaw, J. P. Shi, M. L. Kirwan, R. N. Lipcius, Pattern formation in marsh ecosystems modeled through the interaction of marsh vegetation, mussels and sediment, J. Theor. Biol., 543 (2022), 111102. https://doi.org/10.1016/j.jtbi.2022.111102 doi: 10.1016/j.jtbi.2022.111102
    [26] C. H. Zeng, H. Wang, Noise and large time delay: Accelerated catastrophic regime shifts in ecosystems, Ecol. Model., 233 (2012), 52–58. https://doi.org/10.1016/j.ecolmodel.2012.03.025 doi: 10.1016/j.ecolmodel.2012.03.025
    [27] L.F. Lafuerza, R. Toral, Exact solution of a stochastic protein dynamics model with delayed degradation, Phys. Rev. E, 84 (2011), 051121. https://doi.org/10.1103/PhysRevE.84.051121 doi: 10.1103/PhysRevE.84.051121
    [28] S. L. Pan, Q. M. Zhang, T. Kang, A. Meyer-Baese, X. Li, Finite-time stability of a stochastic tree-grass-water-nitrogen system with impulsive and time-varying delay, Int. J. Biomath., 17 (2023), 2350052. https://doi.org/10.1142/S1793524523500523 doi: 10.1142/S1793524523500523
    [29] Q. L. Han, T. Yang, C. H. Zeng, H. Wang, Z. Liu, Y. Fu, et al., Impact of time delays on stochastic resonance in an ecological system describing vegetation, Physica A, 408 (2014), 96–105. https://doi.org/10.1016/j.physa.2014.04.015 doi: 10.1016/j.physa.2014.04.015
    [30] J. Li, G. Q. Sun, Z. Jin, Interactions of time delay and spatial diffusion induce the periodic oscillation of the vegetation system, Discrete Contin. Dyn. Syst. Ser. B, 27 (2022), 2147–2172. https://doi.org/10.3934/dcdsb.2021127 doi: 10.3934/dcdsb.2021127
    [31] N. F. Britton, Spatial structures and periodic travelling waves in an integro-differential reaction-diffusion population model, SIAM J. Appl. Math., 50 (1990), 1663–1688. https://doi.org/10.1137/0150099 doi: 10.1137/0150099
    [32] H. Wang, On the existence of traveling waves for delayed reaction-diffusion equations, J. Differ. Equ., 247 (2015), 887–905. https://doi.org/10.1016/j.jde.2009.04.002 doi: 10.1016/j.jde.2009.04.002
    [33] G. Y. Lue, M. X. Wang, Stability of planar waves in reaction-diffusion system, Sci. China Math., 54 (2011), 1403–1419. https://doi.org/10.1007/s11425-011-4210-0 doi: 10.1007/s11425-011-4210-0
    [34] M. D. Burlica, D. Rosu, I. I. Vrabie, Abstract reaction-diffusion systems with nonlocal initial conditions, Nonlinear Anal-Theor., 94 (2014), 107–119. https://doi.org/10.1016/j.na.2013.07.033 doi: 10.1016/j.na.2013.07.033
    [35] G. F. Webb, A reaction-diffusion model for a deterministic diffusive epidemic, J. Math. Anal. Appl., 84 (1981), 150–161. https://doi.org/10.1016/0022-247X(81)90156-6 doi: 10.1016/0022-247X(81)90156-6
    [36] S. J. Guo, S. Z. Li, On the stability of reaction-diffusion models with nonlocal delay effect and nonlinear boundary condition, Appl. Math. Lett., 103 (2020), 106197. https://doi.org/10.1016/j.aml.2019.106197 doi: 10.1016/j.aml.2019.106197
    [37] M. Aguerrea, G. Valenzuela, On the minimal speed of traveling waves for a nonlocal delayed reaction-diffusion equation, Nonlinear Oscil., 13 (2010), 1–9. https://doi.org/10.1007/s11072-010-0096-y doi: 10.1007/s11072-010-0096-y
    [38] Q. Xue, G. Q. Sun, C. Liu, Z. G. Guo, Z. Jin, Y. P. Wu, et al., Spatiotemporal dynamics of a vegetation model with nonlocal delay in semi-arid environment, Nonlinear Dyn., 99 (2020), 3407–3420. https://doi.org/10.1007/s11071-020-05486-w doi: 10.1007/s11071-020-05486-w
    [39] C. Liu, F. G. Wang, Q. Xue, L. Li, Z. Wang, Pattern formation of a spatial vegetation system with root hydrotropism, Appl. Math. Comput., 420 (2022), 126913. https://doi.org/10.1016/j.amc.2021.126913 doi: 10.1016/j.amc.2021.126913
    [40] Q. Xue, C. Liu, L. Li, G.Q. Sun, Z. Wang, Interactions of diffusion and nonlocal delay give rise to vegetation patterns in semi-arid environments, Appl. Math. Computation, 399 (2021), 126038. https://doi.org/10.1016/j.amc.2021.126038 doi: 10.1016/j.amc.2021.126038
    [41] S. J. GUO, Stability and bifurcation in a reaction-diffusion model with nonlocal delay effect, J. Differ. Equ., 259 (2015), 1409–1448. https://doi.org/10.1016/j.jde.2015.03.006 doi: 10.1016/j.jde.2015.03.006
    [42] C. H. Wang, H. Wang, S. L. Yuan, Precipitation governing vegetation patterns in an arid or semi-arid environment, J. Math. Biol., 87 (2023), 22. https://doi.org/10.1007/s00285-023-01954-0 doi: 10.1007/s00285-023-01954-0
    [43] L. Eigentler, J. A. Sherratt, An integrodifference model for vegetation patterns in semi-arid environments with seasonality, J. Math. Biol., 81 (2020), 875–904. https://doi.org/10.1007/s00285-020-01530-w doi: 10.1007/s00285-020-01530-w
    [44] Z. L. Zhen, J. D. Wei, J. B. Zhou, L. X. Tian, Wave propagation in a nonlocal diffusion epidemic model with nonlocal delayed effects, Appl. Math. Comput., 339 (2018), 15–37. https://doi.org/10.1016/j.amc.2018.07.007 doi: 10.1016/j.amc.2018.07.007
    [45] S. A. Gourley, S. JWH, Dynamics of a food-limited population model incorporating nonlocal delays on a finite domain, J. Math. Biol., 44 (2002), 49–78. https://doi.org/10.1007/s002850100109 doi: 10.1007/s002850100109
    [46] G. Q. Sun, C. H. Wang, Z. Y. Wu, Pattern dynamics of a Gierer-Meinhardt model with spatial effects, Nonlinear Dyn., 88 (2017), 1385–1396. https://doi.org/10.1007/s11071-016-3317-9 doi: 10.1007/s11071-016-3317-9
    [47] Y. L. Song, R. Yang, G. Q. Sun, Pattern dynamics in a Gierer-Meinhardt model with a saturating term, Appl. Math. Model., 46 (2017), 476–491. https://doi.org/10.1016/j.apm.2017.01.081 doi: 10.1016/j.apm.2017.01.081
    [48] C. Liu, L. L. Chang, Y. Huang, Z. Wang, Turing patterns in a predator-prey model on complex networks, Nonlinear Dyn., 99 (2020), 3313–3322. https://doi.org/10.1007/s11071-019-05460-1 doi: 10.1007/s11071-019-05460-1
    [49] J. Y. Zhou, Y. Ye, A. Arenas, S. Gomz, Y. Zhao, Pattern formation and bifurcation analysis of delay induced fractional-order epidemic spreading on networks, Chaos, Soliton. Fract., 174 (2023), 113805. https://doi.org/10.1016/j.chaos.2023.113805 doi: 10.1016/j.chaos.2023.113805
  • This article has been cited by:

    1. Guo-zheng Quan, Yi-fan Zhao, Qi Deng, Ming-guo Quan, Wei Xiong, Flow Stress Models for 40Cr10Si2Mo Steel and Their Application in Numerical Simulation of Hot Forming, 2024, 1059-9495, 10.1007/s11665-024-10024-5
    2. Feiyun Cong, Zongheng Hong, Xiaojie Lin, Huimin Li, Peak Ratio Iteration-Based Leading-Edge Detection Algorithm in UWB Localization, 2023, 131, 0929-6212, 1663, 10.1007/s11277-023-10517-x
    3. Jia Xu, Teng Xiao, Pin Lv, Zhe Chen, Chao Cai, Yang Zhang, Zehui Xiong, Tracing Human Stress From Physiological Signals Using UWB Radar, 2024, 11, 2327-4662, 32773, 10.1109/JIOT.2024.3410033
    4. Qingxi Zeng, Yu Jin, Haonan Yu, Xia You, A UAV Localization System Based on Double UWB Tags and IMU for Landing Platform, 2023, 23, 1530-437X, 10100, 10.1109/JSEN.2023.3260311
    5. Dansong Yue, Shuqi Shang, Kai Feng, Haiqing Wang, Xiaoning He, Zelong Zhao, Ning Zhang, Baiqiang Zuo, Dongwei Wang, Research on the Model of a Navigation and Positioning Algorithm for Agricultural Machinery Based on the IABC-BP Network, 2023, 13, 2077-0472, 1769, 10.3390/agriculture13091769
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1085) PDF downloads(52) Cited by(2)

Figures and Tables

Figures(10)  /  Tables(3)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog