数学建模课件(3)

2019-08-30 23:57

定义3.4 如果y1(k),y2(k),?,yn(k)是齐次方程

?(y)[En?p1(k)En?1???pn(k)I]y?0 (3.2)

的n个解,c(k)是这些函数的Casorati矩阵,且如果对于某个m?N,det[c(m)]?0,则?y1(k)?Ey1(k)这n个解在N上是线性相关的 c(k)?????n?1y1(k)?Ey2(k)Ey2(k)?En?1????y2(k)??Eyn(k)? ???n?1Eyn(k)?yn(k)证 如果对某个m?N,det[c(m)]?0这就意味着c(m)中的列向量(或行向量)是线性相关的,因此存在一个非零向量 L?(c1,c2,?,cn)T ?y1(m)?Ey1(m)使得 c(m)L?????n?1y1(m)?Ey2(m)Ey2(m)?En?1????y2(m)??c1????Eyn(m)c???2??? ????????n?1Eyn(m)??cn?yn(m)?表示零向量 c(m)L?? (3.8)

令 z(k)?c1y1(k)?c2y2(k)???cnyn(k) 由于yi(k)(i?1,2,?,n)是齐次方程(3.2)的解,所以它们的线性组合z(k)仍旧是(3.2)

n?1的一个解.再则由(3.8)直接推出 z(m)?Ez(m)???Ez(m)?0

因此利用定理3.3知z(k)是一个零函数.即 z(k)?0 k?N. 故导出 c1y1(k)?c2y2(k)???cnyn(k)?? k?N.

因为c1,c2,?cn是不全为零的n个数,所以解y1(k),y2(k),?,yn(k)是线性相关的.

系3.2 令y1(k),y2(k),?,yn(k)是齐次线性差分方程(3.2)的n个解,则 (i)这n个解是线性相关的充分必要条件是:在N上det[c(k)]?0(k?N);

(ii) 这n个解是线性无关的充分必要条件是:对所有的(k?N),det[c(k)]?0. 从这个系的结论知,引进函数的Casorati矩阵之目.当这组函数是齐次线性差分方程的解时,则根据它们的Casorati矩阵行列式在N上是否为0,即可判断这组解函数的线性相关性和线

性无关性.而对任意的n个函数(不受方程约束)而言,这个系的结论就不一定成立.

H是表示齐次程(3.2)的所有解构成之集合.前面已指出H是向量空间V?的一个子

空间.下面就来确定H的维数.

定义3.5 n 阶齐次差分方程(3.2)的任一解可以表示为方程(3.2)的n个线性无关解的一个线性组合.

证 令yi(k)(i?1,2,?,n)是方程(3.2)满足如下的初始条件:

?y1(0)?1,yi(0)?0,y1(2)???y1(n?1)?0??y2(0)?1,y2(1)?0,y2(2)???y2(n?1)?0 (3.9) ???????y(0)?y(1)???y(n?2)?0,y(n?1)?1nnn?n的唯一解.初始条件(3.9)可简单地写成 Ei?1yi(0)?1 i?1,其他初值均为零.显然,当k?0时,这n个解的Casorati矩阵是单位矩阵

?y1(0)?y1(1)c(0)??????y1(n?1)y2(0)y2(1)?y2(n?1)????yn(0)??1??0yn(1)?????????yn(n?1)??00102,?,n

?????00??0? ???1?∴ det[c(k)]?1?0 这就表示了yi(k)(i?1,2,?,n)(k?N)是线性无关的.所以方程(3.2)就n个线性

无关的解.令y(k)是(3.2)的一个解,并考虑

z(k)?y(0)y1(k)?y(1)y2(k)???y(n?1)yn(k)

因为z(k)只是上面所考考的n个线性无关解一个线性组合.所以z(k)也是(3.2)的一个解.而且有 z(0)?y(0)y1(0)?y(1)y2(0)???y(n?1)yn(0)?y(0)

z(1)?y(0)y1(1)?y(1)y2(1)???y(n?1)yn(1)?y(1) ??z(n?1)?y(0)y1(n?1)?y(1)y2(n?1)??

?y(n?1)yn(n?1)?y(n?1)

故对于m?0,1,2,?,n?1时 Ez(0)?z(m)?y(m)

由唯一性定理知,在(k?N)上y(k)?z(k),因此方程(3.2)的任一个解都可以表示为方程(3.2)的n个线性无的解的一个线性组合.

定义3.3 n阶线性齐次差分方程

?(y)?[Enm?p1(k)En?1???pn?1(k)I]y?0 (3.2)

的任意n个线性无关解的一个集合,都称为(3.2)的一个解基.

例3 考虑议程 (E2?3E?2I)y(k)?0 (3.10)

即可写成 y(k?2)?3y(k?1)?2y(k)?0 (3.10)

容易验证y1(k)?1,及y2(k)?2k是此方程的两个线性无关的解.所以这两个函数

ky1(k)?1与y2(k)?2构成了方程(3.10)的一个解基,并且

y(k)?c1?7?c2?2 (3.11)

是二阶线性齐次差分方程(3.10)的含有两个参数c1、c2的解族,实质上就是和在常微分方程教程中所说的通解意义一样,因此我们今天不防称由(3.11)的表示的解为方程(3.10)的通解.

与线性算子 ??En?p1(k)En?1???pn?1(k)E?pn(k)I

有关的是齐次方程 ??[y(k)]?? (3.2) 以及形如 ??[y(k)]?f(k) (3.2) 的非齐次方程,这里f(k)??.

方程(3.2)与(3.12)之间的关系密切,它们的解之间有下面的结论.

