线性常微分方程组求解 代码

Review 常系数齐次线性ODE的特征解法x(n)n

+ a1 x

( n 1)

λ + a1λ特征根 重数

n 1

+ L + an 1 x′ + an x = 0

+ L + an 1λ + an = 0.线性无关解 λt

λ (实) λ (实)

1kλt αt

e

e ,te , , t Lαt αt αt

λt

k 1 λt

e

α ± iβ

1k

e cos β t , e sin β t e cos β t , te cos β t ,L , t e cos β t , eα t sin β t , teα t sin β t ,L , t k 1eα t sin β tk 1 α t

α ± iβ

常系数非齐次线性ODE的待定系数法

x ( n ) + a1 x ( n 1) + L + an 1 x′ + an x = f (t ) f (t ) special solution x(t )

q (t )t k eλt , q real polynomial, p (t )e , λ ∈ deg(q ) ≤ deg( p ), p real polynomial, k = multiplicity of λ as an

第五章 线性微分方程组

[教学目标]

1. 理解线性微分方程组解的存在唯一性定理，掌握一阶齐（非齐）线性微分方程组解

的性质与结构，

2. 理解n 阶线性微分方程与一阶线性微分方程组的关系。 3. 掌握非齐次线性微分方程组的常数变易法，

4. 理解常系数齐线性微分方程组基解矩阵的概念，掌握求基解矩阵的方法。 5. 掌握常系数线性微分方程组的Laplce变换法。 [教学中难点]求解常系数非齐次线性微分方程组 [教学方法] 讲授，实践。 [教学时间] 16学时

[教学内容] n 阶线性微分方程与一阶线性微分方程组的关系，一阶线性微分方程组解的存在唯一性定理；齐（非齐）线性微分方程组解的性质与结构，求解非齐次线性微分方程组的常数变易法；常系数齐线性微分方程组的基解矩阵及求基解矩阵的方法；求常系数线性微分方程组的Laplce变换法。 [考核目标]

1.线性微分方程组解的性质与结构。 2.能够求解常系数线性微分方程组。

§5.1 存在唯一性定理

5.1.1记号和定义 考察形如

??a11(t)x1?a12(t)x2???a1n(t)xn?f1(t)?x1?x??a(t)x?a(t)x???a(t)x?f(t)?22112222nn2

Jordan标准型在求解线性微分方程组的应用

姓名：XX远 学号：20092426 班级：2009121

摘要：该应用在于矩阵不能对角化的时候利用Jordan标准型这种最简化的结果来做题。证明关于一般方阵（不能保证对角化）的某些命题，需要用到Jordan标准型，比如用Jordan标准型求解线性微分方程组。

关键字：Jordan标准型，线性微分方程组，特征值

矩阵内容，是大学学习中必须学习的知识点！其广泛的应用性，还有在处理数据上的优越性，矩阵是学习很多知识体系的支柱，在数据结构，自动控制原理，常微分计算等等上都是基础！

矩阵的对角化用处很大，因为对角化后，对矩阵加乘等等运算都可以简单很多，尤其在涉及特征值的方面！但是许多时候矩阵不能对角化。这时候相似变换的最好结果就是Jordan

标准型的形式，因为矩阵的Jordan标准型是最间的！

Jordan标准型的用处就在于矩阵不能对角化的时候利用Jordan标准型这种最简化的结果来做题。证明关于一般方阵（不能保证对角化）的某些命题，需要用到Jordan标准型，

比如用Jordan标准型求解线性微分方程组。

第三篇 第十章 常微分方程（组）求解

Matlab常微分方程（组）求解 一、 求微分方程的解

（一） 相关函数（命令）及简介

1， dsolve('equ1','equ2',…):Matlab求微分方程的解析解。

equ1,equ2,…为方程（或条件）。写方程（或条件）时用Dy表示y关于自变量的一阶导数，用D2y表示y关于自变量的二阶导数，依次类推。

