Re: [問題]要如何將行列式使用拉普拉斯展開?

看板Mathematica作者 (養花種魚數月亮賞星星)時間12年前 (2012/12/12 00:22), 編輯推噓4(402)
留言6則, 4人參與, 最新討論串2/5 (看更多)
※ 引述《candy88257 (阿泰斯)》之銘言: : 在此說的拉普拉斯是指: http://tinyurl.com/bhhlt4q : 上面網址是維基百科 : 有辦法做到這點嗎? : 想把10幾階的矩陣使用拉普拉斯展開,不想用Det算 : 因為用Det怪怪的,行列式值會有分母....,明明所有數值都沒分數的說... : 求高手幫忙!! 感謝!!! temp=12; A:=N@Table[Sin[m+2n]+Cos[m+2n]z,{m,1,temp},{n,1,temp}]; myminors[mymatrix_?MatrixQ,i_,j_]:=Block[{dim}, dim=Length@mymatrix; mymatrix[[All,Drop[Range[dim],{j}]]][[Drop[Range[dim],{i}]]]] mydet[mymatrix_?MatrixQ]:= Total[Table[Det@myminors[mymatrix,1,i]*mymatrix[[1,i]]*If[OddQ[i],1,-1], {i,Length@A}]] {temp=#;Det[A/.z->1]-mydet[A/.z->1]}&/@Range[2,12] mydet[A]//Simplify -- 養花種魚數月亮賞星星 http://chungyuandye.twbbs.org -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 106.1.135.51

12/12 01:47, , 1F
感謝!!!!!! 真是非常感謝您!
12/12 01:47, 1F

12/13 09:18, , 2F
您好,這樣做的話似乎只能用在12的情況?更高階便必須做
12/13 09:18, 2F

12/13 09:20, , 3F
更多的展開@@...共會有n!/11!項ORZ
12/13 09:20, 3F

12/13 10:37, , 4F
是的,更高階的話要將每個minor繼續分解
12/13 10:37, 4F

12/13 17:00, , 5F
但是怎會有誤差@@? 不曉得到底哪個是對的...
12/13 17:00, 5F

12/13 18:47, , 6F
矩陣的元如果不是純實數而帶未知數也不知該怎麼pivot
12/13 18:47, 6F
文章代碼(AID): #1GnrtXwq (Mathematica)
討論串 (同標題文章)
文章代碼(AID): #1GnrtXwq (Mathematica)