一阶线性微分方程求解

浙江省精品课程--高等数学AⅠ教案（同济六版）2013----------宁波工程学院

补讲2 常数变易法、可降阶方程

1、主要教学目标

1、一阶线性微分方程的标准形式及其解法；

2、三种可降阶微分方程的解法；

2、重点内容

1、一阶线性微分方程的解法及解的结构； 2、常数变易法；

3、三种可降阶微分方程的解法。 3、难点分析

1、用变量代换将伯努利方程转化为线性方程并求解； 2、常数变易法、用变量代换法求解微分方程。 4、对教材的处理及其教学提示

微分方程求解重在掌握思想方法，积分运算不宜过难，淡化伯努利(Bernoulli)方程的标准形式及其解法

5、作业布置P315-1（1）； 2（1）；3； P323-1（1、5、7）；4

一、线性方程

?P(x)dx. 1、通解公式 y?Ce?2、非齐次线性方程的解法----常数变易法

实质: 未知函数的变量代换。新未知函数u(x)?原未知函数y(x),

?P(x)dx?P(x)dxP(x)dx?u(x)[?P(x)]e?, 作变换y?u(x)e?，求导 y??u?(x)e??P(x)dxP(x)dx?Q(x),积分得 u(x)??Q(x)e?将y和y?代入原方程得u?(x)e?dx?C,

3、

一阶线性微分方程教学中的一点体会

摘要：通过对一个一阶线性微分方程组的求解，既让学生能够掌握简单的一阶线性微分方程组求解方法，又可以让学生较好地体会到《线性代数》课程的重要性。

关键词：一阶线性微分方程组；特征值；特征向量；线性变换 中图分类号：g642.1 文献标志码：a 文章编号：1674-9324（2013）19-0168-01

本学期，由于课程设置的调整，部分新生第一学期就开始学习《线性代数》这门课程了。在和他们交谈的过程中，部分学生反映这门课程没有什么作用，内容上主要是一些具体的运算，比如：计算行列式的值，计算矩阵的和、积、逆，求方阵的特征值、特征向量以及将方阵对角化等等。好像与实际应用一点关联也没有。 事实上，《线性代数》是一门非常重要的课程，其在很多专业课程中都有广泛的应用，只是学生没有认识到这一点。因此，我们有必要选择一些较为合适的例题，通过这些例题的讲解，既能够让学生易于接受，又可以让学生认识到《线性代数》课程的重要性，从而更好地激发他们的学习热情。为此，在一阶线性微分方程的教学当中，可以通过下面这个例题的讲解来达到我们的目的。 例：求一阶线性微分方程组

■=-x1（t）+2x2（t）■=3x1（t）-2x2（t） （1）

第二节

一阶线性偏微分方程的解法

一、线性偏微分方程 1、线性算子 算子是一种数学法则，把它作用在一个函数上便产 生了另外一个函数。 2 2 3 2 例如，L 3 及M 2 x 2 2 x x y y x y都是偏微分算子。 u 2u 3u 将其作用于函数u便有：L[u ] 3 x x y y2 2u u 2 M [u ] x 2015/10/13 x 2 y 2

u 2u 3u 于是偏微分方程 3 f ( x, y)便可简单 x x y y记为L[u ] f 或Lu f .

算子L若满足:L[au bv] aL[u] bL[v] 其中，a, b为常数；u, v为函数，则称L为线性算子。

2015/10/13

2.线性微分方程解的叠加原理

定理1:若u1 , u2 ,..., un是某个线性齐次微分方程L[u ]=0 的解，则 ci u i 也为此方程的解。(ci 为任意常数)i 1 n

定理2:若ui 是L[u ] fi (i 1, 2,...)的解，且 ciui收敛，i 1

