對於householder.doc

Similar documents
普通高等学校本科专业设置管理规定

AU = U λ c 2 c 3 c n C C n,, n U 2 U2 C U 2 = B = b 22 b 23 b 2n b 33 b 3n b nn U = U ( U 2, U AU = = = ( ( U 2 U 2 U AU ( U2 λ λ d 2 d 3 d n b 22 b 2

國立屏東教育大學碩士班研究生共同修業要點

填 写 要 求 一 以 word 文 档 格 式 如 实 填 写 各 项 二 表 格 文 本 中 外 文 名 词 第 一 次 出 现 时, 要 写 清 全 称 和 缩 写, 再 次 出 现 时 可 以 使 用 缩 写 三 有 可 能 涉 密 和 不 宜 大 范 围 公 开 的 内 容 不 可 作 为

untitled

Welch & Bishop, [Kalman60] [Maybeck79] [Sorenson70] [Gelb74, Grewal93, Maybeck79, Lewis86, Brown92, Jacobs93] x R n x k = Ax k 1 + Bu k 1 + w

目 录 ONTENTS 2012 年 第 4 期 ( 总 第 9 期 ) 刊 首 语 陈 玉 兰 图 书 馆 情 缘 图 书 馆 帮 助 我 走 进 文 学 殿 堂 胡 尹 强 1 心 里 种 下 一 亩 田 叶 小 芳 3 馆 藏 撷 萃 国 学 宝 典 要 要 叶 四 库 全 书 曳 杨 云 芳

附3

1

No Petrus Ramus 3 John Donne Thomas Adams John Cotton Thomas Hooker Anthony van Dyck 5 Grand Tour p See Perry

C J. C. Caldwell 訛 輯 輥 訛 輰 輥 Victor Nee 1 輥 輱 訛 ~

: 29 : n ( ),,. T, T +,. y ij i =, 2,, n, j =, 2,, T, y ij y ij = β + jβ 2 + α i + ɛ ij i =, 2,, n, j =, 2,, T, (.) β, β 2,. jβ 2,. β, β 2, α i i, ɛ i

2009 年第 54 卷第 12 期 : 1779 ~ csb.scichina.com SCIENCE IN CHINA PRESS,,, ;, , 200


关于印发广东省“互联网+”现代农业行动计划( 年)的通知

( 总 第 1124 期 ) 浙 江 省 人 民 政 府 主 办 2016 年 7 月 18 日 出 版 目 省 政 府 令 浙 江 省 行 政 执 法 证 件 管 理 办 法 ( 浙 江 省 人 民 政 府 令 第 346 号 ) (3) 省 政 府 办 公 厅 文 件 浙 江 省 人 民 政 府

開 創 科 學 計 算 的 研 究 與 職 場 生 涯 13 候 模 型 與 預 測 天 文 以 及 數 位 內 容 產 業 等 等, 這 麼 多 與 我 們 生 活 息 息 相 關 的 產 業, 背 後 有 沒 有 任 何 的 共 通 點? 數 學, 又 在 這 些 產 業 中 扮 演 了 任 何


<4D F736F F D F B0E6B8DFB1BBD2FDD6B8CAFDC7B0D1D42E646F63>

Experimental design: 實驗設計:


LaDefense Arch Petronas Towers 2009 CCTV MOMA Newmark Hahn Liu 8 Heredia - Zavoni Barranco 9 Heredia - Zavoni Leyva

C35N32.dvi

[9] R Ã : (1) x 0 R A(x 0 ) = 1; (2) α [0 1] Ã α = {x A(x) α} = [A α A α ]. A(x) Ã. R R. Ã 1 m x m α x m α > 0; α A(x) = 1 x m m x m +

EL ECTR IC MACH IN ES AND CON TROL Vol113 No11 Jan. 2009,, (, ) :, X 2Y,,,,,,, P ID P ID P ID,, : ; ; ; P ID : TM33 : A : X

标题

Vol. 36 ( 2016 ) No. 6 J. of Math. (PRC) HS, (, ) :. HS,. HS. : ; HS ; ; Nesterov MR(2010) : 90C05; 65K05 : O221.1 : A : (2016)

<4D F736F F D20B0D7D0A1C1AEBDE0D7A8BFAF C5C5B0E E646F63>

微 分 方 程 是 经 典 数 学 的 一 个 重 要 分 支, 常 用 来 描 述 随 时 间 变 化 的 动 态 系 统, 被 广 泛 应 用 于 物 理 学 工 程 数 学 和 经 济 学 等 领 域. 实 际 上, 系 统 在 随 时 间 的 变 化 过 程 中, 经 常 会 受 到 一 些

#4 ~ #5 12 m m m 1. 5 m # m mm m Z4 Z5

Microsoft Word 博士招生简章.doc

关于2007年硕士研究生培养方案修订几点要求的说明

微积分 授课讲义

說 等 目 前 可 以 看 到 的 是 政 府 擬 訂 了 推 動 計 畫 與 目 標, 也 有 發 展 的 技 術 條 件, 但 電 子 書 的 內 容 卻 非 常 匱 乏 數 位 出 版 產 業 發 展 策 略 與 行 動 計 畫 這 項 計 畫 目 標 是 五 年 內 提 供 十 萬 本 電

Microsoft Word - concept-culture.doc

1. 引 言 1.1 職 業 訓 練 局 的 高 峰 進 修 學 院 致 力 為 所 有 金 融 服 務 業 從 業 員 或 有 意 晉 身 該 行 業 的 人 士 提 供 優 質 專 業 培 訓 課 程, 以 助 香 港 維 持 其 領 先 國 際 金 融 中 心 的 地 位. 1.2 香 港 証

中国政法大学2011年本科教学质量报告

<4D F736F F D C6D5CDA8B8DFB5C8D1A7D0A3B1BEBFC6D7A8D2B5C9E8D6C3C9EAC7EBB1ED>

鋼構造論文集第 20 巻第 79 号 (2013 年 9 月 ) AN EVALUATION METHOD FOR ULTIMATE COMPRESSIVE STRENGTH OF STAINLESS STEEL PLATES BASED ON STRESS-STRAIN DIAGRAM * **

<4D F736F F D DB3AFC3FDA4E52DBDD7A4E5A5FEA4E52E646F63>

當前台灣制憲修憲問題之試析─以中華民國憲法之病理病徵為中心

rate extent (2-3) physiological availability systematic availability bio-phase availability biological availability bioavailability 1977 FDA Bio

WL100014ZW.PDF

勤 學 * 卓 越 * 快 樂 成 長 本 校 在 老 師 群 策 群 力 共 同 討 論 下, 型 塑 了 學 校 願 景 : 勤 學 卓 越 快 樂 成 長 ( 一 ) 勤 學 運 用 真 的 力 量 培 養 勤 學, 以 語 文 教 為 基 礎 紮 根 ( 二 ) 卓 越 利 用 美 的 感

标题

穨ld-sms-c

警察人員教考用合一制度之研究

校字〔2007〕号

(domestic crime) (World Trade Organization, WTO) ( 94 89) (Washington Post) ( )

Microsoft Word - 146專題-12 陳振盛.李麗雲.doc

2 199 Navier-Stokes { u t - v 2 u + u u + 1 p = f 1 ρ u = 0 ux y Case25CTA z= { u LC6 60% T n = 0 2 CAS Case25pre ux y z= 0 Case25post u = u x u

24 系 统 与 控 制 纵 横 2015 年 第 1 期 定 系 统 及 导 航 系 统 (1961) 等 这 些 教 材 对 我 们 在 北 大 新 学 科 方 向 的 建 设 帮 助 很 大 二 师 从 林 先 生 1962 年 我 国 对 国 民 经 济 实 行 调 整 巩 固 充 实 提

Microsoft Word - 07.docx

untitled

Microsoft Word - medi_ethics_chi.doc

邱顯比、林清珮、楊宗庭

( ) : ( ) (CIP) /.. :,003. () ISBN O4 44 CIP (00) : : 7 : 7007 : (09 ) : : :850 mm 68 mm / 3 :0.5 :60 :00 0

现代天文学7.ppt

\\Pc6\f\湖工大(自然科学)\自科

successful and it testified the validity of the designing and construction of the excavation engineering in soft soil. Key words subway tunnel

208-教育-文化

Lecture #4: Several notes 1. Recommend this book, see Chap and 3 for the basics about Matlab. [1] S. C. Chapra, Applied Numerical Methods with MATLAB

m m m

No

Microsoft Word ±q¾Ôª§¹q¼v.doc

MAXQ BA ( ) / 20

ebook105-12

上市公司股东大会投票信息公告( )

上市公司股东大会投票信息公告( )

49274h1.pdf

(p.29). (a) F Qq r 2 ()() N (b) Q 2 r 2 F ( 2 )() Q 0 5 C 2. (a) F (b) F 3. 7 (p.42). (a) T (b) F (c) T 2. (a) A (b) (c) 4. (a) 4 (b) (

编 委 会 主 任 : 樊 明 武 副 主 任 : 王 建 枝 叶 朝 辉 刘 献 君 罗 俊 林 萍 华 郑 楚 光 黄 光 英 程 时 杰 编 委 : 王 少 平 王 红 卫 王 国 斌 王 擎 文 志 雄 邓 晓 芒 史 铁 林 吴 毅 李 高 李 亮 李 保 峰 李 维 嘉 陈 安 民 肖

Microsoft PowerPoint 年招生宣传.ppt

北京电影学院

五 研 究 方 向 及 学 制 研 究 方 向 导 师 : 见 附 件 1 学 制 : 三 年 招 生 名 额 : 预 计 招 收 7 名 ( 以 国 家 当 年 正 式 下 达 的 计 划 为 准 ) 2

458 (25),. [1 4], [5, 6].,, ( ).,,, ;,,,. Xie Li (28),,. [9] HJB,,,, Legendre [7, 8],.,. 2. ( ), x = x x = x x x2 n x = (x 1, x 2,..., x

的 ( 四 ) 非 法 集 资 有 哪 些 主 要 表 现 形 式? 非 法 集 资 活 动 涉 及 内 容 广 泛, 表 现 形 式 多 样, 主 要 有 以 下 几 种 : 1 不 具 有 房 产 销 售 的 真 实 内 容 或 者 不 以 房 产 销 售 为 主 要 目 的, 以 返 本 销

<4D F736F F D DBED4B2A4A9D2A9DBA5CDB8EAB054A7B9BEE3AAA92DB3D5A468AF5A2E646F63>

(creep) 500~ ~ 30 MPa 9Cr-1Mo 100 MPa ( ) 9Cr-1Mo Threshold Stress Larson-Miller Manson-Haferd Ω

正式发文

A Thesis in Applied Mathematics Surface of second degree in 3-Minkowski space b Sun Yan Supervisor: Professor Liu Huili Northeastern Universit Decembe

11 : 1341,( r) g { t ( r) exp [i ( t - ) ]} - ( r) 5 { t ( r) exp [i ( t - ) ]} 5 t =, (7) = n ( r) (5) (4), 5 a 5 t ( r) exp [i ( t - ) ] + 5 a 5 ( -

Introduction to Hamilton-Jacobi Equations and Periodic Homogenization

14 43 Cholesky θ jm ω H jm ω N Cholesky n n 2 / Yang 8 FF 1 FF Cholesky V j pδt = Re j h jm qδt exp i mδω { [ ( ) pδt m = 1 n ]


行政院大陸委員會.民國九十四年年報(95.4)

国盛证券投资报告

untitled

投影片 1

VLBI2010 [2] 1 mm EOP VLBI VLBI [3 5] VLBI h [6 11] VLBI VLBI VLBI VLBI VLBI GPS GPS ( ) [12] VLBI 10 m VLBI 65 m [13,14] (referen

<4D F736F F D20CAFDD1A7D1A7D4BABACFB2A2B3F6C6ACCEC4BCFE2E646F63>

2 职 业 核 心 知 识 ( 66 学 分 ) 经 营 环 境 分 析 财 务 管 理 实 务 管 理 学 基 础 统 计 学 基 础 初 级 会 计 中 级 会 计 企 业 制 度 设 计 出 纳 实 务 审 计 基 础 与 实 务 职 业 资 格 培 训 ( 第 三 五 学 期 ) 0 7.5

<4D F736F F D20A661B0ECAE74B2A7AABAA4A4B0EAB941A7F8B0F2C2A6B1D0A87C2E646F63>

運動觀光及休閒產業管理學術研討會論文格式

Trigonometric identities

<453A5C6E6F F6F6B5CBFECC5CC5CB9A4D4BA5C31322D3133D1A7C4EA5CB5DA32D1A7C6DA5CC6E4CBFB5CBEABC6B7BFCE5C C4EAB6C8C8ABCAA1B8DFB5C8D6B0D2B5D1A7D0A3BEABC6B7BFCEB3CCC9EAB1A82DB9A4B3CCC1A6D1A7342E646F63>


11 25 stable state. These conclusions were basically consistent with the analysis results of the multi - stage landslide in loess area with the Monte


0 9 Dryden 80 Waverley l Charles Firth A Commentary on Macaulay s History of England London Macmillan and Co. 98. William R Thayer Macaulay Fif

Transcription:

1. 2. 3. Householder 4. 5. Exponential (e At ) 6. e At 7. (a) Euler [1] (b) [2,3,4,5,6,7] (c) [8] (d) [9] 8. -2-

(1) A (A T = A) Cayley-Hamilton [10] [11] e At t e At -3-

(2) 1. Householder n R n n n H H=I 2(nn T /n n) Householder Householder (Householder reflection) (H=H T ) (H T H=I) (H 2 Householder H ν (u)=u 2(nn T /n n)u H ν =I 2(nn T /n T n) I H ν = I 2(nn T /n T n) Householder H ν (u)=u 2(nn T /n T n){u} n H ν (a) H ν (u)= u 2 (nn T /n n)u H ν = I 2(nn T /n n) (b) (H ν ) 2 =I H ν -1 =H ν (c) H (A) H T AH=D (2.1.1) D -4-

(d) n=2 ν R 2 cosα n= sin α H ν =I- 2 ν T ν cos2α sin2α ν 2 = sin2α cos2α (2.1.2) u H n n H n (u) 1 X X 2 u n H n (u) X 1 X 1 H n (n) (1) ν (X 1 cos+x 2 sin=0) p q r n=3 n= (p 2 +q 2 +r 2 =1) H= 1 0 0 0 1 0 0 0 1 p r -2 q { p q r} = 2 1 2 p 2 pq 2 pr 2 pq 1 2q 2 2qr 2 pr 2qr 2 1 2r -5-

2. e At x & = Ax, x(0)=x o (2.2.1) x x& x(t) o A n n A (2.2.1) x=e At x o (2.2.2) A P n n + P n-1 n-1 + + P 1 + P o =0 (2.2.3) Cayley-Hamilton [5] P n n + P n-1 n-1 + + P 1 + P o I=0 (2.2.4) -6-

3. A R n (or C n ) A T = A [10] 4. A AC=CD D C A A=CDC -1 (2.4.1) f(x) (x-a) f(a) f(x)=(x a)q(x)+f(a) (2.4.2) Q(x) f(a) f(x)=(ax 2 +bx+c)q(x)+px+q (2.4.3) Q(x) px+q x f(a)=(aa 2 +ba+ci)q(a)+pa+qi (2.4.4) a,b,c,p,q I [11] -7-

5. A 3 3 A= 0 c b c 0 a b a 0 (2.5.1) 2 2 2 A 0a + b + c i 2 2 a + b + f(a)=e At c 2 i f(a) =[a(t)a 3 +b(t)a 2 +c(t)a+d(t)i]s(t)+ [p(t)a 2 +q(t)a+r(t)i] (2.5.2) s(t) p(t)a 2 +q(t)a+ri e At ( (similar transform)[5] Q B=Q -1 AQ B A ) e At 1 cos ω t = ω 2 A 2 1 cos ω t + ω 2 A+I (2.5.3) ω= a + 2 2 2 + b c 6. (2.5.2) e At t Q=e At QQ T =I (2.6.1) -8-

a. t=0 e At I( ) b. t=π/2 e At?? c. t=π e At?? ( Householder Givens QR Jacobi ) e At t 2 e At, t R t =? Householder ( H n =I-2(nn T /n n) ) ν n=(p,q,r)=? A A T =A n(2) -9-

Cayley-Hamilton [10] t (1) A 2 2 3 3,, n n (2) 2 2 4 4 (3) (Euler) (4) t (mapping) (5) [9] -10-

1. 2. Cayley-Hamilton 3. t 4. t 5. e At 6. -11-

1. H. Goldstein, Classical Mechanics, 2 nd ed., Addision Wesley, London, 1985. 2. I. E. Leonard, The Matrix Exponential, SIAM Rev., 38(1996), pp.507-512. 3. Thomas J. R. Hughes and James Winget, Finite Rotation Effects in Numerical Integration of Rate Constitutive Equations Arising in Large Deformation Analysis, Short Comm., 1979. 4. R. Rubinstein and S. N. Atluri, Objectivity of Incremental Constitutive Relations over Finite Time Steps in Computational Finite Deformation Analysis, Comp. Meth. App. Mech. Eng., Vol.36, pp.277-290, 1983. 5. D. B. Metzger and R. N. Dubey, Corotational Rates in Constitutive Modeling of Elastic-Plastic Deformation, J. Plasticity, Vol.4, pp.341-368,1987. 6. Juan C. Simo and Jerrold E. Marsden, On the Rotated Stress Tensor and Material Version of the Doyle-Ericksen Formula, 1983. 7. D. P. Flanagan and L. M. Taylor, An Accurate Numerical Algorithin for Stress Integration with Finite Rotations, Comp. Meth. App. Mech. Eng., Vol. 62, pp.305-320, 1987. -12-

8. G. H. Golub and C. F. Van Loan, Matrix Computation, 2 nd ed., The Johns Hopkins University Press, Baltimore, 1987. 9.,, NSC-88-2211-E-019-005,,, 1998. 10. C. Ray Wylie and Louis C. Barrett, Advanced Engineering Mathematics, 6 th, McGraw-Hill, 1995. 11., -,, 15 4, 1987, pp261-265. 12. D. J. Winter, Matrix Algebra, Macmillin Publishing Company, Macmillan, Republic of Singpore, 1992. 13. R. A. Horn and C. R. Johnson, Matrix Analysis, Taipei Publications Frading Company, 1990. 14. A Jennings and J. J. Mckeown, Matrix Computation 2 nd Edition, John Wiley & Sons Ltd., Baffins Lane, Chichester, West Sussex, England, 1992. -13-

15. R. Bellman, Introduction to Matrix Analysis, Society Industrial and Applied Mathematics, Philadelphia, 1995. 16. M. M. Mehrabadi, S. C. Cowin and J. Jaric, Six-Dimensional Orthogonal Tensor Representation of The Rotation about an Axis in Three Dimensions, Int. J. Solids Structures, 32(1995), pp.439-449. 17. Y. Saad, Iterative Methods for Sparse Linear Systems, PWS, 1996. 18. P. Lancaster and M. Tismenetsky, The Theory of Matrices, 2 nd ed., Academic Press, 1985. 19. C. Morler and C. V. Loan, Nineteen Dubious Ways to Compute the Exponential of a Matrix, SIAM Review, Vol.20, No.4, pp.801-836, 1978. 20. M. Vidyasagar, A Novel Method of Evaluating e At in Closed Form, IEEE Transactions on Automatic Control, pp.600-601, 1970. 21. R. C. Ward, Numerical Computation of the Matrix Exponential with Accuracy Estimate, SIAM J. Numer Anal., Vol.14, No.4, pp.600-610, 1977. -14-

-15- ( NSC-88-2211-E-019-005) (1) ( ) (2) (Mathematica Matlab) (3)