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