第四节 一阶常系数线性差分方程
一阶线性常系数差分方程及其应用
r=[.09;.09;-.1;-.1;-1.9;-1.9;-2.09;-2.09]; % 增长率
x=[15;-15;85;-85;85;-85;15;-15];
% 初始值
for n=1:20
x(:,n+1)=(1+r).*x(:,n);
% 迭代计算
end s{1}='单调增趋于正无穷大,r>0,x_0>0';
3.2.4 按揭贷款
1. 问题提出
购买商品房,首付至少两成,余款做按揭贷款, 如何设计合适的按揭计划.
2. 问题分析
个人住房按揭贷款通常有两种分期还本付息方 式,一种是等额本息还款法,每月还款计算公式为:
每月还款额=贷款本金×月利率× (1+月利率)还款月数/[(1+月利率)还款月数-1]
3.2.4 按揭贷款
解答(续) 结论 (1)在中等和较差的自然环境下,由于 1 r 0 ,且 x0 0 ,所以 xk 单调衰减趋于 0,即沙 丘鹤将濒于灭绝;在 1 r 0 范围内,r 的绝对值越 大, xk 单调衰减得越快. (2)在较好的自然环境下,由于 r 0 ,且 x0 0 , 所以 xk 单调增趋于无穷大,即沙丘鹤数量将无限增长.
100
100
0
0
-100 0
5 10 15 20
单 调 增 趋 于 正 无 穷 大 ,r>0,x0>0
-100 0
5 10 15 20
单 调 减 趋 于 负 无 穷 大 ,r>0,x0<0
100
100
0
0
-100 0
5 10 15 20
单 调 减 趋 于 0,-1<r<0,x0>0
一阶常系数线性差分方程
目录
• 引言 • 差分方程的基本理论 • 一阶常系数线性差分方程的求解方法 • 一阶常系数线性差分方程的应用 • 一阶常系数线性差分方程的数值解法 • 一阶常系数线性差分方程的变种与扩展
01 引言
差分方程的概念
差分方程是描述离散 时间系统动态行为的 数学模型。
差分方程在经济学、 物理学、生物学等领 域有广泛应用。
分析经济周期
通过差分方程模型分析经济变量的 周期性变化,为政策制定提供参考。
评估政策效果
模拟不同政策对经济系统的影响, 评估政策的实施效果。
在信号处理中的应用
滤波处理
利用一阶常系数线性差分 方程构建数字滤波器,对 信号进行滤波处理,去除 噪声和干扰。
信号预测
基于差分方程模型对信号 未来走势进行预测,实现 信号的实时跟踪和监控。
与常系数线性差分方程不同,变 系数线性差分方程的系数可以随 时间变化,这使得方程的求解更
加复杂。
求解方法
变系数线性差分方程通常无法通 过简单的代数方法求解,而需要 使用迭代法、变换法或数值方法
等更复杂的求解方法。
应用领域
变系数线性差分方程在经济学、 金融学、信号处理等领域有广泛 应用,如描述股票价格、利率、
05 一阶常系数线性差分方程 的数值解法
欧拉法
基本思想
利用泰勒级数展开式,忽略高阶项, 得到差分方程的近似解。
迭代公式
通过给定的初始值,利用迭代公式逐 步求解差分方程的解。
误差分析
欧拉法是一种显式方法,其局部截断 误差与步长成正比,全局误差随步长 减小而减小。
稳定性分析
对于某些问题,欧拉法可能不稳定, 需要采用其他方法。
01
第一、二节差分方程的基本概念 一阶常系数线性差分方程
二阶线性常系数非齐次差分方程
2 yt + 3 − 3 yt + 2 + 4 yt +1 − 5 yt = 0
t t t t
三阶线性齐次差分方程
五.线性差分方程解的基本定理 线性差分方程解的基本定理 定理10.1 定理 如果 y1 ( t ), y2 ( t ),L , ym ( t ) 是齐次线性差分方程 的 m 个解 则它们的线性组合 个解,则它们的线性组合
2 2
解 ∆yt = f ( t + 1) − f ( t )
= [( t + 1) 2 + 2( t + 1)] − ( t 2 + 2t )
= 2t + 3
∆ yt = f ( t + 2) − 2 f ( t + 1) + f ( t )
2
= [( t + 2) + 2( t + 2)] − 2[( t + 1) + 2( t + 1)]
F ( t , y t , ∆y t , ∆2 y t , ∆3 y t , L , ∆n y t ) = 0
定义10.2 定义
含有自变量 t 和两个或两个以上
的函数值 yt , yt +1 ,L , yt + n的方程 称为差分方程 的方程,称为差分方程 称为差分方程. 出现在差分方程中的未知函数下标的最大差, 出现在差分方程中的未知函数下标的最大差 称为差分方程的阶. 称为差分方程的阶
F ( t , yt , yt +1 , yt + 2 ,L , yt + n ) = 0
注 两个定义不完全等价 例如
∆ y t + ∆y t = 0
差分方程方法
第四章 差分方程方法在实际中,许多问题所研究的变量都是离散的形式,所建立的数学模型也是离散的,譬如,像政治、经济和社会等领域中的实际问题。
有些时候,即使所建立的数学模型是连续形式,例如像常见的微分方程模型、积分方程模型等等,但是,往往都需要用计算机求数值解。
这就需要将连续变量在一定条件下进行离散化,从而将连续型模型转化为离散型模型,因此,最后都归结为求解离散形式的差分方程解的问题。
关于差分方程理论和求解方法在数学建模和解决实际问题的过程中起着重要作用。
下面就不同类型的差分方程进行讨论。
所谓的差分方程是指:对于一个数列{}n x ,把数列中的前1+n 项()n i x i ,2,1,0=关联起来所得到的方程。
4.1常系数线性差分方程4.1.1 常系数线性齐次差分方程 常系数线性齐次差分方程的一般形式为02211=+⋯+++---k n k n n n x a x a x a x (4.1)其中k 为差分方程的阶数,()k i a i ,,2,1 =为差分方程的系数,且()n k a k ≤≠0。
对应的代数方程02211=++++--k k k k a a a λλλ (4.2) 称为差分方程的(4.1)的特征方程,其特征方程的根称为特征根。
常系数线性齐次差分方程的解主要是由相应的特征根的不同情况有不同的形式。
下面分别就特征根为单根、重根和复根的情况给出差分方程解的形式。
1. 特征根为单根设差分方程(4.1)有k 个单特征根 k λλλλ,,,,321 ,则差分方程(4.1)的通解为nk k n n n c c c x λλλ+++= 2211,其中k c c c ,,,21 为任意常数,且当给定初始条件()0 i i x λ= ()k i ,,2,1 = (4.3)时,可以唯一确定一个特解。
2. 特征根为重根设差分方程(4.1)有l 个相异的特征根()k l l ≤≤1,,,,321λλλλ 重数分别为l m m m ,,,21 且k m li i =∑=1则差分方程(4.1)的通解为n l i m i li n i m i i n i m i i n n c n c n c x lλλλ112112111121-=-=-=∑∑∑+++=同样的,由给定的初始条件(4.3)可以唯一确定一个特解。
一阶线性差分方程的解法分析
高中数学新课标选修内容“一阶线性差分方程”的解法分析江西省高中数学课程标准研究组舒昌勇(341200)在高中数学新课标选修系列4的“数列与差分”专题中,一阶常系数线性差分方程x n+1=kx n+b (1)是讨论的重点,其一般形式为x n+1=kx n+f(n) (2)其中k为已知的非零常数,f(n)为n的已知函数.当f(n)≠0时,方程(2)称为非齐次的,f(n)=0时,方程x n+1=kx n(3)称为齐次的,并称(3)为(2)相应的齐次方程.方程(1)是方程(2)当f(n)为常数的情况,是方程(2)能用待定系数法求特解时所具有的几种特殊形式里最简单的一种.我们来讨论方程(1)和(3)通解的求法.1 求一阶齐次差分方程x n+1=kx n的通解用迭代法,给定初始值为x0,则一阶齐次差分方程x n+1=kx n的通解为x1 = kx0,x2=kx1=k2x0,x3=kx2=k3x0,…,一般地,有x n= kx0-1= k(k n-1x0)= k n x0,n = 1,2,…,由于x0表示初始值,可任意给定,所以可视其为任意常数,不妨用c来表示.又根据差分方程通解的定义:如果差分方程的解中含有与方程的阶数相同个数的相互独立的任意常数,则为其通解,故一阶线性齐次方程x n+1=kx n的通解可表为x n=k n c(c为任意常数).对于每一个任意给定的初始值x0,都能得到方程相应于该初始值的一个特解.而求特解只要将给定的初始值x0代入通解求出待定常数c即可.2 求一阶非齐次差分方程x n+1=kx n+b的通解2.1探索一阶非齐次差分方程x n+1=kx n+b通解的结构设数列﹛y n﹜,﹛z n﹜为方程(3)的任意两个解,则y n+1=k y n +b (4)z n+1= k z n +b (5)(4)-(5) 得y n +1-z n +1=k(y n- z n )这意味着一阶非齐次线性差分方程任意两个解的差为相应齐次差分方程的解.从而,若a n为非齐次方程(3)的任意一个解,b n为非齐次方程(3)的一个特解,则a n-b n就为相应齐次方程的一个解.为了探索一阶非齐次差分方程通解的结构,我们对它的任意一个解a n 作适当变形:a n=a n+b n- b n= b n +( a n - b n)这表明,一阶非齐次差分方程的任意一个解可表示为它的一个特解与相应齐次方程一个解的和的形式.从而非齐次方程的通解等于其一个特解加上相应齐次方程的通解.2.2 求一阶非齐次差分方程(3)的通解①用迭代法,设给定的初始值为x0,依次将n=0,1,2,…代入(3),有x1=kx0+bx2=kx1+b=k(kx0+b)+b =k2x0+b(1+k)x3=kx2+b= k[k2x0+b(1+k)]+b= k3x0+b(1+k+k2)……x n=k n x0+b(1+k+k2+…+k n-1)ⅰ)当k ≠1时, 1+k+k 2+…+k n-1 =k k n --11 此时x n =k n x 0+k k b n --1)1(=k n (x 0-k b-1)+k b -1由于x 0表示初始值,可任意给定,故可设其为任意常数,从而x 0-k b-1 也为任意常数.令x 0-kb-1=c ,则(3)的通解可表为 x n =k n c+k b -1 (c 为任意常数)ⅱ)当k=1时,1+k+k 2+…+k n-1=n此时x n =x 0+nb由于x 0可任意给定,即其可为任意常数,故(3)的通解可写为x n =c+nb (c 为任意常数)②待定系数法与求解常微分方程类似,待定系数法也是求非齐次线性差分方程一个特解的一种较为简便、常用的方法.其基本思想是:根据方程的非齐次项f(n)的特点,用与f(n)形式相同但系数为待定的函数,作为方程的特解(称为试解函数),然后将该试解函数代入方程,以确定试解函数(特解)中的待定系数,从而求出方程的一个特解.ⅰ)当k ≠1时,设方程(3)有一特解x n =A ,其中A 为待定常数,将其代入(3),有 A=kA+b , A=k b -1 , 即x n =kb -1 知此时方程(3)的通解为x n = k n c+kb -1 (c 为任意常数) ⅱ)当k=1时,方程(3)为x n+1=x n +b ,知其解数列的一阶差分为常数,可设其有形如x n =An 的特解,代入(3),有A(n+1)=An+b , 得A=b , 即x n =bn知此时方程(3)的通解为x n = k n c+bn= c+bn (c 为任意常数)例1 求差分方程2y t+1+5y t =0的通解,并求满足y 0=2的特解.解 将原方程改写成y t+1=(-25)y t , 故其通解为y t =(-25)t c , c 为任意常数.用y 0=2代入通解:2=(-25)0c , 得 c = 2 .满足初值y 0=2的特解为y t =2(-25)t. 例2 求下列差分方程的通解(1)x n+1=x n +4(2)x n+1+x n =4解(1)方程中有k=1,b=4 .其通解为x n =c+4n ,(c 为任意常数).(2)原方程可化为 x n+1= -x n +4 ,方程中k=-1,b=4 ,其通解为 x n = (-1)n c+)1(14--= (-1)n c+2 ,(c 为任意常数).例3 某学术报告厅的座位是这样的安排的:每一排比前一排多2个座位.已知第一排有30个座位,(1)若用y n 表示第n 排的座位数,试写出用y n 表示y n+1的公式. (2)第10排的座位是多少个?(3)若用S n 表示前n 排的座位数,试写出用S n 表示S n+1的公式. (4)若该报告厅共有20排,那么一共有多少个座位?解 (1)y n+1= y n +2 n =1,2,…(2)解上述差分方程,其中k=1,b=2 ,通解为 y n =2n+c ,c 为任意常数 .由已知y 1=30,代入,得c = 28 .特解为y n =2n+28 , y 10=2×10+28=48(个) .(3)S n+1=S n +y n+1=S n +[2(n+1)+28]可得表达式为 S n+1=S n +2n+30 , n=1,2,…(4)先解上述差分方程,由S n+1-S n =2n+30 ,即△S n =2n+30,知S n 的表达式为n 的二次函数,设S n =An 2+Bn+C ,则△S n =A (n+1)2+B (n+1)+C -An 2-Bn -C=2A n+ A+B = 2n+30 .可得 A=1, B=29 . 又由初始条件 y 1= 30= S 1,有30 =A+B+C ,故C=0 .因此本问题的特解S n = n 2+29n , n =1,2,…S 20= 202+29×20=980(个).注意:在本例小题(1)中每排座位数的表达式y n+1=y n +2 y n+1-y n =2,与小题(2)中前n+1排座位数表达式S n+1=S n +2n+30即S n+1-S n =2n+30都属一阶非齐次线性差分方程x n+1=kx n +f(n)类型,但前者属f(n)为常数的情况,而后者属f(n) 为n 的一次函数的情况,利用差分有关知识,知S n 的表达式是关于n 的二次函数.参考文献[1] 教育部.普通高中数学课程标准(实验)[S].北京:人民教育出版社,2003.83-85.[2] 严士健,张奠宙,王尚志. 普通高中数学课程标准(实验)解读[M].南京:江苏教育出版社,2004.218-228.[3] 张银生,安建业.微积分[M].北京:中国人民大学出版社,2004.431,448-460.[4] 黄立宏,戴斌祥.大学数学(一)[M]. 北京:高等教育出版社,2002.380-389 .(本文刊于中学数学教学(合肥),2006,6.)。
一阶常系数线性差分方程.ppt
(10 13) (10 14)
方程 (10 14) 变形后改写为
yn1 ayn , n 0, 1, 2, 这是等比数列所满足的关系式, 由等比数列通项公式
可以得到 yn (a)n y0 , n 0, 1, 2,
从而得到方程 (10 14) 的通解
y C(a)n , n 0, 1, 2,
(10 15)
其中C 为任意常数.
二、非齐次方程的特解与通解
1. f (n) Pm (n), Pm(n) 为m 次多项式, 则方程(10 13) 为
yn1 ayn Pm (n)
(10 16)
例1 求差分方程 yn1 yn n 3 的通解. 解 因a 1, 对应齐次方程的通解为
y C 1n C
设 y(n) a0n2 a1n, 代入原方程, 有
a0(n 1)2 a1(n 1) a0n2 a1n n 3
比较系数得
a0
1, 2
其中a0,a1,,am为待定系数, 代入方程后, 比较同幂次系数, 可以解代数方程确定待定系数. 若 a 1, 要使方程恒等, 则应设
y(n) nQm (n) a0nm1 a1nm am1n2 amn. 代入方程, 比较同幂次系数, 可以解出式中的待定系数 a0 , a1 , ,am .
y(n) (a)n1 f (0) (a)n2 f (1) (a) f (n 2)
f (n 1)
n1
(a)i f (n i 1)
3 4
2n
由
y0
4,
差分方程2
致价格下跌,农民觉得养猪赔钱,于是转而经营其
它农副产品 . 过一段时间猪肉上市量减少,供不应求 导致价格上涨,原来的饲养户觉得有利可图,又重 操旧业,这样下一个时期会重新出现供大于求 , 价格 下跌的局面. 在没有外界干预的条件下,这种现象将
一直循环下,在完全自由竞争的市场体系中,这种
现象是永远不可避免的 . 由于商品的价格主要由需求
yt a t y0 ( t 0,1,2,).
yt 1 ayt 0
(4)
(2) 特征方程法求解:设
Y t ( 0)
是方程 (4) 的解,代入(4),得
t 1 a t 0 ( 0),
化简得:
a 0,
即
a.
t 1 a t 0
利率,按年复利计息,则 S t 与 r 有如下关系式:
S t 1 S t rS t (1 r ) S t ,
t 0,1,2,,
这是关于 S t 的一个一阶常系数齐次线性差分方程,
其通解为
S t (1 r ) t S 0 ,
t 0,1,2,,
其中 S 0为初始存款总额.
二、 一阶常系数线性差分方程 一阶常系数线性差分方程的一般形式为
yt 1 ayt f ( t ),
其中 a 0 为常数,f ( t ) 为已知函数. 当 f (t ) 0 时,称方程
yt 1 ayt 0
(3)
(a 0)
(4)
为一阶常系数齐次线性差分方程. 若 f (t ) 0, 则 (3) 称为一阶常系数非齐次线性 差分方程.
2
a r 12n I (1 ) I 2 12 (1 r 12n ) 1 12 x, r 12
10.2 一阶常系数线性差分方程
yt 1 ayt 0
(10.2.2)
微积分 第10章 差分方程初步
10.2一阶常系数线性差分方程
方程(10.2.4)称为差分方程(10.2.2) 的特征方程.
a
为特征方程的根,简称特征根. 求出特征根即得方程(10.2.2)的
t y ( a ) 一个特解 t ,再由§10.1的定理2可得(10.2.2)的通解为
B1 cos t B2 sin t
比较上式两边同类型三角函数的系数,得
b1 (a cos ) b2 sin B1 b1 sin b2 (a cos ) B2
因为
(10.2.10)
k
,故方程组(10.2.10)的系数行列式
微积分 第10章 差分方程初步
代入原方程,有
b0 (t 1)3t 1 3b0t 3t 3t
1 1 t * t 1 b y t 3 t 3 解得 0 ,故即原方程的一个特解为 t 3 3
从而原方程通解为
yt C(3)t t 3t 1
( C 为任意常数).
微积分 第10章 差分方程初步
10.2一阶常系数线性差分方程
yt 1 ayt 0
(10.2.2)
称方程(10.2.2)为方程(10.2..1)所对应的一阶常系数齐次线性 差分方程.
微积分 第10章 差分方程初步
10.2一阶常系数线性差分方程
由§10.1的定理3可知,若求得(10.2.1)所对应齐次方程(10.2.2)
* y y (t ) ,则方程 的通解 yt y (t ) 及方程(10.2.1)的一个特解 t
t
微积分 第10章 差分方程初步
差分方程基本知识
3. 常系数线性差分方程及解的性质
定义4 形如
ytn a1 ytn1 an1 yt1 an yt f ( x)
(1)
的差分方程称为n 阶常系数线性差分方程,其中
a1 , a2 , , an 为常数,且 an 0, f (t )为已知函数.
当 f (t) 0时,差分方程(1)称为齐次的,
例如,
yt2 2 yt1 yt 3t
是一个二阶差分方程, 可以化为
yt 2 yt1 yt2 3t2.
如果将原方程的左边写为
( yt2 yt1 ) ( yt1 yt ) yt1 yt 2 yt ,
则原方程还可化为 2 yt 3t.
若 f (t) 0, 则 (3) 称为一阶常系数非齐次线性
差分方程.
1. 常系数齐次线性差分方程的通解 对于一阶常系数齐次线性差分方程
yt1 ayt 0
(4)
通常有如下两种解法.
(1) 迭代法求解: 设 y0 已知,则
yn ayn1 a(ayn2 ) a2 yn2 an1 y1 an y0 ,
由于 a 1 , b 5 , a b,
22
故可设其特解为: yt* kbt .
代入方程,解得:k c 1 ,
ba 2
故原差分方程通解为:
yt
Y
yt*
A
1 t 2
1 2
5 t 2
.
(三) f (t) ctn (c为常数), 则差分方程为
2
于是原方程的通解为
其中C为任意常数.
yt
C(
1)t , 2
差分方程简介
它的通解是 y x Cx A ( A 是任何实常数). ( 3) y x Pn1 ( x ) ( n 1次多项式) 通解 y x Pn ( x ) ( n次多项式 )
4 n y x 0
通解 y x是n 1次多项式.
二、一阶常系数线性差分方程
形如: y x 1 ay x f ( x ) 齐次方程: y x 1 ay x 0
y x ( x 2 ) ( x 1)2 x 2 2 x 1
2 y x 2 ( x 2 ) (2 x 1) 2( x 1) 1 ( 2 x 1) 2
3 y x ( 2 y x ) ( 2) 2 2 0
x n x( x 1)( x 2)( x n 1) , x 0 1
例2 设 求 x n
解: x n ( x 1)n x n
( x 1) x( x 1)( x 1 n 1) x( x 1)( x n 1) [( x 1) ( x n 1)]x( x 1)( x n 2) nx n1
2. 差分方程 有某种商品 t 时期的供给量St与需求 一个例子: 量Dt都是这一时期价格Pt 的线性函数:
St a bPt (a , b 0) , Dt c dPt (c, d 0)
设 t 时期的价格Pt由 t –1时期的价格 Pt 1与供给量 及需求量之差 St 1 Dt 1 按如下关系确定.
Pt Pt 1 ( St 1 Dt 1 )
( 为常数),
即
Pt [1 (b d )]Pt 1 (a c )
这样的方程就是差分方程.
一阶常系数线性差分方程
微积分Calculus一阶常系数线性差分方程一一阶常系数线性差分方程概念1一般形式:1()x x y py f x +−=其中为不等于零的常数,为已知函数。
p ()f x ()f x 若不恒等于零,称以上方程为一阶常系数非齐次线性差分方程。
()f x 若恒等于零,称以上方程为一阶常系数齐次线性差分方程。
齐次线性差分方程的解法1yx =pyx−1=p ∙py x−2=p ∙p ∙py x−3=⋯=p x y 010x x y py +−=一阶齐次线性差分方程:将上述方程变形为:则有:记得一阶齐次线性差分方程的通解:0C y =xx y Cp = (为任意常数)C 二一阶常系数线性差分方程的解法y x+1=py x求差分方程130x x y y ++=的通解。
因为,将其代入通解公式得:3p =−(3)x x y C =− (为任意常数)C 13x xy y +=−将原方程变形为:例解一阶非齐次线性差分方程:1()x x y py f x +−=下面介绍对的三种特殊形式求非齐次差分方程特解的方法。
()f x 非齐次线性差分方程的解法2(1)(为常数,)()f x k =k 0k ≠差分方程变为:1x xy py k +−= 设其特解形式为:s x y Ax *=(其中为待定常数),A s1,p ≠①取即:0s =x y A*=1,p =②取即:1s =x y Ax*=x y A *=将代入差分方程求得A将代入差分方程求得Ax y Ax *=21716x x y y +++=求差分方程的通解.对应齐次差分方程:的通解为:217x x y y +++=0(7)xx y C =− (为任意常数)C p =−7≠1,设特解为y x ∗=A代入原方程得:2A =故原差分方程通解为:2(7)x x y C =+−(为任意常数)C 例解(2)(其中为常数,且)()xf x ka =k a ,0a >0a ≠非齐次差分方程变为:1x x x y py ka +−= 设特解形式为:x sx y Aa x*=①时,取即p a ≠0s =x x y Aa *=②p a =1s =x x y Axa *=时,取即求差分方程的通解11242x x x y y ++−=原方程化简为122xx x y y +−=对应齐次差分方程通解为2xx y C = (为任意常数)C 2p a ==由于,所以原方程得特解形式为:2xx y Ax =代入原方程得:1(1)2222x x xA x Ax ++−=12A =例解原方程特解为:11222x x x y x x *−==所以原方程通解为:12(2)x x x y x C −=+(为任意常数)C。
一阶线性差分方程组的公式解法
一阶线性差分方程组的公式解法
一阶线性差分方程组的公式可以用来描述这样的方程组:\begin{equation} \frac{d x(t)}{d t} = A x(t) + b, \end{equation}其中A是n×n的实矩阵,b是n维实列向量,x(t)是n维实列向量,表示在t的时刻的状态。
其中A矩阵中的每一个元素都代表了系统中变量之间的关系,b向量中的每一个元素则表示了外部影响系统的因素。
一阶线性差分方程组的解法有两种:
1. 直接解法:给定A和b,可以通过求解线性方程组来求解x(t)。
2. 通过求解积分方程的解法:可以将一阶线性差分方程组转化为积分方程,通过求解积分方程来求解x(t)。
一阶线性差分方程组在实际应用中有着广泛的应用,它可以用来描述物理系统,也可以用来描述计算机系统中的控制,通信,信息处理等系统。
它具有解析解和数值解两种解法,并且解法相对简单,使用起来也比较方便。
因此,一阶线性差分方程组在工程和科学中具有重要的意义,广泛地应用于控制理论,系统建模,动力学,信号处理等领域。
★差分方程
当 为常数时, yx = x和它的各阶差商有倍数关系, 所以可设 yx = x为方程(11)的解.
代如方程(11)得
x+2 + ax+1 + bx = 0,
2 + a + b = 0,
(12)
方程(12)称为齐次差分方程(11)的特征方程.
由特征方程的根的情况可得齐次方程的通解:
特征方程的解
两个不相等的实根 1, 2 两个相等实根 1 = 2
定义3 含有未知函数几个时期值的符号的方程, 称 为差分方程.
其一般形式为
G(x, yx, yx+1, , yx+n) = 0.
(2)
定义3中要求 x, yx, yx+1, , yx+n不少于两个.
例如, yx+2 + yx+1 = 0为差分方程, yx = x不是差分方 程.
差分方程式(2)中, 未知函数下标的最大差数为 n, 则 称差分方程为n 阶差分方程.
例6 求差分方程 yx+1 yx = x +1 的通解.
解 对应的齐次方程 yx+1 yx = 0的通解为 y*x C.
这里 a = 1, 设 yx x(B0 B1x), 代入差分方程, 得
(x+1)[B0+B1(x+1)] x(B0+B1x) = x +1. 整理, 得
2B1 x + B0 + B1 = x +1.
解 对应的齐次方程的特征方程为
2 3 + 2 = 0.
方程的根为
1 = 1, 2 = 2,
因为 q = 2 =2, 设特解为 y x Bx2x ,
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
lim Pt = P ;
19
δ t Pt = ( P0 P )( ) + P β
δ (3) 若初始价格不等于平衡价格,且 = 1 ,这时现行价 若初始价格不等于平衡价格, β 称为临界情形; 格 Pt 交替取两个值 P0 和 2 P P0 , 称为临界情形 ;
δ (4) 若初始价格不等于平衡价格,且 > 1 ,这时现行价 若初始价格不等于平衡价格, β 格 Pt 围绕平衡价格上下波动的振幅将随着 t 增大而无限
增大, 增大 , 从而产生价格的大波动 .
20
练习: 练习:
P394 习题九
21
谢谢大家! 谢谢大家!
22
�
为一阶常系数齐次线性差分方程, 一阶常系数齐次线性差分方程, 齐次线性差分方程
(2)
否则,称为一阶常系数非齐次线性差分方程. 否则,称为一阶常系数非齐次线性差分方程. 一阶常系数非齐次线性差分方程 (2)称为 对应的齐次线性差分方程. 称为(1)对应的齐次线性差分方程. 称为 对应的齐次线性差分方程
y t+ 1 y t = 2( At + A + B At B )2 t
= 2A2 ≠ t 2 ,
t t
不存在这样的特解. 不存在这样的特解.
10
t 的通解. 例4 求一阶常系数线性差分方程 y t + 1 2 y t = t 2 的通解.
解 设特解 y t = t ( At + B )2 , 代入方程得
14
t
π
π
例6 求差分方程 2 y t + 1 + 10 y t 5t = 0 的通解 .
5 解 首先把差分方程改写为标准形式 y t +1 + 5 y t = t . 2 a = 5 ≠ 1 , 设特解 y t = A t + B ,代入方程得
y t+ 1 + 5 y t = A( t + 1) + B + 5( At + B )
δ α +γ β Pt + δ Pt 1 = α + γ Pt + Pt 1 = ( t = 0, 1, ) β β 求出方程的通解是
δ α +γ Pt = C + β β +δ
t
17
δ α +γ Pt = C + β β +δ α +γ 令P = , 称 为平 衡 价 格(或均衡价格). β +δ
δ t 则由于 Pt 中包含因子 ( P0 P )( ) , 现行价格 Pt 将始终 β
δ (2) 若初始价格不等于平衡价格,且 < 1 , 这时现行价 若初始价格不等于平衡价格, β 格 Pt 围绕平衡价格上下波动的振幅将随着 t 增大而逐渐减
少 ,且
t → +∞
围绕平衡价格上下波动; 围绕平衡价格上下波动 ;
yt = C(a)t + yt
其中 y
(1)的一个特解, t 是
t = 0, 1, 2, ,C 是任意常数.
3
当 f(x)是多项式,指数函数,正弦函数,余弦函 ( )是多项式,指数函数,正弦函数, 数以及它们的和差或乘积时,一般可用待定系数法 待定系数法求 数以及它们的和差或乘积时,一般可用待定系数法求 (2)的一个特解. 的一个特解. 的一个特解 讨论三种情形: 讨论三种情形: 情形1 情形1 情形2 情形2 情形3 情形3
t
π
= ( B 2 A) cos
π
2
2
t ( A + 2 B ) sin
π
2
t ≡ 5 cos
π
2
t
A = 2, B = 1 ,
得特解为 y = 2 cos
t
t
π
2
t + sin
π
2
t,
通解为 y t = C 2 2 cos
π
2
t + sin
π
2
t , C为任意常数. 为任意常数. 为任意常数
表中 Q m (t ) 是待定系数的 m 次多项式.
8
求一阶常系数线性差分方程 y t + 1 y t = t 2 t 的通解. 通解. 例3
t 解 设特解 yt = ( At + B )2 ,代入方程得
y
t +1
y = ( 2 At + 2 A + 2 B At B )2 t
t
= ( At + 2 A + B )2 t = t 2 t , A = 1, B = 2 ,
第四节
1
一阶常系数线性差分方程标准形式为
yt +1 + ayt = f (t)
时有定义. 时有定义.
(1)
函数 f (t ) 当 t = 0, 1, 2, 其中 t = 0, 1, 2, ,常数 a ≠ 0 ,
如果当 t = 0, 1, 2, 时有 f ( t ) ≡ 0 , 则称方程
yt+1 +ayt = 0
Qm (t ) d t
t Qm (t ) d t
表中 Q m (t ) 是待定系数的 m 次多项式.
12
例5
求线性差分方程 yt + 1 2 y t = 5 cos
t
π
2
t 的通解. 的通解.
解 设特解 y = A cosபைடு நூலகம்
π
2
t + B sin
π
2
π
2
t , 代入方程得
2 2
y
t +1
2 y = A sin t + B cos t 2( A cos π t + B sin π t )
Q s , t = S ( Pt 1 ) ;
(2) 现 期 该 产品 的 销 售量 Q d , t 由 现 行 价 格 Pt 确 定 , 即
Q d , t = D( Pt ) ;
(3) 现期该产品的供需平衡,即 Q d , t = Q s , t . 现期该产品的供需平衡,
常见的供给函数S与需求函数 均为线性函数 常见的供给函数 与需求函数D均为线性函数,于 与需求函数 均为线性函数, 是得到方程组 Q s , t = γ + δPt 1 ( 3) ( 4) Q d , t = α β Pt ( 5) 16 Q d , t = Q s , t
利用初始价格 P0 确定常数 C = P0 P ,
t
可得现行价格 Pt 的公式
δ Pt = ( P0 P ) + P β
18
t
δ t Pt = ( P0 P )( ) + P β
由此可得以下简单结论: 由此可得以下简单结论: (1) 若初始价格 P0 等于平衡价格 P ,则现行价格 Pt 将始 终等于平衡价格;否则,若初始价格不等于平衡价格, 终等于平衡价格 ;否则,若初始价格不等于平衡价格,
11
, 当 f ( t ) = Pm ( t ) d t 时, 一般 可按下表设定非齐次 差分方程 y t + 1 + ay t = f ( t ) 的一个特解 y t :
f (t )
d 与系数 a 的关系
t
特解 y
t 的形式
Pm ( t ) d Pm ( t ) d
t
a+d ≠0 a+d =0
7
一般, 当 f (t ) 是多项式 Pm (t ) 时 ,可按下表设定 非齐次差分方程 y t + 1 + ay t = f ( t ) 的一个特解 y t :
f (t )
系数 a 的取值
特解 y
t 的形式
Pm (t ) Pm (t )
a ≠ 1 a = 1
Q m (t )
t Q m (t )
13
一般, 一般 , 当 f ( t ) = M cos ω t + N sin ω t , 其中 是常数, 可以设 M , N , ω 是常数,且 0 < ω < 2π , ω ≠ π ,可以设 特解为
y = A cos t + B sin t 2 2
其中 A, B 是两个待定常数 .
如果所给差分方程不是标准形式的, 如果所给差分方程不是标准形式的,必须首先把 它化为标准形式才能应用上面给出的通解公式和选取 它化为标准形式才能应用上面给出的通解公式和选取 标准形式 特解的有关结论. 特解的有关结论.
y
t +1
y = A( t + 1 ) + B ( A t + B ) = A ≠ 3 t 2 ,
t
没有这样的特解. 没有这样的特解.
6
例2 求一阶常系数线性差分方程 yt +1 yt = 3t 2 的通解. 的通解. 解 设特解 y t = t ( At + B ) = At 2 + B t , 代入方程得
t
y
t +1
y = 2[ A( t + 1) 2 + B ( t + 1) At 2 B t ] 2 t
t
t
t
= 2( 2 At + A + B )2 ≡ t2
1 1 A= ,B= , 4 4 1 y t = t ( t 1)2 t , 得特解为 4 t t2 t 为任意常数. 为任意常数 从而通解为 y t = (C + )2 , C为任意常数. 4 4
Q s , t = γ + δPt 1 Q d , t = α β Pt Q d , t = Q s , t
行价格 Pt , 并研究其变化规律 .