资源描述：
《牛顿-拉夫逊法潮流计算matlab程序（newton raphson method of power flow calculation matlab program）》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、牛顿-拉夫逊法潮流计算matlab程序（NewtonRaphsonmethodofpowerflowcalculationmatlabprogram）主程序”powerflow_nr。”[bus_res功能，s_res]=powerflow_nr_2%牛顿-拉夫逊法解潮流方程的主程序【巴士线]=opdf_；%打开数据文件的子程序，返回总线（节点数据）和线（线路数据）回主程序[NB，MB=大小（总线）；【NL，毫升]=大小（线）；%计算总线和线矩阵的行数和列数【巴士线，NPQ，NPV，nodenum]=num_（公共汽车、线）；%对节点重新排
2、序的子程序Y=y_（公共汽车、线）；%计算节点导纳矩阵的子程序该=fopen（'result。M’，’w'）；fprintf（Diego，“---------------通过longdinhohe---------------http://blog.sina.com.cn/longdinhohe---------------NN’）；关闭文件（该）；%在当前目录下生成”的结果。”文件，写入节点导纳矩阵长格式每股收益=1.0e-10；%设定误差精度T=100%开始迭代计算，设定最大迭代次数为100，以便不收敛情况下及时跳出【DP、DQ]=dp
3、q_（Y，总线，NPQ，NPV）；%计算功率偏差DP和DQ的子程序J=jac_（公共汽车、Y、NPQ）；%计算雅克比矩阵的子程序UD=0（NPQ，NPQ）；我=1：水UD（我，我）=总线（I，2）；%生成电压对角矩阵结束丹谷=J【DPDQ]；当=丹谷（1：nb-1,1）；%计算相角修正量杜=UD*（丹谷（注：铌npq-1,1））；%计算电压修正量总线（1：NPQ，2）=总线（1：NPQ，2）-杜；%修正电压总线（1：nb-1,3）=总线（1：nb-1,3）-党；%修正相角如果（max（ABS））4、否满足精度误差，如满足则跳出，否则返回继续迭代结束总线=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，“---------------牛顿－拉夫逊法
5、潮流计算结果----------NN节点计算结果：N节点节点电压节点相角（角度）节点注入功率N’）；对于i=1：NBfprintf（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)));EndFprintf(MYF,'n线路计算结果:n节点i节点j线路功率s(
6、I,j)线路功率s(J,i)线路损耗ds(I,j)n');ForI=1:NLFprintf(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
7、),(_imagsres(I,5);EndFclose(MYF);%迭代结束后继续在"result.m"写入节点计算结果和线路计算结果程序结束"子程序1ompF_.M"作用为打开数据文件Function[]=ompF_buslineDfilepathname[],uigetfile(m='*','Selectdatafile');%数据文件类型为m文件,窗口标题为"Selectdatafile"Ifpathname==0Error("YoumustselectavalidDatafile'%如果没有选择有效文件,则出现错误提示)ElseLf
8、ile=length(dfile去除文件后缀%);(dfile(1:lfile-2);%打开文件End子程序2"num.m"作用为对节点重排序,并修改相应的线路数据Functio
