解释这行代码 P[i] = U[i]*np.sum(U*(G[i,:]*np.cos(Angle[i]-Angle) + B[i,:]*np.sin(Angle[i]-Angle)))

这行代码的目的是将一个复杂的矩阵运算结果赋值给P数组的第i个元素。具体来说，它包括以下几个步骤：

1. G和B是两个矩阵，分别代表导纳矩阵中的实部和虚部。
2. U是一个向量，代表每个节点的电压幅值。
3. Angle是一个向量，代表每个节点的电压相角。
4. np.cos(Angle[i]-Angle)和np.sin(Angle[i]-Angle)分别计算了每个节点与第i个节点的电压相角之差的余弦和正弦值。
5. G[i,:]和B[i,:]是G和B矩阵的第i行，代表第i个节点与其他节点之间的实部和虚部导纳值。
6. 将G[i,:]与np.cos(Angle[i]-Angle)逐元素相乘，得到一个与节点数相同的向量。
7. 将B[i,:]与np.sin(Angle[i]-Angle)逐元素相乘，得到一个与节点数相同的向量。
8. 将U与上述两个向量分别逐元素相乘，得到两个长度为节点数的向量。
9. 将这两个向量相加，并对结果求和，得到一个标量。
10. 将该标量赋值给P数组的第i个元素。

综上所述，这行代码的作用是计算节点i的功率注入值，其中包括与其他节点之间的导纳值、电压幅值和相角等多个因素的综合影响。