当前位置:文档之家› 奇摄动Volterra型积分微分方程非线性边值问题.pdf

奇摄动Volterra型积分微分方程非线性边值问题.pdf

奇摄动Volterra型积分微分方程非线性边值问题.pdf
奇摄动Volterra型积分微分方程非线性边值问题.pdf

数值积分与微分方程

2.3 数值积分 2.3.1 一元函数的数值积分 函数1 quad 、quadl 、quad8 功能 数值定积分,自适应Simpleson 积分法。 格式 q = quad(fun,a,b) %近似地从a 到b 计算函数fun 的数值积分,误差为10-6。 若给fun 输入向量x ,应返回向量y ,即fun 是一单值函数。 q = quad(fun,a,b,tol) %用指定的绝对误差tol 代替缺省误差。tol 越大,函数计 算的次数越少,速度越快,但结果精度变小。 q = quad(fun,a,b,tol,trace,p1,p2,…) %将可选参数p1,p2,…等传递给函数 fun(x,p1,p2,…),再作数值积分。若tol=[]或trace=[],则用缺省值进行计算。 [q,n] = quad(fun,a,b,…) %同时返回函数计算的次数n … = quadl(fun,a,b,…) %用高精度进行计算,效率可能比quad 更好。 … = quad8(fun,a,b,…) %该命令是将废弃的命令,用quadl 代替。 例2-40 >>fun = inline(‘3*x.^2./(x.^3-2*x.^2+3)’); equivalent to: function y=funn(x) y=3*x.^2./(x.^3-2*x.^2+3); >>Q1 = quad(fun,0,2) >>Q2 = quadl(fun,0,2) 计算结果为: Q1 = 3.7224 Q2 = 3.7224 补充:复化simpson 积分法程序 程序名称 Simpson.m 调用格式 I=Simpson('f_name',a,b,n) 程序功能 用复化Simpson 公式求定积分值 输入变量 f_name 为用户自己编写给定函数()y f x 的M 函数而命名的程序文件名 a 为积分下限 b 为积分上限 n 为积分区间[,]a b 划分成小区间的等份数 输出变量 I 为定积分值 程序 function I=simpson(f_name,a,b,n) h=(b-a)/n; x=a+(0:n)*h; f=feval(f_name,x); N=length(f)-1;

微分方程的积分因子求解法

