函数矩阵与矩阵微分方程讲稿.ppt

《函数矩阵与矩阵微分方程讲稿.ppt》由会员分享，可在线阅读，更多相关《函数矩阵与矩阵微分方程讲稿.ppt（46页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、第一章第一章 第一节第一节 函数函数关于函数矩阵与矩阵微分方程北京理工大学高数教研室*第一页，讲稿共四十六页哦称为函数矩阵，其中所有的元素称为函数矩阵，其中所有的元素都是定义在闭区间都是定义在闭区间 上的实函数。上的实函数。函数矩阵与数字矩阵一样也有加法，数乘，乘法，函数矩阵与数字矩阵一样也有加法，数乘，乘法，转置等几种运算，并且运算法则完全相同。转置等几种运算，并且运算法则完全相同。例：例：已知已知(),1,2,;1,2,ijaxim jn,a b北京理工大学高数教研室*第二页，讲稿共四十六页哦1sin1cos,11xxxxxxABexex计算计算定义：定义：设设 为一个为一个 阶函数矩阵，

2、如果存阶函数矩阵，如果存在在 阶函数矩阵阶函数矩阵 使得对于任何使得对于任何 都有都有那么我们称那么我们称 在区间在区间 是是可逆的可逆的。,2()TxAB AB AABn()B x,xa bn()()()()A x B xB x A xI()A x,a b()A x北京理工大学高数教研室*第三页，讲稿共四十六页哦称称 是是 的逆矩阵，一般记为的逆矩阵，一般记为例例：已知已知那么那么 在区间在区间 上是可逆的，其逆为上是可逆的，其逆为()B x()A x1()Ax11()0 xxA xe()A x1()10 xxxxeAxe3,5北京理工大学高数教研室*第四页，讲稿共四十六页哦函数矩阵可逆的充

3、分必要条件函数矩阵可逆的充分必要条件定理定理:阶矩阵阶矩阵 在区间在区间 上可逆的充分必上可逆的充分必要条件是要条件是 在在 上处处不为零，并且上处处不为零，并且其中其中 为矩阵为矩阵 的伴随矩阵。的伴随矩阵。定义：定义：区间区间 上的上的 型矩阵函数不恒型矩阵函数不恒等于零的子式的最高阶数称为等于零的子式的最高阶数称为 的的秩秩。mn()A x,a b()A x,a b1*1()()()AxA xA x*()A x()A x,a b()A x北京理工大学高数教研室*第五页，讲稿共四十六页哦特别地，设特别地，设 为区间为区间 上的上的 阶矩阵函数，阶矩阵函数，如果如果 的秩为的秩为 ，则称，则

4、称 一个一个满秩矩阵满秩矩阵。注意：对于阶矩阵函数而言，满秩与可逆不是等价注意：对于阶矩阵函数而言，满秩与可逆不是等价的。即：可逆的一定是满秩的，但是满秩的却不一的。即：可逆的一定是满秩的，但是满秩的却不一定是可逆的。定是可逆的。例例 ：已知已知()A x,a bn()A xn()A x201()A xxx北京理工大学高数教研室*第六页，讲稿共四十六页哦那么那么 。于是。于是 在任何区间在任何区间 上的秩都是上的秩都是2。即。即 是满秩的。但是是满秩的。但是 在在 上是否可逆，完全依赖于上是否可逆，完全依赖于 的取值。当区间的取值。当区间 包含有原点时，包含有原点时，在在 上有零点，从而上有零

5、点，从而 是不可逆的是不可逆的。定义：定义：如果如果 的所有各元素的所有各元素 在在 处有极限，即处有极限，即()A xx()A x,a b()A x()A x,a b,a b,a b()A x,a b()A x()()ijm nA xax()ijax0 xx北京理工大学高数教研室*第七页，讲稿共四十六页哦0lim()(1,;1,)ijijxxaxaim jn其中其中 为固定常数。则称为固定常数。则称 在在 处有处有极限极限，且记为且记为其中其中ija0 xx0lim()xxA xA111212122212nnmmmnaaaaaaAaaa()A x北京理工大学高数教研室*第八页，讲稿共四十六页

