线性微分方程通解和特解的直接解法

一类常系数非齐次线性微分方程通解和特解的直接解法

作者：温大伟， 陈莉， 王红芳， 魏瑾， WEN Da-wei， CHEN Li， WANG Hong-fang， WEI Jin 作者单位：兰州城市学院数学学院,甘肃兰州,730070

刊名：

甘肃高师学报

英文刊名：JOURNAL OF GANSU NORMAL COLLEGES

年，卷(期)：2010，15(2)

被引用次数：0次

参考文献(7条)

1.朱灵用升阶法求常系数非齐次线性微分方程的特解[期刊论文]-高等数学研究 2002(2)

2.张菁一类二阶常系数非齐次线性微分方程通解的求解方法[期刊论文]-高等数学研究 2008(3)

3.徐千里一类常系数非齐次线性微分方程的通解 2000(4)

4.余智君二阶常系数非齐次线性微分方程通解的简易求解法 2008(8)

5.王李常系数非齐次线性微分方程的公式解法[期刊论文]-海南大学学报(自然科学版) 1999(3)

6.东北师范大学数学系微分方程教研室常微分方程 1982

7.王高雄.周之铭.朱思铭.王寿松常微分方程 1978

相似文献(10条)

1.期刊论文唐生强.唐清干n阶常系数非齐次线性微分方程的通解-湖南农业大学学报（自然科学版）2004,30(5) 为研究n阶常系数非齐次线性常微分方程解的问题,求证了n阶常系数非齐次线性常微分方程的通解和特解的积分表达式.利用韦达定理和一个变量替换,对n阶常系数非齐次线性微分方程进行降阶,导出该方程的一个用积分表示的通解公式,并根据特征根的不同情形给出了通解的各种形式及相应的通解和特解公式.

2.期刊论文纪华霞.童宏胜一类n阶非齐次线性微分方程的通解-重庆职业技术学院学报2008,17(2)

一般n阶非齐次线性微分方程的解法是比较困难的,通过一阶非齐次线性微分方程的常数变易法推广到n阶非齐次线性微分方程,得出求通解的方法. 3.期刊论文张菁一类二阶常系数非齐次线性微分方程通解的求解方法-高等数学研究2008,11(3)

根据一类二阶常系数非齐次线性微分方程系数的特点,利用降阶法,给出了求其通解的一种简便方法.当方程的系数满足新方法的要求时,非齐次项的选择范围较大,不局限于通常的两类型.

4.期刊论文佘智君.SHE Zhi-jun二阶常系数非齐次线性微分方程通解的简易求解法-重庆工学院学报（自然科学版）2008,22(8)

介绍了求解二阶常系数非齐次线性微分方程的2种简易方法--降阶法和积分法,扩大了可求解二阶常系数非齐次线性微分方程的范围,并举例说明了它们的应用.

5.期刊论文王焕求二阶和三阶常系数非齐次线性微分方程特解的一个公式-高等数学研究2006,9(3)

基于微分算子分裂的思想,受到一阶线性方程求解公式的启发,运用多重积分交换积分顺序的技巧,得到求二阶和三阶常系数非齐次线性微分方程特解的一般性公式.

6.期刊论文郭世贞.张继红.GUO Shi-zhen.ZHANG Ji-hong常系数非齐次线性微分方程特解的一种公式化解法--高等数学中微分方程教学方法的一种新尝试-大学数学2005,21(5)

给出了常系数非齐次线性微分方程特解的一种新的公式化求解方法.它有助于学生全面了解方程的解法,便于记忆和应用,并且扩大了可求解方程的范围.

7.期刊论文邢春峰.袁安锋.王朝旺.XING Chun-feng.YUAN An-feng.WANG Chao-wang求二阶线性常系数非齐次微分方程通解的一种新方法-北京联合大学学报（自然科学版）2009,23(3)

为了更多地得到理论上和应用上占有重要地位的二阶常系数线性非齐次微分方程的通解,这里使用常数变易法,在先求得二阶常系数线性齐次微分方程一个特解的情况下,将二阶常系数线性非齐次微分方程转化为可降阶的微分方程,从而给出了一种运算量较小的二阶常系数线性非齐次微分方程通解的一般公式,并且将通解公式进行了推广,实例证明该方法是可行的.

8.期刊论文文淑慧n阶常系数非齐次线性微分方程通解的一种显式表达式-职大学报2001,""(4)

本文讨论某类特定的n阶常微分方程的求解问题,并给出一种简单可行的解法.

9.期刊论文鞠晶.Ju Jing常数变易法求解一类四阶常系数非齐次线性微分方程-中国科技信息2007,""(7)

利用常数变易法求解具有实特征根的四阶常系数非齐次线性微分方程,在无需求其特解及基本解组的情况下给出其通解公式,井举例验证公式的适用性.

10.期刊论文米翠兰.MI Cui-lan常数变易法求解三阶常系数非齐次线性微分方程-唐山师范学院学报2005,27(2) 用解一阶微分方程的常数变易法求解三阶常系数非齐次线性微分方程y″＇+py″+qy′+sy=f(x),其优点是无需求特解,无须求基本解组,但可求通解,并且给出了一个通用的公式.

本文链接：https://www.wendangku.net/doc/9810494261.html,/Periodical_gsgsxb201002002.aspx

授权使用：中共汕尾市委党校(zgsw)，授权号：f23a3c0c-45d9-41a4-a408-9dcf014c27d6

下载时间：2010年8月11日

二阶变系数线性微分方程的特解

