求解非线性偏微分方程的傅里叶变换方法

来源：网络收集时间：2026-05-21

导读：第３１卷第３期２０１１年６月广东第二师范学院学报ＪｏｕｒｎａｌｏｆＧｕａｎｄｏｎＵｎｉｖｅｒｓｉｔｏｆＥｄｕｃａｔｉｏｎｇｇｙＶｏｌ．３１Ｎｏ．３Ｊｕｎ．２０１１求解非线性偏微分方程的傅里叶变换方法李华刚１，石智伟２（１．广东第二

第３１卷　第３期

２０１１年６月广东第二师范学院学报ＪｏｕｒｎａｌｏｆＧｕａｎｄｏｎＵｎｉｖｅｒｓｉｔｏｆＥｄｕｃａｔｉｏｎ　　　　ｇｇｙ　　Ｖｏｌ．３１　Ｎｏ．３Ｊｕｎ．２０１１

求解非线性偏微分方程的傅里叶变换方法

李华刚１，石智伟２

（１．广东第二师范学院物理系，广东广州５１０３０３；

）２．广东工业大学信息工程学院，广东广州５１０００６

非线性偏微分方程数值求解在物理和数学上是一项基础工作．通过应用傅立叶变换得　　摘要：

收敛快速的迭代方法．这种迭代方法易于学生掌握和使用，能应用在ｍ到一种原理简单、ａｔｌａｂ程序

设计、数值分析、计算机辅助教学等课程教学中，有助于学生初步掌握非线性偏微分方程迭代求解

方法的学习．

关键词：非线性偏微分方程；傅立叶变换；数值解

（）中图分类号：Ｏ１７５　文献标识码：Ａ　文章编号：１００７７５４２０１１０３０４４４　－８－０－０

随着计算机性能的提高，许多非线性领域本身得到快速发展，非线性问题在许多领域成为数学、物理及

］１－４其应用上的研究热点［非线性偏微分方程一般不存在解析解，因此数值求解成为有力的研究工具．非线．

性偏微分方程的求解方法有许多，比如多重网格法、牛顿迭代法和牛顿法的改进方法等．多重网格法原理简

但是需要仔细调整试探解．牛顿迭代法和牛顿法的改进方法收敛速度比较快，但是需要具有线性代数的单，

理论基础．通过应用傅立叶变换，我们得到一种原理简单，收敛速度快的迭代方法．对于初学的学生，这是一个非常好的迭代方法，易于被学生掌握和应用．这种迭代方法能应用在ｍ数值分析、计算机ａｔｌａｂ程序设计、

辅助教学等课程教学中．

１　方法原理

傅立叶变换可将平方可积的函数表示成复指数函数的积分或级数形式：

Ｆ（ω）＝

傅立叶变换的逆变换为：∫

∞∞－∞ｉｘ－ωｘ）ｅｄｘ，ｆ（（）１ｘ）＝ｆ（

）维可写为：以非线性薛定谔方程为例，非线性薛定谔方程在（１＋１

２　２　，（）ｉｕ｜ｕ３＝｜２＋ｄｚ２ｘ

其中ｕ表示波函数，边界条件是ｕ（在无穷远处为零．假设其解的形式为ｚ为传输方向、ｘ为横向坐标．ｘ）∫－∞ｉｘωＦ（ｅｄω）ω．（）２

２－作者简介：李华刚，男，河北吴桥人，广东第二师范学院物理系讲师．

等：求解非线性偏微分方程的傅里叶变换方法　　　　　　　　　　　　　２０１１年第３期　　　　　　　　　　李华刚，·４５·

ｉｚ５］β，，）可化为即空间光孤子解［方程（ｕ（ｘ，ｚ）＝ＮＵ（ｘ）ｅ３

２　２２　（）Ｕ｜Ｕ＝０，４－βＵ＋２＋Ｎ｜２ｘ

）其中β为空间光孤子的传播常数，的振幅为１．方程（应用傅立叶变换的微分Ｎ为空间光孤子振幅和Ｕ（ｘ）４

可得性质，

Ｆ（ω）＝Ｎ２∞－∞２ｉｘ－ω（Ｕ）ｅｄｘ｜Ｕ｜　，２β＋２（）５

２２ｉｘｉｘ－ωω（］Ｕ）ｅｄｘ－（ｅｄω）ω｜Ｕ｜　

－∞－∞（），（）ｂｘ＝６Ｕ

其中Ｆ（因为傅立叶变换默认的边ｘ）的傅立叶变换，ｂ（ｘ）是Ｕ（ｘ）在横向上各点的传播常数．ω）为Ｕ（

）的边界条件相同，界条件是Ｕ（所以方程的边界条件可以不用设置．试探解的ｘ）在无穷远处为零与方程（４∞［（Ｎ２∞

迭代的过程就是使ｂ（同时空间光孤子的ｂ（ｘ）不是常数即横向各点不相同，ｘ）为一常数即函数是一条直线，

珔，珔为ｂ（波形Ｕ（所以我们在迭代过程中令β＝ｂｘ）趋向稳定解．ｂｘ）的平均值．

２　数值模拟

ｘ／２－设试探解为高斯函数Ｕ＝Ｎ首先模拟Ｎ＝１的解，即基本空间光孤子解．这种情况的解析解为ｅ．

（，／图１为解析解和数值解的对照图，从图上可以明显看出两条曲线完全重合，即数值Ｕ（ｘ）ｅｃｈｘ）２．＝ｓβ＝１２

数值模拟程序见附录）

解为正确的解（．

）解析解和数值解对照图；（）横向上各点的传播常数ｂ（；（）迭代次数与误差变化图１　（ａｂｘ）ｃ

