Function Fafangcheng(canshu() As Double, P() As Double, ByVal m As Integer, N As Integer, L() As Double) Dim zhuanchengB() As Double Dim zhuanchengB1() As Double ReDim zhuanchengB(0 To N - 1, 0 To m - 1) ReDim zhuanchengB1(0 To N - 1, 0 To N - 1) Dim zhuancanShu() As Double ReDim zhuancanShu(0 To N - 1, 0 To m - 1) Dim i, j As Integer For i = 0 To N - 1 For j = 0 To m - 1 zhuancanShu(i, j) = canshu(j, i) Next j Next i Call mmult(zhuancanShu, P, N, m, m) For i = 0 To N - 1 For j = 0 To m - 1 zhuanchengB(i, j) = Results(i, j) Next j Next i '''''BT*P Call mmult(zhuanchengB, canshu, N, m, N) For i = 0 To N - 1 For j = 0 To N - 1 zhuanchengB1(i, j) = Results(i, j) Next j Next i '''''BTPB Call niJuzhen(zhuanchengB1, N) Dim nichengB() As Double ReDim nichengB(0 To N - 1, 0 To N - 1) For i = 0 To N - 1 For j = 0 To N - 1 nichengB(i, j) = Results(i, j) Next j Next i '''''BTPB逆 Call mmult(nichengB, zhuancanShu, N, N, m) Dim zaizhuanB() As Double ReDim zaizhuanB(0 To N - 1, 0 To m - 1) For i = 0 To N - 1 For j = 0 To m - 1 zaizhuanB(i, j) = Results(i, j) Next j Next i '''''BTPB逆BT Dim zaizhuanB1() As Double ReDim zaizhuanB1(0 To N - 1, 0 To m - 1) Call mmult(zaizhuanB, P, N, m, m) For i = 0 To N - 1 For j = 0 To m - 1 zaizhuanB1(i, j) = Results(i, j) Next j Next i ReDim re(0 To N - 1) Call mmult(zaizhuanB1, L, N, m, 1) For i = 0 To N - 1 re(i) = Results(i, 0) Next i End Function i = 0 To N - 1 For j = 0 To m - 1 zhuanchengB(i, j) = Results(i, j) Next j Next i '''''BT*P Call mmult(zhuanchengB, canshu, N, m, N) For i = 0 To N - 1 For j = 0 To N - 1 zhuanchengB1(i, j) = Results(i, j) Next j Next i '''''BTPB Call niJuzhen(zhuanchengB1, N) Dim nichengB() As Double ReDim nichengB(0 To N - 1, 0 To N - 1) For i = 0 To N - 1 For j = 0 To N - 1 nichengB(i, j) = Results(i, j) Next j Next i '''''BTPB逆 Call mmult(nichengB, zhuancanShu, N, N, m) Dim zaizhuanB() As Double ReDim zaizhuanB(0 To N - 1, 0 To m - 1) For i = 0 To N - 1 For j = 0 To m - 1 zaizhuanB(i, j) = Results(i, j) Next j Next i '''''BTPB逆BT Dim zaizhuanB1() As Double ReDim zaizhuanB1(0 To N - 1, 0 To m - 1) Call mmult(zaizhuanB, P, N, m, m) For i = 0 To N - 1 For j = 0 To m - 1 zaizhuanB1(i, j) = Results(i, j) Next j Next i ReDim re(0 To N - 1) Call mmult(zaizhuanB1, L, N, m, 1) For i = 0 To N - 1 re(i) = Results(i, 0) Next i End Function