差分方程及matlab求解
称为n阶差分方程的通解.在通解中给任意常数 C1,C2,…,Cn 以确定的值所得的解,称为n阶差分方程的 特解.
例如,函数yt=at+C(a为已知常数,C为任意常数)是差
分方程yt+1-yt=a的通解.而函数yt=at,yt=at-1,…均是这个 差分方程的特解. 由差分方程的通解来确定它的特解,需要给出确 定特解的定解条件. n阶差分方程F(t,yt,yt+1,…,yt+n)=0常见的定解条件为 初始条件.
方程的通解
1 t 1 1 t -1 2 t 1 t 1 t +1 yt = A( ) + ( ) ( 2 - 1) = A ( ) + 2 2 3 2 2 3 2 A = A - 为任意常数 . 3
2.待定系数法求特解
情形Ⅰ f(t)为常数. 方程变为yt+1+ayt=b, a,b均为非零常数. 试以 yt = (为待定常数)形式的特解代入方程得 +a =(1+a) =b. 当a≠-1时,可求得特解
三、 差分方程的解
定义4 如果将已知函数yt=j(t)代入方程F(t,yt,yt+1,…, yt+n)=0, 使 其 对 t=…,-2,-1,0,1,2,… 成 为 恒 等 式 , 则 称 yt=j(t)为方程的解.含有n个任意(独立)常数C1,C2,…,Cn 的解
yt=j(t,C1,C2,…,Cn)
解,那么,非齐次线性差分方程的通解为: y(t)=yA(t)+ y (t)
即
y(t)=A1y1(t)+A2y2(t)+…+Anyn(t)+ y(t), 这里A1,A2,…,An为n个任意(独立)常数.
第二节 一阶常系数线性差分方程
一阶常系数线性差分方程的一般形式为 yt+1+ayt=f(t) 和 yt+1+ayt=0, 其中f(t)为t的已知函数,a≠0为常数.分别称为一阶常 系数非齐次线性差分方程和其对应的齐次差分方程.
差分方程初步
第一节 差分方程的基本概念
一、 差分的概念 定义1 设函数yt=f(t)在t=…,-2,-1,0,1,2,…处有定义,对 应的函数值为…,y-2,y-1,y0,y1,y2,…,则函数yt=f(t)在时间 t的一阶差分定义为 Dyt=yt+1-yt=f(t+1)-f(t). 依此定义类推,有 Dyt+1=yt+2-yt+1=f(t+2)-f(t+1), Dyt+2=yt+3-yt+2=f(t+3)-f(t+2), ………………
一阶差分的性质
(1) 若yt=C(C为常数),则Dyt=0; (2) 对于任意常数k,D(kyt)=kDyt; (3) D(yt+zt)=Dyt+Dzt.
定义2 函数yt=f(t)在时刻t的二阶差分定义为一阶差分 的差分,即 D2yt= D (D yt)= D yt+1- D yt =(yt+2-yt+1)-(yt+1-yt)=yt+2-2yt+1+yt. 依此定义类推,有 D2yt+1= Dyt+2- Dyt+1=yt+3-2yt+2+yt+1, D2yt+2= Dyt+3- Dyt+2=yt+4-2yt+3+yt+2, ……………… 类推,计算两个相继的二阶差分之差,便得到三阶差分 D3yt= D2yt+1- D2yt=yt+3-3yt+2+3yt+1-yt, D3yt+1= D2yt+2- D2yt+1=yt+4-3yt+3+3yt+2-yt+1, ………………
一般地,k阶差分(k为正整数)定义为
Dk yt = D ( Dk -1 yt ) =D
k -1
yt +1 - D
k -1
yt ( k = 1,2,3, )
i = ( -1) i C k yt + k - i i =0
k
这里
k! C = i! ( k - i )!
i k
二、 差分方程 定义3 含有未知函数yt=f(t)以及yt的差分Dyt, D2yt,…的 函数方程,称为常差分方程(简称差分方程);出现在差 分方程中的差分的最高阶数,称为差分方程的阶. n阶差分方程的一般形式为 F(t,yt, Dyt,…, Dnyt)=0, 其中F是t,yt, Dyt,…, Dnyt的已知函数,且Dnyt一定要在方
3、在法国著名的 Lascaux 洞穴中保留着古 代人类遗留下来的壁画。从洞穴中取出的 木炭在 1950 年做过检测, 测得碳 14 的衰减 系数为每克每分钟 0.97 个, 已知碳 14 的半 衰期为 5568 年,试求这些壁画的年龄(精 确到百年) 。
4、牛顿发现在温差不太大的情况下,物 体冷却的速度与温差成正比。现设正常 体温为 36.5 ,法医在测量某受害者尸 体时测得体温约为 32 度, 一小时后再次 测量,测的体温约为 30.5 度,试推测该 受害者的受害时间。
四、 线性差分方程及其基本定理
形如
yt+n+a1(t)yt+n-1+a2(t)yt+n-2+…+an-1(t)yt+1+an(t)yt=f(t)
的差分方程,称为n阶非齐次线性差分方程.其中 a1(t),a2(t),…,an-1(t),an(t) 和 f(t) 都 是 t 的 已 知 函 数 , 且 an(t)≠0,f(t)≠0.而形如 yt+n+a1(t)yt+n-1+…+an-1(t)yt+1+an(t)yt=0
例 解
求差分方程 t +1 - 2 yt = 5的通解. y
a = -2 -1, b = 5
yt = A 2t - 5, A为任意常数 .
情形Ⅱ f(t)为t的多项式. 不妨设f(t)=b0+b1t(t的一次多项式),即 yt+1+ayt=b0+b1t, t=1,2,…,
其中a,b0,b1均为常数,且a≠0,b1≠0.
yt=(-a)ty0+(-a)t-1f(0)+(-a)t-2f(1)+…+f(t-1)
=(-a)ty0+ y t , (t=0,1,2,…),
其 中 y t = ( - a ) t -1 f ( 0) + ( - a ) t - 2 f (1) + f ( t - 1) = ( - a ) i f ( t - i - 1)
试以特解 y t=a+bt,(a,b为待定系数)代入方程得
a+b (t+1)+a(a+bt)=b0+b1t,
上式对一切t值均成立,其充分必要条件是: (1 + a )a + b = b0 (1 + a )b = b1
当1+a≠0时,即a≠-1时,
b0 b1 a= 1 + a (1 + a )2 b1 b= 1+ a
其中A1,A2,…,An为n个任意(独立)常数.
定理4(非齐次线性差分方程通解结构定理) 如果 y (t)是非齐次线性方程yt+n+a1(t)yt+n-1+a2(t)yt+n-2
+…+an-1(t)yt+1+an(t)yt=f(t)的一个特解,yA(t)是其对应的齐
次线性方程yt+n+a1yt+n-1 +a2yt+n-2 +…+an-1yt+1+anyt=0的通
程中出现.
定义3′ 含有两个或两个以上函数值yt,yt+1,…的函数方 程,称为(常)差分方程,出现在差分方程中未知函数下 标的最大差,称为差分方程的阶. n阶差分方程的一般形式为
F(t,yt,yt+1,…,yt+n)=0,
其中F为t,yt,yt+1,…,yt+n的已知函数,且yt和yt+n一定要 在差分方程中出现.
1. 一个半球状雪堆,其体积融化的速率与半 球面面积S成正比,比例系数k > 0。设融化中 雪堆始终保持半球状,初始半径为R且3小时 中融化了总体积的7/8,问雪堆全部融化还需 要多长时间?
2、 早期肿瘤的体积增长满足 Malthus 模型 ( =λ V,其中λ 为常数) (1)求肿瘤的增倍时 , 间σ 。根据统计资料,一般有σ (7,465)(单 位为天) ,肺部恶性肿瘤的增倍时间大多大于 70 天而小于 465 天(发展太快与太慢一般都不 是恶性肿瘤) 故σ 是确定肿瘤性质的重要参数 , 之一(2)为方便起见,医生通常用肿瘤直径来 表示肿瘤的大小,试推出医生用来预测病人肿 瘤直径增大速度的公式 D=
是方程的解,其中A1,A2,…,Am为任意常数.
n阶齐次线性差分方程yt+n+a1yt+n-1 +a2yt+n-2
定理2
+…+an-1yt+1+anyt=0一定存在n个线性无关的特解.
定理3(齐次线性差分方程通解结构定理) 如 果 y1(t),y2(t),…,yn(t) 是 齐 次 线 性 差 分 方 程 yt+n+a1yt+n-1 +a2yt+n-2 +…+an-1yt+1+anyt=0的n个线性无关 的特解,则方程的通解为: yA(t)=A1y1(t)+A2y2(t)+…+Anyn(t),