二阶变系数线性微分方程的特解 张金战 ( 陇南师范高等专科学校, 甘肃成县 742500) 摘要: 在已知二阶变系数齐次微分方程的一个非零特解的条件下, 可以得到 该齐次微分方程和与它对应的非齐次微分方程的通解, 本文给出了在二阶变系数齐次微分方程的系数满足一定条件下的特解形式. 关键词: 线性微分方程; 特解; 通解 中图分类号: O 175.1 文献标识码: A 文章编号: 1008- 9020( 2007) 02- 014- 02 1 、引言对于方程( 2) 的特解的确定, 有以下结论: 2二阶变系数线性微分方程是指定理 1 若存在实数 a,使 a+ap(x)+q(x)=0, 则方程( 2) 有特 ax 解 y=e. 1y"+p(x)y'+q(x)y=f(x) ( 1) 2axax2ax 证明 : 设 a+ap(x)+q(x)=0, 将 y=e,y'=ae, y"=ae代入方 111y"+p(x)y'+q(x)y=0 ( 2) 2axaxaxax 2程( 2) 的左端得 : ae+aep (x)+eq (x)=e[a+ap (x)+q (x)]=0, 即其中 p( x) ,q(x),f(x)都是关于 x 的连续函数, 方程( 1) 称为 ax y=e是方程( 2) 的特解. 1二阶变系数非齐次线性微分方程, 方程( 2) 称为方程( 1) 对应 x推论1 若 q(x)+p(x)+1=0,则方程( 2) 有特解 y=e. 1的齐次微分方程. 在已知方程( 2) 的一个非零特解的条件下, - x推论 2 若 q(x)- p(x)+1=0,则方程( 2) 有 特解 y=e. 1文[1]给出了求方程( 2) 的通解的刘维尔公式, 文[2]、文[3]给出 推论 3 若 q(x)=0,则方程( 2) 有特解 y=1. 1了方程( 1) 的一个通解公式.这样将求解方程( 1) 和( 2) 的问题 2 定理 2 若 k?1 且 k(k- 1)+kxp(x)+xq(x)=0,则方程( 2) 有特就转化成了找出方程( 2) 的一个非零特解的问题 , 但求方程 k解 y=x. 1( 2) 的特解没有一般方法, 通常用观察法, 多数情况下难以操 2kk- 1证明 : 设 k (k- 1)+kxp (x)+xq (x)=0, 将 y=x,y'=kx,y"=k

偏微分方程数值解法试题与答案

