※ 引述《BumbEgg (该期待什麽呢?)》之铭言： : Find all eigenvalues and a basis for the corresponding eigenspace for the : ／ 3 1 0 ＼ 100 : matrix . A =∣ 0 1 0 ∣ Use your answer to compute A B : ＼ 4 2 1 ／ : T : where B = 〔 2 2 8 〕 . : -----------------------------------------以上为题目 : 书上解得特徵根为 1 跟 3 : ／ 0 ＼ ／ -1 ＼ : V(1) = ker(A - I) = span( v1 = ∣ 0 ∣ , v2 =∣ 2 ∣ ) : ＼ 1 ／ ＼ 0 ／ : ／ 1 ＼ : V(3) = ker(A - 3I) = span( v3 = ∣ 0 ∣) : ＼ 2 ／ : -------到这步都还能理解,以下就想不出所以然------- : 因为 B = 2(v1) + (v2) + 3(v3) : 100 100 : A B = A 〔2(v1) + (v2) + 3(v3)〕 : 100 100 100 : = 2 A (v1) + A (v2) + 3 A (v3) : 100 100 100 : = 2× 1 (v1) + 1 (v2) +3×3 (v3) : =... : 这是书上的解法,部分打出 : 1.我不了解为何 B 要用, ker 1 跟 3 的 base 去取代 ? 有哪一章节的定理提 : 到这方法吗 ? 不太能理解为何要这样处理 : 2.另ㄧ个不懂得点是, A = PD(P^-1) , P与D皆可求出,但书上的解法是直接将 : A的100次方带入, 而并不是 (P)乘(D的100次)乘(P的反矩阵) 再乘 B 矩阵.. : 3. 若想不出书上的解法式不是只能用( A^100 ) = (P)(D^100)(P^-1)处理? : 有其他的方法吗 ? 这题计算还蛮繁复的.. 我是指其实解答用的就是你知道的方法 只是你看不出来，大概是不熟悉矩阵运算跟一般算式间的对应 假设A = [v1 v2 v3] 是个矩阵，则 A * [a b c]^t = a*v1 + b*v2 + c*v3 就是说，A乘上一个向量，其实就是以那个向量作系数，产生一组A的行向量的线性组合 回到上面的题目，左边是原解答，右边是用矩阵列式 接续算完kernel space後 | A = PD(P^-1)，其中P = [v1 v2 v3] | 因为 B = 2(v1) + (v2) + 3(v3) | 因为 (P^-1)B = [2 1 3]^t | 所以 B = P [2 1 3]^t | (A^100)B = (A^100)(2v1 + v2 + 3v3) | (A^100)B = (A^100)P [2 1 3]^t = 2(A^100)v1 + (A^100)v2 + 3(A^100)v3 | = PD^100(P^-1) P [2 1 3]^t = 2(1^100)v1 + (1^100)v2 + 3(3^100)v3 | = PD^100 [2 1 3]^t 到这步可以看一下，因为D^100 = [1^100 0 0 ] [0 1^100 0 ] [0 0 3^100] 右边的算式把D^100乘到後面的向量就是 P [2(1^100) (1^100) 3(3^100)]^t 再套用上面说的矩阵乘向量展开法，跟左式根本是一样的 回到你上面的问题: 1. 因为这其实就是对角化计算背後的观念 利用Ax = λx，把(A^100)B 简化成 P(λ^100)B' 只是照观念一步一步做和单纯矩阵运算的差别 2. 其实算出[2 1 3]^t，就是把P(D^100)(P^-1)B中的(P^-1)B乘开得来的 3. 当然不一定要用先算出P^-1，然後再乘上B的计算方式 可以用 [P|B] 运算到 [I|(P^-1)B] 或设x = [x1 x2 x3]^t，解[v1 v2 v3]x = B，就是 -x2 + x3 = 2 2x2 = 2 x1 + 2x3 = 8 这题很容易就能解出x = [2 1 3]^t (因为求得的eigenvector里面0很多) 所以会有解答的算法比较快的感觉，其实都是一样的作法 就看你觉得哪个计算方式你比较习惯就好 --

※ 发信站: 批踢踢实业坊(ptt.cc) ◆ From: 218.166.98.161 ※ 编辑: ssccg 来自: 218.166.98.161 (08/30 00:47) ※ 编辑: ssccg 来自: 218.166.98.161 (08/30 15:31)