Re: [工數] 為什麼重根時要多乘一個t
※ 引述《eolith123 (deer)》之銘言:
: (4) 2 " 4
: y - 2k y + k y = 0 where k is constant
: at
: Let y = e
: 4 2 2 4
: => a - 2k a + k = 0, a = ±k (k,k,-k,-k)
: kt kt -kt -kt
: y = Ae + Bte + Ce + Dte
: 請問為什麼在重根時需要再乘上一個t?
可以用線性方程組的觀點來看
x1 = y [x1]
' , x =[x2]
x2 = y [x3]
'' [x4]
x3 = y
(3)
x4 = y
(4) 2 '' 4 (4) 2 4
因為 y - 2k y + k y = 0 <=> y = 2k y'' -k y
[ 0 1 0 0 ]
=> x' = [ 0 0 1 0 ] x
[ 0 0 0 1 ]
4 2
[-k 0 2k 0 ]
令A = [ 0 1 0 0 ]
[ 0 0 1 0 ]
[ 0 0 0 1 ]
4 2
[-k 0 2k 0 ]
對A對角化
4 2 2 4 2 2 2 2 2
det(A-cI) = 0 => c -2k c + k = (c - k ) = (c+k) (c-k) = 0
=> c = -k, -k, k, k
[ 1 ] [1 ]
(A+kI)v = 0 => v = a[-k ] , (A-kI)v = 0 => v = b [k ]
2 2
[ k ] [k ]
3 3
[-k ] [k ]
所以兩個eigenvalue的eigenspace dimension都小於重根數 => 不可對角化
但如果找A的Jordan form,則中間的那個J matrix應該是像
[k 1 0 0]
J = [0 k 0 0]
[0 0 -k 1]
[0 0 0 -k]
所以原式變為
-1 -1 -1
x' = QJQ x <=> Q x' = JQ x
-1 -1 '
令 z = Q x => z' = Q x
[z1]
=> z' = Jz , z = [z2]
[z3]
[z4]
=> z1' = kz1 + z2
z2' = kz2
z3' = -kz3 + z4
z4' = -kz4
kt kt kt -kt -kt -kt
=> z2 = c1e , z1 = c1 te + c2 e , z4 = c3 e , z3 = c3 te + c4 e
所以
kt kt
[c1 te + c2 e ]
x = Q kt
[ c1 e ]
-kt -kt
[c3 te + c4 e ]
-kt
[ c3 e ]
所以
kt kt kt -kt -kt -kt
y = x1 =q11(c1 te + c2 e ) + q12(c1 e ) + q13(c3 te + c4 e )+q14(c3 e )
kt kt -kt -kt
= Ae + B te + C e + D te
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 220.129.27.62
※ 文章網址: http://www.ptt.cc/bbs/Math/M.1398192056.A.8CD.html
※ 編輯: yueayase (220.129.27.62), 04/23/2014 02:42:09
推
04/23 03:40, , 1F
04/23 03:40, 1F
→
04/23 06:09, , 2F
04/23 06:09, 2F
→
04/24 00:46, , 3F
04/24 00:46, 3F
→
04/24 00:46, , 4F
04/24 00:46, 4F
→
04/24 00:50, , 5F
04/24 00:50, 5F
→
04/24 00:50, , 6F
04/24 00:50, 6F
→
04/24 00:50, , 7F
04/24 00:50, 7F
推
04/24 02:17, , 8F
04/24 02:17, 8F
推
04/24 02:19, , 9F
04/24 02:19, 9F
推
04/24 15:13, , 10F
04/24 15:13, 10F
→
01/02 15:44,
5年前
, 11F
01/02 15:44, 11F
→
07/07 12:03,
4年前
, 12F
07/07 12:03, 12F
討論串 (同標題文章)
本文引述了以下文章的的內容:
完整討論串 (本文為第 2 之 4 篇):