电力系统潮流计算Matlab实际算例
答案:2 悬赏:40 手机版
解决时间 2021-02-18 12:28
- 提问者网友：贪了杯
- 2021-02-17 19:13
谁有这个的牛顿拉夫逊法和PQ节点的实际算例，谢谢急用，有的话发到我邮箱shijie.47@163.com
最佳答案
- 五星知识达人网友：一秋
- 2021-02-17 19:26
网上搜 PSAT，里面有潮流算例
全部回答
- 1楼网友：你可爱的野爹
- 2021-02-17 20:16
% 配电网前推回代潮流计算程序
% 使用ieee 33节点配电系统作为算例，可实现弱环网情况下的潮流计算
countnum=0;
branchdata = [
1 2 0.0922 0.0470;
2 3 0.4930 0.2511;
3 4 0.3660 0.1864;
4 5 0.3811 0.1941;
5 6 0.8190 0.7070;
6 7 0.1872 0.6188;
7 8 0.7114 0.2351;
8 9 1.0300 0.7400;
9 10 1.0440 0.7400;
10 11 0.1966 0.0650;
11 12 0.3744 0.1238;
12 13 1.4680 1.1550;
13 14 0.5416 0.7129;
14 15 0.5910 0.5260;
15 16 0.7463 0.5450;
16 17 1.2890 1.7210;
17 18 0.7320 0.5740;
2 19 0.1640 0.1565;
19 20 1.5042 1.3554;
20 21 0.4095 0.4784;
21 22 0.7089 0.9373;
3 23 0.4512 0.3083;
23 24 0.8980 0.7091;
24 25 0.8960 0.7011;
6 26 0.2030 0.1034;
26 27 0.2842 0.1447;
27 28 1.0590 0.9337;
28 29 0.8042 0.7006;
29 30 0.5075 0.2585;
30 31 0.9744 0.9630;
31 32 0.3105 0.3619;
32 33 0.3410 0.5302;
]; % 支路，阻抗
nodedata = [
2 100.00 60.00;
3 90.00 40.00;
4 120.00 80.00;
5 60.00 30.00;
6 60.00 20.00;
7 200.00 100.00;
8 200.00 100.00;
9 60.00 20.00;
10 60.00 20.00;
11 45.00 30.00;
12 60.00 35.00;
13 60.00 35.00;
14 120.00 80.00;
15 60.00 10.00;
16 60.00 20.00;
17 60.00 20.00;
18 90.00 40.00;
19 90.00 40.00;
20 90.00 40.00;
21 90.00 40.00;
22 90.00 40.00;
23 90.00 50.00;
24 420.00 200.00;
25 420.00 200.00;
26 60.00 25.00;
27 60.00 25.00;
28 60.00 20.00;
29 120.00 70.00;
30 200.00 600.00;
31 150.00 70.00;
32 210.00 100.00;
33 60.00 40.00;
]; % 节点，负荷
ub = 12.66; % 电压基准 kv
sb = 10; % 功率基准 mva
zb = ub^2/sb; % 阻抗基准 ohm
branchdata(:,[3,4]) = branchdata(:,[3,4]) / zb; % 阻抗标幺化
nodedata(:,[2,3]) = nodedata(:,[2,3]) / sb / 1000;% 功率标幺化
nn = 33; % 节点数
a0 = zeros(nn);
for n = 1:nn-1
a0(branchdata(n,1),branchdata(n,2)) = 1;
end % 形成 a0
associatedmatrix=0;
for n=2:nn-1
associatedmatrix(n,n)=1;
temp=branchdata(n-1,1);
associatedmatrix(n,1:n-1)=associatedmatrix(temp,1:n-1);
end
a0t = a0'; % 形成 a0 的转置
s = [0;-nodedata(:,2) - i*nodedata(:,3)]; % 形成 s
zl = [0;branchdata(:,3) + i*branchdata(:,4)]; % 形成 zl
v = ones(nn,1);
v(1) = 1; % 各个节点电压赋初值
il(nn,1) = -conj(s(nn) / v(nn)); % 最末支路电流赋初值
delta = 1; % 收敛判据赋初值
tempv = v; % 赋初值，用于记忆上次迭代结果
while delta > 1e-8
countnum=countnum+1;
in = conj(s ./ v); % 节点注入电流
for n = 1:nn-1
il(nn-n) = a0(nn-n,nn-n+1:end) * il(nn-n+1:end) - in(nn-n);
end % 电流回代过程
for n = 2:nn
v(n) = a0t(n,1:n-1) * v(1:n-1) - zl(n) * il(n);
end % 电压前推过程
delta = max(abs(v-tempv)); % 更新收敛判据
tempv = v; % 记忆迭代结果
end
vangle(:,1)=abs(v);
vangle(:,2)=angle(v)/3.1415*180;
for i=1:nn-1
st=branchdata(i,1);
en=branchdata(i,2);
sij(i,1)=v(st)*conj((v(st)-v(en))/zl(i+1));
sji(i,1)=v(en)*conj((v(en)-v(st))/zl(i+1));
end
希望你能用到
我要举报
如以上问答信息为低俗、色情、不良、暴力、侵权、涉及违法等信息，可以点下面链接进行举报！
大家都在看
推荐资讯