BPA生成导纳阵

sunek

回复 1# goodluckyang


, I: w" z# ^/ D; F8 s9 a    感谢浩哥分享，这个导纳阵算潮流的时候怎么输出我挺好奇的

这个导纳阵我是这么理解的
" Z, G: c( r5 q+ P( x
ROW=1
. n' E/ v4 G* Z5 m# M) t+ \导纳阵第一行 i=1  代表与节点1 相连的节点的导纳
  JBUS 6 Gij 0 Bij 4181BFCE
+ Z1 P+ X- V0 [  `# f6 K, [  导纳阵第六列 j=6 代表节点6与节点1的互导纳 互电导Gij=0 互电纳Bij=4181BFCE
  同理
8 M8 w, |/ a, c  JBUS 1 Gij 0 Bij C1D348FE
  导纳阵第1列 节点1 自电导Gij=0 自电纳Bij=C1D348FE


关于4181BFCE

计算机存储数据用的是补码 这里应该是十六进制补码
5 n! M0 B- V* T/ m7 [( ~导纳阵 互电纳是正的 自电纳是负的 互电纳Bij=4181BFCE正数 自电纳Bij=C1D348FE负数 符合

% C& `( U# N6 `- x% Z; Z. a4181BFCE=（0100 0001 1000 0001 1011 1111 1100 1110）二进制补码=（0100 0001 1000 0001 1011 1111 1100 1110）二进制原码=（1099022286）十进制
C1D348FE=（1100 0001 1101 0011 0100 1000 1111 1110）二进制补码=（1011 1110 0010 1100 1011 0111 0000 0010）二进制原码=（-1043117826）十进制
9 ^9 L" h* E6 j
4 }, ?0 h' z1 Q. w至于小数点在哪一位 我觉得BPA可能采用定点数运算 还真得拿实际的数据比较下才知道
' n. N4 Z3 R( U% n0 ?# h5 B
* r7 A  @& d' C! j* ?0 U0 N/ P另：BPA里线路、变压器参数格式 是按6.5的格式


希望对你有帮助
sunek

redplum

顶楼主一下

goodluckyang

回复 11# sunek
) [, o8 O% B- `) d- G
$ T/ R8 ^9 P8 X/ L5 m
2 _' s( J/ J2 V! D: y# q* z   sunek版主， 这个导纳阵在潮流计算时，可能无法输出，在稳定计算时通过DEBUG予以输出；    嗯，关于版主的说明，我按照你的思路，验证下是不是与真实导纳阵结果一样，验证完之后向版主汇报

goodluckyang

回复 11# sunek
) U8 z/ ]% }" U6 @0 t4 T$ u/ }! Y      提供的思路不错，应该是这样的，但是验证结果不太理想。
如下：
2 ]% ]6 `1 z) A! b, w/ Y1，首先  ROW=1时，Bij 4181BFCE  与   Bij C1D348FE  之和应为零，而结果好像不是为零。
/ L! X, |) z9 Y# i3 Y2，本例中导纳阵的结果应为，如图：

0 [" c  H" i9 d3，对应的173611的十六进制（源码反码补码一致）2A62B，去对应相应的十六进制，未果。

感谢版主回答，继续挖。

sunek

回复 14# goodluckyang


! t& f) ~6 F$ _3 W    在稳定计算的时候怎么通过DEBUG输出？
    这个导纳阵我看不太懂，是下三角吗？节点4、5、6、7、8的自导纳是啥？写成矩阵可能好看一些。

goodluckyang

回复 15# sunek
! }( J8 V' q* x/ ~8 f1 n DEBUG输出如下：
# A, I  X3 r/ x2 k# {7 E  @
, k( C9 E# C, q; Y: g9节点的导纳阵如下：

可以与BPA的对应下，但是对应不起来

linhaiming

BPA能输出导纳矩阵么？希望牛人解答，顶你！

周沭铭

你搞得 太深了