• 沒有找到結果。

高維度離散型週期性雙穩定型系統的波之傳遞

N/A
N/A
Protected

Academic year: 2021

Share "高維度離散型週期性雙穩定型系統的波之傳遞"

Copied!
4
0
0

加載中.... (立即查看全文)

全文

(1)

1. Introduction

We consider the following lattice dynamical system for unknown u = {uj}j∈Z:

(1.1) u0

j = dj+1uj+1+ djuj−1− (dj+1+ dj)uj + fj(uj), j ∈ Z,

where fj ∈ C1+α[0, 1] for some α ∈ (0, 1) for j ∈ Z, fj+N = fj and dj+N = dj > 0 for all j ∈ Z for some positive integer N. The equation (1.1) can be regarded as a spatial discrete

version of the following reaction-diffusion equation

ut= (d(x)ux)x+ f (x, u),

where d(x) and f (x, u) are periodic in x. In biology, let uj denote the density of a certain species in a periodic patchy environment. Assuming the species at site j can only interact with those at the nearby sites, then the equation (1.1) describes the rate of change of density of this species at each site j. It is equal to the sum of the source fj(uj) at site j and the fluxes qj±1 from sites j ± 1 to site j:

qj+1 := dj+1[uj+1− uj], qj−1 := dj[uj−1− uj],

where dj, dj+1 are the diffusion constants. See [8, 15, 16] for more references and details. It is trivial that for a given initial data {uj(0)} ∈ [0, 1] there exists a unique solution u to (1.1) for t ≥ 0 such that 0 ≤ uj(t) ≤ 1 for all t ≥ 0 and j ∈ Z. We are interested in the wave propagation phenomenon. In particular, we are interested in special solutions U of (1.1) for

t ∈ R satisfying the following conditions:

Uj(t + N/c) = Uj−N(t), t ∈ R, j ∈ Z, (1.2)

Uj(t) → 1 as j → −∞, Uj(t) → 0 as j → +∞, locally in t ∈ R, (1.3)

for some nonzero constant c. We shall call a solution (c, U) of (1.1)-(1.3) as a traveling wave solution. The constant c is the wave speed and U is the profile. In this paper, we shall always assume that

(1.4) fj(0) = fj(1) = 0 ∀ j ∈ Z.

The study of traveling wave for lattice dynamical system has attracted a lot of attention for past years, see, e.g., the works [1, 2, 3, 4, 6, 7, 10, 11, 12, 13, 14, 17, 18, 19, 20]. The main concerns are existence, uniqueness, and stability of traveling waves. Typically, there are two different nonlinearities, namely, monostable and bistable cases. In the monostable case, we have

(1.5) f0

j(1) < 0 < fj0(0) ∀ j ∈ Z, fj(s) > 0 ∀ s ∈ (0, 1), j ∈ Z. For the bistable case, we have f0

j(0) < 0 and fj0(1) < 0 for all j ∈ Z. If N = 1, then fj+1 = fj and dj+1 = dj for all j. This is the so-called homogeneous media case. In general, if N > 1, then it is called the periodic case.

In this paper, we shall focus on the periodic monostable case. We refer the reader to the work [5] and the references cited therein for the periodic bistable case. In [5], the existence, uniqueness and stability of traveling waves for periodic bistable case are studied in details.

The existence of traveling waves for monostable case in periodic media was first obtained by Hudson and Zinner [11, 12] under the extra assumption

(1.6) f0

j(0)s − Ms1+α ≤ fj(s) ≤ fj0(0)s, ∀ s ∈ [0, 1], j ∈ Z,

(2)

2

for some constants M > 0 and α ∈ (0, 1). Recently, one of the authors and Hamel [10] gave a different approach to prove the existence of traveling waves for all speeds c ≥ c∗ for some positive minimal speed c∗. Moreover, it is also shown in [10] that the condition c ≥ c is not only a sufficient condition but also a necessary condition for the existence of traveling waves. For reader’s convenience, we recall some properties of traveling wave from [10]. Let (c, U ) be a traveling wave solution of (1.1)-(1.3) with c 6= 0. Then we have 0 < Uj(t) < 1 for all (j, t) ∈ Z × R; Uj(t) → 0 as t → −∞; Uj(t) → 1 as t → ∞; Uj0(t) > 0 for all t ∈ R and

