由买买提看人间百态

boards

本页内容为未名空间相应帖子的节选和存档,一周内的贴子最多显示50字,超过一周显示500字 访问原贴
Computation版 - 如何快速求一个5000*5000的对称矩阵的逆
相关主题
Fast solver of linear equation systemMATLAB积分求教
dense matrix、sparse matrix终于完成一模拟matlab基本矩阵运算的C库
在线等:紧急求助Matlab高手(ODE问题)求一个特殊矩阵的特征值
[合集] 请教:关于用arpack求教PCA
Binary integer programming in Matlab?怎么用Matlab把一个3x3的矩阵diagonalize?
请教一个计算速度的问题求助:奇怪的Matlab运算错误
BLAS 矩阵相乘问题请问怎样让一m长的向量和一mxn的矩阵相乘,仍为mxn
calculation of inverse of "almost-singular" matrix如何normalize矩阵
相关话题的讨论汇总
话题: matlab话题: 5000话题: 矩阵话题: 对称话题: lu
进入Computation版参与讨论
1 (共1页)
t**********n
发帖数: 68
1
我需要 解 一个线性方程 A*X = D, 其中A 是一个 5000×5000 的对称矩阵,
X = Inv(A)*D.
我试过用 LU 的方法,不过速度很慢。 不只各位能不能提供一些比较快速算法。
多谢了
p*****e
发帖数: 310
2
5000*5000不大啊,LU还慢?

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

s****h
发帖数: 921
3
放到MATLAB里,让它自己找算法求x. x=A\D;
这个根本不需要求逆.
求逆是特别费时间和内存的任务,不到万不得已不要求.
l******n
发帖数: 9344
4
直接用gauss消去

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

p*****e
发帖数: 310
5
这个就是LU分解啊,不知道lz为什么觉得慢

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

b***e
发帖数: 38
6
cholesky
h***z
发帖数: 233
7
How fast do you need it to be?

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

M**********o
发帖数: 203
8
I prefer Matlab, which is highly optimized for solving equations, x = A / D
.
l*****i
发帖数: 3929
9
你这个解一维的方程的确非常快

D

【在 M**********o 的大作中提到】
: I prefer Matlab, which is highly optimized for solving equations, x = A / D
: .

M**********o
发帖数: 203
10
楼主不就是要解方程么,那您有啥高见

【在 l*****i 的大作中提到】
: 你这个解一维的方程的确非常快
:
: D

l*****i
发帖数: 3929
11
你的方法对2维及2维以上不管用,呵呵

【在 M**********o 的大作中提到】
: 楼主不就是要解方程么,那您有啥高见
h***z
发帖数: 233
12
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*****i
发帖数: 3929
13
我说的就是这个/不对

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.

1 (共1页)
进入Computation版参与讨论
相关主题
如何normalize矩阵Binary integer programming in Matlab?
请问C里面有没有现成的matrix computation as matlab请教一个计算速度的问题
大家是怎么在MATLAB里读入矩阵的?BLAS 矩阵相乘问题
[help]should I switch from Matlab to C?calculation of inverse of "almost-singular" matrix
Fast solver of linear equation systemMATLAB积分求教
dense matrix、sparse matrix终于完成一模拟matlab基本矩阵运算的C库
在线等:紧急求助Matlab高手(ODE问题)求一个特殊矩阵的特征值
[合集] 请教:关于用arpack求教PCA
相关话题的讨论汇总
话题: matlab话题: 5000话题: 矩阵话题: 对称话题: lu