常微分方程的积分因子求解法 内容摘要:本文给出了几类特殊形式的积分因子的求解方法,并推广到较一般的形式。 关键词: 全微分方程,积分因子。 一、 基本知识 定义1.1 对于形如 0),(),(=+dy y x N dx y x M (1.1) 的微分方程,如果方程的左端恰是x ,y 的一个可微函数),(y x U 的全微分,即d ),(y x U = dy y x N dx y x M ),(),(+,则称(1.1)为全微分方程. 易知,上述全微分方程的通解为 ),(y x U =C , (C 为任意常数). 定理1.1 (全微分方程的判别法)设),(y x M ,),(y x N 在x ,y 平面上的单连通区域G 内具有连续的一阶偏导数,则(1.1)是全微分方程的充要条件为 x y x N y y x M ??=??),(),( (1.2) 证明见参考文献[1]. 定义1.2 对于微分方程(1.1),如果存在可微函数),(y x μ,使得方程 ),(y x μ0),(),(),(=+dy y x N y x dx y x M μ (1.3) 是全微分方程,则称),(y x μ为微分方程(1.1)的积分因子. 定理1.2 可微函数),(y x μ为微分方程(1.1)的积分因子的充要条件为 x y x y x N ??),(ln ),(μ-y y x y x M ??),(ln ),(μ=x y x N y y x M ??-??),(),( (1.4) 证明:由定理1.1得,),(y x μ为微分方程(1.1)的积分因子的充要条件为 x y x N y x y y x M y x ??=??)),(),(()),(),((μμ, 展开即得:

全微分方程及积分因子

全微分方程及积分因子

全微分方程及积分因子 内容:凑微分法,全微分方程的判别式,全微分方程的公式解,积分因子的微分方程,只含一个变量的积分因子和其他特殊形式的积分因子。由于有数学分析多元微积分的基础,本节的定理1可以简化处理。对课本中第三块知识即全微分方程的物理背景可以留到后面处理,对第四块知识增解和失解的情况要分散在本章各小节,每次都要重视这个问题。关于初等积分法的局限性可归到学习近似解法时一起讲解。 重点:全微分方程的公式解和积分因子的计算,难点为凑微分法和积分因子的计算。 习题1(1,3,5),2,3 思考题:讨论其他特殊形式的积分因子。 方程:0),(),(=+dy y x N dx y x M 判定:全微分?x N y M ??≡?? 解法:C dy y x N dx y x M y y x x =+??00),(),(0 初值问题0=C 积分因子:x N y M y M x N ??-??=? ???????-??μμμ1

)(x μ: N x N y M dx d ?? -??=μμ1 )(y μ: M x N y M dy d ??- ??-=μμ1 1.解下列方程: 1)0)(222=-+dy y x xydx 解:x N y M ?? ≡??=x 2 ??=-+x y C dy y xydx 002 )0(2既 C y y x =-3/32 2)0)2(=+---dy xe y dx e y y 解:x N y M ??≡??=y e -- ??=-+-y x y C dy y dx e 00)2(既C y xe y =--2 3)0)1(222=---+dy y x dx y x x 解:x N y M ??≡??=y x --221 ??=---+x y C dy y dx y x x 002)1(2 C y y y x x =-+---+23 232322)(32 )(32 )(32 既C y x x =-+23 2 2)(32 4)0)ln (3 =++dy x y dx x y

微分积分公式(全集)

高中大学数学微分与积分公式(全集) (高中大学数学) 一、001 01101lim 0 n n n m m x m a n m b a x a x a n m b x b x b n m --→∞?=??+++? =??? L L (系数不为0的情况) 二、重要公式(1)0sin lim 1x x x →= (2)()1 0lim 1x x x e →+= (3 ))1n a o >= (4 )1n = (5)limarctan 2 x x π →∞ = (6)lim tan 2 x arc x π →-∞ =- (7)limarccot 0x x →∞ = (8)lim arccot x x π→-∞ = (9)lim 0x x e →-∞ = (10)lim x x e →+∞ =∞ (11)0 lim 1x x x + →= 三、下列常用等价无穷小关系(0x →) sin x x : tan x x : arcsin x x : arctan x x : 2 11cos 2 x x -: ()ln 1x x +: 1x e x -: 1ln x a x a -: ()11x x ? +-?: 四、导数的四则运算法则 ()u v u v '''±=± ()uv u v uv '''=+ 2u u v uv v v '''-??= ??? 五、基本导数公式

⑴()0c '= ⑵1x x μμμ-= ⑶()sin cos x x '= ⑷()cos sin x x '=- ⑸()2tan sec x x '= ⑹()2cot csc x x '=- ⑺()sec sec tan x x x '=? ⑻()csc csc cot x x x '=-? ⑼()x x e e '= ⑽()ln x x a a a '= ⑾()1 ln x x '= ⑿()1 log ln x a x a '= ⒀( )arcsin x '= ⒁( )arccos x '= ⒂()21arctan 1x x '= + ⒃() 2 1arccot 1x x '=-+⒄()1x '= ⒅ '= 六、高阶导数的运算法则 (1)()()() () () ()()n n n u x v x u x v x ±=±???? (2)()() ()()n n cu x cu x =???? (3)()()() ()n n n u ax b a u ax b +=+???? (4)()()() () ()()()0 n n n k k k n k u x v x c u x v x -=?=????∑ 七、基本初等函数的n 阶导数公式 (1)() () !n n x n = (2)() () n ax b n ax b e a e ++=? (3)() () ln n x x n a a a = (4)()() sin sin 2n n ax b a ax b n π? ?+=++??? ??? ? ? (5) ()() cos cos 2n n ax b a ax b n π? ?+=++??? ??? ? ?

非线性薛定谔方程数值解的MATLAB仿真

[键入作者姓名] [键入文档标题] ——利用分步快速傅里叶变换对光纤中光信号的传输方程进行数值求解

1、非线性薛定谔方程 非线性薛定谔方程(nonlinear Schrodinger equation ,NLSE)是奥地利物理学家薛定谔于1926 年提出的,应用在量子力学系统中。由于量子力学主要研究粒子的动力学运动状态,所以不能运用牛顿力学公式来表示。通常在量子力学中,研究系统的状态一般通过波函数(x ,t)来表示。而对波函数的研究主要是求解非线性薛定谔方程。本文主要研究光脉冲在光纤中传输状态下的演变。 一般情况下,光脉冲信号在光纤中传输时,同时受到光纤的色散和非线性效应的影响。通过Maxwell 方程,考虑到光纤的色散和非线性效应,可以推导出光信号在光纤中的传输方程,即非线性薛定谔方程。NLSE 是非线性偏微分方程,一般很难直接求出解析解,于是通过数值方法进行求解。具体分为两大类:(1)分布有限差分法(split-step finite differencemethod ,SSFD);(2)分步傅里叶变换法(split-step Fourier transform method ,SSFT)。一般情况,在达到相同精度,由于分步傅里叶变换法采用运算速度快的快速傅里叶变换,所以相比较有限差分法运算速度快一到两个数量级。于是本文介绍分步傅里叶变换法来对光纤中光信号的传输方程,即非线性薛定谔方程进行数值求解。并通过MATLAB 软件对结果数值仿真。 非线性薛定谔方程的基本形式为: 22||t xx iu u u u =+ 其中u 是未知的复值函数. 目前,采用分步傅立叶算法(Split step Fourier Method)求解非线性薛定谔方程的数值解应用比较多。分步傅立叶方法最早是在1937年开始应用的,这种方法己经被证明是相同精度下数值求解非线性薛定愕方程最快的方法,部分原因是它采用了快速傅立叶变换算法(Fast Fourier Transform Algorithm)。基于MATLAB 科学计算软件以及MATLAB 强大的符号计算功能,完全可以实现分步傅立叶数值算法来对脉冲形状和频谱进行仿真。 一般情况下,光脉冲沿光纤传播时受到色散和非线性效应的共同作用,假设当传输距离 很小的时候,两者相互独立作用,那么,根据这种思想可建立如下分步傅立叶数值算法的数 学模型: 把待求解的非线性薛定谔方程写成以下形式: ??()U D N U z ?=+? (I ) (II )

非线性薛定谔方程的孤子解和怪波解

非线性薛定谔方程的孤子解和怪波解 摘要:光纤中光波的传输模型一直是当前研究的热点理论模型之一,从非线性薛定谔方程到金格堡-朗道方程,都试图对其进行更好的阐释,其次对于非线性动力学系统中,非线性薛定谔方程的解有呈现出非常多有趣的特征,对于其中特定解的研究能够让我们了解脉冲演化的本质,所以本文主要从孤子解的传输入手,并且简单介绍了怪波解的解形式。 薛定谔方程又称薛定谔波动方程,是量子力学的一个基本方程,同时又是量子力学的基本假设之一,由奥地利物理学家薛定谔1926年在《量子化就是本征值问题》中提出的,它在量子力学中的地位非常重要,相当于牛顿定律对于经典力学一样。 随着人们对世界的不断探索,非线性现象逐渐走进人们的视野,这种现象一般大都用非线性偏微分方程的数学模型来描述,显然线性方程已经不能满足人们的需求。 1973年,Hasegawa从含有非线性项的色散方程中推导出了非线性薛定谔方程。非线性薛定谔方程(NLS)是普适性很强的一个基本方程,最简单的形式是: 其中为常数。因为这个方程在几乎所有的物理分支及其他科学领域得到了广泛的应用,如超导,光孤子在光纤中传播,光波导,等离子体中的Langnui波等,所以许多学者对此方程的研究投入了很大的热情,至今还在生机勃勃的向前发展着。 1 分步傅里叶法计算演化过程 对于处理非线性性薛定谔方程,常用的数值仿真方式为分步傅里叶方法,为了简单起见,只考虑二阶色散和自相位调制,不考虑高阶色散、自陡以及四波混频等高阶非线性效应。上述方程中做 2 β为二阶色散,γ表示Kerr效应系数,g和α分别代表光纤中的增益和损耗。对上述方程转化到频域,先不考虑增益和损耗。可以得到 2 k k k k k dA i A i a a dz βγ =?+F. 其中2 2 2 k i β β ?=Ω 令() exp k k A B i z β =?可以得到 () 2exp k k k k dB i a a i z dz γβ =-? F 以上方程可以用四阶龙格库塔直接求解,但是速度较慢,所以我们需要做差分处理。 ()() ()()() 2 exp k k k k k B z z B z i a z a z i z z γβ +?- =-? ? F 再利用() exp k k A B i z β =?可以得到 ()()()() ()()() 2 2 exp exp exp k k k k k k k k A z z A i a z a z z i z a z i a z z i z γβ γβ ?? +?=+??? ?? ?? ?? ≈????? ?? F F 然后做傅里叶反变换就可以得到最终的结果 ()()()() 2 1exp exp - k k k k a z z a z i a z z i z γβ ?? +?=????? ?? F F

积分微分方程word版

西南交通大学数值分析题库 用复化梯形公式计算积分 1 ()f x dx ?,要把区间[0,1]一般要等分 41 份才能保 证满足误差小于0.00005的要求(这里(2) () 1f x ∞ ≤) ;如果知道(2) ()0f x >,则 用复化梯形公式计算积分1 ()f x dx ? 此实际值 大 (大,小)。 在以1 0((),())()(),(),()[0,1]g x f x xf x g x dx f x g x C = ∈?为内积的空间C[0,1] 中,与非零常数正交的最高项系数为1的一次多项式是 2 3 x 3. (15分)导出用Euler 法求解 (0)1y y y λ'=??=? 的公式, 并证明它收敛于初值问题的精确 解 解 Euler 公式 1 1,1, ,,k k k x y y h y k n h n λ -----------(5分) 1 011k k k y h y h y λλ ------------------- (10分) () 11(0)n n x n x y h e h n λλλ??=+=+→→ ?? ? 若用复化梯形求积公式计算积分1 x I e dx = ? 区间[0,1]应分 2129 等分,即要 计算个 2130 点的函数值才能使截断误差不超过 71 102 -?;若改用复化Simpson 公式,要达到同样精度区间[0,1]应分12 等分,即要计算个 25 点的函数值 1.用Romberg 法计算积分 2 3 2 x e dx -? 解 []02()()2b a T f a f b -= += 9.219524346410430E-003 10221()222 b a a b T T f -+=+= 5.574989241319070E-003 10 022243 T T S -= = 4.360144206288616E-003 22T = 4.499817148069681E-003 21 122243 T T S -= = 4.141426*********E-003

求解非线性薛定谔方程的一类数值解法

求解非线性薛定谔方程的一类数值解法 张艳敏,刘明鼎 (青岛理工大学琴岛学院,山东青岛266106) 摘要:利用非标准有限差分方法构造了求解非线性薛定谔方程的两个非标准有限差分格式。对于离散后的差分格式,把关于时间和空间的步长函数作为分母逼近导数项。对于非线性项,通过非局部的离散方法计算了这两个非标准有限差分格式的局部截断误差。数值实验结果验证了非标准有限差分格式的有效性。关键词:非线性薛定谔方程;局部截断误差;数值解法中图分类号:O241.82 文献标识码:A 文章编号:2095-7726(2019)03-0008-03 薛定谔方程是物理学中量子力学的一个重要方程,可以用于研究深水波浪理论。柱(球)非线性薛定谔方程常用于描述单色波的一维自调适、光学的自陷现象、固体中的热脉冲传播和等离子体中的Langnui 波[1–5],因此对于此类方程的研究具有非常重要的意义。 薛定谔方程有线性和非线性两种,在本文中,我们研究的是非线性薛定谔方程。非线性薛定谔方程解的解析表达式是很难得到的,因此求解此类方程最常用的就是数值解法。求非线性薛定谔方程数值解的方法主要有差分方法、配置谱方法[6]、有限元方法[7]和平均离散梯度方法[8]等。在本文中,我们利用非标准有限差分方法研究了非线性薛定谔方程的数值解,这种方法已在求解偏微分方程中得到了广泛的应用[9],其优点是对非线性项作非局部离散,对导数项作离散后用步长函数作分母,这样不仅能保持差分方程的数值解与原方程的解析解具有相同的正性,而且能保持较好的数值稳定性。 1非标准有限差分格式的构造 现在我们利用文献[10-12]给出的方法构造非线 性薛定谔方程的两种非标准有限差分格式,要考虑的非线性薛定谔方程为 (1)相应的初边值条件为 其中:为虚数单位;、、和均为连续函数;和均为正数。 为了得到非线性薛定谔方程的差分格式,需要对式 (1)进行离散。首先,需要利用网格对区域进行分割,取空间步长时间步长其 次,在网格点处,定义数值解其中,且下面将分别构造式(1)的两种非标准有限差分格式。 1.1第一种非标准有限差分格式的构造 为了构造式(1)的第一种非标准有限差分格式,我们利用R.E.Mickens 提出的构造非标准有限差分格式的原理[10]和文献[13-14]中提到的方法,并利用给定的记号,对式(1)进行离散。离散后的差分方程为 其中,和为分母函数,且,且分母是通过步长函数逼近得到的。 从式(4)可以看出,和分别取代了和分母函数的选择依据了薛定谔方程解的性质[4]。 记对式(4)进行整理,可得第36卷第3期Vol.36No.3 新乡学院学报 Journal of Xinxiang University 2019年3月Mar.2019 收稿日期:2018-12-21 基金项目:山东省高校科技计划项目(J17KB053);青岛理工大学琴岛学院教育教学研究重点项目(2018003A)作者简介:张艳敏(1981—),女,山东东营人,副教授,硕士,研究方向:偏微分方程数值分析。通信作者:刘明鼎(1982—),男,辽宁大连人,副教授,硕士,研究方向:偏微分方程数值分析。 222 (,)(,) (,)(,)(,),i u x t u x t u x t u x t g x t t x ??=++??(,0)(), u x f x =(2) 01(0,)(), (,)()u t p t u L t p t =ìí =?。 (3)0,0;x H t T <£<£i (,)g x t ()f x 0()p x 1()p x H T [0,][0,]H T ′,h H M =Δt T N =。(,)m n x t (,),n m m n u u x t =(0,1,2,,),m x mh m M ==L Δ(0,1,2,,),n t n t n N ==L ,M N ++??Z Z 。111212 2(),i n n n n n m m m m m n n n m m m u u u u u u u g D D ++---+=++(4) 1D 2D 12exp(Δ)1,D t D =-=24sin ()2 h 1D 2D Δt 2,h 11122 ,,D R D R D ==

非线性薛定谔方程数值解的MATLAB仿真

admin [非线性薛定谔方程数值解的MATLAB仿真]——利用分步快速傅里叶变换对光纤中光信号的传输方程进行数值求解

1、非线性薛定谔方程 非线性薛定谔方程(nonlinear Schrodinger equation ,NLSE)是奥地利物理学家薛定谔于1926 年提出的,应用在量子力学系统中。由于量子力学主要研究粒子的动力学运动状态,所以不能运用牛顿力学公式来表示。通常在量子力学中,研究系统的状态一般通过波函数(x ,t)来表示。而对波函数的研究主要是求解非线性薛定谔方程。本文主要研究光脉冲在光纤中传输状态下的演变。 一般情况下,光脉冲信号在光纤中传输时,同时受到光纤的色散和非线性效应的影响。通过Maxwell 方程,考虑到光纤的色散和非线性效应,可以推导出光信号在光纤中的传输方程,即非线性薛定谔方程。NLSE 是非线性偏微分方程,一般很难直接求出解析解,于是通过数值方法进行求解。具体分为两大类:(1)分布有限差分法(split-step finite differencemethod ,SSFD);(2)分步傅里叶变换法(split-step Fourier transform method ,SSFT)。一般情况,在达到相同精度,由于分步傅里叶变换法采用运算速度快的快速傅里叶变换,所以相比较有限差分法运算速度快一到两个数量级。于是本文介绍分步傅里叶变换法来对光纤中光信号的传输方程,即非线性薛定谔方程进行数值求解。并通过MATLAB 软件对结果数值仿真。 非线性薛定谔方程的基本形式为: 22||t xx iu u u u =+ 其中u 是未知的复值函数. 目前,采用分步傅立叶算法(Split step Fourier Method)求解非线性薛定谔方程的数值解应用比较多。分步傅立叶方法最早是在1937年开始应用的,这种方法己经被证明是相同精度下数值求解非线性薛定愕方程最快的方法,部分原因是它采用了快速傅立叶变换算法(FastFourier Transform Algorithm)。基于MATLAB 科学计算软件以及MATLAB 强大的符号计算功能,完全可以实现分步傅立叶数值算法来对脉冲形状和频谱进行仿真。 一般情况下,光脉冲沿光纤传播时受到色散和非线性效应的共同作用,假设当传输距离 很小的时候,两者相互独立作用,那么,根据这种思想可建立如下分步傅立叶数值算法的数 学模型: 把待求解的非线性薛定谔方程写成以下形式: ??()U D N U z ?=+? (I ) (II )

非线性薛定谔方程求解

CHAPTER IV NUMERICAL SOLUTIONS TO THE NONLINEAR SCHR?DINGER EQUATION 4.1Introduction In general,analytical solutions to the full Maxwell wave equation for a nonlinear optical system do not exist.Even numerical solutions to the wave equation are extremely difficult to implement due to the dimensionality of the problem.The vector form of the wave equation is a four-dimensional(three spatial,one temporal),second-order partial differential equation.Thus,approximations based on propagation conditions and experimental results are needed in order to solve an approximate scalar form of the wave equation,i.e.the nonlinear Schr?dinger equation.However,the approximations listed in the previous chapter do limit the generality and validity of the solutions.For example, the condition extreme nonlinearity,as for the case in supercontinuum generation,is a propagation regime where slowly varying envelope approximation may be violated. The purpose of this chapter is to provide an introduction to a very powerful method in numerically solving the NLSE,known as the split-step Fourier method (SSFM)[15].The chapter will begin with a list pointing the advantages of the SSFM

微分积分公式(全集)

高中大学数学微分与积分公式(全集) (高中大学数学) 一、0 101101 lim 0n n n m m x m a n m b a x a x a n m b x b x b n m --→∞?=??++ +? =?? ? (系数不为 的情况) 二、重要公式( )0sin lim 1x x x →= ( )()1 0lim 1x x x e →+= ( ))1n a o >= ( )1n = ( )limarctan 2 x x π →∞ = ( )lim tan 2 x arc x π →-∞ =- ( )limarccot 0x x →∞ = ( )lim arccot x x π→-∞ = ( )lim 0x x e →-∞ = ( )lim x x e →+∞ =∞ ( )0lim 1x x x + →= 三、下列常用等价无穷小关系(0x →) sin x x tan x x arcsin x x arctan x x 2 11cos 2 x x - () ln 1x x + 1x e x - 1ln x a x a - ()11x x ? +-? 四、导数的四则运算法则 ()u v u v '''±=± ()uv u v uv '''=+ 2u u v uv v v '''-??= ??? 五、基本导数公式 ⑴()0c '= ⑵1x x μμμ-= ⑶()sin cos x x '=

