※ 引述《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)