arxiv:solv-int/ v1 4 Aug 1997

Similar documents
Section 6.1 INTRO to LAPLACE TRANSFORMS

AN INEQUALITY OF OSTROWSKI TYPE AND ITS APPLICATIONS FOR SIMPSON S RULE AND SPECIAL MEANS. I. Fedotov and S. S. Dragomir

Coalgebra, Lecture 15: Equations for Deterministic Automata

Review of Gaussian Quadrature method

WENJUN LIU AND QUÔ C ANH NGÔ

Calculus of variations with fractional derivatives and fractional integrals

INEQUALITIES FOR TWO SPECIFIC CLASSES OF FUNCTIONS USING CHEBYSHEV FUNCTIONAL. Mohammad Masjed-Jamei

QUADRATURE is an old-fashioned word that refers to

Some Hermite-Hadamard type inequalities for functions whose exponentials are convex

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

Czechoslovak Mathematical Journal, 55 (130) (2005), , Abbotsford. 1. Introduction

The area under the graph of f and above the x-axis between a and b is denoted by. f(x) dx. π O

1 The Riemann Integral

Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004

SOME INTEGRAL INEQUALITIES OF GRÜSS TYPE

Definite Integrals. The area under a curve can be approximated by adding up the areas of rectangles = 1 1 +

Numerical Analysis: Trapezoidal and Simpson s Rule

The Form of Hanging Slinky

Section 6.1 INTRO to LAPLACE TRANSFORMS

A basic logarithmic inequality, and the logarithmic mean

DEFINITION The inner product of two functions f 1 and f 2 on an interval [a, b] is the number. ( f 1, f 2 ) b DEFINITION 11.1.

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Farey Fractions. Rickard Fernström. U.U.D.M. Project Report 2017:24. Department of Mathematics Uppsala University

Quadratic Forms. Quadratic Forms

Lecture 3: Curves in Calculus. Table of contents

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1

The Regulated and Riemann Integrals

