同伦分析方法(火腿),得到象征性的近似解的非线性耦合方程参数的导数。这些非线性耦合方程与参数微分包含许多重要的数学物理方程和反应扩散方程。通过选择不同的参数值一般正式的数值解,因此,快速收敛级数解。验证该方法的效率和精度,使用两个著名的例子:耦合的汉堡和mKdV方程。结果表明,同伦摄动方法同伦分析方法的一个特例。
分数微分方程获得了重要性和受欢迎程度在过去三十年左右,主要是由于其演示应用程序在许多看似不同的科学和工程领域。例如,地震可以被建模的非线性振荡部分衍生品和流体动力交通模型和部分金融衍生品可以消除缺陷引起的连续交通流的假设。与分数阶微分方程最近被证明是有价值的工具,许多物理现象的建模(<一个href="#B5">1一个>,<一个href="#B21">2一个>]。这是由于这一事实的现实建模一个物理现象并不仅仅取决于即时时间,还在以前的历史时间也可以成功地通过使用分数微积分。大多数非线性分式方程没有确切的分析解决方案,所以必须使用近似和数值技术。Adomain分解方法(<一个href="#B18">3一个>),同伦摄动法(<一个href="#B19">4一个>),变分迭代法(<一个href="#B20">5一个>),和其他方法已经被用来提供分析近似线性和非线性问题。然而,相应的结果的收敛区域相当小。1992年,辽采用同伦在拓扑的基本思想提出一个通用分析方法对于非线性问题,即同伦分析方法(<一个href="#B12">6一个>- - - - - -<一个href="#B16">10一个>]。该方法已成功应用于解决科学和工程的许多类型的非线性问题,如非牛顿流体的粘性流动<一个href="#B10">11一个>],KdV-type方程[<一个href="#B1">12一个>],财政问题[<一个href="#B25">13一个>),部分洛伦兹系统<一个href="#B2">14一个>(],延迟微分方程<一个href="#B3">15一个>]。火腿含有一定的辅助参数<年代vg height="10.95" id="M1" style="vertical-align:-0.1254pt;width:9.1374998px;" version="1.1" viewbox="0 0 9.1374998 10.95" width="9.1374998" xmlns="http://www.w3.org/2000/svg">
ℎ
这为我们提供了一个简单的方法来调整和控制级数解的收敛速度和收敛区域。
火腿提供一定的优势超过常规的数值方法。使用离散化数值方法产生舍入错误造成损失的准确性,需要大型计算机内存和时间。这种计算方法收益率分析解决方案和标准数值方法具有一定的优势。火腿的方法更好,因为它不涉及离散化变量,因此是免费的四舍五入错误和不需要大型计算机内存和时间。
在本文中,我们扩展火腿讨论的应用程序显式nonlinear-coupled方程的数值解的类型与参数导数形式:<年代pan class="equation" id="EEq1">
在哪里<年代vg height="14.3625" id="M3" style="vertical-align:-3.2316pt;width:15.2125px;" version="1.1" viewbox="0 0 15.2125 14.3625" width="15.2125" xmlns="http://www.w3.org/2000/svg">
和<年代vg height="14.3625" id="M4" style="vertical-align:-3.2316pt;width:17.450001px;" version="1.1" viewbox="0 0 17.450001 14.3625" width="17.450001" xmlns="http://www.w3.org/2000/svg">
(<年代vg height="12.8875" id="M5" style="vertical-align:-1.76814pt;width:46.512501px;" version="1.1" viewbox="0 0 46.512501 12.8875" width="46.512501" xmlns="http://www.w3.org/2000/svg">
=
1
,
2
)的线性和非线性函数<年代vg height="7.1624999" id="M6" style="vertical-align:-0.11285pt;width:7.5374999px;" version="1.1" viewbox="0 0 7.5374999 7.1624999" width="7.5374999" xmlns="http://www.w3.org/2000/svg">
和<年代vg height="7.4250002" id="M7" style="vertical-align:-0.11285pt;width:8.0375004px;" version="1.1" viewbox="0 0 8.0375004 7.4250002" width="8.0375004" xmlns="http://www.w3.org/2000/svg">
分别<年代vg height="7.1750002" id="M8" style="vertical-align:-0.1254pt;width:8.9375px;" version="1.1" viewbox="0 0 8.9375 7.1750002" width="8.9375" xmlns="http://www.w3.org/2000/svg">
和<年代vg height="13.425" id="M9" style="vertical-align:-2.29482pt;width:8.8500004px;" version="1.1" viewbox="0 0 8.8500004 13.425" width="8.8500004" xmlns="http://www.w3.org/2000/svg">
的参数描述的顺序的导数。时可以获得不同的非线性耦合系统的参数之一<年代vg height="7.1750002" id="M10" style="vertical-align:-0.1254pt;width:8.9375px;" version="1.1" viewbox="0 0 8.9375 7.1750002" width="8.9375" xmlns="http://www.w3.org/2000/svg">
或<年代vg height="13.425" id="M11" style="vertical-align:-2.29482pt;width:8.8500004px;" version="1.1" viewbox="0 0 8.8500004 13.425" width="8.8500004" xmlns="http://www.w3.org/2000/svg">
各不相同。的研究(<一个href="#EEq1">1.1一个>)是非常必要的和重要的因为当<年代vg height="7.1750002" id="M12" style="vertical-align:-0.1254pt;width:8.9375px;" version="1.1" viewbox="0 0 8.9375 7.1750002" width="8.9375" xmlns="http://www.w3.org/2000/svg">
和<年代vg height="13.425" id="M13" style="vertical-align:-2.29482pt;width:8.8500004px;" version="1.1" viewbox="0 0 8.8500004 13.425" width="8.8500004" xmlns="http://www.w3.org/2000/svg">
是整数,它包含了许多重要的数学物理方程。
本文组织如下。符号和基本定义给出了部分<一个href="#sec2">2一个>。节<一个href="#sec3">3一个>同伦分析方法。节<一个href="#sec4">4一个>著名火腿申请两个耦合的例子:汉堡和mKdV方程。提出了讨论和结论部分<一个href="#sec5">5一个>。
描述的基本想法火腿,我们考虑的算子形式(<一个href="#EEq1">1.1一个>):<年代pan class="equation" id="EEq6">
在哪里<年代vg height="10.325" id="M43" style="vertical-align:-0.0pt;width:14.8375px;" version="1.1" viewbox="0 0 14.8375 10.325" width="14.8375" xmlns="http://www.w3.org/2000/svg">
非线性算子,<年代vg height="15.4875" id="M44" style="vertical-align:-3.4383pt;width:19.5px;" version="1.1" viewbox="0 0 19.5 15.4875" width="19.5" xmlns="http://www.w3.org/2000/svg">
和<年代vg height="19.862499" id="M45" style="vertical-align:-3.85211pt;width:19.4375px;" version="1.1" viewbox="0 0 19.4375 19.862499" width="19.4375" xmlns="http://www.w3.org/2000/svg">
站的分数阶导数和被定义为(<一个href="#EEq3">2.3一个>),<年代vg height="9.125" id="M46" style="vertical-align:-0.11285pt;width:5.0124998px;" version="1.1" viewbox="0 0 5.0124998 9.125" width="5.0124998" xmlns="http://www.w3.org/2000/svg">
表示一个独立的运营商<年代vg height="13.45" id="M47" style="vertical-align:-2.21957pt;width:37.9375px;" version="1.1" viewbox="0 0 37.9375 13.45" width="37.9375" xmlns="http://www.w3.org/2000/svg">
(
,
)
,<年代vg height="13.45" id="M48" style="vertical-align:-2.21957pt;width:38.4375px;" version="1.1" viewbox="0 0 38.4375 13.45" width="38.4375" xmlns="http://www.w3.org/2000/svg">
(
,
)
是未知函数。
通过概括传统的同伦方法,廖(<一个href="#B12">6一个>构造所谓的零变形方程:<年代pan class="equation" id="EEq7">
在哪里<年代vg height="13.55" id="M51" style="vertical-align:-2.29482pt;width:59.224998px;" version="1.1" viewbox="0 0 59.224998 13.55" width="59.224998" xmlns="http://www.w3.org/2000/svg">
∈
(
0
,
1
]
是嵌入参数,<年代vg height="13.2375" id="M52" style="vertical-align:-1.95624pt;width:32.875px;" version="1.1" viewbox="0 0 32.875 13.2375" width="32.875" xmlns="http://www.w3.org/2000/svg">
ℎ
≠
0
是一个非零辅助参数,<年代vg height="13.45" id="M53" style="vertical-align:-2.21957pt;width:52.799999px;" version="1.1" viewbox="0 0 52.799999 13.45" width="52.799999" xmlns="http://www.w3.org/2000/svg">
(
)
≠
0
是一个辅助函数,<年代vg height="10.325" id="M54" style="vertical-align:-0.0pt;width:11.1875px;" version="1.1" viewbox="0 0 11.1875 10.325" width="11.1875" xmlns="http://www.w3.org/2000/svg">
是一个辅助线性算子,<年代vg height="14.75" id="M55" style="vertical-align:-3.25793pt;width:44.037498px;" version="1.1" viewbox="0 0 44.037498 14.75" width="44.037498" xmlns="http://www.w3.org/2000/svg">
0
(
,
)
,<年代vg height="14.75" id="M56" style="vertical-align:-3.25793pt;width:44.537498px;" version="1.1" viewbox="0 0 44.537498 14.75" width="44.537498" xmlns="http://www.w3.org/2000/svg">
0
(
,
)
是最初的猜测<年代vg height="13.45" id="M57" style="vertical-align:-2.21957pt;width:37.9375px;" version="1.1" viewbox="0 0 37.9375 13.45" width="37.9375" xmlns="http://www.w3.org/2000/svg">
(
,
)
,<年代vg height="13.45" id="M58" style="vertical-align:-2.21957pt;width:38.4375px;" version="1.1" viewbox="0 0 38.4375 13.45" width="38.4375" xmlns="http://www.w3.org/2000/svg">
(
,
)
和<年代vg height="14.6" id="M59" style="vertical-align:-3.13504pt;width:60.375px;" version="1.1" viewbox="0 0 60.375 14.6" width="60.375" xmlns="http://www.w3.org/2000/svg">
1
(
,
,
)
,<年代vg height="14.6" id="M60" style="vertical-align:-3.13504pt;width:60.375px;" version="1.1" viewbox="0 0 60.375 14.6" width="60.375" xmlns="http://www.w3.org/2000/svg">
2
(
,
,
)
分别是两个未知函数。是很重要的,一个伟大的自由选择辅助在火腿。显然,当<年代vg height="13.55" id="M61" style="vertical-align:-2.29482pt;width:34.487499px;" version="1.1" viewbox="0 0 34.487499 13.55" width="34.487499" xmlns="http://www.w3.org/2000/svg">
=
0
和<年代vg height="13.55" id="M62" style="vertical-align:-2.29482pt;width:34.487499px;" version="1.1" viewbox="0 0 34.487499 13.55" width="34.487499" xmlns="http://www.w3.org/2000/svg">
=
1
,以下是适用的:<年代pan class="equation" id="eq2">
分别。因此,正如<年代vg height="9.875" id="M64" style="vertical-align:-2.29482pt;width:7.9124999px;" version="1.1" viewbox="0 0 7.9124999 9.875" width="7.9124999" xmlns="http://www.w3.org/2000/svg">
增加从0到1,解决方案<年代vg height="14.6" id="M65" style="vertical-align:-3.13504pt;width:60.375px;" version="1.1" viewbox="0 0 60.375 14.6" width="60.375" xmlns="http://www.w3.org/2000/svg">
1
(
,
,
)
,<年代vg height="14.6" id="M66" style="vertical-align:-3.13504pt;width:60.375px;" version="1.1" viewbox="0 0 60.375 14.6" width="60.375" xmlns="http://www.w3.org/2000/svg">
2
(
,
,
)
不同的初始猜测<年代vg height="14.75" id="M67" style="vertical-align:-3.25793pt;width:44.037498px;" version="1.1" viewbox="0 0 44.037498 14.75" width="44.037498" xmlns="http://www.w3.org/2000/svg">
0
(
,
)
,<年代vg height="14.75" id="M68" style="vertical-align:-3.25793pt;width:44.537498px;" version="1.1" viewbox="0 0 44.537498 14.75" width="44.537498" xmlns="http://www.w3.org/2000/svg">
0
(
,
)
解决方案<年代vg height="13.45" id="M69" style="vertical-align:-2.21957pt;width:37.9375px;" version="1.1" viewbox="0 0 37.9375 13.45" width="37.9375" xmlns="http://www.w3.org/2000/svg">
(
,
)
,<年代vg height="13.45" id="M70" style="vertical-align:-2.21957pt;width:38.4375px;" version="1.1" viewbox="0 0 38.4375 13.45" width="38.4375" xmlns="http://www.w3.org/2000/svg">
(
,
)
。扩大<年代vg height="14.6" id="M71" style="vertical-align:-3.13504pt;width:60.375px;" version="1.1" viewbox="0 0 60.375 14.6" width="60.375" xmlns="http://www.w3.org/2000/svg">
1
(
,
,
)
,<年代vg height="14.6" id="M72" style="vertical-align:-3.13504pt;width:60.375px;" version="1.1" viewbox="0 0 60.375 14.6" width="60.375" xmlns="http://www.w3.org/2000/svg">
2
(
,
,
)
在泰勒级数对<年代vg height="9.875" id="M73" style="vertical-align:-2.29482pt;width:7.9124999px;" version="1.1" viewbox="0 0 7.9124999 9.875" width="7.9124999" xmlns="http://www.w3.org/2000/svg">
,我们有<年代pan class="equation" id="EEq9">
在哪里<年代pan class="equation" id="EEq10">
如果辅助线性算子,最初的猜测,辅助参数<年代vg height="10.95" id="M76" style="vertical-align:-0.1254pt;width:9.1374998px;" version="1.1" viewbox="0 0 9.1374998 10.95" width="9.1374998" xmlns="http://www.w3.org/2000/svg">
ℎ
辅助函数是如此的正确选择,该系列(<一个href="#EEq9">3.5一个>)是收敛的<年代vg height="13.55" id="M77" style="vertical-align:-2.29482pt;width:34.487499px;" version="1.1" viewbox="0 0 34.487499 13.55" width="34.487499" xmlns="http://www.w3.org/2000/svg">
=
1
,那么我们就有<年代pan class="equation" id="EEq11">
必须原始非线性方程的解决方案之一,廖证明(<一个href="#B14">8一个>]。作为<年代vg height="10.95" id="M79" style="vertical-align:-0.1254pt;width:47.049999px;" version="1.1" viewbox="0 0 47.049999 10.95" width="47.049999" xmlns="http://www.w3.org/2000/svg">
ℎ
=
−
1
和<年代vg height="13.45" id="M80" style="vertical-align:-2.21957pt;width:55.724998px;" version="1.1" viewbox="0 0 55.724998 13.45" width="55.724998" xmlns="http://www.w3.org/2000/svg">
(
t
)
=
1
,(<一个href="#EEq7">3.2一个>)和(<一个href="#EEq7">3.3一个>)成为<年代pan class="equation" id="EEq12">
主要用于在同伦摄动方法<一个href="#B11">20.一个>),而直接获得的解决方案,不使用泰勒级数。根据定义(<一个href="#EEq10">3.6一个>),零变形的控制方程可以推导出方程(<一个href="#EEq7">3.2一个>)。定义了向量<年代pan class="equation" id="eq3">
微分方程(<一个href="#EEq7">3.2一个>)和(<一个href="#EEq7">3.3一个>)<年代vg height="7.1374998" id="M83" style="vertical-align:-0.10033pt;width:11.225px;" version="1.1" viewbox="0 0 11.225 7.1374998" width="11.225" xmlns="http://www.w3.org/2000/svg">
次对嵌入参数<年代vg height="9.875" id="M84" style="vertical-align:-2.29482pt;width:7.9124999px;" version="1.1" viewbox="0 0 7.9124999 9.875" width="7.9124999" xmlns="http://www.w3.org/2000/svg">
然后设置<年代vg height="13.55" id="M85" style="vertical-align:-2.29482pt;width:34.487499px;" version="1.1" viewbox="0 0 34.487499 13.55" width="34.487499" xmlns="http://www.w3.org/2000/svg">
=
0
最后他们除以<年代vg height="10.825" id="M86" style="vertical-align:-0.11285pt;width:16.424999px;" version="1.1" viewbox="0 0 16.424999 10.825" width="16.424999" xmlns="http://www.w3.org/2000/svg">
!
,我们有所谓的<年代vg height="7.1374998" id="M87" style="vertical-align:-0.10033pt;width:11.225px;" version="1.1" viewbox="0 0 11.225 7.1374998" width="11.225" xmlns="http://www.w3.org/2000/svg">
阶变形方程:<年代pan class="equation" id="EEq13">
在哪里<年代pan class="equation" id="EEq15">
应用Riemann-Liouville积分算子<年代vg height="15.8875" id="M90" style="vertical-align:-1.76814pt;width:40.787498px;" version="1.1" viewbox="0 0 40.787498 15.8875" width="40.787498" xmlns="http://www.w3.org/2000/svg">
,
两侧的建筑(<一个href="#EEq13">3.10一个>),我们有<年代pan class="equation" id="EEq16">
应该强调<年代vg height="14.6875" id="M92" style="vertical-align:-3.20526pt;width:46.337502px;" version="1.1" viewbox="0 0 46.337502 14.6875" width="46.337502" xmlns="http://www.w3.org/2000/svg">
(
,
)
,<年代vg height="14.6875" id="M93" style="vertical-align:-3.20526pt;width:46.837502px;" version="1.1" viewbox="0 0 46.837502 14.6875" width="46.837502" xmlns="http://www.w3.org/2000/svg">
(
,
)
为<年代vg height="12.825" id="M94" style="vertical-align:-1.71797pt;width:38.424999px;" version="1.1" viewbox="0 0 38.424999 12.825" width="38.424999" xmlns="http://www.w3.org/2000/svg">
⩾
1
是由线性方程(<一个href="#EEq13">3.10一个>线性边界条件下),来自最初的问题,这可以很容易地解决符号计算软件MATLAB等。对于上述方法的收敛我们参考读者廖的工作。廖(<一个href="#B13">7一个>]证明了,只要一个级数解的同伦分析方法收敛,它必须精确的解决方案之一。因此,重要的是要确保解决方案系列是收敛的。注意,解决方案系列包含辅助参数<年代vg height="10.95" id="M95" style="vertical-align:-0.1254pt;width:9.1374998px;" version="1.1" viewbox="0 0 9.1374998 10.95" width="9.1374998" xmlns="http://www.w3.org/2000/svg">
ℎ
,我们可以通过绘制所谓的正确选择<年代vg height="10.95" id="M96" style="vertical-align:-0.1254pt;width:9.1374998px;" version="1.1" viewbox="0 0 9.1374998 10.95" width="9.1374998" xmlns="http://www.w3.org/2000/svg">
ℎ
级数收敛曲线,以确保解决方案。
为了说明该方法上面所讨论的,我们考虑参数导数的非线性耦合的汉堡方程操作员形式:<年代pan class="equation" id="EEq17">
在哪里<年代vg height="11.0625" id="M103" style="vertical-align:-0.30096pt;width:32.224998px;" version="1.1" viewbox="0 0 32.224998 11.0625" width="32.224998" xmlns="http://www.w3.org/2000/svg">
>
0
,<年代vg height="14.7" id="M104" style="vertical-align:-3.21404pt;width:66.137497px;" version="1.1" viewbox="0 0 66.137497 14.7" width="66.137497" xmlns="http://www.w3.org/2000/svg">
=
/
和部分运营商<年代vg height="15.4875" id="M105" style="vertical-align:-3.4383pt;width:19.5px;" version="1.1" viewbox="0 0 19.5 15.4875" width="19.5" xmlns="http://www.w3.org/2000/svg">
和<年代vg height="19.862499" id="M106" style="vertical-align:-3.85211pt;width:19.4375px;" version="1.1" viewbox="0 0 19.4375 19.862499" width="19.4375" xmlns="http://www.w3.org/2000/svg">
被定义为(<一个href="#EEq3">2.3一个>)。假设初始值<年代pan class="equation" id="EEq18">
确切的解决方案(<一个href="#EEq17">4.1一个>对特殊情况):<年代vg height="13.55" id="M108" style="vertical-align:-2.29482pt;width:64.25px;" version="1.1" viewbox="0 0 64.25 13.55" width="64.25" xmlns="http://www.w3.org/2000/svg">
=
=
1
是<年代pan class="equation" id="EEq19">
应用同伦分析方法,针对(<一个href="#EEq17">4.1一个>)和给定的初始条件(<一个href="#EEq18">4.2一个>),它是方便的选择<年代pan class="equation" id="EEq20">
作为初始近似(<一个href="#EEq17">4.1一个>)。我们选择线性算子<年代pan class="equation" id="EEq21">
随着房地产<年代vg height="13.45" id="M112" style="vertical-align:-2.21957pt;width:55.287498px;" version="1.1" viewbox="0 0 55.287498 13.45" width="55.287498" xmlns="http://www.w3.org/2000/svg">
(
)
=
0
在哪里<年代vg height="7.1875" id="M113" style="vertical-align:-0.13794pt;width:7.3874998px;" version="1.1" viewbox="0 0 7.3874998 7.1875" width="7.3874998" xmlns="http://www.w3.org/2000/svg">
是积分常数。此外,我们定义了一个系统的非线性算子<年代pan class="equation" id="eq4">
我们构造零级和<年代vg height="7.1374998" id="M115" style="vertical-align:-0.10033pt;width:11.225px;" version="1.1" viewbox="0 0 11.225 7.1374998" width="11.225" xmlns="http://www.w3.org/2000/svg">
阶变形方程,<年代pan class="equation" id="eq5">
我们从一个初始近似<年代vg height="13.5625" id="M117" style="vertical-align:-2.21957pt;width:100.1px;" version="1.1" viewbox="0 0 100.1 13.5625" width="100.1" xmlns="http://www.w3.org/2000/svg">
(
,
0
)
=
年代
我
n
(
)
,<年代vg height="13.5625" id="M118" style="vertical-align:-2.21957pt;width:100.6px;" version="1.1" viewbox="0 0 100.6 13.5625" width="100.6" xmlns="http://www.w3.org/2000/svg">
(
,
0
)
=
年代
我
n
(
)
,因此我们可以获得直接的其他组件<年代pan class="equation" id="eq6">
6阶火腿和绝对误差的精确解<年代vg height="10.95" id="M120" style="vertical-align:-0.1254pt;width:47.049999px;" version="1.1" viewbox="0 0 47.049999 10.95" width="47.049999" xmlns="http://www.w3.org/2000/svg">
ℎ
=
−
1
如图<一个href="//www.newsama.com/journals/ijde/2011/545607/fig1/" target="_blank">1一个>。也绝对错误<年代vg height="14.75" id="M121" style="vertical-align:-3.25793pt;width:78.212502px;" version="1.1" viewbox="0 0 78.212502 14.75" width="78.212502" xmlns="http://www.w3.org/2000/svg">
|
(
)
−
6
(
)
|
计算了在桌子上吗<一个href="//www.newsama.com/journals/ijde/2011/545607/tab1/" target="_blank">1一个>。图<一个href="//www.newsama.com/journals/ijde/2011/545607/fig2/" target="_blank">2一个>显示了数值解的导数的非线性耦合的汉堡方程参数<年代vg height="10.95" id="M122" style="vertical-align:-0.1254pt;width:47.049999px;" version="1.1" viewbox="0 0 47.049999 10.95" width="47.049999" xmlns="http://www.w3.org/2000/svg">
ℎ
=
−
1
,<年代vg height="13.55" id="M123" style="vertical-align:-2.29482pt;width:64.25px;" version="1.1" viewbox="0 0 64.25 13.55" width="64.25" xmlns="http://www.w3.org/2000/svg">
=
=
1
。图<一个href="//www.newsama.com/journals/ijde/2011/545607/fig3/" target="_blank">3一个>显示了显式数值解<年代vg height="10.95" id="M124" style="vertical-align:-0.1254pt;width:47.049999px;" version="1.1" viewbox="0 0 47.049999 10.95" width="47.049999" xmlns="http://www.w3.org/2000/svg">
ℎ
=
−
1
,<年代vg height="10.9125" id="M125" style="vertical-align:-0.17555pt;width:48.299999px;" version="1.1" viewbox="0 0 48.299999 10.9125" width="48.299999" xmlns="http://www.w3.org/2000/svg">
=
1
/
4
,<年代vg height="13.55" id="M126" style="vertical-align:-2.29482pt;width:48.212502px;" version="1.1" viewbox="0 0 48.212502 13.55" width="48.212502" xmlns="http://www.w3.org/2000/svg">
=
1
/
3
在<年代vg height="10.9125" id="M127" style="vertical-align:-0.17555pt;width:51.75px;" version="1.1" viewbox="0 0 51.75 10.9125" width="51.75" xmlns="http://www.w3.org/2000/svg">
=
0
。
0
2
。
(一)年代trong>
(b)年代trong>
(一)年代trong>
(b)年代trong>
(一)年代trong>
(b)年代trong>
所建议的廖(<一个href="#B13">7一个>),适当的地区<年代vg height="10.95" id="M148" style="vertical-align:-0.1254pt;width:9.1374998px;" version="1.1" viewbox="0 0 9.1374998 10.95" width="9.1374998" xmlns="http://www.w3.org/2000/svg">
ℎ
是一个水平线段。我们可以调查的影响<年代vg height="10.95" id="M149" style="vertical-align:-0.1254pt;width:9.1374998px;" version="1.1" viewbox="0 0 9.1374998 10.95" width="9.1374998" xmlns="http://www.w3.org/2000/svg">
ℎ
解决方案的收敛性系列gevin火腿,通过绘制曲线对比<我>h我>,如图<一个href="//www.newsama.com/journals/ijde/2011/545607/fig4/" target="_blank">4一个>。
4.1的话。我>年代pan>这个例子使用同伦摄动方法已经得到解决(<一个href="#B6">21一个>]。所绘制的图和表<年代vg height="10.95" id="M151" style="vertical-align:-0.1254pt;width:47.049999px;" version="1.1" viewbox="0 0 47.049999 10.95" width="47.049999" xmlns="http://www.w3.org/2000/svg">
ℎ
=
−
1
是在良好的协议与HPM图绘制。年代pan>
4.2。的非线性耦合mKdV方程参数的导数
为了说明该方法上面所讨论的,我们考虑参数导数的非线性耦合mKdV方程操作员形式:<年代pan class="equation" id="EEq22">
与初始条件,<年代pan class="equation" id="EEq23">
正如我们所知,当<年代vg height="13.55" id="M154" style="vertical-align:-2.29482pt;width:64.25px;" version="1.1" viewbox="0 0 64.25 13.55" width="64.25" xmlns="http://www.w3.org/2000/svg">
=
=
1
(<一个href="#EEq22">4.9一个>)kink-type孤子解<年代pan class="equation" id="EEq24">
由风机(<一个href="#B7">22一个>),<年代vg height="16.7125" id="M156" style="vertical-align:-2.29482pt;width:277.89999px;" version="1.1" viewbox="0 0 277.89999 16.7125" width="277.89999" xmlns="http://www.w3.org/2000/svg">
=
+
(
1
/
4
)
(
−
4
2
−
6
+
6
/
+
3
2
/
2
)
,<年代vg height="13.125" id="M157" style="vertical-align:-1.95624pt;width:32.325001px;" version="1.1" viewbox="0 0 32.325001 13.125" width="32.325001" xmlns="http://www.w3.org/2000/svg">
≠
0
,<年代vg height="13.125" id="M158" style="vertical-align:-1.95624pt;width:31.200001px;" version="1.1" viewbox="0 0 31.200001 13.125" width="31.200001" xmlns="http://www.w3.org/2000/svg">
≠
0
。
应用同伦分析方法,针对(<一个href="#EEq22">4.9一个>)和给定的初始条件(<一个href="#EEq23">4.10一个>),它在方便选择<年代pan class="equation" id="EEq25">
作为初始近似(<一个href="#EEq23">4.10一个>)。我们选择线性算子<年代pan class="equation" id="EEq26">
随着房地产<年代vg height="13.45" id="M161" style="vertical-align:-2.21957pt;width:55.287498px;" version="1.1" viewbox="0 0 55.287498 13.45" width="55.287498" xmlns="http://www.w3.org/2000/svg">
(
)
=
0
在哪里<年代vg height="7.1875" id="M162" style="vertical-align:-0.13794pt;width:7.3874998px;" version="1.1" viewbox="0 0 7.3874998 7.1875" width="7.3874998" xmlns="http://www.w3.org/2000/svg">
是积分常数。此外,我们定义了一个系统的非线性算子<年代pan class="equation" id="eq7">
我们构造零级和<年代vg height="7.1374998" id="M164" style="vertical-align:-0.10033pt;width:11.225px;" version="1.1" viewbox="0 0 11.225 7.1374998" width="11.225" xmlns="http://www.w3.org/2000/svg">
阶变形方程,<年代pan class="equation" id="eq8">
我们从一个初始近似<年代vg height="13.5625" id="M166" style="vertical-align:-2.21957pt;width:181.75px;" version="1.1" viewbox="0 0 181.75 13.5625" width="181.75" xmlns="http://www.w3.org/2000/svg">
(
,
0
)
=
(
/
2
)
+
t
一个
n
h
(
)
,<年代vg height="13.5625" id="M167" style="vertical-align:-2.21957pt;width:242.8125px;" version="1.1" viewbox="0 0 242.8125 13.5625" width="242.8125" xmlns="http://www.w3.org/2000/svg">
(
,
0
)
=
(
/
2
)
(
1
+
(
/
)
)
+
t
一个
n
h
(
)
,<年代vg height="10.9125" id="M168" style="vertical-align:-0.17555pt;width:47.525002px;" version="1.1" viewbox="0 0 47.525002 10.9125" width="47.525002" xmlns="http://www.w3.org/2000/svg">
=
0
。
1
,<年代vg height="10.8625" id="M169" style="vertical-align:-0.13794pt;width:34.674999px;" version="1.1" viewbox="0 0 34.674999 10.8625" width="34.674999" xmlns="http://www.w3.org/2000/svg">
=
1
,<年代vg height="10.9125" id="M170" style="vertical-align:-0.17555pt;width:53.162498px;" version="1.1" viewbox="0 0 53.162498 10.9125" width="53.162498" xmlns="http://www.w3.org/2000/svg">
=
1
/
3
,因此我们可以直接获得其它组件如下:<年代pan class="equation" id="eq9">
6阶火腿和绝对误差的精确解<年代vg height="10.95" id="M172" style="vertical-align:-0.1254pt;width:47.049999px;" version="1.1" viewbox="0 0 47.049999 10.95" width="47.049999" xmlns="http://www.w3.org/2000/svg">
ℎ
=
−
1
如图<一个href="//www.newsama.com/journals/ijde/2011/545607/fig5/" target="_blank">5一个>。也绝对错误<年代vg height="14.75" id="M173" style="vertical-align:-3.25793pt;width:78.212502px;" version="1.1" viewbox="0 0 78.212502 14.75" width="78.212502" xmlns="http://www.w3.org/2000/svg">
|
(
)
−
6
(
)
|
计算了在表吗<一个href="//www.newsama.com/journals/ijde/2011/545607/tab2/" target="_blank">2一个>。图<一个href="//www.newsama.com/journals/ijde/2011/545607/fig6/" target="_blank">6一个>显示了数值解的导数的非线性耦合的汉堡方程参数<年代vg height="10.95" id="M174" style="vertical-align:-0.1254pt;width:47.049999px;" version="1.1" viewbox="0 0 47.049999 10.95" width="47.049999" xmlns="http://www.w3.org/2000/svg">
ℎ
=
−
1
,<年代vg height="13.55" id="M175" style="vertical-align:-2.29482pt;width:64.25px;" version="1.1" viewbox="0 0 64.25 13.55" width="64.25" xmlns="http://www.w3.org/2000/svg">
=
=
1
。图<一个href="//www.newsama.com/journals/ijde/2011/545607/fig7/" target="_blank">7一个>显示了显式数值解<年代vg height="10.95" id="M176" style="vertical-align:-0.1254pt;width:47.049999px;" version="1.1" viewbox="0 0 47.049999 10.95" width="47.049999" xmlns="http://www.w3.org/2000/svg">
ℎ
=
−
1
,<年代vg height="10.9125" id="M177" style="vertical-align:-0.17555pt;width:48.299999px;" version="1.1" viewbox="0 0 48.299999 10.9125" width="48.299999" xmlns="http://www.w3.org/2000/svg">
=
1
/
2
,<年代vg height="13.55" id="M178" style="vertical-align:-2.29482pt;width:48.212502px;" version="1.1" viewbox="0 0 48.212502 13.55" width="48.212502" xmlns="http://www.w3.org/2000/svg">
=
2
/
3
在<年代vg height="10.9125" id="M179" style="vertical-align:-0.17555pt;width:59.575001px;" version="1.1" viewbox="0 0 59.575001 10.9125" width="59.575001" xmlns="http://www.w3.org/2000/svg">
=
0
。
0
0
2
。