Let{Z_(n)}_(n)≥0 be a critical or subcritical d-dimensional branching random walk started from a Poisson random measure whose intensity measure is the Lebesugue measure onℝ^(d).Denote by R_(n):=sup{u>0:Z_(n)({x∈ℝ...Let{Z_(n)}_(n)≥0 be a critical or subcritical d-dimensional branching random walk started from a Poisson random measure whose intensity measure is the Lebesugue measure onℝ^(d).Denote by R_(n):=sup{u>0:Z_(n)({x∈ℝ^(d):∣x∣<u})=0}the radius of the largest empty ball centered at the origin of Zn.In this work,we prove that after suitable renormalization,Rn converges in law to some non-degenerate distribution as n→∈.Furthermore,our work shows that the renormalization scales depend on the offspring law and the dimension of the branching random walk.This completes the results of Révész[13]for the critical binary branching Wiener process.展开更多
We study a counterbalanced random walkS_(n)=X_(1)+…+X_(n),which is a discrete time non-Markovian process andX_(n) are given recursively as follows.For n≥2,X_(n) is a new independent sample from some fixed law̸=0 wit...We study a counterbalanced random walkS_(n)=X_(1)+…+X_(n),which is a discrete time non-Markovian process andX_(n) are given recursively as follows.For n≥2,X_(n) is a new independent sample from some fixed law̸=0 with a fixed probability p,andX_(n)=−X_(v(n))with probability 1−p,where v(n)is a uniform random variable on{1;…;n−1}.We apply martingale method to obtain a strong invariance principle forS_(n).展开更多
We primarily provide several estimates for the heat kernel defined on the 2-dimensional simple random walk. Additionally, we offer an estimate for the heat kernel on high-dimensional random walks, demonstrating that t...We primarily provide several estimates for the heat kernel defined on the 2-dimensional simple random walk. Additionally, we offer an estimate for the heat kernel on high-dimensional random walks, demonstrating that the heat kernel in higher dimensions converges rapidly. We also compute the constants involved in the estimate for the 1-dimensional heat kernel. Furthermore, we discuss the general case of on-diagonal estimates for the heat kernel.展开更多
In this paper we define direct product of graphs and give a recipe for obtaining probability of observing particle on vertices in the continuous-time classical and quantum random walk. In the recipe, the probability o...In this paper we define direct product of graphs and give a recipe for obtaining probability of observing particle on vertices in the continuous-time classical and quantum random walk. In the recipe, the probability of observing particle on direct product of graph is obtained by multiplication of probability on the corresponding to sub-graphs, where this method is useful to determining probability of walk on compficated graphs. Using this method, we calculate the probability of Continuous-time classical and quantum random walks on many of finite direct product Cayley graphs (complete cycle, complete Kn, charter and n-cube). Also, we inquire that the classical state the stationary uniform distribution is reached as t→∞ but for quantum state is not always satisfied.展开更多
Based on the random walk and the intentional random walk, we propose two types of immunization strategies which require only local connectivity information. On several typical scale-free networks, we demonstrate that ...Based on the random walk and the intentional random walk, we propose two types of immunization strategies which require only local connectivity information. On several typical scale-free networks, we demonstrate that these strategies can lead to the eradication of the epidemic by immunizing a small fraction of the nodes in the networks. Particularly, the immunization strategy based on the intentional random walk is extremely efficient for the assortatively mixed networks.展开更多
Unstructured P2P has power-law link distribution, and the random walk in power-law networks is analyzed. The analysis results show that the probability that a random walker walks through the high degree nodes is high ...Unstructured P2P has power-law link distribution, and the random walk in power-law networks is analyzed. The analysis results show that the probability that a random walker walks through the high degree nodes is high in the power-law network, and the information on the high degree nodes can be easily found through random walk. Random walk spread and random walk search method (RWSS) is proposed based on the analysis result. Simulation results show that RWSS achieves high success rates at low cost and is robust to high degree node failure.展开更多
We study the mixing rate of non-backtracking random walks on graphs by looking at non-backtracking walks as walks on the directed edges of a graph. A result known as Ihara’s Theorem relates the adjacency matrix of a ...We study the mixing rate of non-backtracking random walks on graphs by looking at non-backtracking walks as walks on the directed edges of a graph. A result known as Ihara’s Theorem relates the adjacency matrix of a graph to a matrix related to non-backtracking walks on the directed edges. We prove a weighted version of Ihara’s Theorem which relates the transition probability matrix of a non-backtracking walk to the transition matrix for the usual random walk. This allows us to determine the spectrum of the transition probability matrix of a non-backtracking random walk in the case of regular graphs and biregular graphs. As a corollary, we obtain a result of Alon et al. in [1] that in most cases, a non-backtracking random walk on a regular graph has a faster mixing rate than the usual random walk. In addition, we obtain an analogous result for biregular graphs.展开更多
We consider a branching random walk with a random environment m time, in which the offspring distribution of a particle of generation n and the distribution of the displacements of its children depend on an environmen...We consider a branching random walk with a random environment m time, in which the offspring distribution of a particle of generation n and the distribution of the displacements of its children depend on an environment indexed by the time n. The envi- ronment is supposed to be independent and identically distributed. For A C R, let Zn(A) be the number of particles of generation n located in A. We show central limit theorems for the counting measure Zn (-) with appropriate normalization.展开更多
In this paper, considering both cluster heads and sensor nodes, we propose a novel evolving a network model based on a random walk to study the fault tolerance decrease of wireless sensor networks (WSNs) due to node...In this paper, considering both cluster heads and sensor nodes, we propose a novel evolving a network model based on a random walk to study the fault tolerance decrease of wireless sensor networks (WSNs) due to node failure, and discuss the spreading dynamic behavior of viruses in the evolution model. A theoretical analysis shows that the WSN generated by such an evolution model not only has a strong fault tolerance, but also can dynamically balance the energy loss of the entire network. It is also found that although the increase of the density of cluster heads in the network reduces the network efficiency, it can effectively inhibit the spread of viruses. In addition, the heterogeneity of the network improves the network efficiency and enhances the virus prevalence. We confirm all the theoretical results with sufficient numerical simulations.展开更多
A general formulation of the stochastic model for random walk in time-random environment and an equivalent definition is established in this paper. Moreover, some basic probability relations similar to the classical c...A general formulation of the stochastic model for random walk in time-random environment and an equivalent definition is established in this paper. Moreover, some basic probability relations similar to the classical case which are very useful in the corresponding research of fractal properties are given. At the end, a typical example is provided to show the recurrence and transience. Key words random environment - random walk in timerandom environment - skew product Markov chain CLC number O 211.6 Foudation item: Supported by the National Natural Science Foundation of China (10371092) and Foundation of Wuhan University.Biography: Zhang Xiao-min (1977-), male, Ph. D candidate, research direction: stochastic processes and random fractal.展开更多
In this paper, some experimental studies on the impact of effluent from an exhaust tower of an underground tunnel with special construction are reported. By measuring the flow field downstream of the tower in NJU mete...In this paper, some experimental studies on the impact of effluent from an exhaust tower of an underground tunnel with special construction are reported. By measuring the flow field downstream of the tower in NJU meteorological wind tunnel, some flow characteristics in the make area were established. Based on these, an advanced random\|walk dispersion model was set up and applied successfully to the simulation of dispersion in the wake area. The modelling results were in accordance with wind tunnel measurements. The computed maximum of ground surface concentration in the building case was a factor of 3-4 higher than that in the flat case and appeared much closer to the source. The simulation indicated that random walk modelling is an effective and practical tool for the wake stream impact assessment.展开更多
We consider a random walk on Z in random environment with possible jumps {-L,…, -1, 1}, in the case that the environment {ωi : i ∈ Z} are i.i.d.. We establish the renewal theorem for the Markov chain of "the envi...We consider a random walk on Z in random environment with possible jumps {-L,…, -1, 1}, in the case that the environment {ωi : i ∈ Z} are i.i.d.. We establish the renewal theorem for the Markov chain of "the environment viewed from the particle" in both annealed probability and quenched probability, which generalize partially the results of Kesten (1977) and Lalley (1986) for the nearest random walk in random environment on Z, respectively. Our method is based on (L, 1)-RWRE formulated in Hong and Wang the intrinsic branching structure within the (2013).展开更多
In this paper, we study the scaling for the mean first-passage time (MFPT) of the random walks on a generalized Koch network with a trap. Through the network construction, where the initial state is transformed from...In this paper, we study the scaling for the mean first-passage time (MFPT) of the random walks on a generalized Koch network with a trap. Through the network construction, where the initial state is transformed from a triangle to a polygon, we obtain the exact scaling for the MFPT. We show that the MFPT grows linearly with the number of nodes and the dimensions of the polygon in the large limit of the network order. In addition, we determine the exponents of scaling efficiency characterizing the random walks. Our results are the generalizations of those derived for the Koch network, which shed light on the analysis of random walks over various fractal networks.展开更多
A numerical model has been developed to simulate the transport and fate of oil spilled at sea. The model combines the transport and fate processes of spilled oil with the random walk technique. Oil movement under th...A numerical model has been developed to simulate the transport and fate of oil spilled at sea. The model combines the transport and fate processes of spilled oil with the random walk technique. Oil movement under the influence of tidal currents, wind driven currents, and turbulent eddies is simulated by the PLUME RW dispersion model developed by HR Wallingford. The weathering processes in the model represent physical and chemical changes of soil slicks with time, and comprise mechanical spreading, dispersion, evaporation and emulsification. Shoreline stranding is determined approximately using a capacity method for different shoreline types. This paper presents details of the model, and describe the results of various sensitivity tests. The model is suitable for oil spill contingency planning.展开更多
In this article, the dependent steps of a negative drift random walk are modelled as a two-sided linear process. Xn=-u+∑j=-∞^∞ φn-jεj, where { ε, εn; -∞〈n〈+∞} is a sequence of independent, identically di...In this article, the dependent steps of a negative drift random walk are modelled as a two-sided linear process. Xn=-u+∑j=-∞^∞ φn-jεj, where { ε, εn; -∞〈n〈+∞} is a sequence of independent, identically distributed random variables with zero mean, u 〉 0 is a constant and the coefficients {φi; -∞〈i〈∞} satisfy 0〈 ∑j=-∞^∞ |jφj|〈 ∞ . Under the conditions that the distribution function of |ε| has dominated variation and ε satisfies certain tail balance conditions, the asymptotic behavior of P{sup n≥0 (-qu+∑j=-∞^∞ εj βnj)〉x} is discussed. Then the result is applied to ultimate ruin probability.展开更多
The water exchange matrix is an efficient tool to study the water exchange among the sub-areas in large-scale bays. The application of the random walk method to calculate the water exchange matrix is studied. Compared...The water exchange matrix is an efficient tool to study the water exchange among the sub-areas in large-scale bays. The application of the random walk method to calculate the water exchange matrix is studied. Compared with the advection-diffusion model, the random walk model is more flexible to calculate the water exchange matrix. The forecast matrix suggested by Thompson et al. is used to evaluate the water exchange characteristics among the sub-areas fast. According to the theoretic analysis, it is found that the precision of the predicted results is mainly affected by three factors, namely, the particle number, the generated time of the forecast matrix, and the number of the sub-areas. The impact of the above factors is analyzed based on the results of a series of numerical tests. The results show that the precision of the forecast matrix increases with the increase of the generated time of the forecast matrix and the number of the particles. If there are enough particles in each sub-area, the precision of the forecast matrix will increase with the number of the sub-areas. Moreover, if the particles in each sub-area are not enough, the excessive number of the sub-areas can result in the decrease of the precision of the forecast matrix.展开更多
Consider a branching random walk, where the underlying branching mechanism is governed by a Galton-Watson process and the moving law of particles by a discrete random variable on the integer lattice Z. Denote by Zn(z...Consider a branching random walk, where the underlying branching mechanism is governed by a Galton-Watson process and the moving law of particles by a discrete random variable on the integer lattice Z. Denote by Zn(z) the number of particles in the n-th generation in the model for each z ∈ Z. We derive the exact convergence rate in the local limit theorem for Zn(z) assuming a condition like "EN(logN)1+λ 〈 ∞" for the offspring distribution and a finite moment condition on the motion law. This complements the known results for the strongly non-lattice branching random walk on the real line and for the simple symmetric branching random walk on the integer lattice.展开更多
Significant and persistent trajectory-to-trajectory variance are commonly observed in particle tracking experiments,which have become a major challenge for the experimental data analysis.In this theoretical paper we i...Significant and persistent trajectory-to-trajectory variance are commonly observed in particle tracking experiments,which have become a major challenge for the experimental data analysis.In this theoretical paper we investigate the ergodicity recovery behavior,which helps clarify the origin and the convergence of trajectory-to-trajectory fluctuation in various heterogeneous disordered media.The concepts of self-averaging and ergodicity are revisited in the context of trajectory analysis.The slow ergodicity recovery and the non-Gaussian diffusion in the annealed disordered media are shown as the consequences of the central limit theorem in different situations.The strange ergodicity recovery behavior is reported in the quenched disordered case,which arises from a localization mechanism.The first-passage approach is introduced to the ergodicity analysis for this case,of which the central limit theorem can be employed and the ergodicity is recovered in the length scale of diffusivity correlation.展开更多
An approach of solving the finite difference equations with Monte Carlo method was presented. The accuracy of algorithm is guaranteed by the Central Limit Theorem. The computation of the value on each single node is i...An approach of solving the finite difference equations with Monte Carlo method was presented. The accuracy of algorithm is guaranteed by the Central Limit Theorem. The computation of the value on each single node is independent of each other, which makes it easy to realize the algorithm in parallel processing and greatly improves the efficiency while dealing with local area computation. As long as different kinds of boundary conditions are statistics modeled,wide applications can the be made where the finite difference method is of competence.展开更多
基金supported by the National Key R&D Program of China(2022YFA1006102).
文摘Let{Z_(n)}_(n)≥0 be a critical or subcritical d-dimensional branching random walk started from a Poisson random measure whose intensity measure is the Lebesugue measure onℝ^(d).Denote by R_(n):=sup{u>0:Z_(n)({x∈ℝ^(d):∣x∣<u})=0}the radius of the largest empty ball centered at the origin of Zn.In this work,we prove that after suitable renormalization,Rn converges in law to some non-degenerate distribution as n→∈.Furthermore,our work shows that the renormalization scales depend on the offspring law and the dimension of the branching random walk.This completes the results of Révész[13]for the critical binary branching Wiener process.
基金Supported by the National Natural Science Foundation of China(11671373).
文摘We study a counterbalanced random walkS_(n)=X_(1)+…+X_(n),which is a discrete time non-Markovian process andX_(n) are given recursively as follows.For n≥2,X_(n) is a new independent sample from some fixed law̸=0 with a fixed probability p,andX_(n)=−X_(v(n))with probability 1−p,where v(n)is a uniform random variable on{1;…;n−1}.We apply martingale method to obtain a strong invariance principle forS_(n).
文摘We primarily provide several estimates for the heat kernel defined on the 2-dimensional simple random walk. Additionally, we offer an estimate for the heat kernel on high-dimensional random walks, demonstrating that the heat kernel in higher dimensions converges rapidly. We also compute the constants involved in the estimate for the 1-dimensional heat kernel. Furthermore, we discuss the general case of on-diagonal estimates for the heat kernel.
文摘In this paper we define direct product of graphs and give a recipe for obtaining probability of observing particle on vertices in the continuous-time classical and quantum random walk. In the recipe, the probability of observing particle on direct product of graph is obtained by multiplication of probability on the corresponding to sub-graphs, where this method is useful to determining probability of walk on compficated graphs. Using this method, we calculate the probability of Continuous-time classical and quantum random walks on many of finite direct product Cayley graphs (complete cycle, complete Kn, charter and n-cube). Also, we inquire that the classical state the stationary uniform distribution is reached as t→∞ but for quantum state is not always satisfied.
文摘Based on the random walk and the intentional random walk, we propose two types of immunization strategies which require only local connectivity information. On several typical scale-free networks, we demonstrate that these strategies can lead to the eradication of the epidemic by immunizing a small fraction of the nodes in the networks. Particularly, the immunization strategy based on the intentional random walk is extremely efficient for the assortatively mixed networks.
文摘Unstructured P2P has power-law link distribution, and the random walk in power-law networks is analyzed. The analysis results show that the probability that a random walker walks through the high degree nodes is high in the power-law network, and the information on the high degree nodes can be easily found through random walk. Random walk spread and random walk search method (RWSS) is proposed based on the analysis result. Simulation results show that RWSS achieves high success rates at low cost and is robust to high degree node failure.
文摘We study the mixing rate of non-backtracking random walks on graphs by looking at non-backtracking walks as walks on the directed edges of a graph. A result known as Ihara’s Theorem relates the adjacency matrix of a graph to a matrix related to non-backtracking walks on the directed edges. We prove a weighted version of Ihara’s Theorem which relates the transition probability matrix of a non-backtracking walk to the transition matrix for the usual random walk. This allows us to determine the spectrum of the transition probability matrix of a non-backtracking random walk in the case of regular graphs and biregular graphs. As a corollary, we obtain a result of Alon et al. in [1] that in most cases, a non-backtracking random walk on a regular graph has a faster mixing rate than the usual random walk. In addition, we obtain an analogous result for biregular graphs.
基金partially supported by the National Natural Science Foundation of China(NSFC,11101039,11171044,11271045)a cooperation program between NSFC and CNRS of France(11311130103)+1 种基金the Fundamental Research Funds for the Central UniversitiesHunan Provincial Natural Science Foundation of China(11JJ2001)
文摘We consider a branching random walk with a random environment m time, in which the offspring distribution of a particle of generation n and the distribution of the displacements of its children depend on an environment indexed by the time n. The envi- ronment is supposed to be independent and identically distributed. For A C R, let Zn(A) be the number of particles of generation n located in A. We show central limit theorems for the counting measure Zn (-) with appropriate normalization.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 61103231 and 61103230)the Innovation Program of Graduate Scientific Research in Institution of Higher Education of Jiangsu Province, China (Grant No. CXZZ11 0401)
文摘In this paper, considering both cluster heads and sensor nodes, we propose a novel evolving a network model based on a random walk to study the fault tolerance decrease of wireless sensor networks (WSNs) due to node failure, and discuss the spreading dynamic behavior of viruses in the evolution model. A theoretical analysis shows that the WSN generated by such an evolution model not only has a strong fault tolerance, but also can dynamically balance the energy loss of the entire network. It is also found that although the increase of the density of cluster heads in the network reduces the network efficiency, it can effectively inhibit the spread of viruses. In addition, the heterogeneity of the network improves the network efficiency and enhances the virus prevalence. We confirm all the theoretical results with sufficient numerical simulations.
文摘A general formulation of the stochastic model for random walk in time-random environment and an equivalent definition is established in this paper. Moreover, some basic probability relations similar to the classical case which are very useful in the corresponding research of fractal properties are given. At the end, a typical example is provided to show the recurrence and transience. Key words random environment - random walk in timerandom environment - skew product Markov chain CLC number O 211.6 Foudation item: Supported by the National Natural Science Foundation of China (10371092) and Foundation of Wuhan University.Biography: Zhang Xiao-min (1977-), male, Ph. D candidate, research direction: stochastic processes and random fractal.
文摘In this paper, some experimental studies on the impact of effluent from an exhaust tower of an underground tunnel with special construction are reported. By measuring the flow field downstream of the tower in NJU meteorological wind tunnel, some flow characteristics in the make area were established. Based on these, an advanced random\|walk dispersion model was set up and applied successfully to the simulation of dispersion in the wake area. The modelling results were in accordance with wind tunnel measurements. The computed maximum of ground surface concentration in the building case was a factor of 3-4 higher than that in the flat case and appeared much closer to the source. The simulation indicated that random walk modelling is an effective and practical tool for the wake stream impact assessment.
文摘We consider a random walk on Z in random environment with possible jumps {-L,…, -1, 1}, in the case that the environment {ωi : i ∈ Z} are i.i.d.. We establish the renewal theorem for the Markov chain of "the environment viewed from the particle" in both annealed probability and quenched probability, which generalize partially the results of Kesten (1977) and Lalley (1986) for the nearest random walk in random environment on Z, respectively. Our method is based on (L, 1)-RWRE formulated in Hong and Wang the intrinsic branching structure within the (2013).
基金Project supported by the Research Foundation of Hangzhou Dianzi University,China (Grant Nos. KYF075610032 andzx100204004-7)the Hong Kong Research Grants Council,China (Grant No. CityU 1114/11E)
文摘In this paper, we study the scaling for the mean first-passage time (MFPT) of the random walks on a generalized Koch network with a trap. Through the network construction, where the initial state is transformed from a triangle to a polygon, we obtain the exact scaling for the MFPT. We show that the MFPT grows linearly with the number of nodes and the dimensions of the polygon in the large limit of the network order. In addition, we determine the exponents of scaling efficiency characterizing the random walks. Our results are the generalizations of those derived for the Koch network, which shed light on the analysis of random walks over various fractal networks.
文摘A numerical model has been developed to simulate the transport and fate of oil spilled at sea. The model combines the transport and fate processes of spilled oil with the random walk technique. Oil movement under the influence of tidal currents, wind driven currents, and turbulent eddies is simulated by the PLUME RW dispersion model developed by HR Wallingford. The weathering processes in the model represent physical and chemical changes of soil slicks with time, and comprise mechanical spreading, dispersion, evaporation and emulsification. Shoreline stranding is determined approximately using a capacity method for different shoreline types. This paper presents details of the model, and describe the results of various sensitivity tests. The model is suitable for oil spill contingency planning.
基金Research supported by National Science Foundation of China (70671018 and 10371117)
文摘In this article, the dependent steps of a negative drift random walk are modelled as a two-sided linear process. Xn=-u+∑j=-∞^∞ φn-jεj, where { ε, εn; -∞〈n〈+∞} is a sequence of independent, identically distributed random variables with zero mean, u 〉 0 is a constant and the coefficients {φi; -∞〈i〈∞} satisfy 0〈 ∑j=-∞^∞ |jφj|〈 ∞ . Under the conditions that the distribution function of |ε| has dominated variation and ε satisfies certain tail balance conditions, the asymptotic behavior of P{sup n≥0 (-qu+∑j=-∞^∞ εj βnj)〉x} is discussed. Then the result is applied to ultimate ruin probability.
基金supported by the National Natural Science Foundation of China(No.10702050)
文摘The water exchange matrix is an efficient tool to study the water exchange among the sub-areas in large-scale bays. The application of the random walk method to calculate the water exchange matrix is studied. Compared with the advection-diffusion model, the random walk model is more flexible to calculate the water exchange matrix. The forecast matrix suggested by Thompson et al. is used to evaluate the water exchange characteristics among the sub-areas fast. According to the theoretic analysis, it is found that the precision of the predicted results is mainly affected by three factors, namely, the particle number, the generated time of the forecast matrix, and the number of the sub-areas. The impact of the above factors is analyzed based on the results of a series of numerical tests. The results show that the precision of the forecast matrix increases with the increase of the generated time of the forecast matrix and the number of the particles. If there are enough particles in each sub-area, the precision of the forecast matrix will increase with the number of the sub-areas. Moreover, if the particles in each sub-area are not enough, the excessive number of the sub-areas can result in the decrease of the precision of the forecast matrix.
文摘Consider a branching random walk, where the underlying branching mechanism is governed by a Galton-Watson process and the moving law of particles by a discrete random variable on the integer lattice Z. Denote by Zn(z) the number of particles in the n-th generation in the model for each z ∈ Z. We derive the exact convergence rate in the local limit theorem for Zn(z) assuming a condition like "EN(logN)1+λ 〈 ∞" for the offspring distribution and a finite moment condition on the motion law. This complements the known results for the strongly non-lattice branching random walk on the real line and for the simple symmetric branching random walk on the integer lattice.
基金Project supported by the National Natural Science Foundation of China(Grants Nos.11705064,11675060,and 91730301).
文摘Significant and persistent trajectory-to-trajectory variance are commonly observed in particle tracking experiments,which have become a major challenge for the experimental data analysis.In this theoretical paper we investigate the ergodicity recovery behavior,which helps clarify the origin and the convergence of trajectory-to-trajectory fluctuation in various heterogeneous disordered media.The concepts of self-averaging and ergodicity are revisited in the context of trajectory analysis.The slow ergodicity recovery and the non-Gaussian diffusion in the annealed disordered media are shown as the consequences of the central limit theorem in different situations.The strange ergodicity recovery behavior is reported in the quenched disordered case,which arises from a localization mechanism.The first-passage approach is introduced to the ergodicity analysis for this case,of which the central limit theorem can be employed and the ergodicity is recovered in the length scale of diffusivity correlation.
文摘An approach of solving the finite difference equations with Monte Carlo method was presented. The accuracy of algorithm is guaranteed by the Central Limit Theorem. The computation of the value on each single node is independent of each other, which makes it easy to realize the algorithm in parallel processing and greatly improves the efficiency while dealing with local area computation. As long as different kinds of boundary conditions are statistics modeled,wide applications can the be made where the finite difference method is of competence.