一．填空（1553=?分） 1．若步长趋于零时，差分方程的截断误差0→lm R ，则差分方程的解lm U 趋近于微分方 程的解lm u . 此结论_______（错或对）； 2．一阶Sobolev 空间{} )(,,),()(21 Ω∈''=ΩL f f f y x f H y x 关于内积=1),( g f _____________________是Hilbert 空间； 3．对非线性（变系数）差分格式，常用 _______系数法讨论差分格式的_______稳定性； 4．写出3 x y =在区间]2,1[上的两个一阶广义导数：_________________________________， ________________________________________； 5．隐式差分格式关于初值是无条件稳定的. 此结论_______（错或对）。 二．（13分）设有椭圆型方程边值问题 用1.0=h 作正方形网格剖分 。 （1）用五点菱形差分格式将微分方程在内点离散化； （2）用截断误差为)(2 h O 的差分法将第三边界条件离散化； （3）整理后的差分方程组为 三．（12）给定初值问题 x u t u ??=?? ， ()10,+=x x u 取时间步长1.0=τ，空间步长2.0=h 。试合理选用一阶偏心差分格式（最简显格式）, 并以此格式求出解函数),(t x u 在2.0,2.0=-=t x 处的近似值。 1．所选用的差分格式是： 2．计算所求近似值： 四．（12分）试讨论差分方程 ()h a h a r u u r u u k l k l k l k l ττ + - = -+=++++11,111 1 逼近微分方程 0=??+??x u a t u 的截断误差阶R 。 思路一：将r 带入到原式，展开后可得格式是在点（l+1/2,k+1/2）展开的。 思路二：差分格式的用到的四个点刚好是矩形区域的四个顶点，可由此构造中心点的差分格 式。

微分方程习题及解答

第十二章 微分方程 §12.1 微分方程基本概念、可分离变量的微分方程、齐次微分方程 一、单项选择题 1. 下列所给方程中，不是微分方程的是( ) ． (A)2xy y '=； (B)222x y C +=； (C)0y y ''+=； (D)(76)d ()d 0x y x x y y -++=． 答(B). 2. 微分方程4(3)520y y xy y '''+-=的阶数是( )． (A)1； (B)2； (C)3； (D)4； 答(C). 3. 下列所给的函数，是微分方程0y y ''+=的通解的是( )． (A)1cos y C x =； (B)2sin y C x =； (C)cos sin y x C x =+； (D)12cos sin y C x C x =+ 答(D). 4. 下列微分方程中，可分离变量的方程是( )． (A)x y y e +'=； (B)xy y x '+=； (C)10y xy '--=； (D)()d ()d 0x y x x y y -++=． 答(A). 5. 下列微分方程中，是齐次方程是微分方程的是( )． (A)x y y e +'=； 2(B)xy y x '+=； (C)0y xy x '--=； (D)()d ()d 0x y x x y y -++=． 答(D). 二、填空题 1．函数25y x =是否是微分方程2xy y '=的解? ． 答：是 .

2．微分方程3d d 0,4x x y y y x =+==的解是 ． 答：2225x y +=． 3．微分方程23550x x y '+-=的通解是． 答：32 52 x x y C =++． 4．微分方程ln 0xy y y '-=的通解是 ． 答： Cx y e =. 5．'的通解是 ． 答：arcsin arcsin y x C =+． 6．微分方程 (ln ln )xy y y y x '-=-的通解是． 答：Cx y e x =． 三、解答题 1．求下列微分方程的通解． (1) 22sec tan d sec tan d 0x y x y x y +=； (2) 2()y xy a y y '''-=+； 解： 解： (3) d 10d x y y x +=； (4) 23d (1)0.d y y x x ++= 解： 解： 2．求下列微分方程满足所给初始条件的特解： (1) 20,0x y x y e y -='==； (2) 2 sin ln ,x y x y y y e π='==； 解： 解： (3) 2d 2d 0,1x x y y x y =+==； (4) d 10d x y y x +=． 解： 解： 3*．设连续函数20()d ln 22x t f x f t ?? =+ ????，求()f x 的非积分表达式． 答：()ln 2x f x e =?．

Maab求解微分方程组及偏微分方程组

第四讲 Matlab 求解微分方程（组） 理论介绍：Matlab 求解微分方程（组）命令 求解实例：Matlab 求解微分方程（组）实例 实际应用问题通过数学建模所归纳得到的方程，绝大多数都是微分方程，真正能得到代数方程的机会很少.另一方面，能够求解的微分方程也是十分有限的，特别是高阶方程和偏微分方程（组）.这就要求我们必须研究微分方程（组）的解法：解析解法和数值解法. 一．相关函数、命令及简介 1.在Matlab 中，用大写字母D 表示导数，Dy 表示y 关于自变量的一阶导数，D2y 表示y 关于自变量的二阶导数，依此类推.函数dsolve 用来解决常微分方程（组）的求解问题，调用格式为： X=dsolve(‘eqn1’,’eqn2’,…) 函数dsolve 用来解符号常微分方程、方程组，如果没有初始条件，则求出通解，如果有初始条件，则求出特解. 注意，系统缺省的自变量为t 2.函数dsolve 求解的是常微分方程的精确解法，也称为常微分方程的符号解.但是，有大量的常微分方程虽然从理论上讲，其解是存在的，但我们却无法求出其解析解，此时，我们需要寻求方程的数值解，在求常微分方程数值解方面，MATLAB 具有丰富的函数，我们将其统称为solver ，其一般格式为： [T,Y]=solver(odefun,tspan,y0) 说明：(1)solver 为命令ode45、ode23、ode113、ode15s 、ode23s 、ode23t 、ode23tb 、ode15i 之一. (2)odefun 是显示微分方程'(,)y f t y =在积分区间tspan 0[,]f t t =上从0t 到f t 用初始条件0y 求解. (3)如果要获得微分方程问题在其他指定时间点012,,,,f t t t t L 上的解，则令tspan 012[,,,]f t t t t =L (要求是单调的). (4)因为没有一种算法可以有效的解决所有的ODE 问题，为此，Matlab 提供了多种求解器solver ，对于不同的ODE 问题，采用不同的solver.

第三章 一阶线性微分方程组 第四讲 常系数线性微分方程组的解法1

第四讲常系数线性微分方程组的解法（4课时） 一、目的与要求: 理解常系数线性微分方程组的特征方程 式, 特征根, 特征向量的概念, 掌握常系数线性微分方程组的基本解组的求法. 二、重点：常系数线性微分方程组的基本解组的求法. 三、难点：常系数线性微分方程组的特征方程式, 特征根, 特征向量的概念. 四、教学方法：讲练结合法、启发式与提问式相结合教学法. 五、教学手段：传统板书与多媒体课件辅助教学相结合. 六、教学过程： 1 新课引入 由定理3.6我们已知道，求线性齐次方程组(3.8)的通解问题，归结到求其基本解组. 但是对于一般的方程组(3.8)，如何求出基本解组，至今尚无一般方法. 然而对于常系数线性齐次方程组 dY AY dx (3.20)

其中A 是n n ?实常数矩阵，借助于线性代数中的约当(Jordan)标准型理论或矩阵指数，可以使这一问题得到彻底解决. 本节将介绍前一种方法，因为它比较直观. 由线性代数知识可知，对于任一n n ?矩阵A ，恒存在非奇异的n n ?矩阵T ，使矩阵1 T AT -成为约当标准型. 为此，对方程组(3.20)引入非奇异线性变换 Y TZ = (3.21) 其中()(,1,2, ,),ij T t i j n == det 0T ≠,将方程组 (3.20)化为 1 dZ T ATZ dx -= (3.22) 我们知道，约当标准型 1 T AT -的形式与矩阵A 的特征方程 11121212221 2 det()0 n n n n nn a a a a a a A E a a a λλλλ ---= =- 的根的情况有关. 上述方程也称为常系数齐次方程组(3.20)的特征方程式.它的根称为矩阵A 的特征根.

微分方程(习题及解答)

第十二章 微分方程 § 微分方程基本概念、可分离变量的微分方程、齐次微分方程 一、单项选择题 1. 下列所给方程中，不是微分方程的是( ) ． (A)2xy y '=； (B)222x y C +=； (C)0y y ''+=； (D)(76)d ()d 0x y x x y y -++=． 答(B). 2. 微分方程4(3)520y y xy y '''+-=的阶数是( )． (A)1； (B)2； (C)3； (D)4； 答(C). 3. 下列所给的函数，是微分方程0y y ''+=的通解的是( )． (A)1cos y C x =； (B)2sin y C x =； (C)cos sin y x C x =+； (D)12cos sin y C x C x =+ 答(D). 4. 下列微分方程中，可分离变量的方程是( )． (A)x y y e +'=； (B)xy y x '+=； (C)10y xy '--=； (D)()d ()d 0x y x x y y -++=． 答(A). 5. 下列微分方程中，是齐次方程是微分方程的是( )． (A)x y y e +'=； 2(B)xy y x '+=； (C)0y xy x '--=； (D)()d ()d 0x y x x y y -++=． 答(D). 二、填空题 1．函数25y x =是否是微分方程2xy y '=的解 ． 答：是 . 2．微分方程 3d d 0,4x x y y y x =+==的解是 ． 答：2225x y +=． 3．微分方程2 3550x x y '+-=的通解是 ． 答：32 52 x x y C =++． 4．微分方程ln 0xy y y '-=的通解是 ． 答： Cx y e =. 5'=的通解是 ． 答：arcsin arcsin y x C =+． 6．微分方程 (ln ln )xy y y y x '-=-的通解是． 答： Cx y e x =． 三、解答题 1．求下列微分方程的通解． (1) 22sec tan d sec tan d 0x y x y x y +=； (2) 2()y xy a y y '''-=+； 解： 解： (3) d 10d x y y x +=； (4) 23d (1)0.d y y x x ++=

高阶线性微分方程常用解法介绍

高阶线性微分方程常用解法简介 关键词：高阶线性微分方程 求解方法 在微分方程的理论中，线性微分方程是非常值得重视的一部分内容，这不仅 因为线性微分方程的一般理论已被研究的十分清楚，而且线性微分方程是研究非线性微分方程的基础，它在物理、力学和工程技术、自然科学中也有着广泛应用。下面对高阶线性微分方程解法做一些简单介绍. 讨论如下n 阶线性微分方程：1111()()()()n n n n n n d x d x dx a t a t a t x f t dt dt dt ---++++= （1），其中()i a t （i=1,2,3,,n ）及f(t)都是区间a t b ≤≤上的连续函数，如果 ()0f t ≡，则方程（1）变为 1111()()()0n n n n n n d x d x dx a t a t a t x dt dt dt ---++++= （2），称为n 阶齐次线性微分方程，而称一般方程（1）为n 阶非齐次线性微分方程，简称非齐次线性微分方程，并且把方程（2）叫做对应于方程（1）的齐次线性微分方程. 1.欧拉待定指数函数法 此方法又叫特征根法，用于求常系数齐次线性微分方程的基本解组。形如 111121[]0,(3),n n n n n n n d x d x dx L x a a a x dt dt dt ---≡++++=其中a a a 为常数，称为n 阶常系数齐次线性微分方程。 111111111111[]()()()n t n t t t t n n n n n n n t t n n n n n n n d e d e de L e a a a e dt dt dt a a a e F e F a a a n λλλλλλλλλλλλλλλλ---------≡++++=++++≡≡++++其中=0(4)是的次多项式. ()F λ为特征方程，它的根为特征根. 1.1特征根是单根的情形 设12,,,n λλλ是特征方程111()0n n n n F a a a λλλλ--≡++++=的n 个彼此不相等的根，则应相应地方程（3）有如下n 个解：12,,,.n t t t e e e λλλ（5）我们指出这n 个解在区间a t b ≤≤上线性无关，从而组成方程的基本解组. 如果(1,2,,)i i n λ=均为实数，则（5）是方程（3）的n 个线性无关的实值 解，而方程（3）的通解可表示为1212,n t t t n x c e c e c e λλλ=+++其中12,,,n c c c 为任意常数. 如果特征方程有复根，则因方程的系数是实常数，复根将称对共轭的出现.设1i λαβ=+是一特征根，则2i λαβ=-也是特征根，因而于这对共轭复根

第三章 一线性微分方程组 第四讲 常系数线性微分方程组的解法(1)

第四讲 常系数线性微分方程组的解法（4课时） 一、目的与要求: 理解常系数线性微分方程组的特征方程式, 特征根, 特征向量的概念, 掌握常系数线性微分方程组的基本解组的求法. 二、重点：常系数线性微分方程组的基本解组的求法. 三、难点：常系数线性微分方程组的特征方程式, 特征根, 特征向量的概念. 四、教学方法：讲练结合法、启发式与提问式相结合教学法. 五、教学手段：传统板书与多媒体课件辅助教学相结合. 六、教学过程： 1 新课引入 由定理3.6我们已知道，求线性齐次方程组(3.8)的通解问题，归结到求其基本解组. 但是对于一般的方程组(3.8)，如何求出基本解组，至今尚无一般方法. 然而对于常系数线性齐次方程组 dY AY dx = (3.20) 其中A 是n n ?实常数矩阵，借助于线性代数中的约当(Jordan)标准型理论或矩阵指数，可以使这一问题得到彻底解决. 本节将介绍前一种方法，因为它比较直观. 由线性代数知识可知，对于任一n n ?矩阵A ，恒存在非奇异的n n ?矩阵T ，使矩阵1T AT -成为约当标准型. 为此，对方程组(3.20)引入非奇异线性变换 Y TZ = (3.21) 其中()(,1,2,,),ij T t i j n ==L det 0T ≠,将方程组(3.20)化为 1dZ T ATZ dx -= (3.22) 我们知道，约当标准型1 T AT -的形式与矩阵A 的特征方程 111212122212det()0n n n n nn a a a a a a A E a a a λλλλ---==-L L M M M L

的根的情况有关. 上述方程也称为常系数齐次方程组(3.20)的特征方程式.它的根称为矩阵A 的特征根. 下面分两种情况讨论. (一) 矩阵A 的特征根均是单根的情形. 设特征根为12,,,,n λλλL 这时 12100 n T AT λλλ-??????=?????? 方程组(3.20)变为 11122200n n n dz dx z dz z dx z dz dx λλλ??????????????????????=???????????????? ?????? M M (3.23) 易见方程组(3.23)有n 个解 1110(),00x Z x e λ????????=????????M 220010(),,()0001n x x n Z x e Z x e λλ????????????????==???????????????? L M M 把这n 个解代回变换(3.21)之中，便得到方程组(3.20)的n 个解 12()i i i i x x i i ni t t Y x e e T t λλ???? ??==?????? M (1,2,,)i n =L

微分方程习题及解答

第十二章 微分方程 § 微分方程基本概念、可分离变量的微分方程、齐次微分方程 一、单项选择题 1. 下列所给方程中，不是微分方程的是( ) ． (A)2xy y '=； (B)222x y C +=； (C)0y y ''+=； (D)(76)d ()d 0x y x x y y -++=． 答(B). 2. 微分方程4(3)520y y xy y '''+-=的阶数是( )． (A)1； (B)2； (C)3； (D)4； 答(C). 3. 下列所给的函数，是微分方程0y y ''+=的通解的是( )． (A)1cos y C x =； (B)2sin y C x =； (C)cos sin y x C x =+； (D)12cos sin y C x C x =+ 答(D). 4. 下列微分方程中，可分离变量的方程是( )． (A)x y y e +'=； (B)xy y x '+=； (C)10y xy '--=； (D)()d ()d 0x y x x y y -++=． 答(A). 5. 下列微分方程中，是齐次方程是微分方程的是( )． (A)x y y e +'=； 2(B)xy y x '+=； (C)0y xy x '--=； (D)()d ()d 0x y x x y y -++=． 答(D). 二、填空题 1．函数25y x =是否是微分方程2xy y '=的解 ． 答：是 . 2．微分方程3d d 0,4x x y y y x =+==的解是 ． 答：2225x y +=． 3．微分方程23550x x y '+-=的通解是． 答：32 52 x x y C =++． 4．微分方程ln 0xy y y '-=的通解是 ． 答： Cx y e =. 5'的通解是 ． 答：arcsin arcsin y x C =+． 6．微分方程 (ln ln )xy y y y x '-=-的通解是 ． 答：Cx y e x =． 三、解答题 1．求下列微分方程的通解． (1) 22sec tan d sec tan d 0x y x y x y +=； (2) 2()y xy a y y '''-=+； 解： 解： (3) d 10d x y y x +=； (4) 23d (1)0.d y y x x ++= 解： 解： 2．求下列微分方程满足所给初始条件的特解： (1) 20,0x y x y e y -='==； (2) 2sin ln ,x y x y y y e π='==； 解： 解： (3) 2d 2d 0,1x x y y x y =+==； (4) d 10d x y y x +=． 解： 解：

线性微分方程组

第五章 线性微分方程组 [教学目标] 1. 理解线性微分方程组解的存在唯一性定理，掌握一阶齐（非齐）线性微分方程组解的性质与结构， 2. 理解n 阶线性微分方程与一阶线性微分方程组的关系。 3. 掌握非齐次线性微分方程组的常数变易法， 4. 理解常系数齐线性微分方程组基解矩阵的概念，掌握求基解矩阵的方法。 5. 掌握常系数线性微分方程组的Laplce 变换法。 [教学中难点]求解常系数非齐次线性微分方程组 [教学方法] 讲授，实践。 [教学时间] 16学时 [教学内容] n 阶线性微分方程与一阶线性微分方程组的关系，一阶线性微分方程组解的存在唯一性定理；齐（非齐）线性微分方程组解的性质与结构，求解非齐次线性微分方程组的常数变易法；常系数齐线性微分方程组的基解矩阵及求基解矩阵的方法；求常系数线性微分方程组的Laplce 变换法。 [考核目标] 1.线性微分方程组解的性质与结构。 2.能够求解常系数线性微分方程组。 §5.1 存在唯一性定理 5.1.1记号和定义 考察形如 1 11112211221122222 1122()()()()()()()()()()()()n n n n n n n nn n n x a t x a t x a t x f t x a t x a t x a t x f t x a t x a t x a t x f t '=++++??'=++++?? ??'=++++? （5.1） 的一阶线性微分方程组，其中已知函数()(,1,2,,)ij a t i j n = 和()(1,2,,)i f t i n = 在区间a t b ≤≤上 上是连续的。方程组（5.1）关于12,,,n x x x 及1 2,,,n x x x ''' 是线性的. 引进下面的记号： 1112121 22 212()() ()()() ()()()() ()n n n n nn a t a t a t a t a t a t A t a t a t a t ??????=?? ? ? ?? （5.2） 这里()A t 是n n ?矩阵，它的元素是2 n 个函数()(,1,2,,)ij a t i j n = . 12()()()()n f t f t f t f t ??????=?????? 12n x x x x ??????=?????? 1 2n x x x x '????'??'=???? '?? （5.3）

(完整版)偏微分方程的MATLAB解法

引言 偏微分方程定解问题有着广泛的应用背景。人们用偏微分方程来描述、解释或者预见各种自然现象，并用于科学和工程技术的各个领域fll。然而，对于广大应用工作者来说，从偏微分方程模型出发，使用有限元法或有限差分法求解都要耗费很大的工作量，才能得到数值解。现在，MATLAB PDEToolbox已实现对于空间二维问题高速、准确的求解过程。 偏微分方程 如果一个微分方程中出现的未知函数只含一个自变量，这个方程叫做常微分方程，也简称微分方程；如果一个微分方程中出现多元函数的偏导数，或者说如果未知函数和几个变量有关，而且方程中出现未知函数对几个变量的导数，那么这种微分方程就是偏微分方程。 常用的方法有变分法和有限差分法。变分法是把定解问题转化成变分问题，再求变分问题的近似解；有限差分法是把定解问题转化成代数方程，然后用计算机进行计算；还有一种更有意义的模拟法，它用另一个物理的问题实验研究来代替所研究某个物理问题的定解。虽然物理现象本质不同，但是抽象地表示在数学上是同一个定解问题，如研究某个不规则形状的物体里的稳定温度分布问题，由于求解比较困难，可作相应的静电场或稳恒电流场实验研究，测定场中各处的电势，从而也解决了所研究的稳定温度场中的温度分布问题。 随着物理科学所研究的现象在广度和深度两方面的扩展，偏微分方程的应用范围更广泛。从数学自身的角度看，偏微分方程的求解促使数学在函数论、变分法、级数展开、常微分方程、代数、微分几何等各方面进行发展。从这个角度说，偏微分方程变成了数学的中心。

一、MATLAB方法简介及应用 1.1 MATLAB简介 MATLAB是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。 1.2 Matlab主要功能 数值分析 数值和符号计算 工程与科学绘图 控制系统的设计与仿真 数字图像处理 数字信号处理 通讯系统设计与仿真 财务与金融工程 1.3 优势特点 1) 高效的数值计算及符号计算功能，能使用户从繁杂的数学运算分析中解脱出来； 2) 具有完备的图形处理功能，实现计算结果和编程的可视化； 3) 友好的用户界面及接近数学表达式的自然化语言，使学者易于学习和掌握； 4) 功能丰富的应用工具箱(如信号处理工具箱、通信工具箱等) ，

