牛顿-拉夫逊法潮流计算matlab程序（Newton Raphson method of power flow calculation matlab program）.docVIP

牛顿-拉夫逊法潮流计算matlab程序（Newton Raphson method of power flow calculation matlab program） 主程序”powerflow_nr。” [ bus_res功能，s_res ] = powerflow_nr_2 %牛顿-拉夫逊法解潮流方程的主程序 【巴士线] = opdf_；%打开数据文件的子程序，返回总线（节点数据）和线（线路数据）回主程序 [ NB，MB =大小（总线）； 【NL，毫升] =大小（线）；%计算总线和线矩阵的行数和列数 【巴士线，NPQ，NPV，nodenum ] = num_（公共汽车、线）；%对节点重新排序的子程序 Y = y_（公共汽车、线）；%计算节点导纳矩阵的子程序 该= fopen（result。M’，’w）； fprintf（Diego，“---------------通过longdinhohe --------------- /longdinhohe --------------- NN’）； 关闭文件（该）；%在当前目录下生成”的结果。”文件，写入节点导纳矩阵 长格式 每股收益= 1.0e-10；%设定误差精度 T = 100 %开始迭代计算，设定最大迭代次数为100，以便不收敛情况下及时跳出 【DP、DQ ] = dpq_（Y，总线，NPQ，NPV）；%计算功率偏差DP和DQ的子程序 J = jac_（公共汽车、Y、NPQ）；%计算雅克比矩阵的子程序 UD = 0（NPQ，NPQ）； 我= 1：水 UD（我，我）=总线（I，2）；%生成电压对角矩阵 结束 丹谷= J【DP DQ ]； 当=丹谷（1：nb-1,1）；%计算相角修正量 杜= UD *（丹谷（注：铌npq-1,1））；%计算电压修正量 总线（1：NPQ，2）=总线（1：NPQ，2）-杜；%修正电压 总线（1：nb-1,3）=总线（1：nb-1,3）-党；%修正相角 如果（max（ABS）） EPS】和（max（abs）） 打破 完成判断是否满足精度误差，如满足则跳出，否则返回继续迭代 结束 总线= pq_（公共汽车、Y、NPQ，NPV）；%计算每个节点的有功和无功注入的子程序 【巴士线] = renum_（巴士线，nodenum）；%对节点恢复编号的子程序 ytym = ytym_（线）；%计算线路的等效YT和YM的子程序，以计算线路潮流 bus_res = bus_res_（总线）；%计算节点数据结果的子程序 s_res = s_res_（公共汽车、线、YtYm）；%计算线路潮流的子程序 该= fopen（result。M，A）； fprintf（Diego，“---------------牛顿－拉夫逊法潮流计算结果---------- NN节点计算结果：N节点节点电压节点相角（角度）节点注入功率N’）； 对于i = 1：NB fprintf（Diego，“% 2.0f ，bus_res（I，1））； fprintf（Diego，“% 10.6f ，bus_res（I，2））； fprintf（Diego，10.6f “%， Bus _ RES (I, 3); Fprintf (MYF, 10.6f % + J% 10.6fn (bus, real _ RES (I, 4)), IMAG (bus _ RES (I, 4))); End Fprintf (MYF, n线路计算结果: n节点i 节点j 线路功率s (I, j) 线路功率s (J, i) 线路损耗ds (I, j) n); For I = 1: NL Fprintf (MYF, 2.0f % s, _ RES (I, 1)); Fprintf (MYF, 2.0f % s, _ RES (I, 2)); Fprintf (MYF, 10.6f % + J% 10.6f, real (S _ RES (I, 3), IMAG (S _ RES (I, 3))); Fprintf (MYF, 10.6f % + J% 10.6f, real (S _ RES (I, 4)), IMAG (S _ RES (I, 4))); Fprintf (MYF, 10.6f % + J% 10.6fn, real (S _ RES (I, 5), (_ imag s res (I, 5); End Fclose (MYF);% 迭代结束后继续在 result.m 写入节点计算结果和线路计算结果 程序结束 子程序1 ompF _. M 作用为打开数据文件 Function [] = ompF _ bus line Dfile pathname [], uigetfile (m = *, Select data fi