关于KdVBurgersKuramoto方程精确亚纯固定解.docVIP

关于KdVBurgersKuramoto方程精确亚纯固定解 　　[摘要]寻找非线性方程的精确解是近年来非线性科学研究的一个基本课题，在这篇文章中，我们研究形如 ωu′+kνuu′+k2αu″+k3μu+k4γu″″=0 的一个常微分方程，作为KdVBurgersKuramoto方程的一个特殊形式，我们将在本文精确地求解出这条方程的精确亚纯解，最后对这些解进行计算机模拟。 　　[关键词]微分方程；精确解；亚纯函数；椭圆函数 　　1。引 言 　　关于KdVBurgersKuramoto方程是[2，3]ut+νuux+μuxxx+αuxx+γuxxxx=0，（1。1） 　　这里的u，γ，υ，α都是常数。这条方程是在许多不同的物理环境中产生的一个重要的数学模型，是用来对许多物理现象的描述，许多方法已经应用于构造KdVBurgersKuramoto方程的精确解。在2010年，Conte R和TuenWai N[5]已经对三阶非线性微分方程的亚纯通解进行的归纳，并且证明了此类方程有亚纯解u，u∈W，W={椭圆函数，有理函数，增长级为1的单周期函数}。这条方程是 　　d0u′′′+d1uu′′+d2u′2+d3u2u′+d4u4+c1u′′+c2uu′+c3u3+c4u′+c5u2+c6u+c7=0。（1。2） 　　所有的三阶方程基本都可以由方程（1。2）变化系数而得到。下面我们对方程（1。1）进行行波变换u=u（η），η=kx+ωt，那么可变为 ωu′+kνuu′+k2αu″+k3μu+k4γu″″=0。（1。3） 　　不失一般性，我们可以假设 kυ=120，k4γ=1，通过变换u→Tu，T=4120γ14υ，积分一次可以得到方程u+k3μu″+k2αu′+ωu′+60u2=c，（1。4） 　　这里的c∈ 　　瘙 ? 。本文将在第二部分介绍一种方法来表示亚纯解以及相关知识，第三部分将求解方程（1。4）的亚纯精确解，最后将对这条方程的解进行计算机模拟。 　　2。方法应用 　　在2011年，Kudryashov N A[4]引入了一个方法来寻找非线性常微分方程E[u（z）]=0的亚纯解，记E[u（z）]是关于u（z）和它的导数的一个多项式，对于方程E[u（z）]=0。极点z=0有相对应的洛朗展开式为：u（z）=∑Pk=1c-kzk+∑∞k=0ckzk，0|z|εi，（2。1） 　　这里P0是极点的重数。方程（2。1）的亚纯解可以表示成下面的引理。 　　引理2。1 方程E[u（z）]=0。的亚纯解的形式为 　　（1）周期为2ω1，2ω2的椭圆函数为：u（z）=∑Pk=2-1kc-kk-1！dk-2dzk-2 （z，2ω1，2ω2）+h0，（2。2） 　　椭圆函数解存在的必要条件为c-1=0。 　　（2）周期为T的单周期函数为： u（z）=πT∑Pk=1-1k-1c-kk-1！dk-1dzk-1cotπzT+h0。（2。3） 　　（3）有理函数解为： u（z）=∑Pk=1c-kzk+c0，（2。4） 　　这里的h0，c0是一个常数。 　　定义2。2 [1]设E（z，u）=0是一个n阶的代数微分方程，假设u（z）=∑+∞n=0un（z-z0）n+p 　　（u0≠0，p0，p∈Ζ）是方程的一个亚纯解，将u（z）代入方程，我们可以得到形如E=∑+∞j=0Ejχj+q=0，q是最小的整数，p和q是由方程的某些项来控制，这些项称为控制项，记为E∧（z，u）。对任意的v，E∧（z，u）关于u的导数 E′∧（z，u）v=limλ→0E′∧（z，u+λv）-E′∧（z，u）λ，方程P（i）=limχ→0χ-i-qE′∧（x，u0χp）χi+p=0的根我们称为方程E（z，u）=0的Fuchs指数。 　　3。亚纯解的表示 　　设u是方程（1。4）的亚纯解，利用定义2。2的方法我们可以计算出p=-3，u0=1，更计算出它的Fuchs指数是P（i）=-1，132-12i71，132+12i71，因不存在非负整数的Fuchs指数，那么它的洛朗展开式可以唯一确定，首先求出单周期函数解，根据根据引理2。1的单周期函数解的表达形式以及利用数学软件maple的各个系数的计算，可以得到方程（1。4）单周期函数解的洛朗展开式可以写成 　　u（z）=-1608πTcotπTz-z0k6μ2+138πTcotπT（z-z0）k2α-18πTcotπTz-z03k3μ+L32cotπT（z-z0）+L32cotπTz-z03-1372960k9μ3+79120k5αμ-1120ω， 　　T=13-3608k6μ2+338k28αμ2-1186640k4α2-1315544960k12μ4π，（3。1） 　　现在构建方