一阶线性微分方程组

第4章 一阶线性微分方程组 一 内容提要 1． 基本概念 一阶微分方程组：形如 ??? ????? ???===) ,,,,( ),,,,(),,,,(2121222111 n n n n n y y y x f dx dy y y y x f dx dy y y y x f dx dy ΛΛΛΛΛ （3.1） 的方程组，（其中n y y y ,,,21Λ是关于x 的未知函数）叫做一阶微分方程组。 若存在一组函数)(,),(),(21x y x y x y n Λ使得在[a,b]上有恒等式 ),,2,1))((,),(),(,() (21n i x y x y x y x f dx x dy n i i ΛΛ==成立，则 )(,),(),(21x y x y x y n Λ称为一阶微分方程组(3.1)的一个解 含有n 任意常数n C C C ,,,21Λ的解 ?????? ?===) ,,,,( ),,,,(),,,,(21321222111n n n n C C C x y C C C x y C C C x y ΛΛΛΛΛ??? 称为(3.1)通解。如果通解满方程组 ???????=Φ=Φ=Φ0 ),,,,,,,,( 0),,,,,,,,(0),,,,,,,,(21212121221211n n n n n n n C C C y y y x C C C y y y x C C C y y y x ΛΛΛΛΛΛΛΛ 则称这个方程组为（3.1）的通积分。 满足初始条件,)(,,)(,)(0020021001n n y x y y x y y x y ===Λ的解，叫做初值问题的解。

二阶线性微分方程的解法

二阶常系数线性微分方程 一、二阶常系数线形微分方程的概念 形如 )(x f qy y p y =+'+'' (1) 的方程称为二阶常系数线性微分方程.其中p 、q 均为实数,)(x f 为已知的连续函数. 如果0)(≡x f ,则方程式 (1)变成 0=+'+''qy y p y (2) 我们把方程(2)叫做二阶常系数齐次线性方程,把方程式(1)叫做二阶常 系数非齐次线性方程. 本节我们将讨论其解法. 二、二阶常系数齐次线性微分方程 1．解的叠加性 定理1 如果函数1y 与2y 是式(2)的两个解, 则2211y C y C y +=也是 式(2)的解,其中21,C C 是任意常数. 证明 因为1y 与2y 是方程(2)的解,所以有 0111 =+'+''qy y p y 0222 =+'+''qy y p y 将2211y C y C y +=代入方程(2)的左边,得 )()()(22112211221 1y C y C q y C y C p y C y C ++'+'+''+'' =0)()(2222111 1=+'+''++'+''qy y p y C qy y p y C 所以2211y C y C y +=是方程(2)的解. 定理1说明齐次线性方程的解具有叠加性. 叠加起来的解从形式看含有21,C C 两个任意常数,但它不一定是方程式(2)的通解. 2.线性相关、线性无关的概念

设,,,,21n y y y 为定义在区间I 内的n 个函数,若存在不全为零的常数 ,,,,21n k k k 使得当在该区间内有02211≡+++n n y k y k y k , 则称这n 个函数在区间I 内线性相关,否则称线性无关. 例如 x x 22sin ,cos ,1在实数范围内是线性相关的,因为 0sin cos 12 2≡--x x 又如2,,1x x 在任何区间(a,b)内是线性无关的,因为在该区间内要使 02321≡++x k x k k 必须0321===k k k . 对两个函数的情形,若=21y y 常数, 则1y ,2y 线性相关,若≠2 1y y 常数, 则1y ,2y 线性无关. 3.二阶常系数齐次微分方程的解法 定理 2 如果1y 与2y 是方程式(2)的两个线性无关的特解,则 212211,(C C y C y C y +=为任意常数)是方程式(2)的通解. 例如, 0=+''y y 是二阶齐次线性方程,x y x y cos ,sin 21==是它的 两个解,且≠=x y y tan 2 1常数,即1y ,2y 线性无关, 所以 x C x C y C y C y cos sin 212211+=+= ( 21,C C 是任意常数)是方程0=+''y y 的通解. 由于指数函数rx e y =(r 为常数)和它的各阶导数都只差一个常数因子, 根据指数函数的这个特点,我们用rx e y =来试着看能否选取适当的常数r , 使rx e y =满足方程(2).

二阶常系数齐次线性微分方程求解方法

第六节 二阶常系数齐次线性微分方程 教学目的：使学生掌握二阶常系数齐次线性微分方程的解法，了解二阶常系数 非齐次线性微分方程的解法 教学重点：二阶常系数齐次线性微分方程的解法 教学过程： 一、二阶常系数齐次线性微分方程 二阶常系数齐次线性微分方程 方程 y py qy 0 称为二阶常系数齐次线性微分方程 其中p 、q 均为常数 如果y 1、y 2是二阶常系数齐次线性微分方程的两个线性无关解 那么y C 1y 1C 2y 2就是它的通解 我们看看 能否适当选取r 使y e rx 满足二阶常系数齐次线性微分方程 为此将y e rx 代入方程 y py qy 0 得 (r 2pr q )e rx 0 由此可见 只要r 满足代数方程r 2pr q 0 函数y e rx 就是微分方程的解 特征方程 方程r 2pr q 0叫做微分方程y py qy 0的特征方程 特征方程的两个根r 1、r 2可用公式 2 422,1q p p r -±+-= 求出 特征方程的根与通解的关系 (1)特征方程有两个不相等的实根r 1、r 2时 函数x r e y 11=、x r e y 22=是方程的两个线性无关的解 这是因为

函数x r e y 11=、x r e y 22=是方程的解 又x r r x r x r e e e y y )(212121-==不是常数 因此方程的通解为 x r x r e C e C y 2121+= (2)特征方程有两个相等的实根r 1r 2时 函数x r e y 11=、x r xe y 12=是二阶常系数齐次线性微分方程的两个线性无关的解 这是因为 x r e y 11=是方程的解 又 x r x r x r x r x r x r qxe e xr p e xr r xe q xe p xe 111111)1()2()()()(1211++++=+'+'' 0)()2(121111 =++++=q pr r xe p r e x r x r 所以x r xe y 12=也是方程的解 且x e xe y y x r x r ==1112不是常数 因此方程的通解为 x r x r xe C e C y 1121+= (3)特征方程有一对共轭复根r 1, 2i 时 函数y e ( i )x 、y e (i )x 是微分方程的两个线性无关的复数形式的解 函数y e x cos x 、y e x sin x 是微分方程的两个线性无关的实数形式的解 函数y 1e (i )x 和y 2e (i )x 都是方程的解 而由欧拉公式 得 y 1e ( i )x e x (cos x i sin x ) y 2e (i )x e x (cos x i sin x ) y 1y 22e x cos x )(21cos 21y y x e x +=βα y 1y 2 2ie x sin x )(21sin 21y y i x e x -=βα 故e x cos x 、y 2e x sin x 也是方程解 可以验证 y 1e x cos x 、y 2e x sin x 是方程的线性无关解 因此方程的通解为 y e x (C 1cos x C 2sin x )

偏微分方程组解法

偏微分方程组解法 某厚度为10cm 平壁原温度为20C ?，现其两侧面分别维持在20C ?和120C ?，试求经过8秒后平壁温度分布，并分析温度分布随时间的变化直至温度分布稳定为止。 22x t a t ??=??τ 式中a 为导温系数，/s m c 2；2=a 。 解： 模型转化为标准形式： 2 21x t t a ??=??τ 初始条件为： ()200,=x t 边界条件为： ()120,0=τt ，()20,1.0=τt 函数： pdefun.m %偏微分方程(一维动态传热) function [c,f,s]=pdefun(x,t,u,dudx) c=1/2e-4;f=dudx;s=0; icbun.m %偏微分方程初始条件(一维动态传热) function u0=icbun(x) u0=20; bcfun.m %偏微分方程边界条件(一维动态传热) function [pl,ql,pr,qr]=bcfun(xl,ul,xr,ur,t) pl=ul-120;ql=0;pr=ur-20;qr=0; 命令： x=linspace(0,10,20)*1e-2; t=linspace(0,15,16); sol=pdepe(0,pdefun,icfun,bcfun,x,t); mesh(x,t,sol(:,:,1)) %温度与时间和空间位置的关系图 %画1、2、4、6、8、15s 时刻温度分布图

plot(x,sol(2,:,1)) 1s时刻，（因为本题sol第一行为0时刻） hold on plot(x,sol(3,:,1)) plot(x,sol(5,:,1)) plot(x,sol(7,:,1)) plot(x,sol(9,:,1)) plot(x,sol(16,:,1)) 计算结果： %第8秒时温度分布 x sol(9,:,1) 经过8秒时的温度分布为： x/cm 0 0.5263 1.0526 1.5789 2.1053 2.6316 3.1579 t/C ?120.0000 112.5520 105.1653 97.8994 90.8100 83.9477 77.3562 x/cm 3.6842 4.2105 4.7368 5.2632 5.7895 6.3158 6.8421 t/C ?71.0714 65.1202 59.5200 54.2784 49.3930 44.8518 40.6338 x/cm 7.3684 7.8947 8.4211 8.9474 9.4737 10.0000 t/C ?36.7095 33.0419 29.5877 26.2982 23.1207 20.0000 或者求第8秒时，x=0,2,4,,6,8,10cm处的温度 [uout,duoutdx]=pdeval(0,x,sol(9,:,:),[0,2,4,6,8,10]*1e-2) 120.0000 92.2279 67.5007 47.5765 32.3511 20.0000

系数非线性常微分方程的特解表达式

万方数据

万方数据

万方数据

万方数据

三类常系数非线性常微分方程的特解表达式 作者：陈友朋， 钱明忠， 黄娟娟 作者单位：江苏省盐城师范学院数学科学学院,江苏盐城,224051 刊名： 高等数学研究 英文刊名：STUDIES IN COLLEGE MATHEMATICS 年，卷(期)：2009，12(4) 被引用次数：0次 参考文献(3条) 1.张建梅.孙志田.崔宁关于y″+py'+qy=Aeαx的特解[期刊论文]-高等数学研究 2005(03) 2.曾菊华.胡小英关于常系数线性微分方程的特解表达式[期刊论文]-高等数学研究 2006(04) 3.Π Э 艾利斯哥尔兹.南开大学数学系编译中队.崔士英微分方程 1959 相似文献(10条) 1.期刊论文刘琳琳非齐次常系数常微分方程特解形式的一个推导-喀什师范学院学报2002,23(3) 考虑n阶非齐次常系数线性常微分方程y(n)+Pn-1y(n-1)+…+p1y1+poy=f(x),当它的右端项f(x)=eλχPm(x)时,给出它的特解形式的推导. 2.期刊论文张学凌.王志伟求一类常微分方程特解的程序化方法-天中学刊2008,23(5) 通过对常微分方程常规解法的进一步探讨,推导出一类三阶常系数非齐次线性微分方程求特解的统一表达式,并通过C++语言编程,利用计算机直接输出结果,提高了求解的速度和准确性. 3.期刊论文沈彻明.SHEN Che-ming求非齐次高阶常系数线性常微分方程的特解的一般公式-数学的实践与认识2000,30(4) 本文提出了高阶常系数线性常微分方程的第二类特征代数方程,并利用它获得了求非齐次方程的特解的一般公式. 4.期刊论文赵苏串一类常系数非齐次常微分方程的特解的求法-上海大学学报(自然科学版)1999,5(6) 讨论了形如u+αu=f(x),u(4)+αu.+βu=f(x),其中f(x)=(sinωx)2k或(cosωx)2k(k∈Z+),ω≠0ε,α,β均为常数的特解的求法. 5.期刊论文龚东山.刘岳巍.贾筱景.GONG Dong-shan.LIU Yue-wei.JIA Xiao-jing计算一类常微分方程特解的新方法-河北北方学院学报(自然科学版)2008,24(6) 目的 计算高阶常微分方程特解的方法有待定系数法、常数变易法、拉普拉斯变换法、积分法等,它们的计算工作量一般较大,为弥补上述方法的不足,有必要探究另一种简便实用的新方法--特征函数法.方法 先定义该类高阶常微分方程的对应齐次方程的特征函数,再利用特征函数的导数,可得到非齐次项为特殊函数情形时方程的一个特解.结果 只需求出特征方程的根,就可得到该类高阶常微分方程的一个特解.结论 利用特征函数法可以得到一类常微分方程的一个特解,该方法使用简单,所得特解形式直观. 6.期刊论文龚东山.刘岳巍.牛富俊.GONG Dong-shan.LIU Yue-wei.NIU Fu-jun特征函数在高阶常微分方程特解计算中的应用-吉林师范大学学报(自然科学版)2008,29(4) 通过借助特征函数的导数,得到了非齐次项为特殊函数情形的一类高阶常微分方程的一个特解的一种新的计算方法.运用该方法,还得到了非齐次项为常见情形时方程的一个特解. 7.期刊论文陈新一一类二阶常微分方程的特解 -高等数学研究2010,13(1) 研究一类二阶实常系数非齐次微分方程y″+py′+q=(a0+a1x)eαxsinβx的解法,应用叠加原理和Euler公式,将其化为二阶线性非齐次方程,并利用对应的特征方程给出了这一类方程特解的一般公式,简化这一类微分方程的求解过程. 8.期刊论文张学凌二阶非齐次线性常微分方程特解的算法模型-许昌学院学报2003,22(2) 用迭代算法求二阶非齐次线性常微分方程y"+py'+qy=pn(x)eax=(AnXn+…+Aixi+…+Ao)eax的特解是一种新的尝试,借助C++BUILDER编译器成功地实现了该算法,较圆满地解决了此类微分方程求特解时实际计算上的问题. 9.期刊论文王欣欣.郑秉文用微分算子求常微分方程特解的注记-吉林师范大学学报(自然科学版)2003,24(3) 本文给出常系数线性微分方程最简特解的定义,论证了常系数线性微分方程最简特解的形式,同时给出了用微分算子求常系数线性微分方程最简特解的方法. 10.期刊论文陈新一.唐文玲.CHEN Xin-yi.TANG Wen-ling一类三阶常微分方程的特解公式-甘肃联合大学学报（自然科学版）2007,21(1) 利用比较系数法,推导出三阶常系数微分方程y"'+py"+qy'+ry=(a0+a1x+a2x2)eλx的特解的一般公式.利用这个公式可直接得到此类微分方程的特解. 本文链接：https://www.wendangku.net/doc/9810494261.html,/Periodical_gdsxyj200904014.aspx 授权使用：中共汕尾市委党校(zgsw)，授权号：0494467a-5728-47be-9cc4-9dcf0154b484 下载时间：2010年8月11日

一维偏微分方程的pdepe(matlab)函数 解法

本文根据matlab帮助进行加工，根据matlab帮助上的例子，帮助更好的理解一维偏微分方程的pdepe函数解法，主要加工在于程序的注释上。 Examples Example 1.This example illustrates the straightforward formulation, computation, and plotting of the solution of a single PDE. This equation holds on an interval for times . The PDE satisfies the initial condition and boundary conditions It is convenient to use subfunctions to place all the functions required by pdepe in a single function. function pdex1 m = 0; x = linspace(0,1,20); %linspace(x1,x2,N)linspace是Matlab中的一个指令，用于产生x1,x2之间的N点行矢量。 %其中x1、x2、N分别为起始值、终止值、元素个数。若缺省N，默认点数为100 t = linspace(0,2,5); sol = pdepe(m,@pdex1pde,@pdex1ic,@pdex1bc,x,t);

% Extract the first solution component as u. u = sol(:,:,1); % A surface plot is often a good way to study a solution. surf(x,t,u) title('Numerical solution computed with 20 mesh points.') xlabel('Distance x') ylabel('Time t') % A solution profile can also be illuminating. figure plot(x,u(end,:)) title('Solution at t = 2') xlabel('Distance x') ylabel('u(x,2)') % -------------------------------------------------------------- function [c,f,s] = pdex1pde(x,t,u,DuDx) c = pi^2; f = DuDx; s = 0; % -------------------------------------------------------------- function u0 = pdex1ic(x) u0 = sin(pi*x); % -------------------------------------------------------------- function [pl,ql,pr,qr] = pdex1bc(xl,ul,xr,ur,t) pl = ul; ql = 0; pr = pi * exp(-t); qr = 1;