［收稿日期］ ２００７ ０３ １２
［作者简介］ 刘亦斌 （１９７０ ） ， 男 ， １９９２ 大学毕业 ， 讲师 ， 现主要从事计算机与高等数学方面的教学与研究工作 。
基 于 Ｍａｐｌｅ 系 统 的 Ｒｉｃｃａｔｉ 方 程 自 动 求 解 的 实 现
刘亦斌
（泉州师范学院高职学院 ，福建 泉州 ３６２０００）
［摘要］ 利用非线性演化方程的精确行波解的手工推导和计算的原理 ， 结合计算机符号计算 ，在计算机代
数系统 Ｍ ａｐｌｅ 上 ，开发了基于非线性代数方程组吴文俊消元法的非线性演化方程精确行波解的自动求解
软件包 ，实现了非线性方程 （组） 的求解过程的完全自动化 。
［关键词］ 非线性演化方程 （组） ；行波解 ； 计算机代数系统 Ｍ ａｐｌｅ ； Ｒｉｃｃａｔｉ 方程
［中图分类号］ ＴＰ３１９ ［文献标识码］ Ａ ［文章编号］ １６７３ １４０９ （２００７） ０２ Ｎ０１１ ０４ Ｎ
目前对基于单个 Ｒｉｃｃａｔｉ 方程 、 耦合的 Ｒｉｃｃａｔｉ 方程和三耦合的 Ｒｉｃａｔｔｉ 方程构造非线性演化方程的精
确行波解的计算 ， 往往采用作行波变换 ， 再通过导数代换 ，利用吴文俊消元法来求出行波解
［１］
。 这种手
工推算的代数方法计算繁冗 ， 使得该方法的应用受到很大的局限性 。
要实现基于 Ｒｉｃｃａｔｉ 方程的非线性演化方程 （组） 精确行波解的计算机自动推导 ， 有 ３ 个关键问题
需待解决 。首先要自动推导出行波解的阶数 ｍ ；其次要有效地求解待定系数 ａ
ｉ
等 ，待定常数 ｋ ，ｃ 和方程
（组） 中的
ｐ
，
ｑ
，ｒ ，ｓ 等的非线性代数方程组 ；最后要保证解集的最小化 ， 即自动判断并删除所有的特解 。
Ｄｍａｉｎ 由以下 ６ 个模块构成 ： ① 主模块 Ｄｍａｉｎ （） ， 完成一些初始化工作以及对其他几个子模块的
调用 ； ② Ｆｄｌｓｔ （） ， 根据方程中函数的自变量产生行波变量 ， 同时将偏微分方程 （组） 转化为常微分方
程 （组） ； ③ ｆｉｎｄ ＿ ｍ ＿ ｎ （） ， 通过平衡方程组中各个方程的线性最高阶导数项和最高阶非线性的阶数 ，
确定行波解的阶数 ； ④ ｃｏｅｆｆ （） ， 导出行波解系数所满足的非线性代数方程组 ； ⑤ ｓｏｌ （） ， 简化和分组
非线性代数方程组 ， 利用吴文俊消元法求解 ， 并根据用户的要求去除所有的特解 ； ⑥ Ｐｒｉｎｔ （） ， 删除特
解及重复解 ， 并依特定的格式输出最终的行波解 。 将 Ｄｍａｉｎ 制成 Ｍ ａｐｌｅ 上的一个共享包 “Ｄｍａｉｎ畅 ｍ” ，
通过调用 Ｄｍａｉｎ （参数 １ ， 参数 ２ ， 参数 ３） 求解非线性演化方程 （组） 。
１ 方程 （组） 参数的输入
主模块 Ｄｍａｉｎ（）的参数 １ 要求以列表的形式给出 ，即方程列表 。 列表的每个方程必须以等式的形式
给出 。 参数 ２ 是用来控制是否保留平凡解（在程序中由 ｃ＿ｋｅｙ 控制） 、是否保留复解（由 ｃ＿
ｐ
ａｒ 控制） ，当参
数 ２ 的值为 １（即 ｃ＿ｋｅｙ ＝ ０ ，ｃ＿
ｐ
ａｒ ＝ ０）时 ，方程组的解必须去除所有平凡解和复解 ；当参数 ２ 的值为 ２（即 ｃ
＿ｋｅｙ ＝ １ ，ｃ＿
ｐ
ａｒ ＝ ０）时 ，方程组的解保留了平凡解 ，去除复解 ；当参数 ２ 的值为 ３（即＿ｋｅｙ ＝ ０ ，ｃ＿
ｐ
ａｒ ＝ １）时 ，
方程组的解保留了复解 ，去除平凡解 ；当参数 ２ 的值为 ４（即 ｃ＿ｋｅｙ ＝ １ ，ｃ＿
ｐ
ａｒ ＝ １）时 ，方程组的解保留了所
有的平凡解和复解 。 至于参数 ３（在程序中由
ｐ
＿ｋｅｙ 控制） ，当其值为 １ 时 ，最终输出的方程（组）的形式是
基于单个 Ｒｉｃｃａｔｉ 方程的 ；值为 ２ 时 ，最终输出的方程组的形式是基于耦合的 Ｒｉｃｃａｔｉ 方程的 ；值为 ３ 时 ，最
终输出的方程组的形式是基于三耦合的 Ｒｉｃｃａｔｉ 方程的 。
１畅１ 参数检查
如果错误的参数输入通过了 Ｍ ａｐｌｅ 并不很严格的参数审查 ，可能导致出现错误的结果而用户毫不知
情 。 参数的类型在变量类型声明后由 Ｍ ａｐｌｅ 系统自行检查 ，参数的数目由程序验证 。 如果出错 ，则给出
相应的错误提示信息 。
１畅２ 处理参数 １ ，即方程列表
首先对方程列表中的每个方程通分 ，并取其分子 。这样做的目的是在保持方程解不变的情况下消除
·１１·
长江大学学报 （自科版） 理工卷
２００７ 年 ６ 月 第 ４ 卷 第 ２ 期
Ｊｏｕｒｎａｌ ｏｆ Ｙａｎｇｔｚｅ Ｕｎｉｖｅｒｓｉｔｙ （Ｎａｔ Ｓｃｉ Ｅｄｉｔ） Ｓｃｉ ＆ Ｅｎｇ Ｖ Ｊｕｎ畅２００７ ， Ｖｏｌ畅４ Ｎｏ畅２
资源评论
