作者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