如何快速求一个5000*5000的对称矩阵的逆 - 未名空间MITBBS历史存档

国际科技财经博客移民网络热点娱乐民生时事公众号

Redian新闻

>未名空间

>Computation - 科学计算

如何快速求一个5000*5000的对称矩阵的逆

如何快速求一个5000*5000的对称矩阵的逆# Computation - 科学计算

m*p2008-05-14 07:05

1 楼

【以下文字转载自 Hardware 讨论区】
发信人: munp (毫微纳皮), 信区: Hardware
标题: ibm本本的无线网卡问题
发信站: BBS 未名空间站 (Fri Mar 28 19:30:43 2008)
ibm thinkpad x40,国内的时候买的,
无线网卡intel minipci 2200bg坏了,原因不明,
打算换上一块从dell laptop上面拆下来的minipci无线网卡,也是intel的,
结果bios自检都没通过,强制通过,进入winxp,结果ibm access connections程序不认,
在网上查了一下,ibm thinkpad做了手脚不让用户随便换无线网卡,
网上倒是有个什么工具可以解除这一限制,没敢用,好像涉及到重新刷新bios,
后来按照着原装的那个minipci无线网卡型号,在网上买了一个同型号的,
装上后bios自检通过,ibm access connections也能识别这个同型号的网卡,上网没问题,
可是新问题来了,
机器无法进入休眠状态,也无法进入待机状态,
无论是合上显示屏,还是通过ibm组合键,
造成的结

t*n2008-05-14 07:05

2 楼

我需要解一个线性方程 A*X = D, 其中A 是一个 5000×5000 的对称矩阵,
X = Inv(A)*D.
我试过用 LU 的方法，不过速度很慢。不只各位能不能提供一些比较快速算法。
多谢了

p*e2008-05-14 07:05

3 楼

5000*5000不大啊，LU还慢？

【在 t**********n 的大作中提到】

: 我需要解一个线性方程 A*X = D, 其中A 是一个 5000×5000 的对称矩阵,
: X = Inv(A)*D.
: 我试过用 LU 的方法，不过速度很慢。不只各位能不能提供一些比较快速算法。
: 多谢了

s*h2008-05-14 07:05

4 楼

放到MATLAB里,让它自己找算法求x. x=A\D;
这个根本不需要求逆.
求逆是特别费时间和内存的任务,不到万不得已不要求.

l*n2008-05-14 07:05

5 楼

直接用gauss消去

【在 t**********n 的大作中提到】

p*e2008-05-14 07:05

6 楼

这个就是LU分解啊，不知道lz为什么觉得慢

【在 s****h 的大作中提到】

: 放到MATLAB里,让它自己找算法求x. x=A\D;
: 这个根本不需要求逆.
: 求逆是特别费时间和内存的任务,不到万不得已不要求.

b*e2008-05-14 07:05

7 楼

cholesky

h*z2008-05-14 07:05

8 楼

How fast do you need it to be?

【在 t**********n 的大作中提到】

M*o2008-05-14 07:05

9 楼

I prefer Matlab, which is highly optimized for solving equations, x = A / D
.

l*i2008-05-14 07:05

10 楼

你这个解一维的方程的确非常快

D

【在 M**********o 的大作中提到】

: I prefer Matlab, which is highly optimized for solving equations, x = A / D
: .

M*o2008-05-14 07:05

11 楼

楼主不就是要解方程么，那您有啥高见

【在 l*****i 的大作中提到】

: 你这个解一维的方程的确非常快
:
: D

l*i2008-05-14 07:05

12 楼

你的方法对2维及2维以上不管用，呵呵

【在 M**********o 的大作中提到】

: 楼主不就是要解方程么，那您有啥高见

h*z2008-05-14 07:05

13 楼

It works in MATLAB (though ManofScorpio typed the wrong slash direction).
To solve Ax = b in MATLAB, where A is a matrix (not necessarily square) and
x, b are vectors, you just type "x = A \ b" and MATLAB will automatically do
the "right" thing: picking a fast solver based on the structure of A,
returning least-squares solution if system is over-determined, etc.

【在 l*****i 的大作中提到】

: 你的方法对2维及2维以上不管用，呵呵

l*i2008-05-14 07:05

14 楼

我说的就是这个/不对

and
do

【在 h***z 的大作中提到】

: It works in MATLAB (though ManofScorpio typed the wrong slash direction).
: To solve Ax = b in MATLAB, where A is a matrix (not necessarily square) and
: x, b are vectors, you just type "x = A \ b" and MATLAB will automatically do
: the "right" thing: picking a fast solver based on the structure of A,
: returning least-squares solution if system is over-determined, etc.