→ a149851571:事實上這題求K的反矩陣很容易 再乘以b就是x了 02/02 12:04
→ a149851571:抱歉 我沒看清楚內文就發文 02/02 12:05
→ a149851571:L^t根L應該是反矩陣 所以只要求D的行列式值就可以 02/02 12:06
→ gnsh:不是要正交才會是L^t=L^-1嗎? 02/02 12:10
推 blazesunny:|K|=|L|*|D|*|L^T| [x]=inv([K])[b] 求inv(L)時很煩XD 02/02 12:13
→ blazesunny:這題L^t≠L^-1唷 02/02 12:14
→ gnsh:我覺得算L^-1沒比直接乘開快耶XDDD 剛試了一下 02/02 12:23
→ gnsh:所以求︱K︱沒有類似我打的那種性質歐QQ 02/02 12:24
→ blazesunny:因為inv([L])太煩 求[x]還是暴力乘開用高斯比較快(吧?) 02/02 12:24
→ blazesunny:有|L|=|L^t|的性質 而|L|的話對第一列做降階很快~ 02/02 12:26
→ gnsh:摁摁 依據我算的 求出K後得x不慢 只是很少再算矩陣XD 好容易 02/02 12:29
→ gnsh:眼花XDD 02/02 12:29
→ blazesunny:恩XD 02/02 12:31
推 squallting:這種題目我肯定不只有單純爆力解開法XD 02/02 12:32
推 jack0711:det(L)=1對第一行降階 02/02 12:32
→ jack0711:det(D)=0.5^3 02/02 12:33
→ jack0711:det(K)=det(D) 02/02 12:33
推 harrypotter2:LDU分解，Kx=b→x=E10..E7E6*D^(-1)*E5E4E3E2E1*b 02/02 12:35
→ harrypotter2:E是基本矩陣，你只要對b做列運算在做個D(-1)就好了 02/02 12:36
推 harrypotter2:修正一下，x=E8..E6E5*D^(-1)*E4E3E2E1*b 02/02 12:39
→ gnsh:不好意思XD沒見過這東西 我連基本矩陣都不知道耶.. 02/02 12:40
推 harrypotter2:LU分解那邊，E2E1A=U→A=E1^(-1)E2^(-1)U=LU 02/02 12:44
→ gnsh:歐歐 我知道了 謝謝XD 02/02 12:45
推 SS327:H大不影電影了阿 02/02 12:49
推 harrypotter2:畢業了改走電子業... 02/02 12:51
推 blazesunny:是 [E][L]=[I] [L]=inv[E] inv[L]=[E] inv[L^t]=[E^t] 02/02 12:54
→ blazesunny:這樣嗎? 02/02 12:55
推 pikachu123:這題應該是考LU分解 你把D乘給L^t K=LU 02/02 13:35
→ pikachu123:因為LDL^t是LU過來的 你令LUx=b 02/02 13:36
→ pikachu123:Ly=b 其中y=Ux 在解一次y=Ux 那個蠻好算的 02/02 13:37
→ gnsh:摁摁QQ 感覺跟直接乘開是差不多的意思XD 02/02 13:42
推 pikachu123:你用LU解linear systme本來就沒有比較快 02/02 13:43
→ pikachu123:那個寫程式比較好寫 這題它都做好LU給你了 02/02 13:44
→ pikachu123:就用LU解比較快一點點就是(就那麼一點點) 02/02 13:44
→ blazesunny:原來如此 謝謝:) 02/02 14:08