6、哦如果如果 的各元素的各元素 在在 处连续，即处连续，即则称则称 在在 处处连续连续，且记为，且记为其中其中()A x()ijax0 xx00lim()()(1,;1,)ijijxxaxaxim jn()A x0 xx00lim()()xxA xA x北京理工大学高数教研室*第九页，讲稿共四十六页哦1101201021022020010200()()()()()()()()()()nnmmmnaxaxaxaxaxaxA xaxaxax容易验证下面的等式是成立的：容易验证下面的等式是成立的：设设则则00lim(),lim()xxxxA xAB xB0(1)lim()()xxA xB xAB北京理

7、工大学高数教研室*第十页，讲稿共四十六页哦00(2)lim()(3)lim()()xxxxkA xkAA x B xAB定义定义：如果如果 的所有各元素的所有各元素 在点在点 处处(或在区间或在区间 上上)可导，便称此函数矩阵可导，便称此函数矩阵 在点在点 处处(或在区间或在区间 上上)可导可导，并且，并且记为记为()()ijm nA xax()(1,;1,)ijax im jn0 xx,a b()A x0 xx,a b北京理工大学高数教研室*第十一页，讲稿共四十六页哦00000110120102102202010200d()()()()limd()()()()()()()()()xx xnn

8、mmmnA xA xxA xA xxxaxaxaxaxaxaxaxaxax 北京理工大学高数教研室*第十二页，讲稿共四十六页哦函数矩阵的导数运算有下列性质：函数矩阵的导数运算有下列性质：是常数矩阵的充分必要条件是是常数矩阵的充分必要条件是 设设(1)均可导，则均可导，则 ()A xd()0dA xx()(),()()ijm nijm nA xaxB xb xdd()d()()()dddA xB xA xB xxxx北京理工大学高数教研室*第十三页，讲稿共四十六页哦dd()d()()()()()dddk xA xk x A xA xk xxxx设设 是是 的纯量函数，的纯量函数，是函数矩是函数矩

9、阵，阵，与与 均可导，则均可导，则特别地，当特别地，当 是常数是常数 时有时有()k xx()A x()k x()A x()k xkdd()()ddA xkA xkxx北京理工大学高数教研室*第十四页，讲稿共四十六页哦(4)设设 均可导，且均可导，且 与与 是是可乘的，则可乘的，则因为矩阵没有交换律，所以因为矩阵没有交换律，所以(),()A x B x()A x()B xdd()d()()()()()dddA xB xA x B xB xA xxxx232dd()()2()dddd()()3()ddA xAxA xxxA xA xAxxx北京理工大学高数教研室*第十五页，讲稿共四十六页哦(5)

10、如果如果 与与 均可导，则均可导，则(6)设设 为矩阵函数，为矩阵函数，是是 的纯量函数，的纯量函数，与与 均可导，则均可导，则()A x1()Ax111d()d()()()ddAxA xAxAxxx()A x()xf tt()A x()f tdd()d()()()()dddA xA xA xftftxxx北京理工大学高数教研室*第十六页，讲稿共四十六页哦定义：定义：如果函数矩阵如果函数矩阵 的所有各的所有各元素元素 在在 上可积，则称上可积，则称 在在 上上可积可积，且，且()()ijm nA xax()(1,;1,)ijax im jn,a b()A x,a b111212122212()

11、d()d()d()d()d()d()d()d()d()dbbbnaaabbbbnaaaabbbmmmnaaaaxxaxxaxxaxxaxxaxxA xxaxxaxxaxx北京理工大学高数教研室*第十七页，讲稿共四十六页哦()d()d()()d()d()dbbaabbbaaakA xxkA xxkRA xB xxA xxB xx函数矩阵的定积分具有如下性质：函数矩阵的定积分具有如下性质：例例1：已知函数矩阵已知函数矩阵试计算试计算21()0 xA xx北京理工大学高数教研室*第十八页，讲稿共四十六页哦23231(1)(),(),()(2)()(3)()dddA xA xA xdxdxdxdA x

12、dxdAxdx证明：证明：02()10 xdA xdx220()00 xdA xdx北京理工大学高数教研室*第十九页，讲稿共四十六页哦由于由于 ，所以，所以下面求下面求 。由伴随矩阵公式可得。由伴随矩阵公式可得 3300()00dA xdx3()A xx 2()3dA xxdx 1()Ax北京理工大学高数教研室*第二十页，讲稿共四十六页哦1*23231()()()1001111AxA xA xxxxxxx 再求再求1()dAxdx北京理工大学高数教研室*第二十一页，讲稿共四十六页哦213410()23dxAxdxx例例2：已知函数矩阵已知函数矩阵23sincossin()10 xxxxxA x