（，其次，数值模拟Ｎ＞１的解，即基态孤子解的其他形式．此时的解析解为Ｕ（ｘ）＝ＮｓｅｃｈＮｘ）β＝

／传播常数也与解析解的传播常数相一致．Ｎ２２，图２显示孤子解与解析解相吻合，

　４·６·广东第二师范学院学报第３１卷

）和（）归一化解析解和数值解对照图，）为Ｎ＝１．图２　（其中（ａｃａ５，

（）为Ｎ＝２．（）和（）横向上各点的传播常数ｂ（ｃ５；ｂｄｘ）

综上所述，我们应用傅立叶变换得到了一种比较简单的迭代方法，数值模拟结果显示其能快速的收敛，与解析解相一致．对于初学的学生，这是一个非常好的一种迭代方法，易于学生掌握和应用．

参考文献：

［］］．（）：ｕａｄｒａｔｉｃｒｏｂｌｅｍ［１ＴＩＳＳＥＵＲＦ，ＭＥＥＲＢＥＲＧＥＮ　Ｋ．ＴｈｅｅｉｅｎｖａｌｕｅＪＳＩＡＭ　Ｒｅｖｉｅｗ，２００１，４３２　　　　ｑｐｇ

２３５２８６．－

［］ＧＯ］．２ＬＵＢＧ　Ｈ，ＶＡＮＤＥＲＶＯＲＳＴ　Ｈ　Ａ．Ｅｉｅｎｖａｌｕｅｃｏｍｕｔａｔｉｏｎｉｎｔｈｅ２０ｔｈｃｅｎｔｕｒＪＪｏｕｒｎａｌｏｆ　　　　　　　　　ｇｐｙ［

，ａｎｄＡｌｉｅｄＭａｔｈｅｍａｔｉｃｓ２０００，１２３：３５６５．Ｃｏｍｕｔａｔｉｏｎａｌ　　　－ｐｐｐ

［］ＫＡＴ［］．Ｏ３ＺＯ，ＬＡＨＩＮＩＹ，ＳＩＬＢＥＲＢＥＲＧＹ．ＭｕｌｔｉｌｅｂｒｅａｋｕｏｆｈｉｈｏｒｄｅｒｓａｔｉａｌｓｏｌｉｔｏｎｓＪｔｉｃｓ　　　　　－　　ｐｐｇｐｐ　

，Ｌｅｔｔｅｒ２００８，３３：２８３０２８３２．－

［］４ＣＯＲＣＯＲＡＮＢ，ＭＯＮＡＴＣ，ＰＵＤＯＤ，ｅｔａｌ．Ｎｏｎｌｉｎｅａｒｌｏｓｓｄｎａｍｉｃｓｉｎａｓｉｌｉｃｏｎｓｌｏｗｌｉｈｔｈｏｔｏｎｉｃ　　　　　　　　　　－　ｙｇｐ

［］．，ｗａｖｅｕｉｄｅＪＯｔｉｃｓＬｅｔｔｅｒ２０１０，３５：１０７３１０７５．ｃｒｓｔａｌ　　－ｇｐｙ

［］５ＸＵＺＹ，ＫＡＲＴＡＳＨＯＶ　Ｖ，ＴＯＲＮＥＲＬ．Ｕｅｒｔｈｒｅｓｈｏｌｄｆｏｒｓｔａｂｉｌｉｔｏｆｍｕｌｔｉｏｌｅｏｄｅｓｏｌｉｔｏｎｓｉｎ　　　　　　　－ｍ　　ｐｐｙｐ　

［］．，ｎｏｎｌｉｎｅａｒｍｅｄｉａＪＯｔｉｃｓＬｅｔｔｅｒ２００５，３０：３１７１３１７３．ｎｏｎｌｏｃａｌ　　　－ｐ

ＴｈｅＭｅｔｈｏｄｏｆＳｏｌｖｉｎｏｎｌｉｎｅａｒＰａｒｔｉａｌＤｉｆｆｅｒｅｎｔｉａｌ　　　－　　ｇＮ

ＥｕａｔｉｏｎｂｌｉｎｏｕｒｉｅｒＴｒａｎｓｆｏｒｍ　　ｑｙＡｐｐｙｇＦ

１２ＬＩＨｕａＳＨＩＺｈｉｅｉａｎ　－ｇ　－ｗ　ｇ，

（，，，１．ＤｅａｒｔｍｅｎｔｏｆＰｈｓｉｃｓＧｕａｎｄｏｎＵｎｉｖｅｒｓｉｔｏｆＥｄｕｃａｔｉｏｎＧｕａｎｚｈｏｕ　　　ｐｙｇｇｙｇ　　

，；，Ｇｕａｎｄｏｎ５１０３０３，Ｐ．Ｒ．Ｃｈｉｎａ２．ＳｃｈｏｏｌｏｆＩｎｆｏｒｍａｔｉｏｎＥｎｉｎｅｅｒｉｎＧｕａｎｄｏｎ　　　ｇｇｇｇｇｇ

，，，）ＵｎｉｖｅｒｓｉｔｏｆＴｅｃｈｎｏｌｏＧｕａｎｚｈｏｕＧｕａｎｄｏｎ５１０００６，Ｐ．Ｒ．Ｃｈｉｎａ　ｇｙｇｇｇｙ　

ｔｉｏｆｍａｔｈｅ …… 此处隐藏：1299字，全部文档内容请下载后查看。喜欢就下载吧 ……

求解非线性偏微分方程的傅里叶变换方法.doc 将本文的Word文档下载到电脑，方便复制、编辑、收藏和打印

下载这篇word文档

本文链接：https://www.jiaowen.net/wenku/106412.html（转载请注明文章来源）

上一篇：3.1图形的平移新北师大版八年级下册
下一篇：2018-2019-2018幼儿园园长助理工作计划-范文模板 (1页)