2， simplify(s):对表达式s使用maple的化简规则进行化简。 例如： syms x

simplify(sin(x)^2+cos(x)^2) ans=1

3,[r,how]=simple(s):由于Matlab提供了多种化简规则，simple命令就是对表达式s用各种规则进行化简，然后用r返回最简形式，how返回形成这种形式所用的规则。 例如： syms x

[r,how]=simple(cos(x)^2-sin(x)^2) r=cos(2*x) how=combine

4,[T,Y]=solver(odefun,tspan,y0),求微分方程的数值解。 （1）其中的solver为命令

ode45,ode23,ode113,ode15s,ode23s,ode23t,ode23tb之一

《数学实验》报告

实验名称 常微分方程的求解 学 院 专业班级 姓 名 学 号

2013年5月

一、 【实验目的】

1. 学习在MATLAB中如何求解微分方程的方法；

2. 掌握基本的微分求解命令，学会结合学过的基础知识求解方程； 3. 熟练运用基本的解法即数值解法解微分方程； 4. 注意不同方法下求得微分方程的优缺点。

二、 【实验任务】

xsinxy?1. 求解微分方程为cosy。

''y2. 用数值方法求解下列微分方程，用不同颜色和线形将y和画在同一个

图形窗口里：

y?ty?y?1?2t初始时间：t0=0;终止时间：tf

三、 【实验程序】 1.

y=dsolve('Dy=x*sinx/cosy','x') 2.

定义的程序：

function xdot=exf(t,x)

xdot=[0 1;1 -t]*x+[0;1]*(1-2*t);

主程序：

2

'''

=?;初始条件：y|t?0?0.1 y

同济大学五版高等数学学习资料

第六章 常微分方程

一． 求解下列微分方程： 1. y' ex y

+ex=0.

解.

dydx=ex(e y 1), dye y 1

=exdx ln1 ey

=ex, 1 ey=cee xc

y=ln(1 ce

e x

).

2. dy dx

=(1 y2

)tanx

y(0)=2

解.

dy

1 y

2

=tanxdx

11+12lncy1 y= lncosx, y(0) = 2, 2lnc1+21 2=0, ln

1+y13+cos2x

3(1 y)=lncos2x, y=3 cos2x

二. 求解下列微分方程：

1. x x

1+ey 1 x

dx+ey

y dy=0 xey

x

1 解. dx y dy

=x

. 1+ey

令

x

y

=u,x=yu.(将y看成自变量) dxdy=u+ydudy

, 所以 u+ydudy=eu(u 1)

1+eu duueu euudy1+eu u= +eu

y=1+eu

c= 1

3

同济大学五版高等数学学习资料

u+eu 1dyd(u+eu)dy1+eu

ln= ln=ln= , = , ydu c yu+euyyu+eu

x

cc1u+euy

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

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

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

#$,&-./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

常 微 分 方 程

试卷（一至十） 试 卷（一）

一、填空题（3′×10＝30′）

1、以y1＝e2x，y2=exsinx，y3=excosx为特解的最低阶常系数齐次线性微分方程是 。

2、微分方程4x3y3dx+3x4y2dy=0的通积分是 。 3、柯西问题

dy?x，y（0）=1的解是 。 dx4、方程ydx-xdy=0的积分因子可取 。

5、证明初值问题的毕卡定理所构造的毕卡序列是 。 6、微分方程F(x，y，p)=0若有奇解y=? (x)，则y=? (x) 满足的P-判别式是 。 7、线性微分方程组

dY，Y2（x）…，Yn（x）?A(x)Y的解组Y1（x）

dx在某区间上线性无头的充分必要条件是 。 8、设A＝

1 0 1 0 0 -1 0 0 2 ，则矩阵指数函数exA＝ 。

9、方程y???y??y?0的通解是 。

10、由方程y????3ay???3ay??y?0的通解是 。 二、解下列各方程（7′×4＝28） 1、求方程

dyx?y?1?的通解： dxx?y?32、 (1+x2)y