U0

j(t) → 0 as t → ±∞.

The aim of this paper is to study the uniqueness and stability of traveling waves in the periodic monostable case. Hence we shall always assume that (1.4), (1.5) and (1.6) hold.

Recall from [10] that for each λ ∈ R there exists a unique v = {vj} with maxj∈Zvj = 1 and vj+N = vj > 0 for all j ∈ Z such that

(1.7) M(λ)vj = dj+1e−λvj+1+ djeλvj−1− (dj+1+ dj)vj + fj0(0)vj

for all j ∈ Z, where M(λ) is the largest eigenvalue of (1.7). Moreover, there exists λ∗ > 0 such that c∗ = M(λ)/λ and the mapping c = M(λ)/λ : (0, λ) 7→ c ∈ (c, ∞) is strictly decreasing.

We shall focus our attention on those traveling waves (c, U ), c > c∗, of (1.1)-(1.3) satisfying

(1.8) lim j−ct→∞ Uj(t) e−λ(j−ct)v j = 1,

for some λ > 0 such that M(λ) = cλ and {vj} is the unique eigenvector of (1.7) corresponding to λ such that maxj∈Zvj = 1 and vj+N = vj > 0 for all j ∈ Z.

We now state our stability theorem as follows.

Theorem 1.1. Suppose that there exists a traveling wave (c, U) with c > c∗ such that (1.8)

holds for some λ ∈ (0, λ∗). Let u be the solution of (1.1) for t ≥ 0 with the initial value

{uj(0)} satisfying 0 ≤ uj(0) ≤ 1, uj(0) ≤ e−λ·jvj ∀j ∈ Z, (1.9) lim inf j→−∞ uj(0) > 0, j→∞lim uj(0) e−λ·jv j = 1. (1.10) Then lim

t→∞supj {|[uj(t)/Uj(t)] − 1|} = 0.

The proof of Theorem 1.1 is based on a method in [3] with some nontrivial modifications. In [3], a lattice dynamical system in homogeneous media is studied. There the proof of stability theorem is through a related continuum equation by extending the spatial variable from j ∈ Z to x ∈ R. But, here we shall only use the original equation (1.1) to prove the stability theorem. Moreover, there is only one wave profile for the homogeneous case in [3]. In our periodic lattice dynamical system, there are N wave profiles. This makes the stability analysis more complicated. To overcome this difficulty, we introduce the following transformation

(1.11) Wj(x) := Uj([j − x]/c), equivalently Uj(t) = Wj(j − ct),

which is very useful in the periodic framework. Indeed, this transformation is reminiscent of a similar transformation in the case of partial differential equation (cf. [9]).

(3)

3

By adapting a method used in [4], we have the following uniqueness theorem.

Theorem 1.2. Suppose that (c, U ) and (c, U) are two traveling wave solutions of (1.1)-(1.3)

such that (1.12) lim j−ct→∞ Uj(t) e−λ(j−ct)v j = h, lim j−ct→∞ Uj(t) e−λ(j−ct)v j = ¯h

for some positive constants λ, h and ¯h such that M(λ) = cλ, where {vj} is the eigenvector

of (1.7) corresponding to λ such that vj = vj+N > 0 for all j and max{vj} = 1. Then there

exists ξ ∈ R such that Uj(t) = Uj(t + ξ) for all j ∈ Z, t ∈ R.

This paper is organized as follows. We shall give the proof of Theorem 1.1 in Section 2. The proof of Theorem 1.2 is given in Section 3. In this paper, we shall use both functions

Uj and Wj define in (1.11) alternatively from time to time.

References

[1] P.W. Bates, Xinfu Chen, A. Chmaj, Traveling waves of bistable dynamics on a lattice, SIAM J. Math. Anal. 35 (2003), 520–546.

[2] Xinfu Chen, S.-C. Fu, J.-S. Guo, Uniqueness and asymptotics of traveling waves of monostable dy-namics on lattices, SIAM J. Math. Anal. 38 (2006), 233–258.

[3] Xinfu Chen, J.-S. Guo, Existence and asymptotic stability of traveling waves of discrete quasilinear monostable equations, J. Differential Equations 184 (2002), 549–569.

