二阶非线性中立型微分方程的Sturm比较定理

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

高阶线性微分方程常用解法简介 关键词：高阶线性微分方程 求解方法 在微分方程的理论中，线性微分方程是非常值得重视的一部分内容，这不仅 因为线性微分方程的一般理论已被研究的十分清楚，而且线性微分方程是研究非线性微分方程的基础，它在物理、力学和工程技术、自然科学中也有着广泛应用。下面对高阶线性微分方程解法做一些简单介绍. 讨论如下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 引言 (1) 2 二阶常系数常微分方程的几种解法 (1) 2.1 特征方程法 (1) 2.1.1 特征根是两个实根的情形 (2) 2.1.2 特征根有重根的情形 (2) 2.2 常数变异法 (4) 2.3 拉普拉斯变化法 (5) 3 常微分方程的简单应用 (6) 3.1 特征方程法 (7) 3.2 常数变异法 (9) 3.3 拉普拉斯变化法 (10) 4 总结及意义 (11) 参考文献 (12)

二阶常微分方程的解法及其应用 摘要：本文通过对特征方程法、常数变易法、拉普拉斯变换法这三种二阶常系数常微分方程解法进行介绍，特别是其中的特征方程法分为特征根是两个实根的情形和特征根有重根的情形这两种情况，分别使用特征值法、常数变异法以及拉普拉斯变换法来求动力学方程，现今对于二阶常微分方程解法的研究已经取得了不少成就，尤其在二阶常系数线性微分方程的求解问题方面卓有成效。应用常微分方程理论已经取得了很大的成就，但是，它的现有理论也还远远不能满足需要，还有待于进一步的发展，使这门学科的理论更加完善。 关键词：二阶常微分方程;特征分析法；常数变异法；拉普拉斯变换

METHODS FOR TWO ORDER ORDINARY DIFFERENTIAL EQUATION AND ITS APPLICATION Abstract:This paper introduces the solution of the characteristic equation method, the method of variation of parameters, the Laplasse transform method the three kind of two order ordinary differential equations with constant coefficients, especially the characteristic equation method which is characteristic of the root is the two of two real roots and characteristics of root root, branch and don't use eigenvalue method, method of variation of constants and Laplasse transform method to obtain the dynamic equation, the current studies on solution of ordinary differential equations of order two has made many achievements, especially in the aspect of solving the problem of two order linear differential equation with constant coefficients very fruitful. Application of the theory of ordinary differential equations has made great achievements, however, the existing theory it is still far from meeting the need, needs further development, to make the discipline theory more perfect. Keywords:second ord er ordinary differential equation; Characteristic analysis; constant variation method; Laplasse transform 1 引言 数学发展的历史告诉我们，300年来数学分析是数学的首要分支，而微分方程

二阶线性微分方程的解法

二阶常系数线性微分方程 一、二阶常系数线形微分方程的概念 形如 )(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).

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

高阶线性微分方程常用解法简介 摘要：本文主要介绍高阶线性微分方程求解方法，主要的内容有高阶线性微分方程求解的常 用方法如。 关键词：高阶线性微分方程 求解方法 在微分方程的理论中，线性微分方程是非常值得重视的一部分内容，这不仅 因为线性微分方程的一般理论已被研究的十分清楚，而且线性微分方程是研究非线性微分方程的基础，它在物理、力学和工程技术、自然科学中也有着广泛应用。下面对高阶线性微分方程解法做一些简单介绍. 讨论如下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 为任意常数. 如果特征方程有复根，则因方程的系数是实常数，复根将称对共轭的出现.

变系数_非线性微分方程的求解

变系数/非线性微分方程的求解：Example1: van der Pol equation Rewrite the van der Pol equation (second-order) The resulting system of first-order ODEs is 见：vdp_solve.m及vdp.mdl vdp_solve.m vdp.mdl

Example2: 2 with x(0) = 4 x (0)=0 5(5)5sin()5 +-+= x t x t x 见：exam2_solve.m及exam2.mdl exam2_solve.m exam2.mdl

Example3: ODEs 函数实现及封装说明[以一阶微分方程为例] 510 w i t h (0)4 dx x x dt +==- 引言： 一步Euler 法求解[相当于Taylor 展开略去高阶项]： 11()k k k k k k k k k k k x x x Ax bu t x x t x x t Ax bu ++-==+??=+??=+??+ 补充说明1：对于任意方程/方程组可化为如下一阶形式[方程组]： x Ax Bu =+ 或者(,)(,)M t x x f t x = 补充说明2：ODEs 的解法不同之处在于 1、时间步长的选取（及导数的求解？）：有无误差控制 变步长； 2、积分方法：选用哪几个时间状态信息。 见：my_ode_rough.m[直接求解] / test_my_ode.m[按Matlab/ODEs 方式封装] my_ode_rough.m

