使用numpy求解Zoepritz 方程矩阵伪逆时报错: SVD did not converge

笔者在使用numpy中的pinv函数求解伪逆时系统报错: SVD did not converge. 奇异值分解不收敛

具体原因不太清楚, 应该是因为函数在求解伪逆的算法在迭代过程中难以收敛导致的.

解决方法: 引入scipy中的求解伪逆的函数 scipy.linalg.pinv

 

Zoepritz equations的python程序在: https://github.com/cui-xiaoang96/Zoepritz-equations