【正文】 end 開始 形成節(jié)點(diǎn)導(dǎo)納矩陣 輸入原始數(shù)據(jù) 設(shè)節(jié)點(diǎn)電壓， i 1,2? ,n,is 置迭代次數(shù) 置節(jié)點(diǎn)號 i 1 按式（ 33），（ 34）計(jì)算雅克比矩陣元素 按式（ 32）計(jì)算節(jié)點(diǎn)的，節(jié)點(diǎn)的， 求解修正方程式，得， 雅克比矩陣是否已全部形成？ 計(jì)算平衡節(jié)點(diǎn)及 PV 節(jié)點(diǎn)功率 求， 迭 代次數(shù) k k+1 i i+1 ？ 潮流計(jì)算完成 計(jì)算各節(jié)點(diǎn)電壓的新值： 。 3*x 2 +x 3 11, 3*x 1 ,x 1 。y 3 ]。 y [y 1 。 y 2 3*x 1 *x 2 +x 1 *x 3 11*x 1 。,i [x dx F x ] %在屏幕上輸出每次的 x i ,dx i ,F x i if abs dx eps amp。 sprintf 39。 dx A\b。 for i 1:1, A Fd x 。 x0 [ ]39。 end format long 1 100。y 2 。 y 3 exp x 2 *x 3 +20*x 3 + 10*pi3 /3。 end function y F x y 1 3*x 1 cos x 2 *x 3 。 2*x 1 , 162* x 2 + ,cos x 3 。,i [x dx F x ] %在屏幕上輸出每次的 x i ,dx i ,F x i if abs dx eps1 amp。 sprintf 39。 dx A\b。 for i 1:1, A Fd x 。 x0 [ ]39。 eps1 。 y 3*sin a 。 end a 0:pi/12:2*pi。 for i 1:n y y+ 1 ^ i1 *1/ 2*i1 。 toc。 % 利用公式計(jì)算接地支路的潮流 end end end end 作業(yè) y pifun 10000 y clear all tic。 else S_res k,4 bus J,2 ^2*conj YtYm k,5 。 % 利用公式計(jì)算非接地支路的潮流 else if J 0 S_res k,3 bus I,2 ^2*conj YtYm k,4 。 S_res k,4 bus J,2 ^2* conj YtYm k,3 +conj YtYm k,5 bus I,2 *bus J,2 * conj cos bus I,3 +j*sin bus I,3 * cos bus J,3 +j*sin bus J,3 *conj YtYm k,3 。 if J~ 0 amp。 % 前兩列為節(jié)點(diǎn)編號 for k 1:nl I S_res k,1 。 S_res zeros nl,5 。 % 相角采用角度制 bus_res :,4 bus :,4 + sqrt 1 *bus :,5 。 % bus_res 矩陣儲(chǔ)存著節(jié)點(diǎn)計(jì)算結(jié)果 bus_res :,1:2 bus :,1:2 。 end end 計(jì)算并返回節(jié)點(diǎn)數(shù)據(jù)結(jié)果 [nb,mb] size bus 。 YtYm k,4 Ym+Yt* 1+K 。 YtYm k,5 Yt* 1K /K/K。 end if K 0 % 變壓器線路 : Zt 和 Ym 為折算到 i 側(cè)的值 ,K 在 j 側(cè) YtYm k,3 Yt/K。 end if K 0 amp。 YtYm k,4 Ym。 if K 0 amp。 end Ym line k,5 +j*line k,6 。 Zt line k,3 +j*line k,4 。 for k 1:nl I line k,1 。 % 對 YtYm 矩陣賦初值 0 YtYm :,1:2 line :,1:2 。 break end end end end % 恢復(fù) line 的編號 作用為計(jì)算線路的等效 Yt 和 Ym，以計(jì)算線路潮流 [nl,ml] size line 。 end end end % 利用 bus 矩陣的首列編號重新對 bus 矩陣排序并存入 bus_temp 矩陣中 bus bus_temp。 for i 1 :nb for j 1 : nb if bus j,1 k bus_temp k,: bus j,: 。 bus_temp zeros nb,mb 。 end % 計(jì)算平衡節(jié)點(diǎn)的無功及有功注入 作用為對節(jié)點(diǎn)和線路數(shù)據(jù)恢復(fù)編號 [nb,mb] size bus 。 end end % 利用公式計(jì)算 PV 節(jié)點(diǎn)的無功注入 for j 1:n bus n,4 bus n,4 +bus n,2 *bus j,2 * real Y n,j *cos bus n,3 bus j,3 +imag Y n,j *sin bus n,3 bus j,3 。 % 生成雅克比矩陣 作用為計(jì)算每個(gè)節(jié)點(diǎn)的功率注入 n nPQ+nPV+1。 end end % 分別計(jì)算 L 矩陣的對角及非對角元素 end end end J [H N。 end % 分別計(jì)算 K 矩陣的對角及非對角元素 if j nPQ+1 if i~ j L i,j bus i,2 *bus j,2 * real Y i,j *sin bus i,3 bus j,3 imag Y i,j *cos bus i,3 bus j,3 。 end end % 分別計(jì)算 N 矩陣的對角及非對角元素 if i nPQ+1 if i~ j K i,j bus i,2 *bus j,2 * real Y i,j *cos bus i,3 bus j,3 +imag Y i,j *sin bus i,3 bus j,3 。 end % 分別計(jì)算 H 矩陣的對角及非對角元素 if j nPQ+1 if i~ j N i,j bus i,2 *bus j,2 * real Y i,j *cos bus i,3 bus j,3 +imag Y i,j *sin bus i,3 bus j,3 。 end end % 初始化并計(jì)算 Qi 和 Pi for i 1:nb1 for j 1:nb1 if i~ j H i,j bus i,2 *bus j,2 * real Y i,j *sin bus i,3 bus j,3 imag Y i,j *cos bus i,3 bus j,3 。 for i 1:nb1 for j 1:nb Pi i,1 Pi i,1 +bus i,2 *bus j,2 * real Y i,j *cos bus i,3 bus j,3 +imag Y i,j *sin bus i,3 bus j,3 。 K L]，并初始化 Qi zeros nb1,1 。 L zeros nPQ,nPQ 。 N zeros nb1,nPQ 。 end end end % 利用循環(huán)計(jì)算求取 dP 和 dQ