WHEN IS A FUNCTION NOT FLAT? 1. Introduction. {e 1 0, x = 0. f(x) =

Introduction to Group Theory

1 The Lagrange interpolation formula

ON COMPANION OF OSTROWSKI INEQUALITY FOR MAPPINGS WHOSE FIRST DERIVATIVES ABSOLUTE VALUE ARE CONVEX WITH APPLICATIONS

S. S. Dragomir. 2, we have the inequality. b a

New Expansion and Infinite Series

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes

Improper Integrals, and Differential Equations

Solution to Fredholm Fuzzy Integral Equations with Degenerate Kernel

Lecture 3. In this lecture, we will discuss algorithms for solving systems of linear equations.

A product convergence theorem for Henstock Kurzweil integrals

Improper Integrals. The First Fundamental Theorem of Calculus, as we ve discussed in class, goes as follows:

THE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM

Lecture 20: Numerical Integration III

Chapter 6 Techniques of Integration

Numerical Linear Algebra Assignment 008

Linear Inequalities. Work Sheet 1

ON ALTERNATING POWER SUMS OF ARITHMETIC PROGRESSIONS

1B40 Practical Skills

The presentation of a new type of quantum calculus

Bases for Vector Spaces

The Shortest Confidence Interval for the Mean of a Normal Distribution

Application of Exp-Function Method to. a Huxley Equation with Variable Coefficient *

Torsion in Groups of Integral Triangles

Parse trees, ambiguity, and Chomsky normal form

1.9 C 2 inner variations

STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS

p-adic Egyptian Fractions

Lecture 3 ( ) (translated and slightly adapted from lecture notes by Martin Klazar)

Some circular summation formulas for theta functions

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

A Generalized Inequality of Ostrowski Type for Twice Differentiable Bounded Mappings and Applications

dx dt dy = G(t, x, y), dt where the functions are defined on I Ω, and are locally Lipschitz w.r.t. variable (x, y) Ω.

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?

Topics Covered AP Calculus AB

Lecture 14: Quadrature

MATH1050 Cauchy-Schwarz Inequality and Triangle Inequality

u(t)dt + i a f(t)dt f(t) dt b f(t) dt (2) With this preliminary step in place, we are ready to define integration on a general curve in C.

Math 554 Integration

ON THE C-INTEGRAL BENEDETTO BONGIORNO

Lecture 09: Myhill-Nerode Theorem

10 Vector Integral Calculus

Theory of Computation Regular Languages. (NTU EE) Regular Languages Fall / 38

Heavy tail and stable distributions

CLOSED EXPRESSIONS FOR COEFFICIENTS IN WEIGHTED NEWTON-COTES QUADRATURES

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

A Companion of Ostrowski Type Integral Inequality Using a 5-Step Kernel with Some Applications

Lecture 1. Functional series. Pointwise and uniform convergence.

GENERALIZED ABSTRACTED MEAN VALUES

P 1 (x 1, y 1 ) is given by,.

S. S. Dragomir. 1. Introduction. In [1], Guessab and Schmeisser have proved among others, the following companion of Ostrowski s inequality:

Math& 152 Section Integration by Parts

ODE: Existence and Uniqueness of a Solution

Polynomials and Division Theory

Matrix Algebra. Matrix Addition, Scalar Multiplication and Transposition. Linear Algebra I 24

Lecture 08: Feb. 08, 2019

CMPSCI 250: Introduction to Computation. Lecture #31: What DFA s Can and Can t Do David Mix Barrington 9 April 2014

Improvements of some Integral Inequalities of H. Gauchman involving Taylor s Remainder

1 Nondeterministic Finite Automata

COSC 3361 Numerical Analysis I Numerical Integration and Differentiation (III) - Gauss Quadrature and Adaptive Quadrature

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

Review on Integration (Secs ) Review: Sec Origins of Calculus. Riemann Sums. New functions from old ones.

Things to Memorize: A Partial List. January 27, 2017

MAA 4212 Improper Integrals

Linearly Similar Polynomials

df dt f () b f () a dt

Research Article On Existence and Uniqueness of Solutions of a Nonlinear Integral Equation

arxiv:nlin/ v1 [nlin.si] 25 Jun 2004

4.1. Probability Density Functions

Theory of Computation Regular Languages

The Riemann Integral

Closure Properties of Regular Languages

Linear Systems with Constant Coefficients

Transcription:

Two-dimensionl soliton cellulr utomton of deutonomized Tod-type A. Ngi 1,2, T. Tokihiro 1, J. Stsum 1, R. Willox 1,3 nd K. Kiwr 4 1 Grdute School of Mthemticl Sciences, University of Tokyo, rxiv:solv-int/9708001v1 4 Aug 1997 Kom 3-8-1, Meguro-ku, Tokyo 153, Jpn 3 Dienst Ten, Vrie Universiteit Brussel, Pleinln 2, 1050 Brussels, Belgium 4 Deprtment of Electricl Engineering, Doshish University, Tne, Kyoto 610-03, Jpn (Received ) Astrct A deutonomized version of the two-dimensionl Tod lttice eqution is presented. Its ultr-discrete nlogue nd soliton solutions re lso discussed. PACS: 03.20.+i; 03.40.Kf; 04.60.Nc Keywords: deutonomized, ultr-discrete, two-dimensionl Tod lttice eqution. 2 E-mil: ngi@st.t.u-tokyo.c.p to pper in Phys. Lett. A.

1 Introduction The importnce of discrete soliton equtions hs een recognized in mny fields such s mthemtics, physics nd engineering, owing to the enormous development of computer science. Quite recently, two new types of discretiztion hve een proposed. One of these is to unequlize lttice intervls [1], which we cll deutonomiztion in this pper. This ide comes from Newton s interpoltion formul [2] nd now mny pplictions in numericl computtion such s in the ρ-lgorithm nd in dptive numericl integrtion re expected. Another type of discretiztion is to discretize dependent s well s independent vriles, known s ultr-discretiztion. One of the most importnt ultr-discrete soliton systems is the so-clled soliton cellulr utomton or SCA for short [3, 4]. In recent ppers [5, 6], generl method to otin the SCA from discrete soliton equtions ws proposed. This procedure hs een pplied to other soliton equtions [7] s well. The im of this pper is to study deutonomized version of the two-dimensionl(2d) Tod lttice eqution nd its ultr-discrete nlogue, which is considered s deutonomized version of the 2D Todtype SCA [7]. In section 2, we review the result y Hirot et l. on the 2D Tod lttice eqution [8] nd derive its deutonomized version. In section 3, we present n ultr-discrete nlogue of the deutonomized 2D Tod lttice eqution, soliton solutions for which re lso discussed. Concluding remrks re given in section 4. 2 Deutonomiztion of the 2D Tod lttice eqution We strt with discrete nlogue of the 2D Tod lttice eqution in iliner form [8], ( + l m τ(l,m,n))τ(l,m,n) ( + l τ(l,m,n)) m τ(l,m,n) = τ(l,m 1,n+1)τ(l +1,m,n 1) τ(l+1,m 1,n)τ(l,m,n), l,m,n Z, (1) where + l nd m stnd for difference opertors defined y + l τ = τ(l +1,m,n) τ(l,m,n) δ Rewriting eq. (1) y using eq. (2), we otin, mτ = τ(l,m,n) τ(l,m 1,n). (2) ǫ (1 δǫ)τ(l +1,m 1,n)τ(l,m,n) τ(l +1,m,n)τ(l,m 1,n) +δǫτ(l,m 1,n+1)τ(l +1,m,n 1) = 0, (3) which possesses prticulr solution expressed s [8] τ(l,m,n) = f (1) (l,m,n) f (1) (l,m,n+1) f (1) (l,m,n+n 1) f (2) (l,m,n) f (2) (l,m,n+1) f (2) (l,m,n+n 1)... f (N) (l,m,n) f (N) (l,m,n+1) f (N) (l,m,n+n 1). (4) In eq. (4) ech entry f (i) (l,m,n) (i = 1,2,,N) stisfies the following dispersion reltions; { + l f(i) (l,m,n) = f (i) (l,m,n+1), m f(i) (l,m,n) = f (i) (l,m,n 1). (5) 1

When we tke f (i) (l,m,n) = p n i (1+δp i ) l (1+ǫp 1 i ) m +q n i (1+δq i ) l (1+ǫq 1 i ) m (6) s solution for eq. (5), the τ-function (4) gives n N-soliton solution. As opposed to the utonomous eq. (1), where there is constnt lttice intervl in ech direction, we now present deutonomized version, where lttice intervl etween neighoring lttice points cn e chosen freely. Let us replce f (i) (l,m,n) in eq. (6) y (1+δ p i ) f (i) (l t,m,n) = p n i (1+ǫ p 1 i ) +qi n (1+δ q i ) (1+ǫ q 1 i ), (7) δ t = l t+1 l t, ǫ = m +1 m, t, Z, (8) where the products represent, for exmple, (1+δ p i ) = (1+δ p i ) t 1, =0 1 t = 0, 1 (1+δ p i ) 1 t 1. =t (9) Then the following dispersion reltions hold; + l t f (i) (l t,m,n) f(i) (l t+1,m,n) f (i) (l t,m,n) δ t m f (i) (l t,m,n) f(i) (l t,m,n) f (i) (l t,m,n) ǫ = f (i) (l t,m,n+1), = f (i) (l t,m,n 1). (10) By using the Plücker identity for determinnts, we see tht the τ-function given y f (1) (l t,m,n) f (1) (l t,m,n+1) f (1) (l t,m,n+n 1) f (2) (l t,m,n) f (2) (l t,m,n+1) f (2) (l t,m,n+n 1) τ(l t,m,n) =, (11)... f (N) (l t,m,n) f (N) (l t,m,n+1) f (N) (l t,m,n+n 1) with ech f (i) solution of eq. (10), stisfies iliner eqution, or equivlently, ( + l t m τ(l t,m,n))τ(l t,m,n) + l t τ(l t,m,n) m τ(l t,m,n) = τ(l t,m,n+1)τ(l t+1,m,n 1) τ(l t+1,m,n)τ(l t,m,n), (12) (1 δ t ǫ )τ(l t+1,m,n)τ(l t,m,n) τ(l t+1,m,n)τ(l t,m,n) +δ t ǫ τ(l t,m,n+1)τ(l t+1,m,n 1) = 0, (13) which is considered s deutonomized 2D Tod lttice eqution (in the sense explined ove). 2

It is interesting to note tht when we replce f (i) (l,m,n) y f (i) (l t,m,n k ) = k 1 c (η c p i ) (1+δ p i ) (1+ǫ p 1 i ) (1+δ q i ) (η c q i ) c (1+ǫ q 1 i ) k 1 +, (14) η c = n c+1 n c, (15) we cn construct new iliner eqution deutonomized with respect to ll independent vriles. However, it reduce to eq. (12) y simple dependent vrile trnsformtion. 3 Ultr-discretiztion of the deutonomized 2D Tod lttice eqution In this section, we present n ultr-discrete nlogue of the deutonomized Tod lttice eq. (13) nd we construct its soliton solutions. Let us denote τ(l t,m,n) s τ t,,n for simplicity. Eqution (13) is rewritten s (1 δ t ǫ )τ t+1,,n τ t,+1,n τ t+1,+1,n τ t,,n +δ t ǫ τ t,,n+1 τ t+1,+1,n 1 = 0. (16) When we introduce the new dependent vrile S t,,n s eq. (16) is equivlent to τ t,,n = exp[s t,,n ], (17) (1 δ t ǫ )+δ t ǫ exp [ e n ( + n + t )( + n + )S t,,n] = exp[ + t + S t,,n], (18) Let us define difference opertor y + t = e t 1, + = e 1, + n = e n 1. (19) = e n ( + n + t )( + n + ) (20) for convenience ske. Tking logrithm nd operting on oth sides of eq. (18), we hve [ log(1 δ t ǫ )+ log 1+ δ ] tǫ exp( S t,,n ) = + t 1 δ t ǫ + S t,,n. (21) We finlly tke n ultr-discrete limit of eq. (21). Setting we otin S t,,n = vε t,,n, δ t = e θt/ε, ǫ = e σ/ε (θ t,σ Z 0 ), (22) ε + t + v t,,n = e n ( + n + t )( + n + )F(v t,,n θ t σ ), (23) F(x) = mx(0,x) (24) in the limit ε +0. We hve rewritten lim ε +0 vε t,,n s v t,,n. Eqution (23) is considered s n ultrdiscrete nlogue of the deutonomized 2D Tod lttice eqution. 3

Next, we discuss soliton solutions for eq. (23). It is nturl to consider tht soliton solutions for eq. (23) re otined y n ultr-discretiztion of those for eq. (12). First of ll, one-soliton solution for eq. (12) is given y τ t,,n = 1+η 1, (25) ) n η 1 = α 1 ( p1 q 1 Introducing new prmeters nd vrile s nd tking the limit ε +0, we otin 1+δ p 1 1+δ q 1 1+ǫ q 1 1 1+ǫ p 1 1. (26) e P1/ε = p 1, e Q1/ε = q 1, e A1/ε = α 1, (27) e θt/ε = δ t, e σ/ε = ǫ, (28) ρ ε t,,n = εlogτ t,,n (29) lim ε +0 ρε t,,n ρ t,,n = mx(0,k 1 ) (30) K 1 = A 1 +n(p 1 Q 1 )+ {F(P 1 θ ) F(Q 1 θ )} = + {F( Q 1 σ ) F( P 1 σ )}, (31) =0 A one-soliton solution for eq. (23) is given y t 1, 0 t = 0, 1 t 1. =t v t,,n = ρ t+1,+1,n 1 +ρ t,,n+1 ρ t+1,,n ρ t,+1,n (33) = mx(0,k 1 P 1 +Q 1 +F(P 1 θ t ) F(Q 1 θ t ) F( P 1 σ )+F( Q 1 σ )) +mx(0,k 1 +P 1 Q 1 ) mx(0,k 1 +F(P 1 θ t ) F(Q 1 θ t )) mx(0,k 1 F( P 1 σ )+F( Q 1 σ )). (34) Secondly, we construct two-soliton solution. Eqution (12) dmits two-soliton solution, τ t,,n = 1+η 1 +η 2 +θ 12 η 1 η 2, (35) ) n η i = α i ( pi q i 1+δ p i 1+δ q i 1+ǫ q 1 i 1+ǫ p 1 i (32) (i = 1,2), (36) θ 12 = (p 2 p 1 )(q 1 q 2 ) (q 1 p 2 )(q 2 p 1 ). (37) In order to tke n ultr-discrete limit for the ove solution, we suppose without loss of generlity tht the inequlity, 0 < p 1 < p 2 < q 2 < q 1 (38) 4

holds. Introducing new prmeters nd vrile s nd tking the limit of smll ε, we otin e Pi/ε = p i, e Qi/ε = q i, e Ai/ε = α i (i = 1,2, P 1 < P 2 < Q 2 < Q 1 ), (39) ρ ε t,,n = εlogτ t,,n (40) lim ε +0 ρε t,,n ρ t,,n = mx(0,k 1,K 2,K 1 +K 2 +P 2 Q 2 ), (41) K i = A i +n(p i Q i )+ {F(P i θ ) F(Q i θ )} + {F( Q i σ ) F( P i σ )} (i = 1,2). (42) Finlly, n N-soliton solution cn lso e constructed, under the ssumption tht the inequlity, holds. Through the sme limiting procedure, we hve N ρ t,,n = mx µ i K i + 0 < p 1 < p 2 < < p N < q N < < q 2 < q 1 (43) µ=0,1 i=1 1 i<i N µ i µ i (P i Q i ), (44) where mx is the mximiztion over ll possile comintions of µ 1 = 0,1, µ 2 = 0,1,, µ N = 0,1. µ=0,1 Since the ounded nd uniform convergenceof the dependent vrile ρ ε t,,n s ε +0 is the cornerstone of the ultr-discretiztion, it should e cler tht eq. (44), which is the ultr-discretiztion of the N- soliton solution for the discrete 2D Tod lttice eqution, is n N-soliton solution for the ultr-discrete eqution. We show one- nd two-soliton solutions for eq. (23) with σ,θ t constnt in Figures 1,3 nd with σ,θ t chosen rndomly in Figures 2,4. It should e noted tht the solutions with σ,θ t constnt re equivlent to those discussed in [7]. We see tht ritrriness of σ nd θ t ffects the vlues of the dependent vrile v n,t,. 5

t=0...451...154......253......352......451......154... n...253......352......451......154......253....352... t=5...154.....253......352......451......154......253... n...352......451......154......253.....352... 451... t=10...253......352......451......154......253......352... n...451......154......253......352....451... 4... t=15...352......451......154......253......352......451... n...154......253......352.....451... 54... t=20...451......154......253......352......451......154... n...253......352......451... 154... 3... Figure 1: One-soliton solution for eq. (23) (P 1 = 5,Q 1 = 4,σ = 1,θ t = 3). The vlue 0 is replced y. (period) for simplicity. 6

t=0...144......154......143......144......154......143... n...154......43......44......44......451....45... t=5...252......142......152......252......142......252... n...142......42......43......54.....53... 43... t=10...55......44......45......55......44......55... n...44......34......341......352....351... 41... t=15...153......142......143......153......142......153... n...42......43......44......44....44... 4... t=20...143......153......253......143......253......143... n...43......44......55.....54... 44... Figure 2: One-soliton solution for eq. (23) (P 1 = 5,Q 1 = 4,σ,θ t : rndomly chosen). 7

t=0...45...154.....253......352......451...123...154...12344321...1235312344321... n...1234434544321....12344321.3521... 4321...451......154......253....352... 51... t=5...4...154....253......352......451......154...123...253...12344321...1352.12344321... n...1234453344321....1234432145321... 4321..154......253.....352... 451... t=10...154...253.....352......451......154......253...123...352...12344321...451..12344321... n...1235412344321....1234433544321... 4321.2531......352....451... 4... t=15...253....352......451......154......253......352...123...451...12344321...154...12344321... n...1253.12344321....1234454344321... 432135321.....451... 54... Figure 3: Two-soliton solution for eq. (23) (P 1 = 5,P 2 = 2,Q 1 = 4,Q 2 = 3,σ = 1,θ t = 3) 8

t=0...4...252....341......263......24... n...34...2223...133..133333111.....23241332.....11233331.522... 3221...141......242......441....262... t=5...5...361....143......352......33......43...2223 n...142..133333111.....23331332.....112333311522... 3221...24......341......153....351... 3... t=10...16...461....253......452......43......153...2224...242..134444222.. n...244313431.....1123334225331... 4332111.34......1441......253... 1451... 3... t=15...163.....145......252......24......34...11113333...43..1232222......1133222321... n 112233322.4311... 211...52......361.....143... 352... Figure 4: Two-soliton solution for eq. (23) (P 1 = 5,P 2 = 2,Q 1 = 4,Q 2 = 3,σ,θ t : rndomly chosen) 9

4 Concluding Remrks We hve presented deutonomized 2D Tod lttice eqution nd its ultr-discrete nlogue. We hve lso found soliton solutions for the ltter cse. The deutonomiztion resulting in eq. (23) cuses ritrriness in the vlues of the soliton solutions, property which is not oserved in the utonomous cse. Deutonomiztion nd ultr-discretiztion of discrete soliton equtions re expected to find more pplictions in the field of computer science such s convergence ccelertion, interpoltion nd sorting lgorithms. Acknowledgements It is plesure to thnk Professor Ryogo Hirot of Wsed University for fruitful discussions, especilly on deutonomized soliton equtions. One of the uthors(r.w.) is postdoctorl fellow of the Fund for Scientific Reserch(F. W. O.), Flnders (Belgium). He would lso like to cknowledge the support of the F. W. O. through moility grnt. The present work is prtilly supported y Grnt-in-Aid from the Jpn Ministry of Eduction, Science nd Culture. References [1] R. Hirot, J. Phys. Soc. Jpn. 67 (1997) 283. [2] L. M. Milne-Thomson, The Clculus of Finite Differences, Chelse, New York, 1933. [3] D. Tkhshi nd J. Stsum, J. Phys. Soc. Jpn. 59 (1990) 3514. [4] D. Tkhshi nd J. Mtsukidir, Phys. Lett. A 209 (1995) 184. [5] T. Tokihiro, D. Tkhshi, J. Mtsukidir nd J. Stsum, Phys. Rev. Lett. 76 (1996) 3247. [6] J. Mtsukidir, J. Stsum, D. Tkhshi, T. Tokihiro nd M. Torii, Phys. Lett. A 225 (1997) 287. [7] S. Moriwki, A. Ngi, J. Stsum, T. Tokihiro, M. Torii, D. Tkhshi nd J. Mtsukidir, to pper in London Mth. Soc. Lecture Notes Series, Cmridge Univ. Press. [8] R. Hirot, M. Ito nd F. Kko, Prog. Theor. Phys. Suppl. 94 (1988) 42. [9] R. Hirot, Direct Method in Soliton Theory (Iwnmi, Tokyo, 1992) [in Jpnese]. 10