则u ci ui 是L[

辽宁工程技术大学

数 学 建 模 课 程 成 绩 评 定 表

赵常新 魏文楷 潘洋 一阶常微分方程模型—人口模型与预测

数学建模

一阶常微分方程模型—人口模型与预测

一．摘要：

二．模型的背景问题描述

三．模型假设

四.分析与建立模型

下表列出了中国1982-1998年的人口统计数据，取1982年为起始年（t 0），

N0 101654万人，Nm 200000万人。

要求：（1）建立中国人口的指数增长模型，并用该模型进行预测，与实际人口数据进行比较。

（2）建立中国人口的Logistic模型，并用该模型进行预测，与实际人口数据进行比较。

（3）利用MATLAB图形，标出中国人口的实际统计数据，并画出两种模型的预测曲线。

赵常新 魏文楷 潘洋 一阶常微分方程模型—人口模型与预测

（4）利用MATLAB图形，画出两种预测模型的误差比较图，并分别标出其误差。

模型一：指数增长模型（马尔萨斯（Malthus）模型）

假设：人口净增长率r是一常数

符号：x(t) t时刻时的人口，可微函数x0 t 0时的人口 则r

x(t t) x(t)

x(t) t

dx

于是x（t）满足如下微分方程： dt rx

x(0) x0解为：x(t) x0ert 模型二：Lo

Some Properties of Solutions of Periodic Second Order

Linear Differential Equations

1. Introduction and main results

In this paper, we shall assume that the reader is familiar with the fundamental results and the stardard notations of the Nevanlinna's value distribution theory of meromorphic functions [12, 14,

(f)and (f)to denote respectively the order 16]. In addition, we will use the notation (f)，

of growth, the lower order of growth and the exponent of convergence of the zeros of a meromorphic function f， e(f)（[see 8]），the

二阶及高阶可降阶微分方程的求解与应用

摘要：根据自己的理解对几类可降阶的微分方程的解题技巧做了一

些总结归纳，并且将这些技巧在应用中得到体现。

关键词：微分方程 可降阶 应用

前言：通过参考大量论文后可以很清楚地发现，高阶微分方程的求

解没有统一的方法，并且几乎所有的论文在介绍高阶微分方程解题方法时均试图用二阶微分方程的求解来类推到高阶方程的求解中.归纳后即根据二阶齐次线性微分方程解的结构总结出求此方程通解的一种方法，再解出非齐次线性微分方程的一个特解就可以得到非齐次微分方程的通解。本篇文章主要是对一些比较特殊而实际应用很强的二阶常系数线性非齐次方程进行研究，从而推导出具有特殊性质的高阶微分方程的解法，用于解决在实际过程中会碰到的问题。

一、三类可降阶的二阶及高阶微分方程

可降阶方程作为一类具有特殊性质的二阶方程，具、有固定的解题模式，经过听取老师上课以及自己课后的整理，总结出三种可降阶类型。

1、形如：y''?f(x) 的方程

个人觉得这种类型方程是所有可降阶方程中最简单的一类，因此最先讨论。 方法：只需令

p?y'?，则p'?y''?积分可得p??f(x)dx?C1，

也就得到了y'??f(x)dx?C1，

二阶椭圆偏微分方程实例求解(附matlab代码)用的是五点差分法。

《微分方程数值解法》期中作业实验报告

二阶椭圆偏微分方程第一边值问题

姓名： 学号： 班级：

2013年11月19

日

二阶椭圆偏微分方程实例求解(附matlab代码)用的是五点差分法。

二阶椭圆偏微分方程第一边值问题

摘要

对于解二阶椭圆偏微分方程第一边值问题，课本上已经给出了相应的差分方程。而留给我的难题就是把差分方程组表示成系数矩阵的形式，以及对系数进行赋值。解决完这个问题之后，我在利用matlab解线性方程组时，又出现“out of memory”的问题。因为99*99阶的矩阵太大，超出了分配给matlab的使用内存。退而求其次，当n=10，h=1/10或n=70，h=1/70时，我都得出了很好的计算结果。然而在解线性方程组时，无论是LU分解法或高斯消去法，还是gauseidel迭代法，都能达到很高的精度。

关键字：二阶椭圆偏微分方程差分方程out of memory LU分解高斯消去法gauseidel迭代法

一、题目重述

解微分方程：

(eyux(x,y))x (exuy(x,y))y (x y)ux(x,y) (x y)uy(x,y) u(x,y) ye xe e y x 1 e

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

《数学实验》报告

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

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

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

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

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

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

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

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

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

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

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