-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathtest.py
91 lines (86 loc) · 3.3 KB
/
test.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
from input_data import input_net_args
from NL_iteration import NL_Iteration
# Tips ：
# 序号PQ节点在前，PV节点在后
# 输入参数为标幺值
# 1.数据输入
Line_arg = [
# 导线首端 导线末端 串联电阻 串联电抗 并联电导 并联电纳
[1 , 2 , 0.00037 , 0.00404 , 0 , 0.12264],
[2 , 3 , 0.00023 , 0.00309 , 0 , 0.16046],
[2 , 4 , 0.000585 , 0.07445 , 0 , 1.61600],
[4 , 5 , 0.00006 , 0.00170 , 0 , 0.53379],
[5 , 3 , 0.00052 , 0.00702 , 0 , 1.46965],
[5 , 10 , 0.00030 , 0.00459 , 0 , 1.00998],
[5 , 9 , 0.00068 , 0.00830 , 0 , 1.63571],
[10 , 14 , 0.00005 , 0.00045 , 0 , 0.09629],
[10 , 9 , 0.00013 , 0.00148 , 0 , 0.30356],
[9 , 11 , 0.00243 , 0.00942 , 0 , 0.03080],
[9 , 13 , 0.00034 , 0.00455 , 0 , 0.40547],
[9 , 8 , 0.00070 , 0.00885 , 0 , 1.82756],
[8 , 12 , 0.00025 , 0.00046 , 0 , 0.97515],
[8 , 7 , 0.00046 , 0.00057 , 0 , 1.04222],
[7 , 15 , 0.00049 , 0.00511 , 0 , 0.24952],
[7 , 6 , 0.00046 , 0.00637 , 0 , 1.25707],
[6 , 16 , 0.00015 , 0.00187 , 0 , 0.09495],
[6 , 3 , 0.00113 , 0.01521 , 0 , 0.71247],
# 变压器
[17 , 14 , 0.00023 , 0.019500 , 0 , 0],
[18 , 4 , 0.00012 , 0.012199 , 0 , 0],
[19 , 15 , 0.00060 , 0.044110 , 0 , 0],
[20 , 16 , 0.00022 , 0.019280 , 0 , 0]
]
Node_args = [
# 节点序号 类型 参数
[1 , "pq" , {"p":3.98600, "q":0.1640} ],
[2 , "pq" , {"P":3.46400, "Q":0.0000} ],
[3 , "pq" , {"P":8.31100, "Q":-0.686} ],
[4 , "pq" , {"P":0.00000, "Q":0.0000} ],
[5 , "pq" , {"P":1.17400, "Q":3.8850} ],
[6 , "pq" , {"p":15.0900, "q":-3.781} ],
[7 , "pq" , {"P":3.18600, "Q":0.7900} ],
[8 , "pq" , {"P":5.72100, "Q":3.6850} ],
[9 , "pq" , {"P":14.9710, "Q":-1.466} ],
[10, "pq" , {"P":9.07100, "Q":0.7590} ],
[11, "pq" , {"P":6.29200, "Q":-3.761} ],
[12, "pq" , {"P":5.47700, "Q":-2.305} ],
[13, "pq" , {"P":12.2850, "Q":1.2770} ],
[14, "pq" , {"P":0, "Q":0} ],
[15, "pq" , {"P":0, "Q":0} ],
[16, "pq" , {"P":0, "Q":0} ],
[17, "pv" , {"P":11.3330, "v":0.9931} ],
[18, "pv" , {"P":55.7030, "v":0.9884} ],
[19, "pv" , {"P":5.98600, "v":0.9815} ],
[20, "slack" , {"V":1.0, "Theta":-0.7828} ]
]
Init_val = [
# 节点序号 参数
[1, {"e":1.02, "f":0}],
[2, {"e":1, "f":0}],
[3, {"e":1, "f":0}],
[4, {"e":0.605533, "f":-0.781194}],
[5, {"e":0.556615, "f":-0.813259}],
[6, {"e":0.63727226, "f":-0.6927696}],
[7, {"e":0.561588141, "f":-0.786478}],
[8, {"e":0.4674631, "f":-0.8591612}],
[9, {"e":0.428363, "f":-0.8940802}],
[10, {"e":0.463429, "f":-0.87664366}],
[11, {"e":0.4069631, "f":-0.9118022}],
[12, {"e":0.4503565, "f":-0.8830939}],
[13, {"e":0.3723258, "f":-0.9114882}],
[14, {"e":0.4080243, "f":-0.9054081}],
[15, {"e":0.5949248, "f":-0.7806451}],
[16, {"e":1.02, "f":0}],
[17, {"e":0.3723258, "f":-0.9114882}],
[18, {"e":0.4080243, "f":-0.9054081}],
[19, {"e":0.5949248, "f":-0.7806451}],
[20, {"e":1.02, "f":0}]
]
# 数据处理
args = input_net_args(Line_arg, Node_args, Init_val)
args.gen_node_admittance_matrix()
args.gen_node_infos()
args.gen_init_values()
# 迭代
nl = NL_Iteration(args)
nl.start_iteration()