python 实现matlab左除
import numpy as np
import scipy
def mldivide_new(A,b):
[mm,nn]=np.shape(A)
ss=np.shape(b)[1]
x=np.zeros((nn,ss))
if np.shape(A)[0]==np.shape(A)[1]:
if A==np.tril(A):
for ii in range(ss):
x[:,ii]=backward_sub(A,b[:,ii])
return
elif A==np.triu(A):
for ii in range(ss):
x[:,ii]=backward_sub(A,b[:,ii])
return
else:
if A==A.T:
[R,p]=np.linalg