基于MATLAB牛顿拉夫逊法进行夏季最小潮流计算
网上找到计算潮流的版本之一如下：看到一个帖子：http://www.ilovematlab.cn/viewthread.php?tid=13338&pid=230344&page=1&extra=page%3D1#pid230344
有谁可以解释一下B2各列参数的意义吗？谢谢!
请输入节点数:n=12
请输入支路数:nl=11
请输入平衡母线节点号:isb=1
请输入误差精度:pr=0.00001
请输入由各支路参数形成的矩阵:B1=
请输入各节点参数形成的矩阵:B2=
导纳矩阵Y=
Columns 1 through 7
31.2739 -83.2494i -2.7123 +14.1874i -7.8023 +16.2549i -9.6108 +19.7021i -6.7873 +14.3288i -4.2052 +13.8772i 0
-2.7123 +14.1874i 2.7123 -13.9634i 0 0 0 0 0
-7.8023 +16.2549i 0 7.9751 -22.2912i 0 0 0 0
-9.6108 +19.7021i 0 0 9.7768 -25.1143i 0 0 0
-6.7873 +14.3288i 0 0 0 6.9534 -19.6955i 0 0
-4.2052 +13.8772i 0 0 0 0 10.0733 -37.9479i -5.7121 +18.6596i
0 0 0 0 0 -5.7121 +18.6596i 5.7121 -18.5071i
0 0 -0.1641 + 5.8760i 0 0 0 0
0 0 0 -0.1578 + 5.2584i 0 0 0
0 0 0 0 -0.1578 + 5.2584i 0 0
0 0 0 0 0 -0.1559 + 5.5822i 0
-0.1686 + 6.0349i 0 0 0 0 0 0
Columns 8 through 12
0 0 0 0 -0.1686 + 6.0349i
0 0 0 0 0
-0.1641 + 5.8760i 0 0 0 0
0 -0.1578 + 5.2584i 0 0 0
0 0 -0.1578 + 5.2584i 0 0
0 0 0 -0.1559 + 5.5822i 0
0 0 0 0 0
0.1559 - 5.5822i 0 0 0 0
0 0.1499 - 4.9955i 0 0 0
0 0 0.1499 - 4.9955i 0 0
0 0 0 0.1559 - 5.5822i 0
0 0 0 0 0.1822 - 6.5242i
迭代次数
4
没有达到精度要求的个数
21 21 14 0
各节点的实际电压标么值E为(节点号从小到大排列):
Columns 1 through 7
1.0500 0.9951 - 0.0992i 1.0293 - 0.0184i 1.0305 - 0.0105i 1.0255 - 0.0160i 1.0671 + 0.0619i 1.0925 + 0.1283i
Columns 8 through 12
1.0225 - 0.0954i 1.0044 - 0.0833i 0.9983 - 0.0890i 1.0110 - 0.0197i 1.0211 + 0.0930i
各节点的电压大小V为(节点号从小到大排列):
1.0500 1.0000 1.0295 1.0305 1.0256 1.0689 1.1000 1.0269 1.0079 1.0023 1.0112 1.0253
各节点的电压相角时shita为(节点号从小到大排列):
0 -5.6958 -1.0221 -0.5827 -0.8925 3.3184 6.6965 -5.3305 -4.7429 -5.0961 -1.1182 5.2027
各节点的功率S为(节点号从小到大排列):
Columns 1 through 7
1.3716 + 0.2627i -1.6000 - 0.5772i -0.0000 + 0.0000i 0.0000 - 0.0000i 0.0000 - 0.0000i -0.0000 - 0.0000i 1.5000 + 0.0966i
Columns 8 through 12
-0.4753 - 0.2947i -0.4074 - 0.3608i -0.4074 - 0.3608i -0.4753 - 0.2947i 0.6000 + 0.3720i
各条支路的首端功率Si为(顺序同您输入B1时一样):
1.6349 + 0.2888i
0.4827 + 0.0380i
0.4139 + 0.1626i
0.4150 + 0.0691i
-0.9770 + 0.0037i
-0.4753 - 0.2947i
-0.4074 - 0.3608i
-0.4074 - 0.3608i
-0.4753 - 0.2947i
-1.4711 - 0.3611i
-0.5980 - 0.2994i
各条支路的末端功率Sj为(顺序同您输入B1时一样):
-1.6000 - 0.5772i
-0.4768 - 0.3478i
-0.4091 - 0.4192i
-0.4092 - 0.4198i
0.9943 + 0.0116i
0.4768 + 0.3478i
0.4091 + 0.4192i
0.4092 + 0.4198i
0.4768 + 0.3495i
1.5000 + 0.0966i
0.6000 + 0.3720i
各条支路的功率损耗DS为(顺序同您输入B1时一样):
0.0349 - 0.2884i
0.0060 - 0.3098i
0.0047 - 0.2566i
0.0058 - 0.3508i
0.0173 + 0.0153i
0.0015 + 0.0531i
0.0017 + 0.0583i
0.0018 + 0.0590i
0.0015 + 0.0548i
0.0289 - 0.2644i
0.0020 + 0.0726i clc
clear
n=input('请输入节点数:n=');
nl=input('请输入支路数:nl=');
isb=input('请输入平衡母线节点号:isb=');
pr=input('请输入误差精度:pr=');
B1=input('请输入由各支路参数形成的矩阵:B1=');
B2=input('请输入各节点参数形成的矩阵:B2=');
Y=zeros(n); e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);
O=zeros(1,n);S1=zeros(nl);
%求导纳矩阵
for i=1:nl
if B1(i,6)==0 %NO.6
p=B1(i,1);q=B1(i,2);
else p=B1(i,2);q=B1(i,1);
end
Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));%非对角元素
Y(q,p)=Y(p,q);%非对角元素
Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;%变压器支路 NO.5变比
Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;
end
disp('导纳矩阵Y=');
disp(Y);
G=real(Y);B=imag(Y);
for i=1:n%n=节点数
e(i)=real(B2(i,3));%NO.3节点电压初值
f(i)=imag(B2(i,3));
V(i)=B2(i,4);%？
end
for i=1:n
S(i)=B2(i,1)-B2(i,2);% NO.1节点发出的功率，NO.2节点负荷大小，S（i）负荷数据
B(i,i)=B(i,i)+B2(i,5);% B=imag(Y)，NO.5？
end
P=real(S);Q=imag(S); B2的第4，5列代表什么含义？
请输入各节点参数形成的矩阵:B2=[
0 0 1.05 0 0 1;
0 1.6+0.992i 1.05 1 0 3;
0 0 1 0 0 2;
0 0 1 0 0 2;
0 0 1 0 0 2;
0 0 1 0 0 2;
0.7708 0 1.1 1.1 0 3;
0 0.539+0.33418i 1 0 0 2;
0 0.66+0.4092i 1 0 0 2;
0 0.66+0.4092i 1 0 0 2;
0 0.77+0.4774i 1 0 0 2;
1+0.62i0.4+0.248i 1 0 0 2
] 谢谢楼主啊,很详细!
页:
[1]
