REDUCEによる有理型差分方程式の保存量の計算.pdfVIP

REDUCEによる有理型差分方程式の保存量の計算 高 　敏 (Min Gao) ? 加藤泰幸 (Yasuyuki Kato) ? 広島大学大学院工学研究科 広島大学大学院工学研究科 伊藤雅明 (Masaaki Ito) ? 広島大学大学院工学研究科 1 はじめに 差分方程式の保存量を求めることは、その方程式自身の可積分性を知る上で重要であるばかりでなく、現 実の現象を記述する微分方程式の差分スキームとしての安定性を議論する上でも欠かせないものである。多 項式型の微分方程式や差分方程式に対しては計算機を用いて具体的な多項式型の保存量を求めるためのア ルゴリズム及びその数式処理プログラムが幾つか開発されている。また、多項式型の微分方程式を差分化 したからといって、その差分方程式や保存量は多項式型になるとは限らず、有理型になることもある。本研 究においては一変数有理型差分方程式の保存量の性質を調べ、その性質を用いて具体的な保存量を求める ための数式処理プログラムを数式処理システム（REDUCE）上にインプリメントした。開発したプログラ ムは、与えられた方程式の保存量を自動的に求めるものである。 2 有理型の保存量 次の非線形振動方程式 d2 x = ?αx ? βx3, 　　α, βは正の定数, 　 (1) dt2 の保存量 H は次式で与えられることが知られている、 1 dx 1 β H = ( )2 + αx2 + x4. (2) 2 dt 2 4 Hirota による差分化を用いると、（1）式は次のように差分化される x(t + δ) + x(t ? δ) 2x(t)A0 2 = 2 2 . (3) 1 ? βδ b1x(t + δ)x(t ? δ) 1 + β b2x(t) 2 2 ここで、A0 = (1 ? 2αδ a2)/(1 + 2αδ a1), a1 + a2 = 1/2, b1 + b2 = 1 である。 x(t ± δ) = xm±1 を用いて（3）式を書き換えると次式が得られる， 2 2 2A0xm ? xm?1 ? b2βδ xmxm?1 xm+1 = . (4) 2 2 2 1 ? 2A0b1βδ xmxm?1 ? b2βδ xm ?gaomin@amath.hiroshima-u.ac.jp ?ykato@amath.hiroshima-u.ac.jp ?ito@amath.hiroshima-u