计算任意矩阵的任意次方

思考题:如何计算如下矩阵的任意次方
( 2 1 0 0 0 2 1 0 0 0 2 1 0 0 0 2 ) \begin{pmatrix}2&1&0&0\\0&2&1&0\\0&0&2&1\\0&0&0&2\end{pmatrix} 2000120001200012
思路:

给出的例子是一个四阶的Jordan块,这是一个提示。我们知道任何方阵都相似于由Jordan阵构成的分块矩阵,也就是 A = T J T − 1 A=TJT^{-1} A=TJT1,其中 J = ( J 1 0 ⋯ 0 0 J 2 ⋯ 0 ⋮ ⋱ ⋮ 0 ⋯ 0 J s ) J=\begin{pmatrix}J_1&0&\cdots&0\\0&J_2&\cdots&0\\\vdots&&\ddots&\vdots\\0&\cdots&0&J_s\end{pmatrix} J=J1000J2000Js是多个循环子空间的直和。所以任意矩阵都可以通过化成Jordan标准型来计算其任意次方:

f ( A ) = A k = ( T ( J 1 0 ⋯ 0 0 J 2 ⋯ 0 ⋮ ⋱ ⋮ 0 ⋯ 0 J s ) T − 1 ) k = T ( J 1 k 0 ⋯ 0 0 J 2 k ⋯ 0 ⋮ ⋱ ⋮ 0 ⋯ 0 J s k ) T − 1 f(A)=A^k=(T\begin{pmatrix}J_1&0&\cdots&0\\0&J_2&\cdots&0\\\vdots&&\ddots&\vdots\\0&\cdots&0&J_s\end{pmatrix}T^{-1})^{k}=T\begin{pmatrix}J_1^{k}&0&\cdots&0\\0&J_2^k&\cdots&0\\\vdots&&\ddots&\vdots\\0&\cdots&0&J_s^k\end{pmatrix}T^{-1} f(A)=Ak=(TJ1000J2000JsT1)k=TJ1k000J2k000JskT1

问题在于如何计算Jordan块的任意次方。