定义3.6 如果u(k)及v(k)是非齐次方程(3.12)的解,则y(k)?u(k)?v(k)是(3.12)的两个解.再根据算子?的线性性 ?(u(k))??(u(k))?f(k)?f(k)?0 即 ?(u(k)??(u(k))?0,则 y(k)?u(k)?u(k)是(3.2)的一个解.

定义3.7 如果z(k)是非齐次方程(3.12)的一个解,且y1(k),y2(k),?,yn(k)是齐次方程(3.2)的一个解基,则非齐次方程(3.12)的每一个解可以表示成

y(k)?c1y1(k)?c2y2(k)???cnyn(k)?z(k) 这里c1,c2?,cn是常量.

k证 由假定y(k)与z(k)都是非齐次方程(3.12)的两个解,根据定理3.6知.它们的差y(k)?z(k)是齐次方程(3.2)的一个解,再根据定理3.5知,y(k)?z(k)可以表示成(3.2)的一个解基(y1(k),y2(k),?,yn(k))的一个线性组合.即

ny(k)?z(k)?c1y1(k)?c2y2(k)???cnyn(k),则y(k)??ci?1i定理证毕. yi(k)?z(k).

例4 考虑非齐次方程 (E2?3E?2I)y(k)?3k (3.11) 容易验证z(k)?(E212?3是方程(3.11)的一个解.事实上,将其代入(3.11)之左端得

12?3k?2k?3E?2I)z(k)?z(k?2)?3z(k?1)?2z(k)??3[k?3?12?3k?1?2?12?3

k92?92?1]?3k?(3.11)的右端

由例3 知y1(k)?7与y2(k)?2k是(3.11)的齐次部分,即方程(3.10)的一组解基,由定3.7和(3.11)的任一解可以表示为 y(k)?c1?7?c2?2k?12?3再考虑由(3.11)与

k?y(k?2)?3y(k?1)?2y(k)?3k初始值y(0)?3,y(1)?1构成的初值问题:?

?y(0)?3,y(1)?11?y(0)?c?7?c??312??2其解答如下所示:由初始条件知 ? ?

?y(1)?c?7?2c?3?112?2?11?c??114 ??c??3?2k故此初值问题的唯一解 y(k)?112?3?2?k12?3

上面例3与例4告诉了我们具有常数的二阶线性齐次差分方程(3.10)和非次差分方程(3.11)存在形如y1(k)?7,y2(k)?2k的解基和特解z(k)?12?3.这些解函数是根据什么规则

k和原理,从方程中把它们找出来呢?这就是我们在下一节将要回答的问题.

1.4 常系数线性差分方程

nn?1???an?1t?an是一个n次的实系数多项定义4.1 设 p(t)?a0t?a1t式.an?0,与p(t)有关的是一个n阶算子多项式P(E),其定义如下:

P(E)?a0En?a1En?1??an?1E?a0I

这样一来 P(E)y?? 是n阶齐次差分方程的简略表示

我们指出,在定义4.1中假定an?0,它就相当于p(0)?0这就对应于n阶算子

??En?p1(k)En?1??pn?1(k)E?pn(k)I 中,假定对所有k?N,p(k)?0时的同

kk?1样意义,其次应当注意到Eb?b?b(b)这类似于微分算子D?kddx Debx?bebx

现在我们容易证明下面的结论:

定义4.1 如果P(E)是一个n阶算子多项式,则 P(E)[bk]?P(b)[bk] 证 P(E)[bk]?(a0En?a1En?1???an?1E?an)bk

?a0E[b]?a1E?a0b?a1bnn?1nkn?1[b]???an?1E[b]?anb

kkkk???an?1b?an)b?P(b)[b]

k系4.1 如果P(E)是一个算子多项式,且p(b)?0,则y(k)?bk是

P(E)y?(a0En?a1En?1???an?1E?anI)y?? (4.1)

的一个解.

证 P(E)(bk)?(a0En?a1En?1???an?1E?anI)bk

?(a0b?a1bnn?1???an?1b?an)b?P(b)(b)??(因p(b)?0)

kk故y(k)?bk是齐次方程(4.1)的一个解.当然,如果我们把限制p(0)?0取掉,这时定理4.1与系4.1的结论仍然成立.因为零函数总是线性齐次差分方程(4.1)的一个解.

Ⅰ. 考虑一阶齐次差分方程 a0Ey(k)?a1y(k)?0 (4.2) 这时P(E)?aE?a1I.我们考虑p(t)?a0t?a1.如果a0a1?0,则t??a1a0ka1a0是多项式

即p(?p(t)的一个零点.)?0.故y(k)?(?a1a0)是一阶齐次常系数线性差分方程的一

k个解.因此 y(k)?(?a1a0)y(0) 表示了方程(5.2)的通解.

Ⅱ. 考虑二阶常系数线性齐次差分方程

a0Ey(k)?a1E(k)?a2E(k)?0 (4.3)

2首先我们写出多项式 p2(t)?a0t?a1t?a2

22根据系4.1,我们应当算出此多项式的零点.即要求解二次代数方程 a0t?a1t?a2?0

的根,亦即解辅助方程 t?c1*t?c2*?0 (4.4) 的根.其中c1*?a1a02,c2*?a2a1为常量,且c2*?0.


数学建模课件(3).doc 将本文的Word文档下载到电脑 下载失败或者文档不完整,请联系客服人员解决!

下一篇:scratch教案 - 变量

相关阅读
本类排行
× 注册会员免费下载(下载后可以自由复制和排版)

马上注册会员

注:下载文档有可能“只有目录或者内容不全”等情况,请下载之前注意辨别,如果您已付费且无法下载或内容有问题,请联系我们协助你处理。
微信: QQ: