%%%此函數(shù)為雅克比矩陣的算法
clcclear
%%電壓賦值
U=input（‘輸入電壓各節(jié)點電壓值‘）;?
U=U‘;??????????????????????????????????????
e=real（U）;
f=imag（U）;
%%初始化
n=length（U）;??????????????????????????%節(jié)點數(shù)量
Bus=[1?2?3?4?5]‘;???????????????????????%節(jié)點號
Vtype=[1?1?1?2?3]‘;?????????????????????%節(jié)點類型??????????????????????????????
swing=find（Vtype==3）;???????????????????%平衡節(jié)點編號
PV=find（Vtype==2）;??????????????????????%PV節(jié)點編號
PQ=find（Vtype==1）;??????????????????????%PQ節(jié)點編號
Vset=[1?1?1?1.1?1.05];??????????????????%電壓設定值
PVPQ=find（Vtype<3）;?????????????????????%PV/PQ節(jié)點編號
%節(jié)點導納矩陣及B，G提取
Y=[6.311-20.4022j?-3.5587+11.3879j?-2.7523+9.1742j?0?0;
????-3.5587+11.3879j?8.5587-31.5272j?-5+15j?5.1948j?0;
????-2.7523+9.1742j?-5+15j?8.5587-31.5272j?0?5.1948j
????0?5.1948j?0?-5.1948j?0;
????0?0?5.1948?0?-5.1948];
G=real（Y）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2140??2018-12-05?16:45??潮流計算\jacobi_matrix.m

?????文件???????4390??2018-12-06?12:02??潮流計算\Load_flow_calculation.m

?????文件???????1736??2018-12-11?09:57??潮流計算\Node_power_unbalance.m

?????文件???????3338??2018-12-12?18:26??潮流計算\zhongjianguocheng.m

?????目錄??????????0??2019-01-09?22:04??潮流計算

-----------?---------??----------?-----??----

????????????????11604????????????????????5