-
1 votesanswersviews
Sympy - 矩阵逆的密集矩阵显示
我正在使用Sympy来运行一个基本的回归示例,希望它能成为分解更复杂的估算器的有效工具 . 但是,我无法弄清楚如何为流程中的每个步骤显示密集矩阵 . 特别是,我试图表示最小二乘估计值: 这是仅有五个观察结果的设置 . from sympy import * y=MatrixSymbol('y',5,1) x=MatrixSymbol('x',5,2) b=MatrixSymbol('b',2,... -
6 votesanswersviews
python中稀疏矩阵的伪逆
我正在使用来自神经成像的数据,并且由于数据量很大,我想为我的代码使用稀疏矩阵(scipy.sparse.lil_matrix或csr_matrix) . 特别是,我需要计算矩阵的伪逆以解决最小二乘问题 . 我找到了sparse.lsqr方法,但效率不高 . 有没有一种方法来计算Moore-Penrose的伪逆(对应于正常矩阵的pinv) . 我的矩阵A的大小约为600'000x2000,并且在矩阵... -
1 votesanswersviews
Matlab中3变量非线性系统的Newton-Raphson方法
我试图在MATLAB中使用Newton-Raphson方法求解3个变量的3个非线性系统 . 以下是3个非线性方程: c * (alpha*I + k_f + k_d + k_n * s + k_p*(1-q))-I *alpha = 0 s * (lambda_b * c* P_C + lambda_r *(1-q))- lambda_b* c * P_C = 0 q * ( ga...