基于matpower，模拟系统发生级联故障的情况，并用牛顿拉夫逊法潮流计算切除的负荷量，实现简单的牛拉法极坐标下的潮流计算。首先将支路导纳矩阵转化为节点导纳矩阵,初始化功率初值,电压初值,接着形成Jacobian矩阵,对矩阵进行LDU分解,前代回代求解修正量,最后得出各个节点的电压幅值与角度。