Matlab 版 (精华区)

发信人: zjliu (秋天的萝卜), 信区: Matlab
标  题: [合集]我的矩阵比较大，怎样求逆？
发信站: 哈工大紫丁香 (2003年08月03日09:40:47 星期天), 站内信件

发信人: zhouxytj (yiyi), 信区: Matlab
标  题: 我的矩阵比较大，怎样求逆？
发信站: 哈工大紫丁香 (Mon Jul 28 09:29:47 2003)

   有一27000*27000的矩阵K，求逆时说内存不够，怎么办？
   我用过:
   Tmp=diag(sparse(ones(size(KK,1),1)));
   ff=KK\Tmp;
   或
   ff=inv(KK);
   都说内存不够，有没有好一些的方法，谢谢！
   

--

※ 来源:．哈工大紫丁香 bbs.hit.edu.cn [FROM: 202.120.182.90]
发信人: nosay (☆纯中纯冰红茶⊙倾派☆), 信区: Matlab
标  题: Re: 我的矩阵比较大，怎样求逆？
发信站: 哈工大紫丁香 (2003年07月28日14:34:11 星期一), 站内信件

我最多只求过27*27的
【 在 zhouxytj (yiyi) 的大作中提到: 】
:    有一27000*27000的矩阵K，求逆时说内存不够，怎么办？
:    我用过:
:    Tmp=diag(sparse(ones(size(KK,1),1)));
:    ff=KK\Tmp;
:    或
:    ff=inv(KK);
:    都说内存不够，有没有好一些的方法，谢谢！
:    


--
 精灵族头号强援 水族最佳同盟

欢迎大家去IQrace版
                版主pujialun&ababyfish欢迎您 
                                      名誉版主：lyjun  
 nosay (☆纯中纯冰红茶⊙倾派☆) 

※ 来源:·哈工大紫丁香 bbs.hit.edu.cn·[FROM: 202.118.228.235]
发信人: Vaccin (Rattle), 信区: Matlab
标  题: Re: 我的矩阵比较大，怎样求逆？
发信站: 哈工大紫丁香 (2003年07月28日15:13:46 星期一), 站内信件

到其他机器上找几条内存条插上

【 在 zhouxytj (yiyi) 的大作中提到: 】
:    有一27000*27000的矩阵K，求逆时说内存不够，怎么办？
:    我用过:
:    Tmp=diag(sparse(ones(size(KK,1),1)));
:    ff=KK\Tmp;
:    或
:    ff=inv(KK);
:    都说内存不够，有没有好一些的方法，谢谢！
:    

--
※ 来源:·哈工大紫丁香 bbs.hit.edu.cn·[FROM: 210.46.77.34]
发信人: aiolos (windforce), 信区: Matlab
标  题: Re: 我的矩阵比较大，怎样求逆？
发信站: 哈工大紫丁香 (2003年07月28日18:08:58 星期一), 站内信件

用最小二乘法逼近，有公式
还有许多逼近的方法，我就知道最小二乘法
但是上次作的时候，还是没有调出来程序，呵呵
能力有限~！~！

【 在 zhouxytj (yiyi) 的大作中提到: 】
:    有一27000*27000的矩阵K，求逆时说内存不够，怎么办？
:    我用过:
:    Tmp=diag(sparse(ones(size(KK,1),1)));
:    ff=KK\Tmp;
:    或
:    ff=inv(KK);
:    都说内存不够，有没有好一些的方法，谢谢！
:    


--
山中的小猫                                                      

※ 来源:·哈工大紫丁香 bbs.hit.edu.cn·[FROM: 202.118.228.179]
发信人: Mathsims (玉子), 信区: Matlab
标  题: Re: 我的矩阵比较大，怎样求逆？
发信站: 哈工大紫丁香 (Mon Jul 28 19:25:17 2003)

这么大规模的矩阵求逆，不光内存不够，数值计算的性能（误差、稳定性等等）也很是问
题。建议先降低它的规模再求其逆。
至于降低规模的方法，用马克思主义的话说，叫做“具体问题具体分析”，恐怕只有你自
己去找了。
【 在 zhouxytj 的大作中提到: 】
:    有一27000*27000的矩阵K，求逆时说内存不够，怎么办？
:    我用过:
:    Tmp=diag(sparse(ones(size(KK,1),1)));
:    ff=KK\Tmp;
:    或
:    ff=inv(KK);
:    都说内存不够，有没有好一些的方法，谢谢！
:    