⑷()cos sin x x '=- ⑸()2tan sec x x '= ⑹()2cot csc x x '=- ⑺()sec sec tan x x x '=? ⑻()csc csc cot x x x '=-? ⑼()x x e e '= ⑽()ln x x a a a '= ⑾()1ln x x '= ⑿()1 log ln x a x a '= ⒀( )arcsin x '= ⒁( )arccos x '= ⒂()21arctan 1x x '= + ⒃() 2 1arccot 1x x '=-+⒄()1x '= ⒅ ' = 六、高阶导数的运算法则 ( )()()()()()()()n n n u x v x u x v x ±=±???? ( )()() () ()n n cu x cu x =???? ( )()() () ()n n n u ax b a u ax b +=+???? ( )()()() ()()()() n n n k k k n k u x v x c u x v x -=?=????∑ 七、基本初等函数的 阶导数公式 ( )()()!n n x n = ( )()()n ax b n ax b e a e ++=? ()() ln n x x n a a a = ()()sin sin 2n n ax b a ax b n π??+=++??? ??? ? ? ()()cos cos 2n n ax b a ax b n π??+=++??? ??? ? ? () () () 1 1! 1n n n n a n ax b ax b +???=- ? +?? + ()() () ()() 1 1! ln 1n n n n a n ax b ax b -?-+=-????+

