作者CRAZYAWIND (怒火烧不尽)
看板Grad-ProbAsk
标题Re: [理工] [工数]-[矩阵]
时间Mon Dec 28 14:43:21 2009
※ 引述《mdpming (ｎＥｗ)》之铭言：
: ※ 引述《JaLunPa (呷懒趴)》之铭言：
: : A=[5 -13]
: : [2 -5] 求 lim A^n
: : n→∞
λ= ±i
┌ 5-i -13 ┐ ┌ 13 ┐
当λ = i │ 2 -5-i │ v1 = │5-i │
└ ┘ └ ┘
┌ 5+i -13 ┐ ┌ 13 ┐
λ = -i │ 2 -5+i │ v2 = │ 5+i │
└ ┘ └ ┘
┌ 13 13 ┐ -1 1 ┌ 5+i -13 ┐
令 s = │5-i 5+i │ s = ── │-5+i 13 │
└ ┘ 26i └ ┘
n
n ┌ (i) 0 n ┐ -1
A = S │ 0 (-i) │ S
└ ┘
┌ n n n n ┐
1 │ 13(5+i)(i) +13(-5+i)(-i) -169(i) +169(-i) │
= ── │ n n n n│
26i │26(i) -26(-i) -13(5-i)(i) +13(5+i)(-i) │
└ ┘
我还是比较喜欢做苦工的方法
: : 题目是92台大大气
: : 谢谢
: (A-λI)X = 0
: [5-λ -13] [x1]
: = 0
: [2 -5-λ] [x2]
: λ1 = i , λ2 = -i
: |A-λI| = (λ-i)(λ+i) = 0
: Caley-Hamiton 定理
: 0 = |A-AI|I = (A-iI)(A+iI)
: 利用带余除法
: n
: x = p(x)(x-i)(x+i) + c1*x + c0
: n
: A = p(A)(A-iI)(A+iI) + c1*A + c0I
: ^^^^^^^^^^^^^^^^^^^
: = 0
: n
: A = c1A + c0I
: 取 x = i
: n
: i = c1*i + c0
: 取 x = -i
: n
: -i = -c1*i + c0
: 得
: c0 = 0
: ...
: 我算到这 接不下去了 有高手能只点一下吗@@
: 求出 C0 C1 就完成 99%了 ><"
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 59.105.159.190
1F：推 shinyhaung:我也比较喜欢这个XD 12/28 14:47
2F：推 mdpming:@@... 12/28 14:50
3F：→ mdpming:我控制 第8章 卡好久了..=.= 12/28 14:50
4F：→ CRAZYAWIND:第8章不难 很制式化= = 我都把他当作文在写了 12/28 14:52
5F：推 mdpming:第八章过了 是不是就海阔天空..=.= 12/28 14:53
6F：→ CRAZYAWIND:不 是另一个痛苦的开始= = 12/28 14:54
7F：推 mdpming:真的甲的...难道老师骗人.. 12/28 14:55
8F：→ CRAZYAWIND:我很讨厌马达系统= = 12/28 14:57
9F：推 JaLunPa:但是n→∞耶 这样要怎麽讨论 12/28 14:58
10F：推 mdpming:难道把 n 通通改成 ∞ @@ 12/28 15:00
11F：→ CRAZYAWIND:这就是n→∞了阿= = 这题特徵值不会收敛 就放着了 12/28 15:00
12F：→ CRAZYAWIND:用Caley-Hamiton 也是这个答案 只不过比较好看一点 12/28 15:01
13F：推 mdpming:Caley-Hamiton 要怎麽继续做下去@@ 12/28 15:05
14F：→ CRAZYAWIND:C1= [(i)^n -(-i)^n]/2i c0=0 然後去乘A吧= = 12/28 15:10
15F：推 mdpming:所以 C0 = 0 12/28 16:09