[4] Xinfu Chen, J.-S. Guo, Uniqueness and existence of traveling waves for discrete quasilinear monos-table dynamics, Math. Ann. 326 (2003), 123–146.

[5] Xinfu Chen, J.-S. Guo, C.-C, Wu, Traveling waves in discrete periodic media for bistable dynamics, Arch. Rational Mech. Anal. (to appear).

[6] S.-N. Chow, J. Mallet-Paret, W. Shen, Stability and bifurcation of traveling wave solution in coupled map lattices, Dynam. Systems Appl. 4 (1995), 1–26.

[7] S.N. Chow, J. Mallet-Paret, W. Shen, Travelling waves in lattice dynamical systems, J. Differential Equations 149 (1998), 249–291.

[8] P.C. Fife, Mathematical aspects of reacting and diffusing systems, Lecture Notes in Biomathematics 28, Springer Verlag, 1979.

[9] P.C. Fife, J.B. McLeod, The approach of solutions of non-linear diffusion equations to traveling front solutions, Arch. Rational Mech. Anal. 65 (1977), 335–361.

[10] J.-S. Guo, F. Hamel, Front propagation for discrete periodic monostable equations, Math. Ann. 38 (2006), 233–258.

[11] W. Hudson, B. Zinner, Existence of traveling waves for a generalized discrete Fisher’s equation, Comm. Appl. Nonlinear Anal. 1 (1994), 23–46.

[12] W. Hudson, B. Zinner, Existence of travelling waves for reaction-dissusion equations of Fisher type in periodic media, In: Boundary Value Problems for Functional-Differential Equations, J. Henderson (ed.), World Scientific, 1995, pp. 187–199.

[13] J. Mallet-Paret, The Fredholm alternative for functional-differential equations of mixed type, J. Dy-nam. Differential Equations 11 (1999), 1–47.

[14] J. Mallet-Paret, The global structure of traveling waves in spatially discrete dynamical systems, J. Dynam. Differential Equations 11 (1999), 49–127.

[15] N. Shigesada, K. Kawasaki, Biological invasions: theory and practice, Oxford Series in Ecology and Evolution, Oxford, Oxford University Press, 1997.

[16] B. Shorrocks, I.R. Swingland, Living in a Patch Environment, Oxford University Press, New York, 1990.

[17] J. Wu, X. Zou, Asymptotic and periodic boundary values problems of mixed PDEs and wave solutions of lattice differential equations, J. Differential Equations 135 (1997), 315–357.

(4)

4

[18] B. Zinner, Stability of traveling wavefronts for the discrete Nagumo equations, SIAM J. Math. Anal. 22 (1991), 1016–1020.

[19] B. Zinner, Existence of traveling wavefront solutions for the discrete Nagumo equation, J. Differential Equations 96 (1992), 1–27.

[20] B. Zinner, G. Harris, W. Hudson, Traveling wavefronts for the discrete Fisher’s equation, J. Differ-ential Equations 105 (1993), 46–62.

參考文獻

相關文件

In the inverse boundary value problems of isotropic elasticity and complex conductivity, we derive estimates for the volume fraction of an inclusion whose physical parameters

Wang, Unique continuation for the elasticity sys- tem and a counterexample for second order elliptic systems, Harmonic Analysis, Partial Differential Equations, Complex Analysis,

The existence of transmission eigenvalues is closely related to the validity of some reconstruction methods for the inverse scattering problems in an inhomogeneous medium such as

Quadratically convergent sequences generally converge much more quickly thank those that converge only linearly.

“Find sufficiently accurate starting approximate solution by using Steepest Descent method” + ”Compute convergent solution by using Newton-based methods”. The method of

denote the successive intervals produced by the bisection algorithm... denote the successive intervals produced by the

化成 reduced echelon form 後由於每一個 row 除了該 row 的 pivot 外, 只剩 free variables (其他的 pivot variable 所在的 entry 皆為 0), 所以可以很快地看出解的形式..

Assume that the partial multiplicities of purely imaginary and unimodular eigenvalues (if any) of the associated Hamiltonian and symplectic pencil, respectively, are all even and