作者aflilfesy (...)
看板NCTU-STAT95G
標題Re: [問題] 統研所網頁@@
時間Sun Jul 30 23:33:51 2006
※ 引述《Oyin (Lazy)》之銘言:
: HW04
: (4)
: 將 Gauss 消去法, Gauss -Jordan消去法, LU分解, Jacobi Mathod, 利用指標寫成副程式和標頭檔.
: */
: 寫出來可以教教我嗎@@
提供我的想法:
Gauss消去法:
老師已經有code了 但是 那是二維的.. 想辦法用指標寫成一維
幾乎是把[i][j] 改成 [i*n+j] 想一下吧...
Gauss -Jordan消去法:
老師有完整的code,但那也是二維,必須用指標改寫成一維,
很不幸地,我看不是很懂他寫什麼>"<
我的做法是...
把Gauss消去法的"消去項步驟(往下減)"搞清楚之後,
再加上 把他寫成往上減,就可以變成對角矩陣的型態囉..
就差不多了.. 再想辦法變成單位矩陣~
LU分解:
在做Gauss消去法時,簡化出來的矩陣,其實就是U(上三角矩陣)
而下三角該如何寫呢??
關鍵應該就在...倍數
例如:第一列*a倍去減第二列, 這a其實就是L(下三角)的某一元素
(我是在翻舊課本時,才想起來啥米是LU分解@@ )
恩~~~~~但這邊我覺得有點怪...
因為我看書的定義,U(上三角矩陣)是沒要求對角項為1,而是L(下三角)對角項是1.
所以我也不知該怎麼辦嚕~...
Jacobia Method:
恩~~ 如果你在google 打看看.. jacobia C 你可以看到code
很不幸地,我又是有看沒有懂,
所以我是照著老師網頁的演算法,寫的哩哩辣辣~!!
感覺還頗粗糙~
後來翻翻我大學的數值分析課本,要用Jacobia Method 其實條件還頗多的!!
難怪我常常輸入方程式 電腦都跑到當~!
痾~~~
剛剛看到前面的文章
請問這次的作業
是不是 不需要改寫成一維指標阿??
--
有幫助到你嘛?? ^^
希望有幫助..
我只能分享我的想法了~~
~ 瞎++ 寫到瞎 ~ =.=a
: ※ 引述《mythology612 (如果能不改變就好了!)》之銘言:
: 是不是線上上不去阿>"<
: 嗚嗚~難得今天要好好努力寫作業的QQ
: 有誰有題目可以PO出來可以我看呢~
: 謝啦^^
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.113.90.16
※ 編輯: aflilfesy 來自: 140.113.90.16 (07/30 23:38)
1F:推 Oyin:^^恩恩~~謝謝你兒~~真厲害~0~ 我幾乎都不會寫@@ 07/31 22:34
2F:推 Y0SHIKI:我是改成一維指標 不過寫一題之後就不想寫了 囧> 08/03 03:05