function?[DPDQ]=culation（VY1dijPQPVIbalanceNode）%電壓幅值?導納矩陣幅值?導納矩陣的相角差?有功?無功?PV節點?平衡節點
%P=P‘;
%Q=Q‘;
dij=sparse（dij）;%轉化為稀疏矩陣

DP=-P+diag（V）*（Y1.*cos（dij））*V;%diag（V）將V轉化為對角陣，對應書本上的公式P65
DQ=-Q+diag（V）*（Y1.*sin（dij））*V;
DP（balanceNode:）=0;

DQ（PVI:）=0;%PV節點沒有無功功率不平和方程式；

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????603??2016-05-28?01:50??FormY.m
?????文件????????1898??2017-01-05?09:16??IEEE30.txt
?????文件????????2755??2017-04-07?15:47??MC.m
?????文件????????2240??2016-05-29?12:45??OutData.m
?????文件????????1219??2016-05-27?23:18??ReadData.m
?????文件??????250620??2016-07-10?19:35??WF1.txt
?????文件??????257609??2016-07-10?19:36??WF2.txt
?????文件?????????354??2016-03-28?19:56??culation.m
?????文件?????????741??2017-03-20?09:51??jacobian.m
?????文件????????1649??2016-06-01?18:18??main.m
?????文件?????????979??2017-03-20?09:55??newLine.m
?????文件?????????364??2017-03-20?09:37??xiuzheng.m