--
不是什么高山，
也不是什么星斗，
我只是一块有信念的石头。

小样儿，敢惹我，看我砸不死你！
※ 来源:．哈工大紫丁香 bbs.hit.edu.cn [FROM: 202.118.226.28]
发信人: ljfcy (毛毛), 信区: Matlab
标  题: Re: 我的矩阵比较大，怎样求逆？
发信站: 哈工大紫丁香 (2003年07月28日20:34:50 星期一), 站内信件

你对这么大的矩阵求逆,就从工程的角度讲都没有任何意义,

【 在 zhouxytj (yiyi) 的大作中提到: 】
:    有一27000*27000的矩阵K，求逆时说内存不够，怎么办？
:    我用过:
:    Tmp=diag(sparse(ones(size(KK,1),1)));
:    ff=KK\Tmp;
:    或
:    ff=inv(KK);
:    都说内存不够，有没有好一些的方法，谢谢！
:    

--
※ 来源:·哈工大紫丁香 bbs.hit.edu.cn·[FROM: 202.118.228.210]
发信人: yrabbits (yrabbits), 信区: Matlab
标  题: Re: 我的矩阵比较大，怎样求逆？
发信站: 哈工大紫丁香 (2003年07月28日22:54:06 星期一), 站内信件

我觉得应该根据你得专业知识，将问题简化，如果这样得矩阵求逆，
对于实际问题来说是没有很多意义的
一定应该可以解耦的。
【 在 zhouxytj (yiyi) 的大作中提到: 】
   有一27000*27000的矩阵K，求逆时说内存不够，怎么办？
   我用过:
   Tmp=diag(sparse(ones(size(KK,1),1)));
   ff=KK\Tmp;
   或
   ff=inv(KK);
   都说内存不够，有没有好一些的方法，谢谢！
   

--

※ 来源:．哈工大紫丁香 bbs.hit.edu.cn [FROM: 202.120.182.90]

--
※ 来源:·哈工大紫丁香 bbs.hit.edu.cn·[FROM: 202.119.14.197]
发信人: able (随心所欲), 信区: Matlab
标  题: Re: 我的矩阵比较大，怎样求逆？
发信站: 哈工大紫丁香 (2003年07月29日14:41:45 星期二), 站内信件


把你的矩阵转化为稀疏矩阵，然后然后进行计算就可以乐。一个对角阵可能有几百K的空间
就可以乐
【 在 zhouxytj (yiyi) 的大作中提到: 】
:    有一27000*27000的矩阵K，求逆时说内存不够，怎么办？
:    我用过:
:    Tmp=diag(sparse(ones(size(KK,1),1)));
:    ff=KK\Tmp;
:    或
:    ff=inv(KK);
:    都说内存不够，有没有好一些的方法，谢谢！
:    

--
※ 来源:·哈工大紫丁香 bbs.hit.edu.cn·[FROM: 202.118.237.10]
发信人: arc (ARC), 信区: Matlab
标  题: Re: 我的矩阵比较大，怎样求逆？
发信站: 哈工大紫丁香 (2003年07月30日18:03:13 星期三), 站内信件


一个原因就是无法通过硬件平台实现，即使实现了运算时间也无法忍受

【 在 ljfcy (毛毛) 的大作中提到: 】
: 你对这么大的矩阵求逆,就从工程的角度讲都没有任何意义,
: 【 在 zhouxytj (yiyi) 的大作中提到: 】
: :    有一27000*27000的矩阵K，求逆时说内存不够，怎么办？
: :    我用过:
: :    Tmp=diag(sparse(ones(size(KK,1),1)));
: :    ff=KK\Tmp;
: :    或
: :    ff=inv(KK);
: :    都说内存不够，有没有好一些的方法，谢谢！

--
※ 来源:·哈工大紫丁香 bbs.hit.edu.cn·[FROM: 202.118.228.214]