二阶非线性常微分方程的打靶法matlab实现

二阶非线性常微分方程的打靶法 1.问题： 试用打靶法求二阶非线性常微分方程亮点边值的数值解： 要求用Matlab 编程计算，请给出一些例子，验证你的算法与程序的正确性。 2.打靶法分析： 非线性打靶法: 非线性打靶法的基本原理是将两点边值问题(1)转化为下面形式的初值问题 令z = y ′,将上述二阶方程降为一阶方程组

3.Matlab 源代码： 创建M 文件： function ys=dbf(f,a,b,alfa,beta,h,eps) ff=@(x,y)[y(2),f(y(1),y(2),x)]; xvalue=a:h:b;%x取值范围 n=length(xvalue) s0=a-0.01;%选取适当的s的初值 x0=[alfa,s0];%迭代初值 flag=0;%用于判断精度 y0=rk4(ff,a,x0,h,a,b); if abs(y0(1,n)-beta)<=eps flag=1; y1=y0; else s1=s0+1; x0=[alfa,s1]; y1=rk4(ff,a,x0,h,a,b); if abs(y1(1,n)-beta)<=eps flag=1; end end if flag~=1 while abs(y1(1,n)-beta)>eps s2=s1-(y1(1,n)-beta)*(s1-s0)/(y1(1,n)-y0(1,n)); x0=[alfa,s2]; y2=rk4(ff,a,x0,h,a,b); s0=s1; s1=s2;

y0=y1; y1=y2; end end xvalue=a:h:b; yvalue=y1(1,:); ys=[xvalue',yvalue']; function x=rk4(f,t0,x0,h,a,b)%rung-kuta法求每个点的近似值（参考大作业一）t=a:h:b;%迭代区间 m=length(t);%区间长度 t(1)=t0; x(:,1)=x0;%迭代初值 for i=1:m-1 L1=f(t(i),x(:,i)); L2=f(t(i)+h/2,x(:,i)'+(h/2)*L1); L3=f(t(i)+h/2,x(:,i)'+(h/2)*L2); L4=f(t(i)+h,x(:,i)'+h*L3); x(:,i+1)=x(:,i)'+(h/6)*(L1+2*L2+2*L3+L4); end 4.举例 求二阶非线性方程的边值问题： 在matlab 控制台中输入： f=@(x,y,z)(x^2+z*x^2); x0l=0; x0u=2*exp(-1); alfa=0; beta=2; h=0.01 dbf(f,x0l,x0u,y0l,y0u,h,1e-6); >> y=ans(:,2); x=ans(:,1); >> plot(x,y,'-r') >> 结果：

二阶非齐次线性微分方程的解法

目 录 待定系数法 常数变异法 幂级数法 特征根法 升阶法 降阶法 关键词：微分方程，特解，通解， 二阶齐次线性微分方程 常系数微分方程 待定系数法 解决常系数齐次线性微分方程[]21220, (1) d x dx L x a a x dt dt ≡++= 12,. a a 这里是常数 特征方程 212()0 F a a λλλ=++= (1.1) (1)特征根是单根的情形 设 12,,,n λλλL 是特征方程的 (1.1)的2个彼此不相等的根，则相应的方程 (1)有如 下2个解： 12,t t e e λλ (1.2) 如果(1,2)i i λ=均为实数，则 (1.2)是方程 (1)的2个线性无关的实值解，而方程 (1)的通解可表示为 1212t t x c e c e λλ=+ 如果方程有复根，则因方程的系数是实系数，复根将成对共轭出现。设 i λαβ=+是一特征根，则i λαβ=-也是特征根，因而与这对共轭复根对应，方程 (1)有两个复值解 (i)t (cos t sin ),t e e i t αβαββ+=+ (i)t (cos t sin ).t e e i t αβαββ-=-

它们的实部和虚部也是方程的解。这样一来，对应于特征方程的一对共轭复根 i λαβ=±，我们可求得方程 (1)的两个实值解 cos ,sin .t t e t e t ααββ （2）特征根有重跟的情形 若10λ=特征方程的 k 重零根，对应于方程 (1)的k 个线性无关的解211,t,t ,k t -L 。 若这个k 重零根1 0,λ≠设特征根为12,,,,m λλλL 其重数为1212,,,k (k 2) m m k k k k ++=L L 。方程 (1)的解为 11112222111,t ,t ;,t ,t ;;,t ,t ;m m m m t t k t t t k t t t k t e e e e e e e e e λλλλλλλλλ---L L L L 对于特征方程有复重根的情况，譬如假设i λαβ=+是k 重特征根，则i λαβ=-: 也是k 重特征根，可以得到方程 (1)的2k 个实值解 2121cos ,cos ,cos ,,cos ,sin ,sin ,sin ,,sin .t t t k t t t t k t e t te t t e t t e t e t te t t e t t e t ααααααααββββββββ--L L 例1 求方程 220d x x dt -=的通解。 解 特征方程 210λ-=的根为121,1λλ==-有两个实根，均是单根，故方程的通 解为 12, t t x c e c e -=+ 这里 12 ,c c 是任意常数。 例2 求解方程 220d x x dt +=的通解。 解 特征方程 210λ+=的根为12,i i λλ==-有两个复根， 均是单根，故方程的通解 为 12sin cos , x c t c t =+ 这里12 ,c c 是任意常数。

非线性微分方程和稳定性

第六章 非线性微分方程和稳定性 6-1 对下列方程求出常数特解，并且画出方程经过()0,0x 的积分曲线的走向，从而判断各驻定解的稳定性；然后作变量替换，使非零驻定解对应于新的方程的零解。 1）+∞<<-∞>>+=02 ,0,0,x B A Bx Ax dt dx 2） ()()0,310≥--=x x x x dt dx 解 1）方程可化为 )( x B A Bx dt dx +=，则其常数特解为 B A x x - ==21,0，即为驻定解。 由于方程为分离变量方程（或迫努利方程），当B A x x - ≠≠,0时，分离变量得 Adt dx B A x x =? ????? ? ?+-11 方程的通解为 At Ce Bx A x =+ 利用初始条件()?? ? ? ?-≠≠=B A x x x x 000,00，得 0 0Bx A x C += ，故得原方程满足初始 条件的解为 ()0)(0≥? ?? ? ??++-= -t e B x A B A t x At （1） 由式（1）和方程右端的表达式，得出 当00>x 时，0>dt dx ，)(t x 递增， 又 B e B x A B B x A At →??? ? ??+->+-00 ,时，+∞→)(t x ， 即)1ln( 10+= →B x A A t t 时，+∞→)(t x 。

当 ???????<-><+>-<>+<0 00,000 00 0 dt dx ,B A x , B x A dt dx ,B A x B x A x 时，有 ()+∞→- →t B A t x )( 所以解（1）的图像如图6-5所示。 图6-5 从解的图像可以看出： 解01=x 不稳定；解B A x -=2稳定。 利用变换B A x y +=，可将原方程化为 2 2 )()(By Ay B A y B B A y A dt dy +-=-+- = 所以原方程的驻定解B A x - =2对应于方程 2 By Ay dt dy +-= 的零解0=y 。 2）由()()031=--x x x ，求得常数解为 3,1,0321===x x x 。 因为()()()31,--=x x x x t f 在全平面上连续可微，故对任意初始点()00,x t ，解唯一存在，当0,0≥≥x t 时有

二阶常微分方程的解法

南京师范大学泰州学院 毕业论文（设计） （一六届） 题目：二阶常微分方程的解法 院（系、部）：数学科学与应用学院 专业：数学与应用数学 姓名：潘陆 学号08120146 指导教师：刘陆军 南京师范大学泰州学院教务处制

摘要：本文主要是介绍了二阶常微分方程众多解法中的三种，分别为特征方程法，拉普拉斯变换法和常数变易法，研究并讨论了二阶常微分方程在特征方程法中特征方程根为实根，复根和重根的情形。我们选用了弹簧振子系统的振子运动，用这三种不同的方法来解决该问题。 关键词：二阶常微分方程；特征根法；常数变易法；拉普拉斯变换 Abstract:The main purpose of this paper is the second-order ordinary many differential equation solution of three, respectively as the characteristic equation method, Laplace transform method and variation of constants method, study and discuss the second-order often differential equation in the characteristic equation of the roots of the characteristic equation for real roots, complex roots and root weight. We choose the spring oscillator the oscillator motion, these three different methods to solve the problem. Keywords: second order ordinary differential equation; Characteristic analysis; constant variation method; Laplasse transform

二阶常系数线性微分方程的解法word版

第八章 8.4讲 第四节 二阶常系数线性微分方程 一、二阶常系数线形微分方程的概念 形如 )(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 2 2sin ,cos ,1在实数范围内是线性相关的,因为 0s i n c o s 12 2≡--x x 又如2,,1x x 在任何区间(a,b)内是线性无关的,因为在该区间内要使 02 321≡++x k x k k 必须0321===k k k . 对两个函数的情形,若 =2 1y 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 ,