|
|
马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
%本程序的功能是用牛顿-拉夫逊法进行潮流计算# u9 j, D9 r7 K& G+ l/ c. A
n=input('请输入节点数:n=');; b3 L2 T+ n" l' l7 N6 Y
nl=input('请输入支路数:nl=');
5 P7 v; u; N K/ rB1=input('请输入由支路参数形成的矩阵:B1=');%变压器侧为1,否则为0
/ i2 `' w' p K7 ]) dX=input('请输入由节点号及其对地阻抗形成的矩阵:X=');
, W( x7 Y2 u9 {5 M" g% }Y=zeros(n);5 z8 p! W |! G) [5 f
for i=1:n
( L" f. N; ]! M% q if X(i,2)~=0;
K7 f6 ?# q" V# v" u) H3 E p=X(i,1);
0 f D, z& w3 A9 A- `; ?2 M Y(p,p)=X(i,2);
: q# D9 y9 s, L3 V1 S+ \ end& A! r' D1 O4 g9 r( I
end! M! r9 b$ B8 V* b. w2 U+ I8 \
for i=1:nl3 s7 `' w+ F2 G
if B1(i,6)==0
7 o6 U$ s/ `/ H+ h: c" M0 P p=B1(i,1);q=B1(i,2);
( Y# X& t( G8 i( x) G4 k else p=B1(i,2);q=B1(i,1);
' [' l }9 A' j# n0 H D end
* A3 K& N% o: F Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));
- P1 r+ w/ m- p Y(q,p)=Y(p,q);
5 s( ?' ]5 t G# T Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;; f, ~5 o$ p0 U4 u
Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;0 [* F% J2 W) V2 b$ Q
end %求导纳矩阵/ w0 s! B% A" f$ b9 ~/ G
我在求导纳矩阵时,关于变压器的处理方面老是弄不清楚怎么折算,请问谁研究过这个程序?导纳矩阵的程序是不是有问题 |
|