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

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

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

摘 要

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

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

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

第!"卷第"期宝鸡文理学院学报#自然科学版$

#$,&-./0’&120&345&’’676&18.9:0/;<=46/=6)09-.0’<=46/=6

(!")&("%&’

(!**">0.

!**"年+月

????????????????????????????????????????????????????????????????

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

杨继明

玉溪师范学院数学系A云南玉溪B#C+"**$

@

摘

要D给出了常系数齐次线性微分方程组初值问题的一个求解公式A并由此推出常系数齐次线性

差分方程组在给定的初始条件下的一个求解公式E

关键词D常系数F线性微分方程F线性差分方程F标准基解矩阵F矩阵的方幂中图分类号D"HC("G

文献标识码D8

文章编号D"**HI"!B"#!**"$*"I**"+I*+

JKLMNOPQRLSLPTUVWVRVQPTQPOUXMLYPUNLK

ZLNL[UWULOSPVWUQM\VKKUMUWRVQP

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

实验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

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

实验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

非线性电路中的微分方程解法

非线性电路理论及应用

周波 电路研-11 2011307080114

非线性电路中的微分方程解法

微分方程数值解法初值： 初值： 所谓初值问题， 所谓初值问题，是函数及其必要的导数在积分的起始 点为已知的一类问题，一般形式为： 点为已知的一类问题，一般形式为：

y

(n)

= f ( x)或y

(n)

= f ( x, y′, y′′,L, y

( n 1)

)

我们先介绍 y′(x) = f (x, y(x)) 简单的一阶问题： 简单的一阶问题：

a≤x≤b

y(a) =y0

(8 1 )

只要f ( x, y )满足(李卜希兹)( Lipschitz条件) : f ( x, y ) f ( x, y ) ≤ L y y

非线性电路中的微分方程解法

由常微分方程的理论可知：上述问题的解唯一存在。 的理论可知：上述问题的解唯一存在。 常微分方程求解求什么？应求一满足初值问题(8—1) 求解求什么？应求一满足初值问题( 的解函数y 如对下列微分方程： 的解函数y = y(x) ,如对下列微分方程：

第八章 序

y ′ = xy 2 dy dy 2 x2 = xy 2 2 = xdx = +c dx y 2

关于数值分析的

常微分方程的数值解法

一、题目 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方法等；难点是理

good

第九节 二阶常系数非齐次线性微分方程

x x

f(x) [P(x)cos x Q(x)sin x]ef(x) P(x)emmm教学目的：掌握自由项为和的二

阶常系数非齐次线性微分方程特解的方法

教学重点：二阶常系数非齐次线性微分方程求特解的待定系数法 教学难点：二阶常系数非齐次线性微分方程求特解的待定系数法 教学内容：

二阶常系数非齐次线性微分方程的形式为：

y py qy f(x)

根据二阶线性微分方程解的结构，要求解二阶常系数非齐次线性微分方程，只需先求得

对应齐次线性微分方程的通解和该非齐次线性微分方程的一个特解即可。而齐次线性微分方程的通解已在上一目得到解决，因此本节将解决非齐次线性微分方程的特解问题。为此，针对自由项的特点，采用如下待定系数法：

根据二阶非齐次线性微分方程解的结构，要求二阶常系数非齐次线性微分方程的通解，

**yy Y就是非齐次方程的通Y只需先求得非齐次方程的特解和对应齐次方程的通解，则

解。而用待定系数法求二阶常系数非齐次线性微分方程y py qy f(x)的特解分两种

情形讨论：

一、f(x) e xPm(x)型

这里 是常数,Pm(x)是m次多项式.

由于指数函数与多项式之积的导数仍是同类型的函数,而现在微分方程右端

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

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

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

§1、 引

言

主要内容

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

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

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

此种关系的数学表达就为

微分方程

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

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

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

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

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