电力系统潮流计算

代码说明：

首先输入原始数据，包括节点个数NB 支路个数NL 发电机个数NG节点数据矩阵 databus(节点类型 节点号 发电机有功 发电机无功 负荷有功 负荷无功 电压幅值 相角，标志，节点功率限值 电压限值）节点类型 1 PQ节点 0 PV节点 -1 平衡节点）支路数据矩阵 dataline(起始节点号 终止节点号 支路电阻 支路电抗 支路变比 支路对地导纳，支路有功潮流（7）,支路无功潮流（8），支路有功逆潮流，支路逆无功潮流，额定电流（线路容量限制），线路每小时故障率，线路停运率，标志，支路标识）发电机组参数矩阵（发电机序号，所在节点号，有功功率（GW），无功功率上限，每小时的故障率，每小时的修复率的倒数,有功下限，无功下限，空载电压，机端电压，发电机电抗）然后使用牛顿-拉夫逊法对电力系统进行潮流计算，输出节点电压幅值和相角