| 查看: 2664 | 回复: 6 | ||
学术老混混金虫 (小有名气)
|
[求助]
【求助】如何判断大型矩阵是否奇异？
|
|
本人在求解非线性方程组（结构动力学问题）的时候，参考Newton相关算法中有一步是判断Jacobi矩阵是否奇异，但是没给出具体的实现。本问题的Jacobi矩阵维数较大（1000×1000），本人采用VC++代码实现，请问有什么较为简便的实现方法？ 用求行列式的方法对于本问题的大型矩阵好像因为数值误差的存在导致即使矩阵奇异，行列式结果也有个较大值。请问是否能够通过对结果的数量级进行大致的判断，是否靠谱？另lapack包里有没有直接求解矩阵行列式的函数，函数名是啥？ 或者有什么别的方法，比如根据特征值判断之类的，总之目的就一个：判断该类大型矩阵的奇异性。最好能给出C++代码的实现，或者能否把Matlab相关函数移植过来也行，呵呵 |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教：矩阵奇异值分解的运算量如何计算？
已经有4人回复
- COMSOL AC/DC模拟出错，显示出现奇异矩阵
已经有3人回复
- 求解奇异矩阵的问题
已经有6人回复
- 求逆矩阵问题、
已经有4人回复
- 近似奇异矩阵 怎么处理
已经有5人回复
- COMSOL 求解出现奇异矩阵
已经有4人回复
- 矩阵相乘，逆向求解，这种方法可行吗？
已经有8人回复
- 根据《雅典宪章》的四个准则：居住、工作、游憩和交通，设一个判断矩阵
已经有3人回复
- 如何得到一个奇异矩阵的线性无关的行号
已经有3人回复
- 用java编程求矩阵的特征值与奇异值
已经有3人回复
学术老混混
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3055.2
- 红花: 2
- 帖子: 144
- 在线: 58.8小时
- 虫号: 2263658
- 注册: 2013-01-27
- 性别: GG
- 专业: 机械结构强度学
2楼2013-04-30 22:59:51
wshaoxin
铁杆木虫 (正式写手)
- 应助: 47 (小学生)
- 金币: 5678.9
- 散金: 1839
- 红花: 14
- 帖子: 929
- 在线: 154.1小时
- 虫号: 1785567
- 注册: 2012-04-29
- 专业: 决策理论与方法
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
学术老混混: 金币+50, ★有帮助, 反幂法的代码你有没有啊？最好是C++的，matlab也行，多谢哈。有的话剩下的50金币也给你 2013-05-02 12:18:42
感谢参与，应助指数 +1
学术老混混: 金币+50, ★有帮助, 反幂法的代码你有没有啊？最好是C++的，matlab也行，多谢哈。有的话剩下的50金币也给你 2013-05-02 12:18:42
|
1000*1000的矩阵用matlab可以求吧，eig（矩阵）就能求出特征根，或者是反幂法迭带求最小特征值，一般的数值分析书上应该都有。当然如果是超高维（＞10000)的矩阵，应该另有它法。 [ 发自手机版 http://muchong.com/3g ] |
3楼2013-05-02 10:53:42
| 顶 |
4楼2013-05-02 12:32:04
wshaoxin
铁杆木虫 (正式写手)
- 应助: 47 (小学生)
- 金币: 5678.9
- 散金: 1839
- 红花: 14
- 帖子: 929
- 在线: 154.1小时
- 虫号: 1785567
- 注册: 2012-04-29
- 专业: 决策理论与方法
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
学术老混混: 金币+50, ★有帮助, 好的，多谢 2013-05-02 17:17:47
学术老混混: 金币+50, ★有帮助, 好的，多谢 2013-05-02 17:17:47
| 我从《数值分析》：何汉林 梅家斌，科学出版社。找了这个，你试试吧？ |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : Matlab程序.docx
2013-05-02 12:50:54, 12.27 K
5楼2013-05-02 12:51:28
学术老混混6楼2013-05-03 08:51:59
学术老混混
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3055.2
- 红花: 2
- 帖子: 144
- 在线: 58.8小时
- 虫号: 2263658
- 注册: 2013-01-27
- 性别: GG
- 专业: 机械结构强度学
7楼2013-05-03 15:34:10
10