对这
个节点进行编号。这种编号方法的根据是导纳矩阵中出线支路数最少的节点所对应的行中非零元素也
2动态地按增加出线支路数最少编号在上述的方法中各节点的出线支路数是按原始网络统计出来的在编号过程中认为固定不变的事实上在节点消去过程中每消去一个节点以后与该节点相连的各节点的出线支路数将发生变化增加减少或保持不变。因此如果每消去一个节点后立即修正尚未编号节点的出线支路数然后选其中支路数最少的一个节点进行编号就可以预期得到更好的效果动
f态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。
3MATLAB编程应用
Matlab是“MatrixLaboratory”的缩写主要包括一般数值分析矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。因为使用Matlab编程运算与人进行科学计算的思路和表达方式完全一致所以不像学习高级语言那样难于掌握而且编程效率和计算效率极高还可在计算机上直接输出结果和精美的图形拷贝所以它的确为一高效的科研助手。
四、设计内容
设计流程图
f1程序
clearclc
重新编号把原题中的节点12345重新依次编号为51234其中14号为PQ节点5号为平衡节点
y0
输入原始数据求节点导纳矩阵
y121006018iy131006018iy141004012i
y151002006i
y231001003iy251008024i
y341008024i
y450
fori15
fforji5
yjiyij
e
d
e
d
Y0
求互导纳
fori15
forj15
ifij
Yijyij
e
d
e
d
e
d
求自导纳
fori15
Yiisumyi
e
d
YY为导纳矩阵
GrealY
BimagY
原始节点功率
S10202i
S2045015i
S304005i
S40601i
S50
PrealS
QimagS
赋初值
Uo
es15U5106
ezeros15
oxo
es81fxo
es81
cou
t0计算迭代次数
whilemaxfx1e5
fori14
forj14
Hij0Nij0Mij0Lij0oPi0oQi0
e
d
e
d
fori14
forj15oPioPiUiUjGijcoseiejBijsi
eiej
oQioQiUiUjGijsi
eiejBijcoseiej
e
d
oPioPiPioQioQiQi
e
d
ffxoPoQ
求雅克比矩阵
当ij时候求HNML如下
fori14
forj14
ifijHijUiUjGijsi
eiejBijcoseiej
NijUiUjGijcoseiejBijsi
eiej
LijHij
MijNij
e
d
e
d
e
d
HNML
当ij时HNML如下
fori14
forj15
ifij
HiiHiiUiUjGijsi
eiejBijcoseiej
NiiNiiUiUjGijcoseiejBijsi
eiej
MiiMiiUiUjGijcoseiejBijsi
eiejLiiLiiUiUjGijsi
eiejBijr