资源描述:
《非线性动力学-5.ppt》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、非线性动力学姚宝恒上海交通大学船舶海洋与建工学院BeyondPerturbationIntroductiontoHomotopyAnalysisMethodOutlineConceptofHomotopyinTopologyBasicideasofHomotopyAnalysismethodExamplesApplicationsofthetheoryinsolvingnonlinearequationsConclusionsReferences“摄动方法”的本质:应用方程中的小(大)物理参数,将一个非线性问题转化为无穷多个线性子问题。优点:物理意义明确;
2、简单、易懂;缺点:(1)依赖小参数,当所研究问题不含小参数时使得摄动展开法面临困难(2)摄动展开解只在参数比较小的情况下能够给出较好的近似,随着“小参数”的增大,近似解精度下降,以致失效。(3)无法确保解的收敛怎样的近似解析方法才是最理想的?不依赖小参数确保解的收敛性,适用于强非线性问题拓扑学中的几个基本概念拓扑和拓扑空间如果对一个非空集合X给予适当的结构,使之能引入微积分中的极限和连续的概念,这样的结构就称为拓扑。具有拓扑结构的空间称为拓扑空间。引入拓扑结构的方法有多种,如邻域系、开集系、闭集系、闭包系、内部系等不同方法。同伦的基本概念两个拓扑空间如果可
3、以通过一系列连续的形变从一个变到另一个,那么就称这两个拓扑空间同伦。同伦的定义设X和Y都是拓扑空间,f和g是X到Y的连续映射,即f:X→Y,g:X→Y,如果存在连续映射H:X×I→Y(这里I=[0,1]),使得对任何x∈X,满足:则称f和g是同伦的,称H是由f到g的一个同伦或伦移,即H(x,0)=f(x),H(x,1)=g(x),同伦是关于映射的等价关系f(x)=H(x,0)g(x)=H(x,1)H(x,q)示意图二、“同伦分析方法”简述拓扑理论传统的同伦概念:其中,q为嵌入变量.易知,q=0时,H(x;0)=f(x);q=1时,H(x;1)=g(x).因
4、此,当嵌入变量q从0增加到1时,函数H(x,q)从f(x)连续变化到g(x).这样,H(x,t)建立起从f(x)到和g(x)之间的联系.在拓扑(topology)〕理论中,这种连续的变化称为同伦(homotopy),表示为Liao提出“广义同伦”之概念:BasicideasofHAME1.非线性代数方程f(x)=0.(构造同伦)设为已知的初始猜测解,嵌入变量为一未知的嵌入变量的函数,我们构造如下的一个单参数的非线性代数方程:(1)当时,上述方程为线性方程即当时,方程(1)变为则,就是原非线性方程f(x)=0的解.因此,当嵌入变量从0变化到1时,从初始猜测解
5、变化到非线性代数方程解,因此方程(1)构造了一个的同伦.设存在无穷阶导数根据Taylor定理,有则如何求?(2)将(1)式对p求一阶导数(3)令得则将(3)式对p再求一次导数(4)(5)令得(6)类似地,可以求得k阶变形导数,则一阶近似公式为(时为牛顿迭代公式)E2.非线性微分方程whereisanonlinearoperator,denotesindependentvariable,isanunknownfunction,respectively.(1)Constructzero-orderdeformationequationWhere∈[0,1]is
6、theembeddingparameter,isanonzeroauxiliaryparameter,isanauxiliaryfunction,isanauxiliarylinearoperator,isaninitialguessof,isaunknownfunction,respectively.(7)Obviously,whenp=0andp=1,itholdsThusasincreasesfrom0to1,thesolutionvariesfromtheinitialguesstothesolution.ExpandinginTaylorseri
7、eswithrespectto,onehaswhere(8)Iftheauxiliarylinearoperator,theinitialguess,theauxiliaryparameter,andtheauxiliaryfunctionaresoproperlychosen,theseries(8)convergesat,onehaswhichmustbeoneofsolutionsoforiginalnonlinearequation.Asand,Eq(7)becomes(9)whichisusedmostlyinthehomotopyanalysi
8、smethod.(2)Constructmth-orderdefo