資源簡介
在matlab上面編寫的,基于原對偶內點法的(跟蹤中心路徑參數內點法)最優潮流,具體編程思想請參見王錫凡《現代電力系統分析》
代碼片段和文件信息
function?[A?dh_dx?dg_dx?H_?d2h_dx_y?d2g_dx_c?d2f_dx?temp?L_Z?U_W]?=?Coeff(num_node?num_PQ?Y?B?X?A2)
[num_branch?~]??=?size(B);
len_x?=?2*(num_node-num_PQ+num_node);?%?14
num_gen?=?num_node-num_PQ;?%?2
num_equa?=?2*num_node;
num_inequa?=?2*num_gen+num_node+num_branch;
z?=?X(1:num_inequa);?l?=?X(num_inequa+1:2*num_inequa);?w?=?X(2*num_inequa+1:3*num_inequa);
u?=?X(3*num_inequa+1:4*num_inequa);?x?=?X(4*num_inequa+1:4*num_inequa+len_x);?y?=?X(4*num_inequa+len_x+1:4*num_inequa+len_x+num_equa);
Xtilde?=?x(2*num_gen+1:len_x);
%%?計算等式約束Jacobian矩陣
dh_dPG?=?zeros(22*num_node);
dh_dQR?=?zeros(22*num_node);
for?k?=?num_PQ+1:num_node
???dh_dPG(k-num_PQk*2-1)?=?1;
???dh_dQR(k-num_PQk*2)?=?1;
end
dh_dXtilde?=?zeros(2*num_node2*num_node);
for?ii?=?1:num_node
????f
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????8814??2017-10-16?09:20??內點法\Coeff.m
?????文件???????1207??2017-10-15?21:54??內點法\dX.m
?????文件???????8554??2019-02-28?15:59??內點法\main.m
?????文件????????855??2017-10-07?10:56??內點法\makeY.m
?????文件???????1097??2017-10-27?14:41??內點法\solution.txt
?????目錄??????????0??2019-02-28?15:56??內點法
-----------?---------??----------?-----??----
????????????????20527????????????????????6
評論
共有 條評論