importnumpyasnpimportmathdefsgn1(x ) : ifx 0: return1elifx==0: return0elifx 0: return-1ai2=NP.mat )欢乐的流沙将dtype=float(n=AI2.shape )0) forIinrange(n-2 ) :c=-1*SGN1 ) AI2 ) I1，I ) *改为ai2[i 1:n，I ) i] l1=l.copy ) L1(1) ) u=NP.vstack () ) a， b ) I=NP.mat(NP.eye(n ) ) h=I-2 * u * u.tai2=h * ai2 * h.terr=1a i3=ai2.copy ) ) fortin range (88 ) I ) ) cV1 ) I1，I )-sV1 ) I1 )=CAI2) V1*AI2I=1sita=math.Atan(AI2[I1，i]/ai2[i，i] ) C=math.cos )