| 查看: 2622 | 回复: 6 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
学术老混混金虫 (小有名气)
|
[求助]
【求助】如何判断大型矩阵是否奇异？
|
||
|
本人在求解非线性方程组（结构动力学问题）的时候，参考Newton相关算法中有一步是判断Jacobi矩阵是否奇异，但是没给出具体的实现。本问题的Jacobi矩阵维数较大（1000×1000），本人采用VC++代码实现，请问有什么较为简便的实现方法？ 用求行列式的方法对于本问题的大型矩阵好像因为数值误差的存在导致即使矩阵奇异，行列式结果也有个较大值。请问是否能够通过对结果的数量级进行大致的判断，是否靠谱？另lapack包里有没有直接求解矩阵行列式的函数，函数名是啥？ 或者有什么别的方法，比如根据特征值判断之类的，总之目的就一个：判断该类大型矩阵的奇异性。最好能给出C++代码的实现，或者能否把Matlab相关函数移植过来也行，呵呵 |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
面上合作单位盖章
已经有5人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有8人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
申请删除本帖
已经有7人回复
-
基金不中，共勉
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有11人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教：矩阵奇异值分解的运算量如何计算？
已经有4人回复
- COMSOL AC/DC模拟出错，显示出现奇异矩阵
已经有3人回复
- 求解奇异矩阵的问题
已经有6人回复
- 求逆矩阵问题、
已经有4人回复
- 近似奇异矩阵 怎么处理
已经有5人回复
- COMSOL 求解出现奇异矩阵
已经有4人回复
- 矩阵相乘，逆向求解，这种方法可行吗？
已经有8人回复
- 根据《雅典宪章》的四个准则：居住、工作、游憩和交通，设一个判断矩阵
已经有3人回复
- 如何得到一个奇异矩阵的线性无关的行号
已经有3人回复
- 用java编程求矩阵的特征值与奇异值
已经有3人回复
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
学术老混混
金虫 (小有名气)
- 应助: 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
学术老混混6楼2013-05-03 08:51:59
10