全微分方程及积分因子

全微分方程及积分因子 内容:凑微分法,全微分方程的判别式,全微分方程的公式解,积分因子的微分方程,只含一个变量的积分因子和其他特殊形式的积分因子。由于有数学分析多元微积分的基础,本节的定理1可以简化处理。对课本中第三块知识即全微分方程的物理背景可以留到后面处理,对第四块知识增解和失解的情况要分散在本章各小节,每次都要重视这个问题。关于初等积分法的局限性可归到学习近似解法时一起讲解。 重点:全微分方程的公式解和积分因子的计算,难点为凑微分法和积分因子的计算。 习题1(1,3,5),2,3 思考题:讨论其他特殊形式的积分因子。 方程:0),(),(=+dy y x N dx y x M 判定:全微分? x N y M ??≡?? 解法:C dy y x N dx y x M y y x x =+??00),(),(0 初值问题0=C 积分因子: x N y M y M x N ??-??=????????-??μμμ1 )(x μ: N x N y M dx d ??-??=μμ1 )(y μ: M x N y M dy d ??-??-=μμ1 1.解下列方程: 1)0)(222=-+dy y x xydx 解: x N y M ??≡??=x 2 ??=-+x y C dy y xydx 002)0(2既 C y y x =-3/32 2)0)2(=+---dy xe y dx e y y

解: x N y M ??≡??=y e -- ??=-+-y x y C dy y dx e 00)2(既C y xe y =--2 3)0)1(222=---+dy y x dx y x x 解: x N y M ??≡??=y x --221 ??=---+x y C dy y dx y x x 002)1(2 C y y y x x =-+---+23232322 )(32)(32)(32 既C y x x =-+2322 )(32 4)0)ln (3=++dy x y dx x y 解: x N y M ??≡??=x 1 C dy y dx x y y x =+??030既C y x y =+4/||ln 4 5)05233 3222=+-+dy y y x dx y y x 解: x N y M ??≡??=326--y x ??=-+-x y C dy y dx y y x 00222253 C y x y x =++-/523 6)02cos )2sin 1(2=-+xdy y dx x y 解: x N y M ??≡??=x y 2sin 2 C ydy dx x y x y =-+??002)2sin 1(

积分和简单的微分方程

第三讲 积分和简单的微分方程 1 对于保守力有p dE F dx =- ,势能极值点就是受力平衡点 2 小量展开能将复杂的表达式简化,用多项式逼近任意函数。重要的公式: 当1x <<时 2 (1)(1)1 (2) n n n x nx x -+=++ + 3 常见的求导公式 1[]'n n x nx -=;[]'x x e e =;[sin ]'cos x x =;[cos ]'sin x x -=;1[ln ]'x x = 积分是变量累计的基本方法。掌握积分之后一方面可以用更为简明的办法处理部分竞赛题,另一方面为同学们自学各种高级课程扫平了障碍。 物理方程常常同时包括某个物理量和这个物理量的导数,这样的方程就叫微分方程。掌握微分方程之后,对于许多问题便可以跳出具体的已知量、未知量的限制,从物理本质的角度,讨论问题的可解性,归纳多题一解的方法。 第一部分 单元函数积分 知识点睛 引入:物理公式分类 物理公式分成:状态方程(初中常见,例如牛二,万有引力)和过程方程(例如动能定理,动量定理)。判定以下方程是状态方程还是过程方程:m V ρ=;F ma =;x vt = 看下面两组方程 U I R =;U IR = q I t = ;q It = 前一组是状态的方程。后一组是过程的方程。当电流是常数的时候,两个式子都是对的。然后电流是变化的时候,前一组方程还成立,后一组得到的就不是电流了,而是电流的平均值。如果还要求结果是瞬时的电流,必须把第二组第一个变成求导数,后一个方程就把乘积变成了对瞬时的电流*时间再求和,也就是我们今天要学的积分。 先看两个例子: 上讲回顾 本讲目标 知识模块

相关主题
文本预览
相关文档 最新文档