Python和科学计算认证群组  - 讨论区

标题:关于numpy高效运算

2012年12月13日 星期四 16:18

 

这是我的两个式子,W为400*4225的矩阵,为稀疏矩阵,但是稀疏系数不大,Dx为4225*4225的稀疏矩阵,稀疏稀疏很大,其中很多为0,DY同Dx,想请问下该怎么用程序表示这两个公式呢?

 

2012年12月14日 星期五 11:35

用scipy的稀疏矩阵库试试看吧。

2012年12月18日 星期二 19:23

对于第一个式子,假设那个需要求逆的矩阵是A,

b =W.T * y

那么第一个式子就是

A * x = b

这个是标准的线性方程组求解问题,不需要对矩阵A求逆。如果A比较稀疏的话,你可以用 scipy.sparse.linalg试试看;如果不稀疏就用linalg.solve试试看。应该都比直接求逆要快。

此外,在优化之前,你应该测量一下每一步运算所需的时间,找到瓶颈,然后有针对性的进行优化。

 

 

如下红色区域有误,请重新填写。

    你的回复:

    请 登录 后回复。还没有在Zeuux哲思注册吗?现在 注册 !

    Zeuux © 2024

    京ICP备05028076号