13、exxx北京理工大学高数教研室*第二十二页，讲稿共四十六页哦试求试求d(2)()dA xx0(1)lim()xA xd(5)()dA xx22d(3)()dA xxd(4)()dA xx例例3：已知函数矩阵已知函数矩阵试求试求证明：证明：sincos()cossinxxA xxx200(),()xxA x dxA x dx北京理工大学高数教研室*第二十三页，讲稿共四十六页哦00000sincos()cossin1cossinsin1cosxxxxxxdxxdxA x dxxdxxdxxxxx同样可以求得同样可以求得222220sincos()2cossinxxxA x dxxxx北京理工大学高

14、数教研室*第二十四页，讲稿共四十六页哦例例4：已知函数矩阵已知函数矩阵试计算试计算22()20300 xxxxexexA xeex3100(),()xA x dxA x dx北京理工大学高数教研室*第二十五页，讲稿共四十六页哦定义定义：设有定义在区间设有定义在区间 上的上的 个连续的函个连续的函数向量数向量如果存在一组不全为零的常实数如果存在一组不全为零的常实数使得对于所有的使得对于所有的 等式等式成立，我们称，在成立，我们称，在 上上 ,a bm12()(),(),()(1,2,)iiiinxax axaxim12,mk kk,xa b1122()()()0mmkxkxkx,a b12(),

15、(),()mxxx线性相关线性相关。北京理工大学高数教研室*第二十六页，讲稿共四十六页哦12(),(),()mxxx否则就说否则就说 线性无关。即如果线性无关。即如果只有在只有在 等式才成立，那么就说等式才成立，那么就说 线性无关。线性无关。定义定义：设设 是是 个定义在个定义在区间区间 上的连续函数向量上的连续函数向量记记120mkkk12(),(),()mxxx12(),(),()mxxxm,a b12()(),(),()(1,2,)iiiinxax axaxim()()d(,1,2,)bTijijagxxxi jm北京理工大学高数教研室*第二十七页，讲稿共四十六页哦以以 为元素的常数矩阵

16、为元素的常数矩阵称为称为 的的Gram矩阵，矩阵，称为称为Gram行列式行列式。定理定理：定义在区间定义在区间 上的连续函数向量上的连续函数向量 线性无关的充要条件是它线性无关的充要条件是它的的Gram矩阵为满秩矩阵。矩阵为满秩矩阵。ijg()ijm nGg12(),(),()mxxxdetG,a b12(),(),()mxxx北京理工大学高数教研室*第二十八页，讲稿共四十六页哦12()(0,),()(,0)xxxx例例：设设则则于是于是 的的Gram矩阵为矩阵为233111221233221d()301d()3babagxxbagggxxba12(),()xx北京理工大学高数教研室*第二十九

17、页，讲稿共四十六页哦33331()0310()3baGba所以所以故当故当 时，时，在在 上是线性无关的。上是线性无关的。33 21det()9Gba12det0,(),()Gxxab,a b北京理工大学高数教研室*第三十页，讲稿共四十六页哦定义：定义：设设 是是 个定义在区间个定义在区间 上的上的 有有 阶导数的函数向量，记阶导数的函数向量，记那么称矩阵那么称矩阵12()(),(),()iiiinxax axax(1,2,)imm1m ,a b11112122122212()()()()()()()()()()()()()nnmmmmnxaxaxaxxaxaxaxA xxaxaxax北京理工

18、大学高数教研室*第三十一页，讲稿共四十六页哦(1)111212122212(1)(1)(1)11121(1)(1)(1)21222(1)(1)12()()()()()()()()()()()()()()()()()()()()()mm mnnnmmmnmmmnmmmnmmmmW xA x A xAxaxaxaxaxaxaxaxaxaxaxaxaxaxaxaxaxaxa(1)()mmnx北京理工大学高数教研室*第三十二页，讲稿共四十六页哦是是 的的Wronski矩阵。矩阵。12(),(),()mxxx(1)(),(),()mA xA xAx其中其中 分别是分别是 的一阶，二阶，的一阶，二阶，阶导

