齐次微分方程的解法

关于数值分析的

常微分方程的数值解法

一、题目 2x y y 求解初值问题 y

y 0 1 0 x 1 ，10等分区间，求节点处的近似值，并对所求结

果与分析解的结果进行比较。

二、方法

欧拉法

三、程序

function E=euler(f,a,b,y0,N)

x=zeros(1,N+1);

y=zeros(1,N+1);

x(1)=a;

y(1)=y0;

h=(b-a)/N;

for n=1:N

x(n+1)=x(n)+h;

y(n+1)=y(n)+h*feval(f,x(n),y(n));

end

T=[x',y']

四、结果

>> format compact

>> euler(inline('y-2*x/y'),0,1,1,10)

T =

0 1.0000

0.1000 1.1000

0.2000 1.1918

0.3000 1.2774

0.4000 1.3582

0.5000 1.4351

0.6000 1.5090

0.7000 1.5803

0.8000 1.6498

0.9000 1.7178

1.0000 1.7848

>> y=[1.0000 1.1000 1.1918 1.2774

第六章 常微分方程的数值解法

§6.0 引言

§6.1 算法构造的主要途径 §6.2 Runge-Kutta Method算法 §6.3 线性多步法

§6.4 线性多步法的一般形式 §6.5 算法的稳定性、收敛性

§6.0 引 言

1． 主要考虑如下的一阶常微分方程初值问题的求解:

?dy?dx?f?x,y????y?x0??y0??

微分方程的解就是求一个函数y=y(x)，使得该函数满足微分方程并且符合初值条件。 2. 例如微分方程:

xy'-2y=4x ；初始条件: y(1)=-3。

于是可得一阶常微分方程的初始问题

???y??2y??x4?y(1)??3。

显然函数y(x)=x2-4x满足以上条件,因而是该初始问题 的微分方程的解。

3. 但是,只有一些特殊类型的微分方程问题能

够得到用解析表达式表示的函数解，而大量的微分方程问题很难得到其解析解，有的甚至无法用解析表达式来表示。因此，只能依赖于数值方法去获得微分方程的数值解。 4． 微分方程的数值解：设微分方程问题的解y(x)的存在区间是[a,b]，初始点x0=a，将[a,b]进行划分得一系列节点x0 , x1 ,...,xn，其中a= x0< x1<…< xn =b。

y(x)的解析表达式不容易得到或根本无

第八章

常微分方程数值解法

摘要：对显式Euler方法来说,当解二阶连续可导时,其局部...(3.10)有解但解不唯一.不论如何选择这八个参数,不可能...算法8.1 经典Runge-Kutta方法本算法用经典Runge-... 关键词：导,论,算法 类别：专题技术

来源：牛档搜索（Niudown.COM）

本文系牛档搜索（Niudown.COM）根据用户的指令自动搜索的结果，文中内涉及到的资料均来自互联网，用于学习交流经验，作品其著作权归原作者所有。不代表牛档搜索（Niudown.COM）赞成本文的内容或立场，牛档搜索（Niudown.COM）不对其付相应的法律责任！

常微分方程数值解法

教学目的 1. 掌握解常微分方程的单步法：Euler方法、Taylor方法和Runge-Kutta方法；2. 掌握解常微分方程的多步法：Adams步法、Simpson方法和Milne方法等；3. 了解单步法的收敛性、相容性与稳定性；多步法的稳定性。

教学重点及难点 重点是解常微分方程的单步法：Euler方法、Taylor方法和Runge-Kutta方法和解常微分方程的多步法：Adams步法、Simpson方法和Milne方法等；难点是理

第九章 常微分方程的数值解法主要内容§1、引言

§2、初值问题的数值解法--单步法§3、龙格-库塔方法

§4、收敛性与稳定性§5、初值问题的数值解法―多步法 §6、方程组和刚性方程 §7、习题和总结

§1、 引

言

主要内容

研究的问题 数值解法的意义

1.什么是微分方程 ? 现实世界中大多数事物 内部联系非常复杂

其状态随着 时间、地点、条件 的不同而不同 找出其状态和状态变化规律之间的相互联系， 也即一个或一些函数与他们的导数之间的关系

此种关系的数学表达就为

微分方程

2.数值求解微分方程的意义如何建立数学模型已在建模课程中得到讨论， 各类微分方程本身和他们的解所具有的特性 已在常微分方程及数学物理方程中得以解释， 本章专门 讨论

如何利用数值方法求解微分方程(组)的问题。

3.什么是微分方程 (组)的解析解?

3.什么是微分方 程(组)的解析解?

一个具有所要求阶连续导数的解析函数，将 它代入微分方程，恰使其所有条件都得到满 足的解称为解析解(或古典解),称为真解或解。 寻找解析解的过程称为求解微分方程。 y f ( x , y ), y ( x 0 ) y0

第七章 常微分方程数值解法

常微分方程中只有一些典型方程能求出初等解(用初等函数表示的解)，大部分的方程是求不出初等解的。另外，有些初值问题虽然有初等解，但由于形式太复杂不便于应用。因此，有必要探讨常微分方程初值问题的数值解法。本章主要介绍一阶常微分方程初值问题的欧拉法、龙格-库塔法、阿达姆斯方法，在此基础上推出一阶微分方程组与高阶方程初值问题的 数值解法；此外，还将简要介绍求解二阶常微分方程值问题的差分方法、试射法。

第一节 欧拉法

求解常微分方程初值问题

?dy??f(x,y) ?dx??y(x0)?y0 （1）

的数值解，就是寻求准确解y(x)在一系列离散节点x0?x1?x2???xn?? 上的近似值 y0,y1,y2,?,yn,?

?yn?称为问题的数值解，数值解所满足的离散方程统称为差分格式，hi称为步长，实用中常取定步长。

?xi?xi?1显然，只有当初值问题(1)的解存在且唯一时，使用数值解法才有意义，这一前提条件由下 面定理保证。

定理 设函数f?x,y?在区域D:a?x?b,???y???

上连续，且在区域D内

第十二章：微分方程

教学目的：

1．了解微分方程及其解、阶、通解，初始条件和特等概念。

2．熟练掌握变量可分离的微分方程及一阶线性微分方程的解法。

3．会解齐次微分方程、伯努利方程和全微分方程，会用简单的变量代换解某些微分方程。

4.会用降阶法解下列微分方程：()()n y f x =， (,)y f x y '''+和(,)y f y y '''=

5.理解线性微分方程解的性质及解的结构定理。

6．掌握二阶常系数齐次线性微分方程的解法，并会解某些高于二阶的常系数齐次线性微分方程。

7.求自由项为多项式、指数函数、余弦函数，以及它们的和与积的二阶常系数非齐次线性微分方程的特解和通解。

8.会解欧拉方程，会解包含两个未知函数的一阶常系数线性微分方程组。

9．会解微分方程组（或方程组）解决一些简单的应用问题。

教学重点：

1、可分离的微分方程及一阶线性微分方程的解法

2、可降阶的高阶微分方程()()n y f x =， (,)y f x y '''+和(,)y f y y '''=

3、二阶常系数齐次线性微分方程；

4、自由项为多项式、指数函数、余弦函数，以及它们的和与积的二阶常系数非齐次线性微分方程；

教学难点：

1、齐次微分方程、伯努利方程和全微分方程；

2、线性微分方程解的性

第十二章：微分方程

教学目的：

1．了解微分方程及其解、阶、通解，初始条件和特等概念。

2．熟练掌握变量可分离的微分方程及一阶线性微分方程的解法。

3．会解齐次微分方程、伯努利方程和全微分方程，会用简单的变量代换解某些微分方程。

4.会用降阶法解下列微分方程：()()n y f x =， (,)y f x y '''+和(,)y f y y '''=

5.理解线性微分方程解的性质及解的结构定理。

6．掌握二阶常系数齐次线性微分方程的解法，并会解某些高于二阶的常系数齐次线性微分方程。

7.求自由项为多项式、指数函数、余弦函数，以及它们的和与积的二阶常系数非齐次线性微分方程的特解和通解。

8.会解欧拉方程，会解包含两个未知函数的一阶常系数线性微分方程组。

9．会解微分方程组（或方程组）解决一些简单的应用问题。

教学重点：

1、可分离的微分方程及一阶线性微分方程的解法

2、可降阶的高阶微分方程()()n y f x =， (,)y f x y '''+和(,)y f y y '''=

3、二阶常系数齐次线性微分方程；

4、自由项为多项式、指数函数、余弦函数，以及它们的和与积的二阶常系数非齐次线性微分方程；

教学难点：

1、齐次微分方程、伯努利方程和全微分方程；

2、线性微分方程解的性

论文常系数线性微分方程的解法

常系数线性微分方程的解法

摘 要

本文主要介绍了常系数线性微分方程的解法。由于在讨论常系数线性微分方程的解法时，需要涉及实变量的复值函数及复值数函数问题，所以文章首先给予介绍。然后，本文又通过构造特征方程运用代数运算分情况给出常系数齐次线性微分方程的具体解法，并引出可以化为此方程形式的欧拉方程。有了前面讨论的结果，对于常系数非齐次线性微分方程可以采用常数变易法，这里没有具体介绍，而是介绍具有某些特殊形式的非齐次线性微分方程的解法，即比较系数法和拉普拉斯变换法。

关键词：复值函数与复指数函数，齐次线性微分方程，欧拉方程，非齐次线性微分方程，比较系数法，拉普拉斯变换法

The Methods of Solving Linear Differential Equation with Constant Coefficients

This paper describes the methods of solving linear differential equation with constant coefficients. First of all, the paper below will describe complex-va

淮北师范大学 2013届学士学位论文

常微分方程数值解法的误差分析

学院、专业 数学科学学院 数学与应用数学 研 究 方 向 计算数学 学 生 姓 名 李 娜 学 号 20091101070 指导教师姓名 陈 昊 指导教师职称 讲 师

年 月 日

常微分方程数值解法的误差分析

李 娜

（淮北师范大学数学科学学院，淮北，235000）

摘 要

自然界与工程技术中的很多现象，往往归结为常微分方程定解问题。许多偏微

分方程问题也可以化为常微分方程问题来近似求解。因此，研究常微分方程的数值解法是有实际应用意义的。数值解法是一种离散化的数学方法，可以求出函数的精确解在自变量一系列离散点处的近似值。随着计算机计算能力的增强以及数值计算方法的发展，常微分方程的数值求解方法越来越多，比较成熟的有Euler法、后退Euler法、梯形方法、Runge

非线性方程和常微分方程的解法

实验8 非线性方程和常微分方程的解法

一、实验目的

学会用MATLAB软件求解非线性方程和常微分方程。

二、实验内容与要求

1. 非线性方程的整值解

（1）最小二乘法

格式：fsolve（’fun’,x0）%求方程fun=0在估计值x附近的近似解。

[例1.72] 求方程x e 0的解。

>>fc=inline(‘x-exp(-x)’);

>>xl=fsolve(fc,0)

xl=

0.5671

问题1.28：求解方程5xsinx-e 0,观察知有多解，如何求之？

先用命令fplot（’5*x^2*sin(x)-exp(-x),0]’,[0,10]）作图1.13，注意5*x^2*sin(x)-exp(-x),0]中的“[ ，0]”是作y=0直线，即x轴。方程在[0，10]区间从图中可看出有4个解，分别在0，3，6，9附近，所以用命令：

>>fun=inline(’5*x^2*sin(x)-exp(-x)’);

>>fsolve(fun,[0,3,6,9],le-6)

得出结果：

ans=

0.5918 3.1407 6.2832 9.4248

【例 1.73】求解方程组x-0.7sinx-0.2co