y(1,1)=(0.00+0.25i);y (1,2)=1/ (0.04+0.25i); y (1,3)=1/ (0.1+0.35i); y (1,4)=0 ;y(1,5)=0; y(2,2)=(0.5i)+(1/1.05^2-1/1.05)/0.015i;y(2,3)=1/ (0.08+0.30i) ;y(2,4)=1/(0.00+0.015i)/1.05;c(1,2)=0.25i; y(3,3)=(0.00+0.25i)+(1/1.05^2-1/1.05)/0.03i;y(3,5)=1/ (0.00+0.03i)/1.05 ; y(4,4)=(1-1/1.05)/0.015i; y(5,5)=(1-1/1.05)/0.03i;的潮流计算

根据题意，给出的是一个五节点的电力系统，其中节点1为平衡节点，其他节点2-5为PQ节点。已知节点的导纳矩阵为：

Y = | 0.00+0.25i   0.25i         1/ (0.1+0.35i) 0              0           |
    | 0.25i        0.5i+(1/1.05^2-1/1.05)/0.015i 1/ (0.08+0.30i) 1/(0.00+0.015i)/1.05 0           |
    | 1/ (0.1+0.35i) 1/ (0.08+0.30i) 0.00+0.25i+(1/1.05^2-1/1.05)/0.03i 0 1/ (0.00+0.03i)/1.05 |
    | 0            1/(0.00+0.015i)/1.05 0           (1-1/1.05)/0.015i 0           |
    | 0            0            1/ (0.00+0.03i)/1.05 0           (1-1/1.05)/0.03i |

根据潮流计算的基本公式：P=|V|*|V|*G+|V|*|V|*B，Q=|V|*|V|*(-G*sinθ+B*cosθ)，其中G和B为导纳矩阵的实部和虚部，θ为节点电压相角，|V|为节点电压幅值。

由于节点1为平衡节点，电压幅值=1，相角=0，所以只需要求解节点2-5的电压幅值和相角即可。

采用牛顿-拉夫逊法进行迭代求解，迭代公式为：


其中，x为节点电压幅值和相角，F(x)为节点潮流方程，J(x)为潮流方程的雅可比矩阵。

采用MATLAB进行求解，代码如下：