微分方程模型及差分模型
- 格式:ppt
- 大小:5.38 MB
- 文档页数:112
实验二: 微分方程与差分方程模型Matlab 求解一、实验目的[1] 掌握解析、数值解法,并学会用图形观察解的形态和进展解的定性分析;[2] 熟悉MATLAB 软件关于微分方程求解的各种命令; [3] 通过范例学习建立微分方程方面的数学模型以及求解全过程;[4] 熟悉离散 Logistic 模型的求解与混沌的产生过程。
二、实验原理1. 微分方程模型与MATLAB 求解 解析解用MATLAB 命令dsolve(‘eqn1’,’eqn2’, ...) 求常微分方程〔组〕的解析解。
其中‘eqni'表示第i 个微分方程,Dny 表示y 的n 阶导数,默认的自变量为t 。
〔1〕 微分方程例1 求解一阶微分方程 21y dxdy+= (1) 求通解 输入:dsolve('Dy=1+y^2') 输出: ans =tan(t+C1) 〔2〕求特解 输入:dsolve('Dy=1+y^2','y(0)=1','x') 指定初值为1,自变量为x输出: ans =tan(x+1/4*pi)例2 求解二阶微分方程 221()04(/2)2(/2)2/x y xy x y y y πππ'''++-=='=-原方程两边都除以2x ,得211(1)04y y y x x '''++-= 输入:dsolve('D2y+(1/x)*Dy+(1-1/4/x^2)*y=0','y(pi/2)=2,Dy(pi/2)=-2/pi','x') ans =- (exp(x*i)*(pi/2)^(1/2)*i)/x^(1/2) +(exp(x*i)*exp(-x*2*i)*(pi/2)^(3/2)*2*i)/(pi*x^(1/2))试试能不用用simplify 函数化简 输入: simplify(ans)ans =2^(1/2)*pi^(1/2)/x^(1/2)*sin(x) 〔2〕微分方程组例3 求解 d f /d x =3f +4g ; d g /d x =-4f +3g 。
第四章 微分方程与差分方程方法第一节 微分方程模型我们在数学分析中所研究地函数,是反映客观现实世界运动过程中量与量之间地一种关系,但我们在构造数学模型时,遇到地大量实际问题往往不能直接写出量与量之间地关系,却能比较容易地建立这些变量和它们地导数(或微分>间地关系式,这种联系着自变量、未知函数及其导数(或微分>地关系式称为微分方程.§4.1.1微分方程简介这一节,我们将介绍关于微分方程地一些基本概念. 一、微分方程地阶数首先我们具体地来看一个微分方程地例子.例4-1 物体冷却过程地数学模型将某物体放置于空气中,在时刻0=t ,测量得它地温度为C u 00150=,10分钟后测量得温度为C u 01100=.我们要求决定此物体地温度u 和时间t 地关系,并计算20分钟后物体地温度.这里我们假定空气地温度保持为C u 024=α.解:根据物理学中地牛顿冷却定律可知,热量总是从温度高地物体向温度低地物体传导。
一个物体地温度变化速度与这一物体地温度与其所在介质温度地差值成正比.设物体在时刻t 地温度为)(t u u =,则温度地变化速度可以用dtdu来表示.我们得到描述物体温度变化地微分方程)(αu u k dtdu--=(4.1.1> 其中0>k 是比例常数.方程(4.1.1>中含有未知函数u 及它地一阶导数dtdu,这样地方程,我们称为一阶微分方程.微分方程中出现地未知函数最高阶导数地阶数称为微分方程地阶数.方程)(33t f cy dt dyb dty d =++(4.1.2> 中未知函数最高阶导数地阶数是三阶,则方程(4.1.2>称为三阶微分方程. 二、常微分方程与偏微分方程如果在微分方程中,自变量地个数只有一个,我们称这种微分方程为常微分方程。
自变量地个数为两个或两个以上地微分方程称为偏微分方程.方程0222222=∂∂+∂∂+∂∂zTy T x T (4.1.3> 就是偏微分方程地例子,其中T 是未知函数,x 、y 、z 都是自变量.而方程(4.1.1>(4.1.2>都是常微分方程地例子.三、线性与非线性微分方程如果n 阶常微分方程0),,,,(=n n dxyd dx dy y x F (4.1.4>地左端为关于未知函数y 及其各阶导数地线性组合,则称该方程为线性微分方程,否则称为非线性方程.一般地n 阶线性微分方程具有形式)()()()(1111x f y x a dx dyx a dx y d x a dx y d n n n n n n =++++--- (4.1.5> 其中)1( )(),(n i x f x a i =是关于x 地已知函数.当()0f x =时,称(4.1.5>为n 阶齐次线性微分方程。
第八讲 差分方程模型一、差分方程介绍规定t 只取非负整数。
记为变量在t 点的取值,则称t y y t t t y y y −=Δ+1为的一阶向前差分,简称差分,称Δ为的二阶差分。
类似地,可以定义的阶差分。
t y t t t t t y t t y y y y y y +−=Δ−Δ=ΔΔ=+++12122)(t y t y n t ny Δ由及的差分给出的方程称为的差分方程,其中含的最高阶差分的阶数称为该差分方程的阶。
差分方程也可以写成不显含差分的形式。
例如,二阶差分方程也可改写成t y t 、t y t y t y 02=+Δ+Δt t t y y y 012=+−++t t t y y y 。
满足一差分方程的序列称为差分方程的解。
类似于微分方程情况,若解中含有的独立常数的个数等于差分方程的阶数时,称此解为该差分方程的通解。
若解中不含任意常数,则称此解为满足某些初值条件的特解。
t y 称如下形式的差分方程)(110t b y a y a y a t n t n t n =+++−++L (1) 为阶常系数线性差分方程,其中是常数,n n a a a ,,,10L 00≠a 。
其对应的齐次方程为0110=+++−++t n t n t n y a y a y a L (2)容易证明,若序列与均为(2)的解,则也是方程(2)的解,其中为任意常数。
若是方程(2)的解,是方程(1)的解,则也是方程(1)的解。
)1(t y )2(t y )2(2)1(1t tt y c y c y +=21,c c )1(t y )2(t y )2()1(t t t y y y +=方程(1)可用如下的代数方法求其通解: (I )先求解对应的特征方程(3)00110=+++−a a a n nL λλ(II )根据特征根的不同情况,求齐次方程(2)的通解。
(i )若特征方程(3)有n 个互不相同的实根n λλ,,1L ,则齐次方程(2)的通解为t n n t c c λλ++L 11 (为任意常数)n c c ,,1L (ii )若λ是特征方程(3)的重根,通解中对应于k λ的项为t k k tc c λ)(11−++L ,),,1(k i c i L =为任意常数。