潮流matlab程序

电力系统潮流计算matlab程序，内含有4、14、30、1024节点的计算

function [ Jacobian] = formJaco（ YcosYsindiagVolPQPVNodebalanceNodenode）
% 形成雅克比矩阵
%  编写时间：2013.1.112013.1.122013.1.13
%  修改时间：2013.1.16（在平衡节点与PV节点置零位置错误）2013.1.17（运行时间过长，原因在于雅克比矩阵的形成每次都需要计算）
% H矩阵的编写：将H矩阵的表达式展开，利用矩阵的左乘右乘来计算;J矩阵与H矩阵的编写类似


VYVsin = sparse（diagVol*Ysin*diagVol）;
VYVcos = sparse（diagVol*Ycos*diagVol）;
sparseQ = sparse（diag（Q））; 
sparseP = sparse（diag（P））; 
%%形成H矩阵
H = sparse（sparseQ-VYVsin）;
H（balanceNode:） = 0;
H（:balanceNode） = 0;
H（balanceNodebalanceNode） = 1;
%%形成J矩阵
J = sparse（VYVcos-sparseP）;
J（PVNode:） = 0;
%J（:PVNode） = 0;
J（balanceNode:） = 0;
J（:balanceNode） = 0;
%%形成N矩阵
%N1 = sparseP;
%N2 = VYVcos;
N = sparse（-sparseP-VYVcos）;
%N（PVNode:） = 0;
N（:PVNode） = 0;
N（balanceNode:） = 0;
N（:balanceNode） = 0;
%%形成L矩阵
%L1 = sparseQ;
%L2 = VYVsin;
L = sparse（sparseQ+VYVsin）;
L = （diag（-1.*ones（node1）））*L;
L（PVNode:） = 0;
L（:PVNode） = 0;
L = L+sparse（PVNodePVNode1nodenode）;
Jacobian =sparse（ [HN;JL]）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1115  2013-01-17 16:59  潮流程序\formJaco.m

     文件       1464  2013-04-16 21:58  潮流程序\formY.asv

     文件       1435  2013-04-16 20:51  潮流程序\formY.m

     文件     185785  2009-04-16 16:10  潮流程序\IEEE1047.dat

     文件      16630  2008-11-16 22:28  潮流程序\IEEE118.dat

     文件       1110  2009-12-26 10:53  潮流程序\IEEE14.dat

     文件       4017  2009-04-16 09:29  潮流程序\IEEE30.dat

     文件      61718  2009-08-14 10:10  潮流程序\IEEE300.dat

     文件        321  2009-12-23 10:10  潮流程序\IEEE4.dat

     文件       1041  2013-01-17 21:23  潮流程序\unbalancePQ.m

     目录          0  2013-04-24 10:31  潮流程序

----------- ---------  ---------- -----  ----

               274636                    11