基于MATLAB的潮流计算源程序代码(共7页).docx

%*电力系统直角坐标系下的牛顿拉夫逊法潮流计算*clearclcload E:dataIEEE014_Node.txtNode=IEEE014_Node;weishu=size(Node);nnum=weishu(1,1); %节点总数 load E:dataIEEE014_Branch.txtbranch=IEEE014_Branch;bwei=size(branch); bnum=bwei(1,1); %支路总数Y=(zeros(nnum);Sj=100;%*节点导纳矩阵*for m=1:bnum; s=branch(m,1); %首节点 e=branch(m,2)