假如是计算 J k , r a n k ( J ) = m J^k,rank(J)=m Jk,rank(J)=m,记 f ( x ) = x k , k ∈ R f(x)=x^k,k\in R f(x)=xk,kR,则
f ( J ) = J k = ( λ 1 0 ⋯ 0 0 λ 1 ⋯ 0 ⋮ ⋱ ⋱ ⋱ ⋮ ⋮ ⋱ λ 1 0 ⋯ ⋯ 0 λ ) k = ( f ( λ ) f ′ ( λ ) / 1 ! f ′ ′ ( λ ) / 2 ! ⋯ f ( m − 1 ) ( λ ) / ( m − 1 ) ! 0 f ( λ ) f ′ ( λ ) / 1 ! ⋯ f ( m − 2 ) ( λ ) / ( m − 2 ) ! ⋮ ⋱ ⋱ ⋱ ⋮ ⋮ ⋱ f ( λ ) f ′ ( λ ) / 1 ! 0 ⋯ ⋯ 0 f ( λ ) ) f(J)=J^k=\begin{pmatrix}\lambda&1&0&\cdots&0\\0&\lambda&1&\cdots&0\\\vdots&\ddots&\ddots&\ddots&\vdots\\\vdots&&\ddots&\lambda&1\\0&\cdots&\cdots&0&\lambda\end{pmatrix}^k= \begin{pmatrix}f(\lambda)&f'(\lambda)/1!&f''(\lambda)/2!&\cdots&f^{(m-1)}(\lambda)/(m-1)!\\0&f(\lambda)&f'(\lambda)/1!&\cdots&f^{(m-2)}(\lambda)/(m-2)!\\\vdots&\ddots&\ddots&\ddots&\vdots\\\vdots&&\ddots&f(\lambda)&f'(\lambda)/1!\\0&\cdots&\cdots&0&f(\lambda)\end{pmatrix} f(J)=Jk=λ001λ01λ0001λk=f(λ)00f(λ)/1!f(λ)f(λ)/2!f(λ)/1!f(λ)0f(m1)(λ)/(m1)!f(m2)(λ)/(m2)!f(λ)/1!f(λ)
在本次思考题中, λ = 2 , m = 4 \lambda=2,m=4 λ=2,m=4,求 k = 0.5 k=0.5 k=0.5次方:
( 2 1 0 0 0 2 1 0 0 0 2 1 0 0 0 2 ) 1 / 2 = ( 2 1 / 2 1 2 ∗ 1 ! ∗ 2 − 1 / 2 − 1 4 ∗ 2 ! ∗ 2 − 3 / 2 3 8 ∗ 3 ! ∗ 2 − 5 / 2 0 2 1 / 2 1 2 ∗ 1 ! ∗ 2 − 1 / 2 − 1 4 ∗ 2 ! ∗ 2 − 3 / 2 0 0 2 1 / 2 1 2 ∗ 1 ! ∗ 2 − 1 / 2 0 0 0 2 1 / 2 ) \begin{pmatrix}2&1&0&0\\0&2&1&0\\0&0&2&1\\0&0&0&2\end{pmatrix}^{1/2}=\begin{pmatrix}2^{1/2}&\frac{1}{2*1!}*2^{-1/2}&-\frac{1}{4*2!}*2^{-3/2}&\frac{3}{8*3!}*2^{-5/2}\\0&2^{1/2}&\frac{1}{2*1!}*2^{-1/2}&-\frac{1}{4*2!}*2^{-3/2}\\0&0&2^{1/2}&\frac{1}{2*1!}*2^{-1/2}\\0&0&0&2^{1/2}\end{pmatrix} 20001200012000121/2=21/200021!121/221/20042!123/221!121/221/2083!325/242!123/221!121/221/2

  • 2
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Quartz是OpenSymphony开源组织在Job scheduling领域又一个开源项目,它可以与J2EE与J2SE应用程序相结合也可以单独使用。Quartz可以用来创建简单或为运行十个,百个,甚至是好几万个Jobs这样复杂的程序。Jobs可以做成标准的Java组件或 EJBs。 Quartz的优势: 1、Quartz是一个任务调度框架(库),它几乎可以集成到任何应用系统中。 2、Quartz是非常灵活的,它让您能够以最“自然”的方式来编写您的项目的代码,实现您所期望的行为 3、Quartz是非常轻量级的,只需要非常少的配置 —— 它实际上可以被跳出框架来使用,如果你的需求是一些相对基本的简单的需求的话。 4、Quartz具有容错机制,并且可以在重启服务的时候持久化(”记忆”)你的定时任务,你的任务也不会丢失。 5、可以通过Quartz,封装成自己的分布式任务调度,实现强大的功能,成为自己的产品。6、有很多的互联网公司也都在使用Quartz。比如美团 Spring是一个很优秀的框架,它无缝的集成了Quartz,简单方便的让企业级应用更好的使用Quartz进行任务的调度。   课程说明:在我们的日常开发中,各种大型系统的开发少不了任务调度,简单的单机任务调度已经满足不了我们的系统需求,复杂的任务会让程序猿头疼, 所以急需一套专门的框架帮助我们去管理定时任务,并且可以在多台机器去执行我们的任务,还要可以管理我们的分布式定时任务。本课程从Quartz框架讲起,由浅到深,从使用到结构分析,再到源码分析,深入解析Quartz、Spring+Quartz,并且会讲解相关原理, 让大家充分的理解这个框架和框架的设计思想。由于互联网的复杂性,为了满足我们特定的需求,需要对Spring+Quartz进行二次开发,整个二次开发过程都会进行讲解。Spring被用在了越来越多的项目中, Quartz也被公认为是比较好用的定时器设置工具,学完这个课程后,不仅仅可以熟练掌握分布式定时任务,还可以深入理解大型框架的设计思想。
[入门数据分析的第一堂课]这是一门为数据分析小白量身打造的课程,你从网络或者公众号收集到很多关于数据分析的知识,但是它们零散不成体系,所以第一堂课首要目标是为你介绍:Ø  什么是数据分析-知其然才知其所以然Ø  为什么要学数据分析-有目标才有动力Ø  数据分析的学习路线-有方向走得更快Ø  数据分析的模型-分析之道,快速形成分析思路Ø  应用案例及场景-分析之术,掌握分析方法[哪些同学适合学习这门课程]想要转行做数据分析师的,零基础亦可工作中需要数据分析技能的,例如运营、产品等对数据分析感兴趣,想要更多了解的[你的收获]n  会为你介绍数据分析的基本情况,为你展现数据分析的全貌。让你清楚知道自己该如何在数据分析地图上行走n  会为你介绍数据分析的分析方法和模型。这部分是讲数据分析的道,只有学会底层逻辑,能够在面对问题时有自己的想法,才能够下一步采取行动n  会为你介绍数据分析的数据处理和常用分析方法。这篇是讲数据分析的术,先有道,后而用术来实现你的想法,得出最终的结论。n  会为你介绍数据分析的应用。学到这里,你对数据分析已经有了初步的认识,并通过一些案例为你展现真实的应用。[专享增值服务]1:一对一答疑         关于课程问题可以通过微信直接询问老师,获得老师的一对一答疑2:转行问题解答         在转行的过程中的相关问题都可以询问老师,可获得一对一咨询机会3:打包资料分享         15本数据分析相关的电子书,一次获得终身学习
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值