【答案】应助回帖

» 本帖已获得的红花（最新10朵）

hk_kong36

11楼2013-07-11 10:25:09

tgvictor99

木虫 (知名作家)

迷糊领域荣誉领主

应助: 40 (小学生)
金币: 2667.4
散金: 1907
红花: 49
沙发: 1
帖子: 5238
在线: 923.8小时
虫号: 2524003
注册: 2013-06-27
专业: 流体力学

【答案】应助回帖

感谢参与，应助指数 +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朵）

hk_kong36

每个人的人生都是时代的产物。每个人的血统、命运和梦想组成了现在这个时代。

12楼2013-07-11 10:51:28

hk_kong36

金虫 (著名写手)

应助: 0 (幼儿园)
金币: 681.7
散金: 1000
红花: 7
帖子: 1066
在线: 262.9小时
虫号: 964842
注册: 2010-03-08
性别: MM
专业: 磁共振成像技术与造影剂

送红花一朵

引用回帖:

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

hk_kong36

金虫 (著名写手)

应助: 0 (幼儿园)
金币: 681.7
散金: 1000
红花: 7
帖子: 1066
在线: 262.9小时
虫号: 964842
注册: 2010-03-08
性别: MM
专业: 磁共振成像技术与造影剂

送红花一朵

引用回帖:

11楼: Originally posted by pippi6 at 2013-07-11 10:25:09
你是就这一个矩阵么？数据传上来，我可以帮你试算一下。

576*576的矩阵，你帮我看看 O(∩_∩)O谢谢

thanks，聆听花开的声音，需要的是沉淀

14楼2013-07-11 14:38:21

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

应助: 413 (硕士)
贵宾: 0.002
金币: 7116.5
散金: 15
红花: 63
帖子: 1639
在线: 798.9小时
虫号: 2469437
注册: 2013-05-14
专业: 计算数学与科学工程计算

引用回帖:

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

矩阵的数据呢？

15楼2013-07-11 14:41:30

hk_kong36

金虫 (著名写手)

应助: 0 (幼儿园)
金币: 681.7
散金: 1000
红花: 7
帖子: 1066
在线: 262.9小时
虫号: 964842
注册: 2010-03-08
性别: MM
专业: 磁共振成像技术与造影剂

送红花一朵

引用回帖:

2楼: Originally posted by pippi6 at 2013-07-10 11:50:47
是稀疏矩阵吗？阶数太高，用直接法有点问题，试着用迭代法。还有就是，非要matlab做吗？

不是稀疏矩阵，不用matlab的话会比较好算吗，如果这样也可以不用的

thanks，聆听花开的声音，需要的是沉淀

16楼2013-07-11 14:45:17

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

应助: 413 (硕士)
贵宾: 0.002
金币: 7116.5
散金: 15
红花: 63
帖子: 1639
在线: 798.9小时
虫号: 2469437
注册: 2013-05-14
专业: 计算数学与科学工程计算

【答案】应助回帖

引用回帖:

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

这都没关系，你就把数据传上来先试试。最好双精度，请告知格式。

17楼2013-07-11 14:49:37

arieskaka

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1342.7
红花: 1
帖子: 22
在线: 51.1小时
虫号: 2488500
注册: 2013-05-30
性别: GG
专业: 电化学

我记得有个数值计算里用的伪逆矩阵可以试试……

18楼2013-07-11 15:15:25

tgvictor99

木虫 (知名作家)

迷糊领域荣誉领主

应助: 40 (小学生)
金币: 2667.4
散金: 1907
红花: 49
沙发: 1
帖子: 5238
在线: 923.8小时
虫号: 2524003
注册: 2013-06-27
专业: 流体力学

引用回帖:

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朵）

hk_kong36

每个人的人生都是时代的产物。每个人的血统、命运和梦想组成了现在这个时代。

19楼2013-07-11 16:27:27

hk_kong36

金虫 (著名写手)

应助: 0 (幼儿园)
金币: 681.7
散金: 1000
红花: 7
帖子: 1066
在线: 262.9小时
虫号: 964842
注册: 2010-03-08
性别: MM
专业: 磁共振成像技术与造影剂

送红花一朵

引用回帖:

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