应甩微分算子理论建立基本解

应甩微分算子理论建立基本解
应甩微分算子理论建立基本解

各种类型的微分方程及其相应解法教程文件

各种类型的微分方程及其相应解法 专业班级:交土01班 姓名:高云 学号:1201110102 微分方程的类型有很多种,解题时先判断微分方程是哪种类型,可以帮助我们更快解题,所以我们有必要归纳整理一下各类型(主要是一阶和二阶)的微分方程及其相应解法。 一、一阶微分方程的解法 1.可分离变量的方程 dx x f dy y g )()(=,或)()(y g x f dx dy = 其特点是可以把变量x 和y 只分别在等式的两边,解法关键是把变量分离后两边积分。 例1.求微分方程ydy dx y xydy dx +=+2的通解. 解 先合并dx 及dy 的各项,得dx y dy x y )1()1(2-=- 设,01,012≠-≠-x y 分离变量得 dx x dy y y 1112-=- 两端积分??-=-dx x dy y y 1112得 ||ln |1|ln |1|ln 2 112C x y +-=- 于是 2212)1(1-±=-x C y 记,21C C ±=则得到题设方程的通解 .)1(122-=-x C y 2.齐次方程 (1))(x y f dx dy = (2) )(c by ax f dx dy ++=(a ,b 均不等于0) 例2求解微分方程.2222xy y dy y xy x dx -=+- 解 原方程变形为=+--=2222y xy x xy y dx dy ,1222?? ? ??+--??? ??x y x y x y x y 令,x y u =则,dx du x u dx dy +=方程化为,1222u u u u dx du x u +--=+ 分离变量得?? ????-+--??? ??--112212121u u u u ,x dx du = 两边积分得 ,ln ln ln 2 1)2ln(23)1ln(C x u u u +=----

随机微分方程

一、一维分岔 考虑一维随机微分方程 ()()()()()()()()() dX = m X dt +X dB t =m X +X X /2dt +X dB t 6.141σσσσ'-???? 生成的连续动态系统 ()()()()()()t t 00t x =x +m s x dx + s x dB s 6.142??σ?-?? () 它是以 x 为初值的(6.1-41)之唯一强解。假定 ()()m 0 = 00 = 0 6.143σ-,() 从而0是?的一个固定点。对此固定点,dB(t)是随机参激。设m(x) 有界,对所有 x 0≠满足椭圆性条件 ()0 6.144x σ≠-() 这保证最多只有一个平稳概率密度。求解与(6.1-41)相应的平稳FPK 方程得平稳概率密度 ()()() () 12 2m u p x C x exp[ ] 6.145u x du σσ-=-? () 于是,上述动态系统有两种可能的平稳状态:不动点(平衡状态)与非平凡平稳运动。前者的不变测度0δ的密度为()x δ,后者的不变测度ν的密度为(6.1-45)。为研究 D-分岔,需计算这两个不变测度的Lyapunov 指数。为此,考虑(6.1-41)的线性化方程

()()()()dV =m X Vdt +X V dB t =[m (X)((X)(X))/2]Vdt VdB t 6.146σσσσ''''''++- () 利用(2.5-6)之解(2.5-11),得(6.1-46)之解 ()()()()()t t V t =V 0exp[(m +/2)X ds +X dB s ] 6.147 σσσ''''-??() 动态系统?关于测度μ的Lyapunov 指数定义为 ()()1 lim ln V t 6.148t t ?λμ→∞=-() (6.1-47)代入(6.1-48),注意()00σ=,得不动点Lyapunov 指数 ()()()()()()()()00 1() lim [ln 000]00 lim 0(6.1-49)?t t t t B t V m ds dB s m m t t ?λδσσ→∞→∞ '''''=++=+=??对以(6.1-45)为密度的不变测度ν,(6.1-47)代入(6.1-48), 假定σ'有界,m /2σσ'''+可积,得Lyapunov 指数 ()0 1 lim (m /2)(X)ds [m (x)(x)(x)/2]p(x)dx 6.150t t R t ?λνσσσσ→∞''''''=+=+-??() 进行分部积分,并利用(6.1-45),最后得 ()2 m(x) -2p(x)dx 0 6.151(x)R ?λνσ?? =<-??? ??() 随机跨临界分岔考虑(6.1-41)的特殊情形 ()()2dX X X dt X dB t 6.152ασ=-+- () 生成的动态系统族α?

微分算子法

高阶常微分方程的微分算子法 摘自《大学数学解题法诠释》 .徐利治,.冯克勤,.方兆本,.徐森林,.1999 高阶方程的求解自然要比一阶方程更为困难,即使是对于线性微分方程。但是有一个例外:常系数线性微分方程。我们可以完整的求出它的通解来,所以常系数线性方程的求解,主要精力是集中在讨论对应的非齐 次方程的特解。本节主要讨论微分算子法。 1.求方程230y y y ''''''--=的通解. 解 记() n n y D y =,将方程写成 32230D y D y Dy --= 或3 2 (23)0D D D y --= 我们熟知,其实首先要解特征方程 32230D D D --= 得0,1,3D =-故知方程有三特解31,,x x e e -,由于此三特解为线性无关,故立得通解 3123x x y C C e C e -=++ 注:本题方程为齐次常系数三阶常微分方程,线性常微分方程的一般形状是 1111()()()()() n n n n n n n d y d y dy L y a x a x dx dx dx a x y f x ---=++++=L 其中系数1(),,()n a x a x L 是某区间(,)a b 上的连续函数,上述方程又可写成 1 1()(()())n n n L y D a x D a x y -≡+++L ()f x = 可以把上面括号整体看作一种运算,常称为线性微分算子。本题中各()i a x 均为实常数,今后也仅对实常系数的情形来进一步发展线性微分算子方法。 2.求解 61160y y y y ''''''-+-= 解 写成 32 (6116)0D D D y -+-= 从特征方程 3 2 06116D D D =-+- (1)(2)(3)D D D =--- 解得 1,2,3D =共三实根,故可立即写成特解 23123x x x y C e C e C e =++ 3.求解 39130y y y y ''''''-++= 解 写成 32 (3913)0D D D y -++= 或 2 (1)(413)0D D D y +-+= 特征方程 2 (1)(413)0D D D +-+=有根 1,23D i =-±,故对应的特解是x e -,2cos3x e x , 2sin 3x e x 从而通解是 22123cos3sin 3x x x y C e C e x C e x -=++ 4.求(4) 45440y y y y y ''''''-+-+=之通解. 解 写成 432 (4544)0D D D D y -+-+= 或 22 (2)(1)0D D y -+= 特征根是2,2,D i =±,对应的特解应是 22,,cos ,sin x x e xe x x ,故写成通解 21234()()cos sin x y x e C C x C x C x =+++ 5.求1 (cos )y y x -''+=的通解 解 本题为非齐次方程,先求出对应的齐次方程 0y y ''+=的通解,写成2 (1)0D y +=,可知特征根为i ±,相应的通解为112cos sin y C x C x =+ 设原方程有特解形为 *12()cos ()sin y C x x C x x =+ 其中12,C C 为待定函数,常数变异告诉我们,应求解下面的方程组 121 12()cos ()sin 0()(cos )()(sin )(cos ) C x x C x x C x x C x x x -?''+=??''''+=?? 或 121 12()cos ()sin 0()sin ()cos (cos ) C x x C x x C x x C x x x -?''+=??''-+=?? (方程组右端为原方程非齐次项1 (cos )x -),解得 1sin ()cos x C x x '=-,2()1C x '= 或 1()ln cos C x x =,2()C x x = 最后得通解为 1*()()()y x y x y x =+

微分算子法典型例题讲解

高阶常微分方程的微分算子法 高阶方程的求解自然要比一阶方程更为困难,即使是对于线性微分方程。但是有一个例外:常系数线性微分方程。我们可以完整的求出它的通解来,所以常系数线性方程的求解,主要精力是集中在讨论对应的非齐 次方程的特解。本节主要讨论微分算子法。 1.求方程230y y y ''''''--=的通解. 解 记()n n y D y =,将方程写成 32230D y D y Dy --= 或32(23)0D D D y --= 我们熟知,其实首先要解特征方程 32230D D D --= 得0,1,3D =-故知方程有三特解31,,x x e e -,由于此三特解为线性无关,故立得通解 3123x x y C C e C e -=++ 注:本题方程为齐次常系数三阶常微分方程,线性常微分方程的一般形状是 1111()()()()() n n n n n n n d y d y dy L y a x a x dx dx dx a x y f x ---=++++= 其中系数1(),,()n a x a x 是某区间(,)a b 上的连续函数,上述方程又可写成 1 1()(()())n n n L y D a x D a x y -≡+++ ()f x = 可以把上面括号整体看作一种运算,常称为线性微分算子。本题中各()i a x 均为实常数,今后也仅对实常系数的情形来进一步发展线性微分算子方法。 2.求解 61160y y y y ''''''-+-= 解 写成 32(6116)0D D D y -+-= 从特征方程 3 2 06116D D D =-+- (1)(2)(3)D D D =--- 解得 1,2,3D =共三实根,故可立即写成特解 23123x x x y C e C e C e =++ 3.求解 39130y y y y ''''''-++= 解 写成 32(3913)0D D D y -++= 或 2(1)(413)0D D D y +-+= 特征方程 2(1)(413)0D D D +-+=有根 1,23D i =-±,故对应的特解是x e -,2cos3x e x , 2sin3x e x 从而通解是 22123cos3sin3x x x y C e C e x C e x -=++ 4.求(4)45440y y y y y ''''''-+-+=之通解. 解 写成 432(4544)0D D D D y -+-+= 或 22(2)(1)0D D y -+= 特征根是2,2,D i =±,对应的特解应是 22,,cos ,sin x x e xe x x ,故写成通解 21234()()cos sin x y x e C C x C x C x =+++ 5.求1(cos )y y x -''+=的通解 解 本题为非齐次方程,先求出对应的齐次方程 0y y ''+=的通解,写成2(1)0D y +=,可知特征根为i ±,相应的通解为112cos sin y C x C x =+ 设原方程有特解形为 *12()cos ()sin y C x x C x x =+ 其中12,C C 为待定函数,常数变异告诉我们,应求解下面的方程组 121 12()cos ()sin 0 ()(cos )()(sin )(cos ) C x x C x x C x x C x x x -?''+=??''''+=?? 或 121 12()cos ()sin 0()sin ()cos (cos ) C x x C x x C x x C x x x -?''+=??''-+=?? (方程组右端为原方程非齐次项1(cos )x -),解得 1s i n ()cos x C x x '=-,2()1C x '= 或 1()ln cos C x x =,2()C x x = 最后得通解为 1*()()()y x y x y x =+ 12cos sin cos ln cos sin C x C x x x x x =+++

偏微分方程理论的归纳与总结

偏微分方程基本理论的归纳与总结 偏微分方程是储存自然信息的载体,自然现象的深层次性质可以通过数学手段从方程中推导出来.最为一种语言,微分方程在表达自然定律方面比文字具有更强的优越性.微分方程是一个庞大的体系,它的基本问题就是解的存在性和唯一性.该学科的主要特征是不存在一种可以统一处理大多数偏微分方程的适定性问题的普适的方法和理论.这是与常微分方程有显著差异的地方.这种特性使得我们将方程分为许多种不同类型,这种分类的依据主要来自数学与自然现象这两个方面.从数学的角度,方程的类型一般总是对应于一些普遍的理论和工具.换句话讲,如果能建立一个普遍性的方法统一处理一大类方程问题,那么这个类型就被划分出来.而从自然现象的角度,我们又可以根据不同的运动类型以及性质将方程进行分类.当然这两种方式常常不能截然区分,通常它们是相互关联的,这就造成方程的概念有许多重叠现象. 根据数学的特征,偏微分方程主要被分为五大类,它们是: (1)线性与拟微分方程,研究这类方程的主要工具是Fourier分析方法; (2)椭圆型方程,它的方法是先验估计+泛函分析手段; (3)抛物型方程,主要是Galerkin方法,算子半群,及正则性估计; (4)双曲型方程,对应于Galerkin方法; (5)一阶偏微分方程,主要工具是数学分析方法. 从自然界的运动类型出发,偏微分方程可分为如下几大类: (1)稳态方程(非时间演化方程); (2)耗散型演化方程,这类方程描述了时间演化过程中伴有能量损耗与补充的自然运动.相变与混沌是它们的主要内容; (3)保守系统,如具有势能的波方程.该系统控制的运动是与外界隔离的,及无能量输入,也无能量损耗.行波现象与周期运动是它们的主要特征; (4)守恒律系统,这类方程是一阶偏微分方程组,它们与保守系统具有类似的性质,可视为物质流的守恒.激波行为是由守恒律系统来控制. 下面具体来介绍三类经典方程: 三类典型方程:椭圆型方程,抛物型方程,双曲型方程,即偏微分方程模型的建立,解问题的解法以及三类典型方程的基本理论. 关于三类典型方程定解问题的解题方法,它们主要是分离变量法、积分变换法、特征线法、球面平均法、降维法和Green 函数方法. 关于三类典型方程的基本理论——极值原理和能量估计,并由此给出了解的唯一性和稳定性的相关结论. 具体来说,关于二阶线性椭圆形方程,我们研究它的古典解和弱解.前者主要介绍了基本解、调和函数的基本性质、Green 函数、极值原理、最大模估计、能量方法和变分原理;而后者的研究则需要知道Sobolev空间的相关知识再加以研究;关于二阶线性抛物型方程,主要研究它的Fourier 变换、特殊的求解方法、基本解、方程式和方程组的最大值原理以及最大模估计、带有非经典边界条件和非局部项的方程式的最大值原理及能量方法;关于二阶线性双曲型方程,主要研究初值问题的求解方法、初值问题的能量不等式与解的适定性、以及混合问题的能量模估计与解的适定性. 椭圆、抛物和双曲这三类线性偏微分方程解的适定性问题,它们分别以拉普拉斯方程、热传导方程和波动方程作为代表.具体地说,对于某些规则的求解区域试图求出满足特定线性偏微分方程和定解条件的具体解,这就决定了存在性问题;再利用方程本身所具有的特殊性质,将证明所求解是唯一的,也就解决了唯一性问题;关于连续依赖性问题,需要在不同函数空

各类微分方程的解法大全

各类微分方程的解法 1.可分离变量的微分方程解法 一般形式:g(y)dy=f(x)dx 直接解得∫g(y)dy=∫f(x)dx 设g(y)及f(x)的原函数依次为G(y)及F(x),则G(y)=F(x)+C为微分方程的隐式通解 2.齐次方程解法 一般形式:dy/dx=φ(y/x) 令u=y/x则y=xu,dy/dx=u+xdu/dx,所以u+xdu/dx=φ(u),即du/[φ(u)-u]=dx/x 两端积分,得∫du/[φ(u)-u]=∫dx/x 最后用y/x代替u,便得所给齐次方程的通解 3.一阶线性微分方程解法 一般形式:dy/dx+P(x)y=Q(x) 先令Q(x)=0则dy/dx+P(x)y=0解得y=Ce- ∫P(x)dx,再令y=u e-∫P(x)dx代入原方程解得u=∫Q(x) e∫P(x)dx dx+C,所以y=e-∫P(x)dx[∫Q(x)e∫P(x)dx dx+C] 即y=Ce-∫P(x)dx +e- ∫P(x)dx∫Q(x)e∫P(x)dx dx为一阶线性微分方程的通解 4.可降阶的高阶微分方程解法 ①y(n)=f(x)型的微分方程 y(n)=f(x) y(n-1)= ∫f(x)dx+C1 y(n-2)= ∫[∫f(x)dx+C1]dx+C2 依次类推,接连积分n次,便得方程y(n)=f(x)的含有n个任意常数的通解②y”=f(x,y’) 型的微分方程 令y’=p则y”=p’,所以p’=f(x,p),再求解得p=φ(x,C1) 即dy/dx=φ(x,C1),所以y=∫φ(x,C1)dx+C2 ③y”=f(y,y’) 型的微分方程

令y’=p则y”=pdp/dy,所以pdp/dy=f(y,p),再求解得p=φ(y,C1) 即dy/dx=φ(y,C1),即dy/φ(y,C1)=dx,所以∫dy/φ(y,C1)=x+C2 5.二阶常系数齐次线性微分方程解法 一般形式:y”+py’+qy=0,特征方程r2+pr+q=0 6.二阶常系数非齐次线性微分方程解法 一般形式: y”+py’+qy=f(x) 先求y”+py’+qy=0的通解y0(x),再求y”+py’+qy=f(x)的一个特解y*(x) 则y(x)=y0(x)+y*(x)即为微分方程y”+py’+qy=f(x)的通解 求y”+py’+qy=f(x)特解的方法: ①f(x)=P m(x)eλx型 令y*=x k Q m(x)eλx[k按λ不是特征方程的根,是特征方程的单根或特征方程的重根依次取0,1或2]再代入原方程,确定Q m(x)的m+1个系数 ②f(x)=eλx[Pl(x)cosωx+P n(x)sinωx]型 令y*=x k eλx[Q m(x)cosωx+R m(x)sinωx][m=max﹛l,n﹜,k按λ+iω不是特征方程的根或是特征方程的单根依次取0或1]再代入原方程,分别确定Q m(x)和R m(x)的m+1个系数

倒向随机微分方程理论

倒向随机微分方程理论的一段往事 (2008-07-18 22:04:36) 转载 分类:数学江湖 标签: 杂谈 转自:https://www.360docs.net/doc/21180051.html,/ 文章是中国金融数学届的狂牛的老头子:彭实戈写的,在这里转给大家欣 赏。按:这个文章回顾了倒向随机微分方程理论产生的一段往事,同样是数学上一个让人愉悦的故事。 当年,我和Pardoux写的关于倒向随机微分方程 简称BSDE理论的那篇文章发表在一个叫《SystemsandControlLetters》的“小杂志”上。那是一个“有心栽花花不开,无意插柳柳成荫”的故事。BSDE的文章发表于1990年,而这项研究的实际完成是在1989年4月。其时我从法国回来,正在复旦大学做博士后 1988年开始。数学系的李训经教授在复旦组织了一个每周一次的控制论讨论班,讨论班的一个重点是随机系统的最优控制问题。当时雍炯敏刚从美国回来,在复旦任副教授,陈叔平在浙大,经常到复旦来参加讨论班。李老师有两个博士生胡瑛和周迅宇 我刚到复旦时,周迅宇还在日本Nisio教授那里,大概属于联合培养,他们都具备了非常好的概率论和随机分析的基础。我说非常好,是相对于我这个刚从法国著名的Pardoux研究团体回来的“洋博士”而言的。当时从国外回来的“洋博士”还不算多,大家都对我们“另眼相待”。回国后看到复旦的这些博士生的基础打得如此之牢固,令我十分佩服。 讨论班的学术气氛很热烈,有两个主攻方向:一是无穷维系统最优控制的最大值原理;一是随机最优控制问题,扩散项含时间的随机控制系统最大值原理是当时大家关心的公开难题之一。那是一个硕果累累的年代,产生了一批令国际同行刮目相看的研究成果,称其为“FudanGroup”。 复旦对于博士后的生活安排得非常周到。我有一个二室一厅的套间,里面是整套全新的家具。胡瑛是这里的常客——几乎每天都来。经常是进门后没说几句话就坐下来,拿出纸和笔来讨论问题,累了就到校园里去散一会儿步,饿了就出去找个饭店或到食堂吃一顿。我们两个合作写了好几篇文章,当时的主攻方向是广义的和无穷维随机系统的最大值原理。李训经和雍炯敏先生也经常来访,我们也经常去李老师家。我们有一些合作的具体题目。休息的时候,也经常谈及几个“大

微分算子法中D的运算

微分算子法中D 的运算 D :微分的意思,如Dx 2=2x , D 3x 2=0 D 1:积分的意思,如D 1x=2x 2 ******************************************************************************* 定理1:)()(F k F e e D kx kx = 注意使用公式时的前后顺序 例: x x x x e e k e e D 22222225)12()1()1(=+=+=+ 推论:) (1)(F 1k F e e D kx kx = (F(k)≠0) 例:x e y y 2=+'' x e y D 22)1(=+ x x x e e e D y 22222*5 1121)1(1=+=+= ****************************************************************************** 定理2:)(sin sin )(F 22a F ax ax D -?= )(cos cos )(F 22a F ax ax D -?= 注意使用公式时的前后顺序 推论:) (1sin sin )(F 122a F ax ax D -?= (F(-a 2) ≠0) 例:x y y 3cos 24=+) ( x y D 3cos 2)1(4 =+ x x x x D x D y 3cos 4113cos 82121)3(13cos 23cos 1)(123cos )1(1222224*=??=+-??=?+?=?+?=遇到sinax,cosax 时,要凑出D 2来。F(D)里有D 2,即可代换为-a 2,代换后继续算F(D)。 ******************************************************************************* 定理3: )()()()(F x v k D F e x v e D kx kx += 注意使用公式时的前后顺序 推论:)() (1)()(F 1x v k D F e x v e D kx kx += 例:x e x y y 22y 44?=+'-''

姚老师最爱的两招:表格法与微分算子法

姚老师最爱的两招:表格法与微分算子法,因为效率高,所以喜欢,仅此而已!录入可是字字辛苦,希望大家珍惜哦! 一、 分部积分的表格法 分部积分主要针对被积函数为两类函数乘积的类型,主要可以归纳为反幂、对幂、幂三、幂指和三指五种,幂可以扩展为多项式函数,三主要指正弦和余弦两类三角函数,基本原则是把其中一类函数拿去凑微分,遵循“反对幂三指”、越往后越先凑微分的原则,前四种称为“终止模式”,最后一种称为“循环模式”。当涉及到幂函数(多项式函数)次数较高时,需多次用到分部积分,计算较繁且易出错,因此介绍一个推广公式: 定理:设(),()u u x v v x ==有1n +阶连续导数,则 (1)()(1)(2)(3)1(1)''''''(1)n n n n n n n uv dx uv u v u v u v u vdx +---++=-+-++-? ?。(此定理及证 明可略,仅告诉大家,我不是瞎编乱造,而是有理论依据的!) 【证:用数学归纳法。 当0n =时,''uv dx uv u vdx =-??。 设1n k =≥时,(1)()(1)(2)(3)1(1)''''''(1)k k k k k k k uv dx uv u v u v u v u vdx +---++=-+-+ +-?? (*) 则当1n k =+时,(2)(1)(1)(1)'k k k k uv dx udv uv u v dx ++++==-???, 将上式的'u (*)式中的u ,则有 (1)()(1)(2)1(2)'''''''(1)k k k k k k u v dx u v u v u v u vdx +--++=-+++-? ?, 从而(2)(1)()(1)(2)2(2)''''''(1)k k k k k k k uv dx uv u v u v u v u vdx ++--++=-+-+ +-??,得证。】 上述式子并不好记,它的一个直观表达就是表格法,如下表。 1))1)2) v v +-- 下面通过例子给予演示: (1)“幂三”型 例1.1 52(325)cos x x x xdx +-+? 解:

分数阶微分方程-课件

分数阶微分方程 第三讲分数阶微分方程基本理论 一、分数阶微分方程的出现背景及研究现状 1、出现背景 分数阶微积分是关于任意阶微分和积分的理论,它与整数阶微积分是统一的,是整数阶微积分的推广。 整数阶微积分作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受,很多问题的数学模型最终都可以归结为整数阶微分方程的定解问题,其无论在理论分析还是数值求解方面都已有较完善的理论。但当人们进入到复杂系统和复杂现象的研究时,经典整数阶微积分方程对这些系统的描述将遇到以下问题: (1)需要构造非线性方程,并引入一些人为的经验参数和与实际不符的假设条件; (2)因材料或外界条件的微小改变就需要构造新的模型; (3)这些非线性模型无论是理论求解还是数值求解都非常繁琐。 基于以上原因,人们迫切期待着有一种可用的数学工具和可依据的基本原理来对这些复杂系统进行建模。分数阶微积分方程非常适合于刻画具有记忆和遗传性质的材料和过程,其对复杂系统的描述具有建模简单、参数物理意义清楚、描述准确等优势,因而成为复杂力学与物理过程数学建模的重要工具之一。 2、研究现状 在近三个世纪里,对分数阶微积分理论的研究主要在数学的纯理论领域里进行,似乎它只对数学家们有用。然而在近几十年来,分数阶微分方程越来越多的被用来描述光学和热学系统、流变学及材料和力学系统、信号处理和系统识别、控制和机器人及其他应用领域中的问题。分数阶微积分理论也受到越来越多的国内外学者的广泛关注,特别是从实际问题抽象出来的分数阶微分方程成为很多数学工作者的研究热点。随着分数阶微分方程在越来越多的科学领域里出现,无论对分数阶微分方程的理论分析还是数值计算的研究都显得尤为迫切。然而由于分数阶微分是拟微分算子,它的保记忆性(非局部性)对现实问题进行了优美刻画的同时,也给我们的分析和计算造成很大困难。 在理论研究方面,几乎所有结果全都假定了满足李氏条件,而且证明方法也和经典微积分方程一样,换句话说,这些工作基本上可以说只是经典微积分方程理论的一个延拓。对分数阶微分方程的定性分析很少有系统性的结果,大多只是给出了一些非常特殊的方程的求解,且常用的求解方法都是具有局限性的。 在数值求解方面,现有分数阶方程数值算法还很不成熟,主要表现为: (1)在数值计算中一些挑战性难题仍未得到彻底解决,如长时间历程的计算和大空间域的计算等; (2)成熟的数值算法比较少,现在研究较多的算法主要集中在有限差分方法与有限单元法; (3)未形成成熟的数值计算软件,严重滞后于应用的需要。

随机微分方程在物理学中的应用

科技大学 本科毕业论文 论文题目:随机微分方程在物理学中的应用院系:物理科学与技术学院 专业:应用物理 姓名:vvv 学号:0700000069 指导教师:xxx

二零一二年三月 摘要 牛顿和莱布尼兹创建了微积分学,为了描述机械动力学、天文学等领域的物理现象,建立了确定性的微分方程。确定性的微分方程在实际问题中有大量的应用。然而在研究实际物理现象的数学模型时,描述一个具体物理现象所用的一组数学方程不会是完全精确的。实际问题中不确定性因素大量存在且往往是问题的关键所在,不可忽视。由于二十世纪中叶大量的含有不确定性的实际问题的出现,以及对模型精确性要求和实际问题复杂性认识的不断提高,不确定性因素越来越多的被考虑到模型的建立中,这就在微分方程的基础上引入了随机因素,促使了随机积分的构建与发展,并在此基础上建立了随机微分方程的相关理论和方法。 随着科技的发展,随机微分方程越来越广泛地应用于模型的建立和分析中。本文针对物理学中存在随机性的特征,提取其中的数学本质,利用数学方法和策略,建立相应的随机微分方程,分析其中数学特征和数学机理,推导相关的公式和性质,通过分析来更好的理解物理学中的随机性问题。 关键词:随机微分方程;布朗运动;matlab模拟;

Abstract. Newton and Leibniz created calculus, in order to describe the mechanical dynamics, astronomy and other fields of physics, the establishment of a deterministic differential equation. Deterministic differential equations large number of practical problems in application. However, the actual physical phenomena in the study mathematical model to describe the physical phenomenon of a specific set of mathematical equations used to not be completely accurate. Practical problems of uncertainties abound and often the crux of the problem can not be ignored. Since the mid-twentieth century, a lot of uncertainty with the actual problems, and the accuracy of the model and actual problems requires understanding the complexity of continuous improvement, more and more uncertainty to the model to be considered in This is the basis of the differential equations introduced random factor

微分算子法

高阶常微分方程的微分算子法 撰写 摘自《大学数学解题法诠释》 .徐利治,.冯克勤,.方兆本,.徐森林,.1999 高阶方程的求解自然要比一阶方程更为困难,即使是对于线性微分方程。但是有一个例外:常系数线性微分方程。我们可以完整的求出它的通解来,所以常系数线性方程的求解,主要精力是集中在讨论对应的非齐 次方程的特解。本节主要讨论微分算子法。 1.求方程230y y y ''''''--=的通解. 解 记() n n y D y =,将方程写成 32230D y D y Dy --= 或32 (23)0D D D y --= 我们熟知,其实首先要解特征方程 32230D D D --= 得0,1,3D =-故知方程有三特解31,,x x e e -,由于此 三特解为线性无关,故立得通解 3123x x y C C e C e -=++ 注:本题方程为齐次常系数三阶常微分方程,线性常微分方程的一般形状是 1111()()()()() n n n n n n n d y d y dy L y a x a x dx dx dx a x y f x ---=++++=L 其中系数1(),,()n a x a x L 是某区间(,)a b 上的连续函数,上述方程又可写成 11()(()())n n n L y D a x D a x y -≡+++L ()f x = 可以把上面括号整体看作一种运算,常称为线性微分算子。本题中各()i a x 均为实常数,今后也仅对实常系数的情形来进一步发展线性微分算子方法。 2.求解 61160y y y y ''''''-+-= 解 写成 32 (6116)0D D D y -+-= 从特征方程 3 2 06116D D D =-+- (1)(2)(3)D D D =--- 解得 1,2,3D =共三实根,故可立即写成特解 23123x x x y C e C e C e =++ 3.求解 39130y y y y ''''''-++= 解 写成 32 (3913)0D D D y -++= 或 2 (1)(413)0D D D y +-+= 特征方程 2 (1)(413)0D D D +-+=有根 1,23D i =-±,故对应的特解是x e -,2cos3x e x , 2sin 3x e x 从而通解是 22123cos3sin 3x x x y C e C e x C e x -=++ 4.求(4) 45440y y y y y ''''''-+-+=之通解. 解 写成 432 (4544)0D D D D y -+-+= 或 22 (2)(1)0D D y -+= 特征根是2,2,D i =±,对应的特解应是 22,,cos ,sin x x e xe x x ,故写成通解 21234()()cos sin x y x e C C x C x C x =+++ 5.求1 (cos )y y x -''+=的通解 解 本题为非齐次方程,先求出对应的齐次方程 0y y ''+=的通解,写成2 (1)0D y +=,可知特征根为i ±,相应的通解为112cos sin y C x C x =+ 设原方程有特解形为 *12()cos ()sin y C x x C x x =+ 其中12,C C 为待定函数,常数变异告诉我们,应求解下面的方程组 121 12()cos ()sin 0()(cos )()(sin )(cos ) C x x C x x C x x C x x x -?''+=??''''+=?? 或 121 12()cos ()sin 0()sin ()cos (cos ) C x x C x x C x x C x x x -?''+=??''-+=?? (方程组右端为原方程非齐次项1 (cos )x -),解得 1sin ()cos x C x x '=-,2()1C x '= 或 1()ln cos C x x =,2()C x x = 最后得通解为 1*()()()y x y x y x =+

各类微分方程的解法大全

创作编号:BG7531400019813488897SX 创作者:别如克* 各类微分方程的解法 1.可分离变量的微分方程解法 一般形式:g(y)dy=f(x)dx 直接解得∫g(y)dy=∫f(x)dx 设g(y)及f(x)的原函数依次为G(y)及F(x),则G(y)=F(x)+C为微分方程的隐 式通解 2.齐次方程解法 一般形式:dy/dx=φ(y/x) 令u=y/x则y=xu,dy/dx=u+xdu/dx,所以u+xdu/dx=φ(u),即du/[φ(u)-u] =dx/x两端积分,得∫du/[φ(u)-u]=∫dx/x 最后用y/x代替u,便得所给齐次方程的通解 3.一阶线性微分方程解法 一般形式:dy/dx+P(x)y=Q(x) 先令Q(x)=0则dy/dx+P(x)y=0解得y=Ce- ∫P(x)dx,再令y=u e-∫P(x)dx代入原方程 解得u=∫Q(x) e∫P(x)dx dx+C,所以y=e-∫P(x)dx[∫Q(x)e∫P(x)dx dx+C] 即y=Ce-∫P(x)dx +e- ∫P(x)dx∫Q(x)e∫P(x)dx dx为一阶线性微分方程的通解 4.可降阶的高阶微分方程解法 ①y(n)=f(x)型的微分方程 y(n)=f(x) y(n-1)= ∫f(x)dx+C1

y(n-2)= ∫[∫f(x)dx+C1]dx+C2 依次类推,接连积分n次,便得方程y(n)=f(x)的含有n个任意常数的通解②y”=f(x,y’) 型的微分方程 令y’=p则y”=p’,所以p’=f(x,p),再求解得p=φ(x,C1) 即dy/dx=φ(x,C1),所以y=∫φ(x,C1)dx+C2 ③y”=f(y,y’) 型的微分方程 令y’=p则y”=pdp/dy,所以pdp/dy=f(y,p),再求解得p=φ(y,C1) 即dy/dx=φ(y,C1),即dy/φ(y,C1)=dx,所以∫dy/φ(y,C1)=x+C2 5.二阶常系数齐次线性微分方程解法 一般形式:y”+py’+qy=0,特征方程r2+pr+q=0 6.二阶常系数非齐次线性微分方程解法 一般形式: y”+py’+qy=f(x) 先求y”+py’+qy=0的通解y0(x),再求y”+py’+qy=f(x)的一个特解y*(x) 则y(x)=y0(x)+y*(x)即为微分方程y”+py’+qy=f(x)的通解 求y”+py’+qy=f(x)特解的方法: ①f(x)=P m(x)eλx型 令y*=x k Q m(x)eλx[k按λ不是特征方程的根,是特征方程的单根或特征方程的重根依次取0,1或2]再代入原方程,确定Q m(x)的m+1个系数 ②f(x)=eλx[Pl(x)cosωx+P n(x)sinωx]型

微分算子法

微分算子法

高阶常微分方程的微分算子法 摘自《大学数学解题法诠释》 .徐利治,.冯克勤,.方兆本,.徐森林,.1999 高阶方程的求解自然要比一阶方程更为困难,即使是对于线性微分方程。但是有一个例外:常系数线性微分方程。我们可以完整的求出它的通解来,所以常系数线性方程的求解,主要精力是集中在讨论对应的非齐 次方程的特解。本节主要讨论微分算子法。 1.求方程230y y y ''''''--=的通解. 解 记() n n y D y =,将方程写成 3 2 230D y D y Dy --= 或3 2(23)0 D D D y --= 我们熟知,其实首先要解特征方程 32230 D D D --= 得0,1,3D =-故知方程有三特解 31,,x x e e -,由于此三特解为线性 无关,故立得通解 31 23x x y C C e C e -=++ 注:本题方程为齐次常系数三阶常微分方程,线性常微分方程的一般形状是 1111()()() ()() n n n n n n n d y d y dy L y a x a x dx dx dx a x y f x ---=++++=L 其中系数1 (),,()n a x a x L 是某区间 (,) a b 上的连续函数,上述方 程又可写成 11()(()())n n n L y D a x D a x y -≡+++L () f x = 可以把上面括号整体看作 一种运算,常称为线性微分 算子。本题中各()i a x 均为实常数,今后也仅对实常系数的情形来进一步发展线性微分算子方法。 2.求解 61160y y y y ''''''-+-= 解 写 成 32(6116)0 D D D y -+-=

各种类型的微分方程及其相应解法

各种类型的微分方程及其相应解法 专业班级:交土01班 姓名:高云 学号:1201110102 微分方程的类型有很多种,解题时先判断微分方程是哪种类型,可以帮助我们更快解题,所以我们有必要归纳整理一下各类型(主要是一阶和二阶)的微分方程及其相应解法。 一、一阶微分方程的解法 1.可分离变量的方程 dx x f dy y g )()(=,或)()(y g x f dx dy = 其特点是可以把变量x 和y 只分别在等式的两边,解法关键是把变量分离后两边积分。 例1.求微分方程ydy dx y xydy dx +=+2的通解. 解 先合并dx 及dy 的各项,得dx y dy x y )1()1(2-=- 设,01,012≠-≠-x y 分离变量得 dx x dy y y 1112-=- 两端积分??-=-dx x dy y y 1112得 ||ln |1|ln |1|ln 2 112C x y +-=- 于是 2212)1(1-±=-x C y 记,21C C ±=则得到题设方程的通解 .)1(122-=-x C y 2.齐次方程 (1))(x y f dx dy = (2) )(c by ax f dx dy ++=(a ,b 均不等于0) 例2求解微分方程.2222xy y dy y xy x dx -=+- 解 原方程变形为=+--=2222y xy x xy y dx dy ,1222?? ? ??+--??? ??x y x y x y x y 令,x y u =则,dx du x u dx dy +=方程化为,1222u u u u dx du x u +--=+ 分离变量得?? ????-+--??? ??--112212121u u u u ,x dx du = 两边积分得 ,ln ln ln 2 1)2ln(23)1ln(C x u u u +=----

微分算子法典型例题讲解

高阶常微分方程的微分算子法 3.求解 39130y y y y ''''''-++= 解 写成 32 (3913)0D D D y -++= 或 2 (1)(413)0D D D y +-+= 特征方程 2 (1)(413)0D D D +-+=有根 1,23D i =-±,故对应的特解是x e -,2cos3x e x , 2sin 3x e x 从而通解是 22123cos3sin 3x x x y C e C e x C e x -=++ 4.求(4) 45440y y y y y ''''''-+-+=之通解. 解 写成 432 (4544)0D D D D y -+-+= 或 22 (2)(1)0D D y -+= 特征根是2,2,D i =±,对应的特解应是 22,,cos ,sin x x e xe x x ,故写成通解 21234()()cos sin x y x e C C x C x C x =+++ 5.求1 (cos )y y x -''+=的通解 解 本题为非齐次方程,先求出对应的齐次方程 0y y ''+=的通解,写成2 (1)0D y +=,可知特征根为i ±,相应的通解为112cos sin y C x C x =+ 设原方程有特解形为 *12()cos ()sin y C x x C x x =+ 其中12,C C 为待定函数,常数变异告诉我们,应求解下面的方程组 121 12()cos ()sin 0()(cos )()(sin )(cos ) C x x C x x C x x C x x x -?''+=??''''+=?? 或 121 12()cos ()sin 0()sin ()cos (cos ) C x x C x x C x x C x x x -?''+=??''-+=?? (方程组右端为原方程非齐次项1 (cos )x -),解得 1sin ()cos x C x x '=-,2()1C x '= 或 1()ln cos C x x =,2()C x x = 最后得通解为 1*()()()y x y x y x =+ 12cos sin cos ln cos sin C x C x x x x x =+++

相关文档
最新文档