请大神帮我看一下 矩阵初等变换的python程序到底怎么错了？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

import numpy as np
def triangleFun(a,b):
row=a.shape[0]
col=a.shape[1]
n=row-1
for i in range(0,n):
  for j in range(i+1,n): #变为上对角元素
    lam =float(a[j,i])/a[i,i]
    a[j,0:n]=a[j,0:n]-lam*a[i,0:n]
    b[j,0:n]=b[j,0:n]-lam*a[i,0:n]  
for p in range(n,-1,-1):
  for q in range(p-1,-1,-1):  #变为下对角元素
    lam=float(a[q,p])/a[p,p]
    a[q,::-1]=a[q,::-1]-lam*a[p,::-1]
    b[q,::-1]=b[q,::-1]-lam*a[p,::-1]
for l in range(0,n+1):
   a[l,l]=a[l,l]/a[l,l]
   b[l,:]=b[l,:]/a[l,l]
print(a,b)

a =np.mat([[1,2,1,4],[2,0,4,3],[4,2,2,1],[-3,1,3,2]])
b=np.eye(4)
print(triangleFun(a,b))

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：python Triangle import range Shape

程序出来了 就是答案不对

def triangleFun(self):
        a=self.a
        b=self.b
        row=a.shape[0]
        col=a.shape[1]
        n=row-1
        for i in range(0,n):
            for j in range(i+1,n+1):
                if a[i,i]==0:
                    continue
                lam =float(a[j,i])/a[i,i]
                a[j,0:n+1]=a[j,0:n+1]-lam*a[i,0:n+1]
                b[j,0:n+1]=b[j,0:n+1]-lam*b[i,0:n+1]
        for p in range(n,-1,-1):
            for q in range(p-1,-1,-1):
                if a[p,p]==0:
                    continue
                lam=float(a[q,p])/a[p,p]
                a[q,0:n+1]=a[q,0:n+1]-lam*a[p,0:n+1]
                b[q,0:n+1]=b[q,0:n+1]-lam*b[p,0:n+1]
        for l in range(0,n+1):
            lam=a[l,l]
            if lam==0: