對(duì)反散射方法所得S-G方程解的數(shù)值模擬

2013-07-07 03:58:52曹建莉

赤峰學(xué)院學(xué)報(bào)·自然科學(xué)版 2013年1期

關(guān)鍵詞：程序方法

史偉，曹建莉

（河南工業(yè)大學(xué) 理學(xué)院，河南鄭州 450001）

對(duì)反散射方法所得S-G方程解的數(shù)值模擬

史偉，曹建莉

（河南工業(yè)大學(xué) 理學(xué)院，河南鄭州 450001）

本文采用反散射方法得到給定初始值的條件下S-G方程的解，再利用maple軟件編程選取適當(dāng)?shù)膮?shù)得到解的數(shù)值模擬.

S-G方程；反散射方法；maple軟件

S-G方程是數(shù)學(xué)物理領(lǐng)域中的熱點(diǎn)方程，反散射方法是孤立子理論中的一種經(jīng)典方法，得到了廣泛的應(yīng)用[1-3]，而與數(shù)學(xué)軟件的結(jié)合卻是近幾年的事[4].本文討論在給定初始條件q(x,t)|t=0=q(x),r(x,t)|t=0=r(x)時(shí)S-G方程的解，首先求出以初始條件中q(x)，r(x)為位勢(shì)的方程（1）的散射數(shù)據(jù)，然后求出當(dāng)q(x,t)，r(x,t)作為S-G方程解的時(shí)候，其相應(yīng)的散射數(shù)據(jù)隨時(shí)間的變化規(guī)律，最后由T-L-M積分方程組，根據(jù)散射數(shù)據(jù)來(lái)恢復(fù)q(x,t)，r(x,t)，最后用maple軟件選取適當(dāng)?shù)膮?shù)給出S-G方程解的數(shù)值模擬.

1 S-G方程

首先考察一般特征值問(wèn)題

其中q=q(x,t),r=r(x,t)，并設(shè)φ1,φ2隨時(shí)間t的演化規(guī)律為

其中A，B，C是x，t，ξ的函數(shù)，并設(shè)ξt=0,由φ1xt=φ1ix,及φ2xt=φ2tx推出A，B，C要滿(mǎn)足以下三個(gè)方程

將A，B，C取為ξ的多項(xiàng)式

代入(3)并比較ξn+1的系數(shù)，得bn=cn=0;比較ξn的系數(shù)，可得anx=0,bnx+2ibn-1=-2qan,cnx-2icn-1=2ran.由此定出an=常數(shù)及bn-1,cn-1.如此反復(fù)，即可定出aj,bj,cj,比較ξ0的系數(shù)，給出非線(xiàn)性演化方程.取n=3,用maple軟件計(jì)算得

（4），（5）中的系數(shù)a0,a1,a2,a3均為常數(shù)，適當(dāng)選取這些常數(shù)，并取q(x)，r(x)建立關(guān)系，方程可約化為：

2 利用反散射方法得到S-G方程的解定義

為S-G特征問(wèn)題的散射數(shù)據(jù).給定散射數(shù)據(jù)（8）后，再恢復(fù)S-G特征值問(wèn)題（1）中的位勢(shì)q(x)，r(x).

積分方程（10）的解未必對(duì)所有的x都是存在唯一的.例如，當(dāng)ρ(ξ)=(ξ)=0，上，下平面各有一個(gè)特征值λ(Imλ>0(Imλ<0)對(duì)應(yīng)的歸一化因子是c，，積分方程化為

用maple解代數(shù)方程得

設(shè)q(x)，r(x)分別是方程（5），（6），（7）的解，q,r及其導(dǎo)數(shù)當(dāng)|x|→∞時(shí)很快趨于零，當(dāng)t=0時(shí)，q=q(x)，r=r(x)，而且以q (x)，r(x)為位勢(shì)的特征值問(wèn)題（1）的散射數(shù)據(jù)為