19、数矩阵。阶导数矩阵。定理：定理：设设 是是 的的Wronski矩阵。如果在区间矩阵。如果在区间 上的某个点上的某个点 ，常数矩阵，常数矩阵 的秩等于的秩等于 ，则，则向量向量 在在 上线性无上线性无关。关。()A x1m()W x12(),(),()mxxx0,xa b0()W xm12(),(),()mxxx,a b,a b北京理工大学高数教研室*第三十三页，讲稿共四十六页哦例例：设设则则因为因为 的秩为的秩为2，所以，所以 与与 线性无线性无关。关。212()(1,),()(,1,)xxx xxex221()1012()011012()101xxxxxxA xexxA xexxxW xex

20、e(0)W1()x2()x北京理工大学高数教研室*第三十四页，讲稿共四十六页哦 1111122112211222221122()()()()()()()()()()()()()()()()()()()()()nnnnnnnnnnndxa t x ta t x ta t x tf tdtdxa t x ta t x ta t x tf tdtdxa t x ta t x ta t x tf tdt 北京理工大学高数教研室*第三十五页，讲稿共四十六页哦()()()()dx tA t x tf tdt1 11 212 12 22121212()()()()()()()()()()()(),(),()

21、()(),(),()nnnnn nTnTnatatatatatatA tatatatx txtxtxtftftftft北京理工大学高数教研室*第三十六页，讲稿共四十六页哦1112121222121212()()()()()()()()()()()(),(),()()(),(),()nnnnnnTnTnatatatatatatA tatatatx tx tx tx tf tf tf tft北京理工大学高数教研室*第三十七页，讲稿共四十六页哦上述方程组的初始条件为上述方程组的初始条件为可以表示成可以表示成定理：定理：设设 是一个是一个 阶常数矩阵，则微分方程组阶常数矩阵，则微分方程组满足初始条件满

22、足初始条件 的解为的解为1010202000(),(),()nnx txx txx tx010200(),Tnx txxxAn()()dx tAx tdt00()x tx北京理工大学高数教研室*第三十八页，讲稿共四十六页哦0()0A t txex定理：定理：设设 是一个是一个 阶常数矩阵，则微分方程阶常数矩阵，则微分方程组组满足初始条件满足初始条件 的解为的解为例例1：设设An()()()dx tAx tf tdt00()x tx000()()0()tA t tA t ttxexefd北京理工大学高数教研室*第三十九页，讲稿共四十六页哦126103114A 求微分方程组求微分方程组 满足初始条

23、满足初始条件件 的解。的解。解：解：首先计算出矩阵函数首先计算出矩阵函数()()dx tAx tdt(0)1 1 1Tx北京理工大学高数教研室*第四十页，讲稿共四十六页哦(1 2)26(1)3(1 3)tttAtttttttt eteteetet etetetet e由前面的定理可知微分方程组由前面的定理可知微分方程组满足初始条件满足初始条件 的解为的解为()()dx tAx tdt(0)1 1 1Tx(0)(1 2)(1)(1)TAttttxe xt et et e北京理工大学高数教研室*第四十一页，讲稿共四十六页哦例例2：设设126103,()0114tteAf te 求微分方程组求微分方

24、程组 满足初始满足初始条件条件 的解。的解。解：解：由上述定理可知满足所给初始条件的微分方程由上述定理可知满足所给初始条件的微分方程组解为组解为()()()dx tAx tf tdt(0)1 1 1Tx北京理工大学高数教研室*第四十二页，讲稿共四十六页哦()0()(0)()tAtA tx te xefd由上面的例题可知由上面的例题可知而而(0)(1 2)(1)(1)TAtttte xt et et e北京理工大学高数教研室*第四十三页，讲稿共四十六页哦()(1 2 2)2()6()()(1)3()()()(1 3 3)tttAtttttttteteteetetetetetete 所以有所以有()1 8()()4()1 4()A tttefett 北京理工大学高数教研室*第四十四页，讲稿共四十六页哦0()000222 18()()4()14()422tttA tttttdefdetdtdttettt222(14)()(12)(122)tttttex tttette 故有故有北京理工大学高数教研室*第四十五页，讲稿共四十六页哦第一章第一章 第一节第一节 函数函数感谢大家观看感谢大家观看9/5/2022第四十六页，讲稿共四十六页哦