24小时热门版块排行榜    

查看: 7627  |  回复: 42

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

你是就这一个矩阵么?数据传上来,我可以帮你试算一下。

» 本帖已获得的红花(最新10朵)

11楼2013-07-11 10:25:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tgvictor99

木虫 (知名作家)

迷糊领域荣誉领主

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
10楼: Originally posted by hk_kong36 at 2013-07-11 09:33:57
是576*576的矩阵,可逆,确实是在解方程组,矩阵元素都特别小,也不对称

求逆的过程会不会很复杂。。。。。...

针对楼主給的情况。我想了一个方法,没测试过,楼主可以试一试。
推论过程:
因为楼主说这个矩阵里的元素都在1e-4的量级,这就导致矩阵的行列式会非常小。所以做处理,假说这个矩阵为A,逆为Ani,令A=1e-4*Abar
因为    A*Ani=I ,那么 1e-4*Abar*Ani=I. 那么Abar*(1e-4*Ani) = I;
所以 inv(Abar) =1e-4*Ani
所以 Ani= 1e4*inv(Abar)
所以Ani =1e4*inv(1e4*A)
这样转换矩阵1e4*A的行列式值就不会接近于0了,matlabb应该就能给出结果了。

» 本帖已获得的红花(最新10朵)

每个人的人生都是时代的产物。每个人的血统、命运和梦想组成了现在这个时代。
12楼2013-07-11 10:51:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hk_kong36

金虫 (著名写手)

送红花一朵
引用回帖:
12楼: Originally posted by tgvictor99 at 2013-07-11 10:51:28
针对楼主給的情况。我想了一个方法,没测试过,楼主可以试一试。
推论过程:
因为楼主说这个矩阵里的元素都在1e-4的量级,这就导致矩阵的行列式会非常小。所以做处理,假说这个矩阵为A,逆为Ani,令A=1e-4*Abar
...

谢谢  我试了下还是一样的提示。。。
thanks,聆听花开的声音,需要的是沉淀
13楼2013-07-11 14:35:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hk_kong36

金虫 (著名写手)

送红花一朵
引用回帖:
11楼: Originally posted by pippi6 at 2013-07-11 10:25:09
你是就这一个矩阵么?数据传上来,我可以帮你试算一下。

576*576的矩阵,你帮我看看  O(∩_∩)O谢谢
thanks,聆听花开的声音,需要的是沉淀
14楼2013-07-11 14:38:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

引用回帖:
14楼: Originally posted by hk_kong36 at 2013-07-11 14:38:21
576*576的矩阵,你帮我看看  O(∩_∩)O谢谢...

矩阵的数据呢?
15楼2013-07-11 14:41:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hk_kong36

金虫 (著名写手)

送红花一朵
引用回帖:
2楼: Originally posted by pippi6 at 2013-07-10 11:50:47
是稀疏矩阵吗?阶数太高,用直接法有点问题,试着用迭代法。还有就是,非要matlab做吗?

不是稀疏矩阵,不用matlab的话会比较好算吗,如果这样也可以不用的
thanks,聆听花开的声音,需要的是沉淀
16楼2013-07-11 14:45:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

引用回帖:
16楼: Originally posted by hk_kong36 at 2013-07-11 14:45:17
不是稀疏矩阵,不用matlab的话会比较好算吗,如果这样也可以不用的...

这都没关系,你就把数据传上来先试试。最好双精度,请告知格式。
17楼2013-07-11 14:49:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

arieskaka

新虫 (初入文坛)

我记得有个数值计算里用的伪逆矩阵可以试试……
18楼2013-07-11 15:15:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tgvictor99

木虫 (知名作家)

迷糊领域荣誉领主

引用回帖:
13楼: Originally posted by hk_kong36 at 2013-07-11 14:35:17
谢谢  我试了下还是一样的提示。。。...

奥,能求出这个矩阵的行列式么。matlab里面好像大于1e300的数就定义为inf就是无穷大了,而小雨1e-300的数就定义为0了。576个1e-4 相乘,结果肯定到了1e-2000了吧,就被自动定义为零了。matalb求逆还是有局限的。没数据部好具体处理啊。

» 本帖已获得的红花(最新10朵)

每个人的人生都是时代的产物。每个人的血统、命运和梦想组成了现在这个时代。
19楼2013-07-11 16:27:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hk_kong36

金虫 (著名写手)

送红花一朵
引用回帖:
19楼: Originally posted by tgvictor99 at 2013-07-11 16:27:27
奥,能求出这个矩阵的行列式么。matlab里面好像大于1e300的数就定义为inf就是无穷大了,而小雨1e-300的数就定义为0了。576个1e-4 相乘,结果肯定到了1e-2000了吧,就被自动定义为零了。matalb求逆还是有局限的。没 ...

算过了  0,所以求逆行不通,我想可能需要别的方法求解方程组。。
thanks,聆听花开的声音,需要的是沉淀
20楼2013-07-11 16:43:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hk_kong36 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 6/300 2026-08-25 00:11 by BZKMTicpDhFj
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +7 scalable 2026-08-24 7/350 2026-08-24 23:01 by anata1209
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] 我面上完蛋了 +9 且听虎啸 2026-08-20 10/500 2026-08-24 20:08 by 苦寒来香
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见