bài toán dưới 533 nút xuất các thông số đều đạt yêu cầu khi thêm EVs; còn lớn hơn sẽ bị lỗi ở phần B và dẵ update hiệu quả hơn 2000bus

[SharkpJIM]

% --- Main IEEE_33bus ---
clear; clc;
define_constants; % Nạp các hằng số của MATPOWER (như PQ, PV, VM, PD, v.v.)

% Tải dữ liệu lưới điện từ file case_numberdata.m
mpc = case_ACTIVSg500S();

% --- PHẦN A: CHẠY LƯỚI CƠ BẢN ---
fprintf('================================================\n');
fprintf(' KẾT QUẢ PHẦN A\n');
fprintf('================================================\n');

% A.1 Chạy phân bố công suất
results_A = runpf(mpc, mpoption('verbose', 0, 'out.all', 0));

if results_A.success
fprintf('-> Chạy phân bố công suất thành công!\n\n');

% A.2 Xuất điện áp mỗi nút 
fprintf('--- 2. Điện áp tại các nút (P.U và kV) ---\n');
fprintf('Nut\t\t|V| (p.u.)\tGoc (do)\tDien ap (kV)\n');
V_base = 12.66; % kV
for i = 1:size(results_A.bus, 1)
    bus_id = results_A.bus(i, BUS_I);
    vm_pu = results_A.bus(i, VM);
    va_deg = results_A.bus(i, VA);
    vm_kv = vm_pu * V_base;
    fprintf('%2d\t\t%.4f\t\t%.2f\t\t%.3f\n', bus_id, vm_pu, va_deg, vm_kv);
end

% A.3 Xuất dòng điện trên mỗi nhánh 
% Tính dòng điện: Iij = Sij / Vi (theo đơn vị có tên hoặc p.u.)
% Sbase = 10 MVA, Vbase = 12.66 kV -> Ibase = S / (sqrt(3)*V)
S_base_MVA = 10;
I_base_kA = S_base_MVA / (sqrt(3) * V_base);
I_base_A = I_base_kA * 1000; % ~456.04 Amps

fprintf('\n--- 3. Dòng điện trên các nhánh ---\n');
fprintf('Tu Nut\tDen Nut\t\tDong (A)\t\tDong (p.u.)\n');

% Lấy công suất dòng chảy (MVA) từ kết quả
S_from = sqrt(results_A.branch(:, PF).^2 + results_A.branch(:, QF).^2);
% Lấy điện áp tại nút đi (From Bus)
Vm_from = results_A.bus(results_A.branch(:, F_BUS), VM);

% Tính dòng I (p.u.) = S (p.u.) / V (p.u.)
% Lưu ý: Trong matpower S tính bằng MVA, cần đổi sang p.u. bằng cách chia baseMVA
I_pu = (S_from / S_base_MVA) ./ Vm_from;
I_abs = I_pu * I_base_A;

for k = 1:size(results_A.branch, 1)
    fbus = results_A.branch(k, F_BUS);
    tbus = results_A.branch(k, T_BUS);
    fprintf('%2d\t\t%2d\t\t\t%.2f\t\t%.4f\n', fbus, tbus, I_abs(k), I_pu(k));
end

% A.4 Tính tổn thất công suất P và Q 
loss = get_losses(results_A);
total_loss_P = sum(real(loss)); % MW
total_loss_Q = sum(imag(loss)); % MVAr
fprintf('\n--- 4. Tổn thất toàn hệ thống ---\n');
fprintf('Tong ton that P: %.4f MW\n', total_loss_P);
fprintf('Tong ton that Q: %.4f MVAr\n', total_loss_Q);

else
fprintf('Phân bố công suất thất bại (Phần A).\n');
end

% --- PHẦN B: TÍNH TOÁN VỚI XE ĐIỆN (EVs) [cite: 47] ---
fprintf('\n\n================================================\n');
fprintf(' KẾT QUẢ PHẦN B (CÓ EVs)\n');
fprintf('================================================\n');

% B.1 Nhập vị trí và công suất trạm sạc
% Yêu cầu: 5 đến 10 trạm, công suất tùy chọn.
% Giả sử ta chọn 6 trạm sạc tại các nút xa nguồn để thấy rõ sụt áp.
mpc_EV = mpc; % Tạo bản sao dữ liệu lưới

% Cấu hình trạm sạc (Bạn có thể thay đổi số liệu ở đây)
% Cột 1: Số thứ tự nút (Bus ID)
% Cột 2: Công suất tác dụng P (kW) - Tự chọn tùy ý
% Cột 3: Công suất phản kháng Q (kVAr) - Thường sạc EV cosphi ~ 1 hoặc 0.95
ev_locations = [
2, 150, 0;
4, 180, 0;
14, 120, 0;
17, 200, 0;
26, 150, 0;
29, 210, 0;
33, 200, 0; % Trạm sạc nhanh lớn
];

fprintf('Đã thêm %d trạm sạc EV vào hệ thống.\n', size(ev_locations, 1));

% Cập nhật tải vào dữ liệu mpc_EV
for i = 1:size(ev_locations, 1)
bus_idx = ev_locations(i, 1);
p_add_kw = ev_locations(i, 2);
q_add_kvar = ev_locations(i, 3);

% Chuyển đổi kW -> MW (Dữ liệu trong runpf dùng MW)
mpc_EV.bus(bus_idx, PD) = mpc_EV.bus(bus_idx, PD) + p_add_kw/1000;
mpc_EV.bus(bus_idx, QD) = mpc_EV.bus(bus_idx, QD) + q_add_kvar/1000;

fprintf(' - Nut %d: Them %.2f kW\n', bus_idx, p_add_kw);

end

% B.2 Tính lại các tham số (Tương tự phần A)

% a) Chạy phân bố công suất [cite: 51]
results_B = runpf(mpc_EV, mpoption('verbose', 0, 'out.all', 0));

if results_B.success
% b) Xuất điện áp [cite: 52]
fprintf('\n--- b) Điện áp nút khi có EVs (So sanh V_A va V_B) ---\n');
fprintf('Nut\t\tV_Cu (p.u)\tV_Moi (p.u)\tChenh lech\n');
for i = 1:size(results_B.bus, 1)
v_old = results_A.bus(i, VM);
v_new = results_B.bus(i, VM);
fprintf('%2d\t\t%.4f\t\t%.4f\t\t%.4f\n', ...
results_B.bus(i, BUS_I), v_old, v_new, v_new - v_old);
end

% c) Xuất dòng điện [cite: 53]
S_from_B = sqrt(results_B.branch(:, PF).^2 + results_B.branch(:, QF).^2);
Vm_from_B = results_B.bus(results_B.branch(:, F_BUS), VM);
I_pu_B = (S_from_B / S_base_MVA) ./ Vm_from_B;
I_abs_B = I_pu_B * I_base_A;

fprintf('\n--- c) Dòng điện nhánh khi có EVs (Hien thi 5 nhanh dau tien) ---\n');
fprintf('(Xem chi tiet trong bien I_abs_B)\n');
fprintf('Tu Nut\tDen Nut\t\tDong Moi (A)\n');
for k = 1:5 
    fprintf('%2d\t\t%2d\t\t\t%.2f\n', ...
        results_B.branch(k, F_BUS), results_B.branch(k, T_BUS), I_abs_B(k));
end

% d) Tính tổn thất [cite: 54]
loss_B = get_losses(results_B);
total_loss_P_B = sum(real(loss_B));
total_loss_Q_B = sum(imag(loss_B));

fprintf('\n--- d) Tổn thất khi có EVs ---\n');
fprintf('Tong ton that P: %.4f MW (Tang them %.4f MW)\n', ...
    total_loss_P_B, total_loss_P_B - total_loss_P);
fprintf('Tong ton that Q: %.4f MVAr\n', total_loss_Q_B);

else
fprintf('Phân bố công suất thất bại (Phần B) - Có thể quá tải hệ thống.\n');
end
-----------------------------------------------------------------------------------------------------------------------------================================================
KẾT QUẢ PHẦN A

-> Chạy phân bố công suất thành công!

--- 2. Điện áp tại các nút (P.U và kV) ---
Nut |V| (p.u.) Goc (do) Dien ap (kV)
1 1.0131 -10.70 12.826
2 1.0125 -10.77 12.818
3 1.0086 -11.99 12.769
4 1.0064 -12.24 12.741
5 1.0195 -1.31 12.907
6 1.0194 -1.32 12.906
7 1.0348 -6.87 13.101
8 1.0296 -8.65 13.034
9 1.0400 -4.07 13.166
10 1.0280 -13.75 13.015
11 1.0280 -13.75 13.015
12 1.0096 -14.45 12.781
13 1.0090 -14.51 12.774
14 1.0321 -3.40 13.066
15 1.0247 -6.01 12.973
16 1.0400 -0.56 13.166
17 1.0400 0.00 13.166
18 1.0400 0.89 13.166
19 1.0205 -14.21 12.920
20 1.0205 -14.21 12.919
21 1.0111 -12.26 12.800
22 1.0090 -12.51 12.773
23 1.0309 -7.88 13.051
24 1.0167 -10.75 12.872
25 1.0167 -10.75 12.871
26 1.0164 -10.79 12.868
27 1.0216 -16.39 12.933
28 1.0205 -16.52 12.920
29 1.0110 -14.50 12.799
30 1.0095 -14.68 12.780
31 1.0348 -16.95 13.101
32 1.0347 -16.97 13.099
33 1.0119 -15.12 12.810
34 1.0115 -15.17 12.806
35 1.0200 -15.79 12.914
36 1.0194 -15.86 12.906
37 1.0241 -0.93 12.965
38 1.0241 -0.93 12.965
39 1.0332 -7.76 13.080
40 1.0268 -8.96 12.999
41 1.0258 -9.08 12.987
42 1.0273 -13.54 13.006
43 1.0273 -13.55 13.005
44 1.0194 -9.71 12.905
45 1.0193 -9.71 12.905
46 1.0082 -13.09 12.764
47 1.0076 -13.17 12.756
48 1.0218 -17.62 12.936
49 1.0400 -17.47 13.166
50 1.0400 -17.43 13.166
51 1.0218 -17.62 12.935
52 1.0214 -17.66 12.931
53 1.0180 -13.17 12.888
54 1.0175 -13.23 12.882
55 1.0260 -13.11 12.990
56 1.0260 -13.11 12.989
57 1.0172 -9.59 12.877
58 1.0106 -12.37 12.794
59 1.0083 -12.66 12.764
60 1.0234 -9.51 12.956
61 1.0234 -9.51 12.956
62 1.0128 -11.84 12.822
63 1.0128 -11.84 12.822
64 1.0122 -11.90 12.815
65 1.0309 -10.60 13.051
66 1.0282 -12.28 13.017
67 1.0281 -12.29 13.016
68 1.0166 -13.44 12.871
69 1.0164 -13.47 12.867
70 1.0258 -14.82 12.986
71 1.0400 -14.74 13.166
72 1.0400 -14.68 13.166
73 1.0400 -14.65 13.166
74 1.0273 -16.79 13.006
75 1.0272 -16.80 13.005
76 1.0226 -0.06 12.946
77 1.0225 -0.07 12.945
78 1.0184 -8.41 12.892
79 1.0183 -8.41 12.892
80 1.0346 -10.13 13.099
81 1.0296 -12.70 13.035
82 1.0400 -4.13 13.166
83 1.0125 -15.51 12.818
84 1.0115 -15.63 12.805
85 1.0284 -15.77 13.019
86 1.0281 -15.81 13.015
87 1.0115 0.67 12.806
88 1.0115 0.66 12.805
89 1.0216 -11.06 12.934
90 1.0215 -11.07 12.933
91 1.0151 -7.55 12.851
92 1.0141 -7.67 12.838
93 1.0216 -6.38 12.933
94 1.0215 -6.39 12.933
95 1.0154 -7.52 12.855
96 1.0149 -7.57 12.849
97 1.0219 -17.99 12.938
98 1.0213 -18.06 12.930
99 1.0260 -9.10 12.990
100 1.0258 -9.13 12.987
101 1.0260 -9.10 12.989
102 1.0173 -16.49 12.879
103 1.0167 -16.57 12.871
104 1.0209 -17.35 12.925
105 1.0196 -17.51 12.908
106 1.0226 -14.69 12.946
107 1.0224 -14.71 12.943
108 1.0189 -10.03 12.900
109 1.0182 -10.11 12.891
110 1.0319 -4.91 13.064
111 1.0310 -5.03 13.052
112 1.0164 -13.68 12.867
113 1.0157 -13.76 12.859
114 1.0216 -16.19 12.934
115 1.0216 -16.20 12.933
116 1.0216 -16.20 12.933
117 1.0240 -13.86 12.964
118 1.0238 -13.88 12.962
119 0.9993 -15.28 12.651
120 0.9989 -15.33 12.646
121 1.0148 -13.89 12.847
122 1.0139 -14.00 12.836
123 1.0322 -3.74 13.068
124 1.0300 -4.44 13.040
125 1.0322 -3.74 13.068
126 1.0322 -3.74 13.068
127 1.0400 0.40 13.166
128 1.0400 -2.38 13.166
129 0.9977 -15.78 12.631
130 0.9961 -15.96 12.611
131 1.0240 -9.38 12.964
132 1.0236 -9.43 12.959
133 1.0162 -7.47 12.865
134 1.0156 -7.54 12.858
135 1.0191 -11.26 12.902
136 1.0187 -11.31 12.897
137 1.0194 -13.94 12.906
138 1.0182 -14.10 12.890
139 1.0200 -10.04 12.913
140 1.0200 -10.04 12.913
141 1.0148 -2.55 12.847
142 1.0122 -2.84 12.815
143 1.0235 6.37 12.957
144 1.0400 18.85 13.166
145 1.0400 9.07 13.166
146 1.0299 -7.82 13.039
147 1.0201 -8.48 12.914
148 1.0191 -8.60 12.902
149 1.0262 -13.63 12.992
150 1.0262 -13.64 12.991
151 1.0319 -7.99 13.063
152 1.0312 -10.37 13.055
153 1.0300 -10.52 13.039
154 1.0143 -10.91 12.841
155 1.0134 -11.01 12.830
156 0.9937 -14.51 12.580
157 0.9921 -14.72 12.559
158 1.0124 -13.45 12.816
159 1.0120 -13.50 12.812
160 1.0217 -14.95 12.935
161 1.0213 -15.00 12.930
162 1.0302 -7.78 13.042
163 1.0263 -8.56 12.993
164 1.0246 -8.76 12.972
165 1.0323 -3.39 13.069
166 1.0272 -4.95 13.004
167 1.0400 -1.17 13.166
168 1.0400 0.27 13.166
169 1.0400 1.89 13.166
170 1.0149 -7.65 12.849
171 1.0147 -7.68 12.847
172 1.0158 -12.70 12.860
173 1.0157 -12.71 12.859
174 1.0157 -12.71 12.858
175 1.0365 -3.52 13.122
176 1.0278 -2.71 13.012
177 1.0276 -2.73 13.010
178 1.0277 -2.72 13.010
179 1.0060 -14.67 12.736
180 1.0058 -14.70 12.734
181 1.0281 -17.78 13.015
182 1.0268 -17.94 12.999
183 1.0291 -6.64 13.028
184 1.0287 -6.68 13.023
185 1.0133 -9.21 12.829
186 1.0127 -9.29 12.821
187 1.0206 -11.45 12.921
188 1.0205 -11.46 12.919
189 1.0301 -17.07 13.042
190 1.0299 -17.10 13.039
191 1.0223 -15.67 12.942
192 1.0220 -15.70 12.939
193 1.0175 -14.16 12.882
194 1.0170 -14.22 12.876
195 1.0301 -7.81 13.041
196 1.0209 -10.05 12.925
197 1.0400 -2.94 13.166
198 1.0400 -5.39 13.166
199 1.0196 -10.20 12.908
200 1.0217 -9.31 12.935
201 1.0216 -9.33 12.933
202 1.0120 -12.20 12.812
203 1.0119 -12.21 12.811
204 1.0113 -12.30 12.803
205 1.0242 -16.58 12.966
206 1.0241 -16.59 12.966
207 1.0238 -16.63 12.962
208 1.0170 -16.30 12.875
209 1.0158 -16.46 12.860
210 1.0075 -12.12 12.754
211 1.0069 -12.19 12.747
212 1.0218 -14.77 12.936
213 1.0203 -14.95 12.917
214 1.0164 -14.77 12.867
215 1.0150 -14.95 12.850
216 1.0167 -15.11 12.871
217 1.0159 -15.20 12.861
218 1.0286 -15.52 13.022
219 1.0286 -15.52 13.022
220 1.0375 -3.85 13.135
221 1.0261 -7.11 12.991
222 1.0400 -3.80 13.166
223 1.0400 -3.56 13.166
224 1.0400 -0.59 13.166
225 1.0400 -2.16 13.166
226 1.0154 -7.65 12.855
227 1.0153 -7.66 12.854
228 1.0053 -10.03 12.727
229 1.0045 -10.13 12.717
230 1.0248 -12.74 12.974
231 1.0400 -12.46 13.166
232 1.0302 -8.68 13.042
233 1.0203 -12.00 12.917
234 1.0194 -12.11 12.906
235 1.0314 -12.52 13.058
236 1.0244 -13.75 12.969
237 1.0244 -13.75 12.969
238 1.0199 -11.72 12.912
239 1.0193 -11.79 12.905
240 1.0253 -16.46 12.981
241 1.0250 -16.50 12.977
242 1.0221 -15.86 12.940
243 1.0209 -16.01 12.925
244 1.0197 -12.92 12.909
245 1.0195 -12.94 12.907
246 1.0351 -7.52 13.105
247 1.0302 -6.84 13.042
248 1.0351 -7.52 13.105
249 1.0351 -7.52 13.105
250 1.0351 -7.52 13.105
251 0.9968 -14.66 12.619
252 0.9947 -14.91 12.593
253 1.0020 -14.69 12.686
254 1.0020 -14.69 12.686
255 1.0020 -14.69 12.686
256 1.0020 -14.69 12.686
257 1.0020 -14.69 12.686
258 1.0400 -14.61 13.166
259 1.0019 -14.71 12.684
260 1.0280 -4.79 13.015
261 1.0277 -4.82 13.011
262 1.0296 -8.96 13.034
263 1.0240 -13.18 12.963
264 1.0296 -8.96 13.034
265 1.0296 -8.96 13.034
266 1.0231 -13.29 12.952
267 1.0182 -18.14 12.891
268 1.0165 -18.36 12.870
269 1.0207 -7.22 12.922
270 1.0206 -7.24 12.920
271 1.0308 -12.71 13.050
272 1.0257 -14.45 12.986
273 1.0256 -14.47 12.984
274 1.0307 -5.70 13.048
275 1.0307 -5.70 13.048
276 1.0186 -13.68 12.895
277 1.0185 -13.68 12.895
278 1.0185 -13.69 12.894
279 1.0069 -11.19 12.747
280 1.0066 -11.22 12.744
281 1.0078 -16.14 12.759
282 1.0060 -16.37 12.736
283 1.0190 -14.21 12.900
284 1.0184 -14.28 12.893
285 1.0243 -14.46 12.967
286 1.0242 -14.47 12.966
287 1.0169 -12.49 12.874
288 1.0160 -12.61 12.863
289 1.0213 -6.47 12.929
290 1.0212 -6.48 12.928
291 1.0091 -14.13 12.775
292 1.0090 -14.14 12.774
293 1.0224 -14.88 12.944
294 1.0221 -14.92 12.939
295 1.0164 -12.01 12.868
296 1.0163 -12.02 12.866
297 1.0260 -10.31 12.989
298 1.0256 -11.75 12.984
299 1.0254 -11.78 12.981
300 1.0241 -4.75 12.965
301 1.0400 -4.65 13.166
302 1.0400 -4.67 13.166
303 1.0226 -4.93 12.947
304 1.0037 -15.51 12.707
305 1.0400 -15.30 13.166
306 1.0400 -15.23 13.166
307 1.0026 -15.65 12.693
308 1.0221 -14.49 12.940
309 1.0220 -14.51 12.938
310 1.0173 -17.14 12.880
311 1.0167 -17.23 12.871
312 1.0216 -16.53 12.933
313 1.0214 -16.56 12.931
314 1.0166 -12.92 12.870
315 1.0166 -12.93 12.870
316 1.0241 -2.84 12.965
317 1.0240 -2.85 12.964
318 1.0279 -3.70 13.014
319 1.0400 -2.00 13.166
320 1.0117 -16.02 12.808
321 1.0093 -16.31 12.778
322 1.0069 -14.16 12.747
323 1.0062 -14.23 12.739
324 0.9995 -15.34 12.654
325 0.9979 -15.56 12.634
326 1.0156 -14.86 12.858
327 1.0139 -15.07 12.836
328 1.0086 -12.77 12.769
329 1.0083 -12.81 12.765
330 1.0189 -10.10 12.899
331 1.0180 -10.22 12.888
332 1.0329 -7.68 13.076
333 1.0243 -9.29 12.968
334 1.0243 -9.29 12.968
335 1.0243 -9.29 12.968
336 1.0236 -5.00 12.959
337 1.0235 -5.01 12.958
338 1.0170 5.08 12.875
339 1.0159 4.93 12.861
340 1.0063 -11.77 12.739
341 1.0059 -11.82 12.734
342 1.0215 -14.50 12.932
343 1.0214 -14.51 12.931
344 1.0214 -14.51 12.931
345 1.0040 -14.88 12.710
346 1.0033 -14.96 12.702
347 0.9977 -15.25 12.631
348 0.9965 -15.40 12.615
349 1.0314 -13.36 13.057
350 1.0400 -12.80 13.166
351 1.0400 -13.12 13.166
352 1.0400 -12.82 13.166
353 1.0400 -10.52 13.166
354 1.0229 -9.13 12.950
355 1.0226 -9.17 12.946
356 1.0140 -16.99 12.837
357 1.0132 -17.09 12.827
358 1.0198 -17.67 12.911
359 1.0196 -17.70 12.908
360 1.0073 -12.61 12.753
361 1.0073 -12.61 12.753
362 1.0073 -12.61 12.753
363 1.0073 -12.61 12.753
364 1.0069 -12.67 12.747
365 1.0096 -14.12 12.781
366 1.0087 -14.23 12.770
367 1.0176 -7.71 12.883
368 1.0175 -7.73 12.881
369 1.0210 -6.53 12.926
370 1.0205 -6.59 12.920
371 1.0292 -8.53 13.030
372 1.0208 -12.08 12.923
373 1.0204 -12.12 12.918
374 1.0068 -14.57 12.746
375 1.0065 -14.61 12.742
376 1.0118 -11.28 12.810
377 1.0106 -11.44 12.794
378 1.0266 -4.04 12.997
379 1.0264 -4.06 12.995
380 1.0165 -10.14 12.869
381 1.0163 -10.16 12.867
382 1.0161 -16.94 12.863
383 1.0156 -17.00 12.858
384 1.0085 -14.75 12.767
385 1.0080 -14.81 12.761
386 1.0295 -4.69 13.033
387 1.0243 -6.46 12.967
388 1.0295 -4.69 13.033
389 1.0295 -4.69 13.033
390 1.0295 -4.69 13.033
391 1.0295 -4.69 13.033
392 1.0295 -4.69 13.033
393 1.0295 -4.69 13.033
394 1.0295 -4.69 13.033
395 1.0295 -4.69 13.033
396 1.0295 -4.69 13.033
397 1.0237 -6.53 12.960
398 1.0240 -6.49 12.964
399 1.0246 -14.71 12.972
400 1.0242 -14.75 12.967
401 1.0126 1.88 12.819
402 1.0113 1.70 12.803
403 1.0246 -14.06 12.972
404 1.0243 -14.09 12.968
405 1.0207 -5.35 12.923
406 1.0207 -5.36 12.922
407 1.0356 -11.07 13.111
408 1.0295 -13.18 13.033
409 1.0356 -11.07 13.111
410 1.0400 -9.33 13.166
411 1.0400 -10.08 13.166
412 1.0400 -9.93 13.166
413 1.0400 -9.78 13.166
414 1.0356 -11.07 13.111
415 1.0117 -13.47 12.808
416 1.0116 -13.48 12.807
417 1.0049 -13.97 12.722
418 1.0039 -14.09 12.710
419 1.0191 -17.04 12.902
420 1.0190 -17.05 12.901
421 1.0141 -10.51 12.839
422 1.0134 -10.59 12.830
423 1.0069 1.04 12.747
424 1.0042 0.71 12.713
425 1.0257 -10.60 12.986
426 1.0250 -12.94 12.977
427 1.0248 -12.96 12.974
428 1.0347 -7.52 13.099
429 1.0245 -9.37 12.971
430 1.0400 -7.49 13.166
431 1.0400 -7.46 13.166
432 1.0400 -7.51 13.166
433 1.0400 -7.20 13.166
434 1.0400 -7.48 13.166
435 1.0347 -7.52 13.099
436 1.0347 -7.52 13.099
437 1.0400 -5.59 13.166
438 1.0400 -5.09 13.166
439 1.0400 -6.00 13.166
440 1.0219 -8.99 12.937
441 1.0131 -13.13 12.826
442 1.0219 -8.99 12.937
443 1.0219 -8.99 12.937
444 1.0219 -8.99 12.937
445 1.0110 -10.64 12.800
446 1.0096 -10.81 12.781
447 1.0165 -10.32 12.869
448 1.0156 -10.44 12.857
449 1.0157 -10.59 12.859
450 1.0157 -10.60 12.859
451 1.0157 -10.60 12.859
452 1.0213 1.20 12.930
453 1.0318 -8.82 13.063
454 1.0165 -13.07 12.869
455 1.0400 -8.25 13.166
456 1.0400 -7.63 13.166
457 1.0268 -9.25 12.999
458 1.0400 -5.02 13.166
459 1.0162 4.60 12.865
460 1.0085 -11.66 12.768
461 1.0062 -10.63 12.738
462 1.0234 -9.97 12.957
463 1.0234 -9.97 12.957
464 1.0320 -7.79 13.065
465 1.0275 -8.98 13.008
466 1.0226 -9.65 12.946
467 1.0183 -7.84 12.892
468 1.0088 -10.59 12.771
469 0.9972 -14.22 12.624
470 1.0195 -7.43 12.907
471 1.0327 -7.47 13.074
472 1.0256 -9.01 12.984
473 1.0196 -6.98 12.908
474 0.9908 -14.56 12.543
475 1.0205 -7.89 12.919
476 1.0186 -13.33 12.895
477 1.0171 -11.50 12.877
478 1.0251 -6.59 12.978
479 1.0341 -7.45 13.092
480 1.0400 -5.58 13.166
481 1.0400 -6.14 13.166
482 1.0400 -6.70 13.166
483 1.0400 -6.17 13.166
484 1.0400 -6.64 13.166
485 1.0137 -13.91 12.834
486 1.0166 -13.05 12.870
487 1.0254 -16.35 12.982
488 0.9966 -13.89 12.617
489 1.0253 -13.75 12.981
490 1.0238 -14.60 12.962
491 1.0182 -14.97 12.891
492 1.0169 -5.99 12.873
493 1.0209 -8.66 12.924
494 1.0209 -8.66 12.924
495 1.0133 -10.97 12.828
496 1.0205 4.85 12.919
497 1.0400 6.92 13.166
498 1.0400 5.96 13.166
499 1.0054 -15.20 12.728
500 1.0238 -13.88 12.961

--- 3. Dòng điện trên các nhánh ---
Tu Nut Den Nut Dong (A) Dong (p.u.)
2 1 1995.05 4.3747
44 1 1663.77 3.6483
421 1 297.29 0.6519
4 3 6222.62 13.6448
62 3 852.92 1.8703
3 479 5534.25 12.1354
6 5 148.66 0.3260
316 5 5355.99 11.7445
401 5 5510.46 12.0832
8 7 11491.42 25.1981
9 7 17039.75 37.3644
9 7 17039.75 37.3644
7 232 10105.01 22.1580
7 262 12543.28 27.5046
8 287 8843.36 19.3915
8 457 2649.77 5.8103
11 10 113.42 0.2487
342 10 2520.35 5.5266
10 349 2619.39 5.7437
13 12 1738.18 3.8114
12 216 1651.10 3.6205
12 304 4184.60 9.1759
12 441 7092.25 15.5517
15 14 19260.66 42.2343
16 14 19710.77 43.2213
17 14 19645.73 43.0787
17 14 19643.48 43.0738
18 14 6953.61 15.2477
386 14 23411.01 51.3351
386 14 23394.46 51.2988
15 91 2852.14 6.2541
15 230 5221.51 11.4496
367 15 7693.73 16.8706
15 369 2675.84 5.8675
15 473 955.74 2.0957
20 19 209.98 0.4604
19 112 1006.66 2.2074
193 19 674.94 1.4800
19 272 1406.11 3.0833
22 21 6956.00 15.2529
21 112 2388.33 5.2371
21 196 9050.79 19.8464
24 23 13719.30 30.0834
151 23 3606.21 7.9076
23 386 10875.01 23.8465
23 464 1170.80 2.5673
23 471 4962.34 10.8813
25 24 105.87 0.2321
26 24 1005.86 2.2056
210 24 5781.40 12.6773
24 360 6845.25 15.0101
28 27 2910.48 6.3820
27 68 4432.67 9.7199
27 382 1646.20 3.6097
30 29 4674.03 10.2491
454 29 5164.87 11.3254
29 491 956.37 2.0971
32 31 445.19 0.9762
31 85 1960.41 4.2987
104 31 1366.79 2.9971
31 189 392.88 0.8615
31 205 707.17 1.5507
31 218 1497.93 3.2846
267 31 2616.53 5.7375
34 33 1166.57 2.5580
83 33 610.69 1.3391
112 33 1732.42 3.7988
36 35 2007.02 4.4010
191 35 361.68 0.7931
35 399 1666.71 3.6547
38 37 5.94 0.0130
76 37 5002.38 10.9691
37 176 4992.73 10.9479
40 39 4046.29 8.8726
39 332 1036.63 2.2731
39 428 3083.94 6.7624
41 40 2668.04 5.8504
40 131 1378.92 3.0237
43 42 204.92 0.4494
81 42 3177.78 6.9682
42 191 2979.67 6.5337
45 44 93.98 0.2061
44 200 1757.68 3.8542
47 46 2033.53 4.4591
46 156 2041.31 4.4761
238 46 4047.76 8.8758
49 48 301.36 0.6608
50 48 516.38 1.1323
51 48 16.13 0.0354
52 48 1080.84 2.3700
358 48 264.25 0.5794
419 48 1092.91 2.3965
54 53 1614.62 3.5405
233 53 4387.26 9.6203
53 276 2863.39 6.2788
56 55 103.50 0.2270
55 66 1640.57 3.5974
55 500 1539.50 3.3758
58 57 12831.17 28.1359
58 57 12831.17 28.1359
57 123 19716.08 43.2330
440 57 5951.69 13.0507
59 58 6597.59 14.4670
58 253 4341.50 9.5199
322 58 7428.65 16.2894
322 58 7426.72 16.2851
61 60 9.19 0.0201
60 131 594.44 1.3035
466 60 591.33 1.2967
63 62 0.00 0.0000
64 62 1621.66 3.5559
287 62 1053.47 2.3100
62 457 3139.06 6.8833
66 65 4477.65 9.8185
80 65 2871.81 6.2972
65 297 2050.04 4.4953
67 66 209.81 0.4601
66 308 2634.79 5.7775
69 68 787.57 1.7270
68 295 2795.05 6.1289
68 314 2376.43 5.2110
71 70 1456.81 3.1945
72 70 454.70 0.9971
73 70 242.64 0.5321
216 70 2010.86 4.4094
70 490 993.14 2.1777
75 74 243.48 0.5339
74 189 514.08 1.1273
240 74 715.23 1.5683
77 76 320.77 0.7034
452 76 5296.94 11.6150
79 78 17.60 0.0386
78 328 7060.30 15.4816
78 387 7077.51 15.5194
81 80 6097.03 13.3694
82 80 15930.71 34.9325
407 80 7036.40 15.4292
149 81 2925.27 6.4145
84 83 2865.55 6.2835
83 342 2313.95 5.0740
86 85 807.53 1.7707
272 85 2686.34 5.8905
88 87 244.98 0.5372
87 141 14592.38 31.9978
87 143 13177.29 28.8949
423 87 2069.03 4.5369
90 89 272.56 0.5977
89 187 1537.15 3.3706
89 457 1803.96 3.9557
92 91 3026.63 6.6367
91 133 182.52 0.4002
94 93 240.73 0.5279
269 93 1885.93 4.1354
336 93 2121.30 4.6515
96 95 1581.18 3.4672
95 170 552.78 1.2121
95 369 2125.77 4.6613
98 97 1713.98 3.7584
181 97 1104.97 2.4230
358 97 1116.72 2.4487
100 99 622.11 1.3641
101 99 40.97 0.0898
99 465 484.40 1.0622
472 99 246.21 0.5399
103 102 1801.30 3.9498
102 308 2992.73 6.5624
419 102 1295.06 2.8398
105 104 3718.41 8.1536
104 191 2862.52 6.2769
107 106 739.11 1.6207
106 191 1854.98 4.0675
236 106 2560.46 5.6145
109 108 2205.45 4.8361
108 472 4497.96 9.8630
108 495 2293.48 5.0291
111 110 2851.61 6.2530
274 110 2167.62 4.7531
110 316 5193.33 11.3878
113 112 2137.07 4.6861
193 112 1419.12 3.1118
112 236 655.52 1.4374
112 283 560.50 1.2290
372 112 6793.97 14.8977
112 499 2550.80 5.5933
115 114 233.79 0.5126
116 114 143.26 0.3141
114 191 795.70 1.7448
114 312 440.29 0.9654
118 117 407.11 0.8927
117 426 2568.88 5.6330
490 117 2179.80 4.7798
120 119 1495.49 3.2793
119 129 1871.34 4.1034
253 119 3347.87 7.3411
122 121 2483.88 5.4466
121 276 986.94 2.1641
121 287 1651.42 3.6212
124 123 4048.78 8.8781
125 123 0.00 0.0000
126 123 0.00 0.0000
127 123 7165.55 15.7124
128 123 7450.44 16.3371
165 123 9219.83 20.2170
124 260 1588.42 3.4831
124 369 2463.14 5.4011
130 129 4891.80 10.7266
304 129 1159.10 2.5417
324 129 2050.98 4.4973
132 131 1161.94 2.5479
131 333 361.80 0.7933
134 133 1982.00 4.3461
133 226 248.27 0.5444
133 369 1410.17 3.0922
473 133 933.91 2.0479
136 135 1239.82 2.7187
135 287 1439.14 3.1557
135 457 2649.85 5.8105
138 137 3648.63 8.0006
137 214 2279.74 4.9990
137 426 2954.30 6.4781
137 426 2954.48 6.4785
140 139 12.04 0.0264
139 376 4118.92 9.0319
429 139 4125.49 9.0463
142 141 7522.93 16.4961
300 141 7450.26 16.3368
144 143 13386.68 29.3540
144 143 13386.68 29.3540
145 143 6933.69 15.2040
145 143 6933.69 15.2040
143 338 5589.63 12.2568
143 401 9375.73 20.5589
143 452 5781.81 12.6782
496 143 6309.47 13.8353
147 146 1267.43 2.7792
146 162 870.86 1.9096
146 195 443.49 0.9725
148 147 2659.17 5.8310
475 147 1576.70 3.4573
150 149 68.99 0.1513
293 149 2433.08 5.3352
489 149 444.41 0.9745
152 151 9628.33 21.1128
428 151 6617.09 14.5098
153 152 4268.55 9.3600
238 152 5854.58 12.8378
155 154 3112.91 6.8259
154 429 6506.99 14.2684
154 488 3395.12 7.4447
157 156 5016.88 11.0009
488 156 2982.12 6.5391
159 158 1065.66 2.3368
158 233 2192.76 4.8082
158 291 1131.29 2.4807
161 160 1199.22 2.6296
160 212 595.31 1.3054
160 399 653.16 1.4322
163 162 3334.81 7.3125
163 162 3334.74 7.3123
162 195 674.31 1.4786
162 220 17623.99 38.6455
162 232 9765.02 21.4125
164 163 4521.08 9.9137
200 163 2161.14 4.7389
166 165 14915.39 32.7061
167 165 8634.56 18.9336
168 165 8567.42 18.7864
169 165 17077.67 37.4475
165 386 10272.39 22.5250
166 228 9446.47 20.7140
260 166 739.12 1.6207
166 289 6217.05 13.6326
171 170 575.81 1.2626
170 226 85.64 0.1878
173 172 145.94 0.3200
174 172 339.27 0.7439
202 172 1128.08 2.4736
172 244 815.21 1.7876
176 175 2007.71 4.4025
175 220 2007.71 4.4025
177 176 403.24 0.8842
178 176 312.81 0.6859
318 176 2452.29 5.3773
180 179 663.52 1.4549
179 291 805.43 1.7661
345 179 428.92 0.9405
179 374 258.47 0.5668
182 181 3797.42 8.3269
181 487 4431.92 9.7182
184 183 1384.75 3.0365
183 247 907.91 1.9908
183 274 2167.12 4.7520
186 185 2113.54 4.6345
445 185 4055.86 8.8936
185 467 6118.67 13.4169
188 187 511.12 1.1208
187 454 1044.19 2.2897
190 189 644.75 1.4138
192 191 759.17 1.6647
191 205 1469.95 3.2233
191 293 1287.54 2.8233
194 193 1331.41 2.9195
196 195 13029.15 28.5700
197 195 12993.31 28.4914
198 195 9329.31 20.4571
371 195 9906.13 21.7219
199 196 3999.10 8.7691
201 200 410.71 0.9006
203 202 218.75 0.4797
204 202 2131.54 4.6740
380 202 3243.96 7.1133
206 205 138.69 0.3041
207 205 1040.49 2.2816
205 312 151.35 0.3319
209 208 3803.01 8.3391
322 208 6470.74 14.1889
356 208 2347.26 5.1470
211 210 1606.88 3.5235
488 210 4186.91 9.1810
213 212 4509.43 9.8882
212 399 443.16 0.9718
408 212 4808.71 10.5444
215 214 3928.26 8.6138
214 403 1714.38 3.7593
217 216 2559.85 5.6132
219 218 15.17 0.0333
218 399 1550.40 3.3997
221 220 9391.84 20.5942
222 220 554.05 1.2149
223 220 866.54 1.9001
224 220 11221.80 24.6069
224 220 11222.86 24.6092
225 220 11324.22 24.8315
225 220 11324.22 24.8315
220 471 20989.75 46.0258
221 447 6897.78 15.1253
221 478 1456.81 3.1945
221 493 3879.80 8.5075
227 226 252.01 0.5526
229 228 2541.50 5.5729
228 468 1020.66 2.2381
228 474 6087.43 13.3484
231 230 387.07 0.8488
285 230 5288.89 11.5974
233 232 13245.09 29.0435
262 232 6287.46 13.7870
232 371 1956.10 4.2893
232 453 3233.65 7.0907
234 233 2776.10 6.0874
233 486 3968.96 8.7030
236 235 3772.23 8.2717
235 271 3782.21 8.2935
235 407 7488.97 16.4216
237 236 124.83 0.2737
283 236 1099.53 2.4110
239 238 1811.89 3.9731
241 240 908.89 1.9930
240 272 1519.30 3.3315
243 242 3442.85 7.5494
242 287 4016.61 8.8075
242 320 1655.30 3.6297
245 244 547.72 1.2010
244 426 739.25 1.6210
247 246 896.61 1.9661
248 246 0.00 0.0000
249 246 0.00 0.0000
250 246 0.00 0.0000
250 246 0.00 0.0000
246 332 1047.13 2.2961
428 246 863.12 1.8926
246 471 231.83 0.5084
252 251 5631.19 12.3479
251 360 4075.28 8.9362
417 251 1592.69 3.4924
254 253 0.00 0.0000
255 253 0.00 0.0000
256 253 0.00 0.0000
257 253 0.00 0.0000
258 253 667.81 1.4644
259 253 384.45 0.8430
253 347 1984.11 4.3507
253 417 1200.71 2.6329
261 260 857.85 1.8811
263 262 18908.74 41.4626
264 262 0.00 0.0000
265 262 0.00 0.0000
266 263 2850.21 6.2499
263 276 5737.92 12.5820
263 326 5820.31 12.7626
263 476 1555.73 3.4114
263 491 4417.69 9.6870
268 267 5439.22 11.9270
312 267 3033.37 6.6515
270 269 314.44 0.6895
475 269 1583.69 3.4727
272 271 3822.68 8.3823
273 272 436.73 0.9576
275 274 0.00 0.0000
277 276 19.01 0.0417
278 276 267.41 0.5864
276 320 6179.37 13.5500
276 485 1169.87 2.5653
280 279 794.30 1.7417
279 340 816.65 1.7907
279 468 1586.23 3.4782
282 281 5503.18 12.0672
310 281 4233.62 9.2834
281 328 5892.29 12.9205
281 384 3001.05 6.5806
284 283 1513.19 3.3181
286 285 381.06 0.8356
285 487 4933.96 10.8191
288 287 2699.52 5.9194
287 415 2072.09 4.5436
287 476 1001.06 2.1951
290 289 170.50 0.3739
468 289 6057.61 13.2830
292 291 356.42 0.7816
294 293 1182.22 2.5923
296 295 361.81 0.7934
314 295 2397.98 5.2582
295 447 4327.91 9.4901
477 295 1222.85 2.6814
298 297 3744.58 8.2110
425 297 5923.09 12.9880
297 440 11582.00 25.3967
299 298 670.59 1.4705
298 426 3104.64 6.8078
301 300 897.86 1.9688
302 300 928.67 2.0364
303 300 4004.29 8.7805
405 300 1561.62 3.4243
300 478 1866.01 4.0917
305 304 310.22 0.6802
306 304 668.41 1.4657
307 304 3579.66 7.8494
309 308 417.06 0.9145
308 500 759.31 1.6650
311 310 2034.81 4.4619
358 310 1623.87 3.5608
313 312 574.88 1.2606
399 312 3206.31 7.0307
315 314 19.91 0.0437
317 316 172.84 0.3790
319 318 1987.57 4.3583
318 336 2123.41 4.6562
378 318 2209.16 4.8442
321 320 7033.82 15.4236
323 322 2123.43 4.6562
324 322 6929.90 15.1957
325 324 4893.16 10.7296
327 326 5635.72 12.3579
491 326 619.75 1.3590
329 328 1218.05 2.6709
331 330 2650.63 5.8122
333 330 4776.71 10.4743
449 330 2139.11 4.6906
333 332 5188.37 11.3769
332 471 5010.18 10.9862
334 333 23.00 0.0504
335 333 38.23 0.0838
337 336 261.44 0.5733
378 336 1620.00 3.5523
492 336 1431.35 3.1386
339 338 3489.27 7.6512
338 459 2157.07 4.7300
341 340 1136.50 2.4921
460 340 370.97 0.8134
343 342 115.99 0.2543
344 342 112.31 0.2463
346 345 2238.35 4.9082
345 415 1822.31 3.9959
348 347 4092.39 8.9737
347 417 2102.69 4.6107
350 349 842.36 1.8471
351 349 553.58 1.2139
352 349 818.15 1.7940
353 349 3283.52 7.2000
349 403 2571.38 5.6385
355 354 921.68 2.0210
421 354 2583.25 5.6645
479 354 3472.66 7.6148
357 356 2596.12 5.6927
356 382 392.27 0.8602
359 358 675.37 1.4809
361 360 0.00 0.0000
362 360 0.00 0.0000
363 360 0.00 0.0000
364 360 1426.88 3.1288
360 417 1395.85 3.0608
366 365 2643.22 5.7960
365 374 1394.69 3.0582
365 454 4037.28 8.8529
368 367 531.10 1.1646
367 461 4113.74 9.0205
367 468 3057.95 6.7054
370 369 1646.45 3.6103
372 371 8105.83 17.7743
373 372 1334.01 2.9252
375 374 1161.10 2.5460
377 376 4221.48 9.2568
376 429 6231.60 13.6645
376 460 1565.43 3.4326
376 469 5396.73 11.8338
376 495 807.13 1.7699
379 378 598.24 1.3118
381 380 470.41 1.0315
493 380 3700.01 8.1133
383 382 1391.88 3.0521
385 384 1714.80 3.7602
441 384 4656.66 10.2110
387 386 8874.28 19.4593
387 386 8874.28 19.4593
388 386 0.00 0.0000
389 386 0.00 0.0000
390 386 0.00 0.0000
391 386 0.00 0.0000
392 386 0.00 0.0000
393 386 0.00 0.0000
394 386 0.00 0.0000
395 386 0.00 0.0000
396 386 0.00 0.0000
440 386 28693.47 62.9184
397 387 1841.84 4.0387
398 387 772.93 1.6949
387 467 3947.54 8.6561
387 470 4144.44 9.0878
400 399 1061.23 2.3270
408 399 8097.92 17.7569
402 401 4123.51 9.0419
404 403 910.92 1.9975
406 405 242.77 0.5323
492 405 1332.86 2.9227
408 407 12895.56 28.2771
409 407 0.00 0.0000
410 407 3456.16 7.5786
411 407 3586.69 7.8648
412 407 3562.82 7.8125
413 407 2687.02 5.8920
414 407 0.00 0.0000
416 415 282.78 0.6201
418 417 2990.92 6.5584
417 468 8269.93 18.1341
417 469 871.90 1.9119
417 474 2375.31 5.2085
420 419 242.19 0.5311
422 421 2290.37 5.0223
424 423 7384.03 16.1915
423 496 8683.04 19.0400
426 425 5923.09 12.9880
427 426 664.61 1.4573
429 428 18248.33 40.0145
430 428 381.97 0.8376
431 428 677.03 1.4846
432 428 762.54 1.6721
433 428 394.88 0.8659
434 428 1081.49 2.3715
435 428 0.00 0.0000
436 428 0.00 0.0000
437 428 8749.68 19.1861
438 428 8699.51 19.0761
439 428 8805.60 19.3087
429 466 1401.98 3.0742
441 440 11747.69 25.7601
442 440 0.00 0.0000
443 440 0.00 0.0000
444 440 0.00 0.0000
446 445 4145.79 9.0908
445 461 667.78 1.4643
448 447 2675.09 5.8659
450 449 52.06 0.1142
451 449 102.70 0.2252
495 449 1990.80 4.3654
454 453 11576.46 25.3846
455 453 3709.96 8.1351
456 453 6277.05 13.7642
485 454 3470.86 7.6108
458 457 4885.52 10.7129
459 496 1495.70 3.2797
463 462 0.00 0.0000
462 465 2637.72 5.7839
462 477 2637.72 5.7839
465 464 3651.44 8.0068
471 464 4221.86 9.2576
470 467 2200.41 4.8250
472 471 6963.11 15.2685
476 486 989.93 2.1707
480 479 2483.71 5.4462
481 479 2029.45 4.4501
482 479 2122.67 4.6545
483 479 1080.13 2.3685
484 479 1313.36 2.8799
485 499 3156.46 6.9214
500 489 307.84 0.6750
494 493 0.00 0.0000
497 496 2964.06 6.4995
498 496 1081.15 2.3707

--- 4. Tổn thất toàn hệ thống ---
Tong ton that P: 91.2224 MW
Tong ton that Q: 1600.8541 MVAr

================================================
KẾT QUẢ PHẦN B (CÓ EVs)

Đã thêm 7 trạm sạc EV vào hệ thống.

• Nut 2: Them 150.00 kW
• Nut 4: Them 180.00 kW
• Nut 14: Them 120.00 kW
• Nut 17: Them 200.00 kW
• Nut 26: Them 150.00 kW
• Nut 29: Them 210.00 kW
• Nut 33: Them 200.00 kW

--- b) Điện áp nút khi có EVs (So sanh V_A va V_B) ---
Nut V_Cu (p.u) V_Moi (p.u) Chenh lech
1 1.0131 1.0131 -0.0000
2 1.0125 1.0125 -0.0000
3 1.0086 1.0086 -0.0000
4 1.0064 1.0064 -0.0000
5 1.0195 1.0195 0.0000
6 1.0194 1.0194 0.0000
7 1.0348 1.0348 -0.0000
8 1.0296 1.0296 -0.0000
9 1.0400 1.0400 0.0000
10 1.0280 1.0280 -0.0000
11 1.0280 1.0280 -0.0000
12 1.0096 1.0096 -0.0000
13 1.0090 1.0090 -0.0000
14 1.0321 1.0321 -0.0000
15 1.0247 1.0247 -0.0000
16 1.0400 1.0400 0.0000
17 1.0400 1.0400 0.0000
18 1.0400 1.0400 0.0000
19 1.0205 1.0205 -0.0000
20 1.0205 1.0205 -0.0000
21 1.0111 1.0111 -0.0000
22 1.0090 1.0090 -0.0000
23 1.0309 1.0309 -0.0000
24 1.0167 1.0167 -0.0000
25 1.0167 1.0167 -0.0000
26 1.0164 1.0164 -0.0000
27 1.0216 1.0216 -0.0000
28 1.0205 1.0205 -0.0000
29 1.0110 1.0110 -0.0000
30 1.0095 1.0095 -0.0000
31 1.0348 1.0348 -0.0000
32 1.0347 1.0347 -0.0000
33 1.0119 1.0119 -0.0000
34 1.0115 1.0115 -0.0000
35 1.0200 1.0200 -0.0000
36 1.0194 1.0194 -0.0000
37 1.0241 1.0241 -0.0000
38 1.0241 1.0241 -0.0000
39 1.0332 1.0332 -0.0000
40 1.0268 1.0268 -0.0000
41 1.0258 1.0258 -0.0000
42 1.0273 1.0273 -0.0000
43 1.0273 1.0273 -0.0000
44 1.0194 1.0193 -0.0000
45 1.0193 1.0193 -0.0000
46 1.0082 1.0082 -0.0000
47 1.0076 1.0076 -0.0000
48 1.0218 1.0218 -0.0000
49 1.0400 1.0400 0.0000
50 1.0400 1.0400 0.0000
51 1.0218 1.0218 -0.0000
52 1.0214 1.0214 -0.0000
53 1.0180 1.0180 -0.0000
54 1.0175 1.0175 -0.0000
55 1.0260 1.0260 -0.0000
56 1.0260 1.0260 -0.0000
57 1.0172 1.0171 -0.0000
58 1.0106 1.0106 -0.0000
59 1.0083 1.0082 -0.0000
60 1.0234 1.0234 -0.0000
61 1.0234 1.0234 -0.0000
62 1.0128 1.0127 -0.0000
63 1.0128 1.0127 -0.0000
64 1.0122 1.0122 -0.0000
65 1.0309 1.0309 -0.0000
66 1.0282 1.0282 -0.0000
67 1.0281 1.0281 -0.0000
68 1.0166 1.0166 0.0000
69 1.0164 1.0164 0.0000
70 1.0258 1.0258 -0.0000
71 1.0400 1.0400 0.0000
72 1.0400 1.0400 0.0000
73 1.0400 1.0400 0.0000
74 1.0273 1.0273 -0.0000
75 1.0272 1.0272 -0.0000
76 1.0226 1.0226 -0.0000
77 1.0225 1.0225 -0.0000
78 1.0184 1.0183 -0.0000
79 1.0183 1.0183 -0.0000
80 1.0346 1.0346 -0.0000
81 1.0296 1.0296 -0.0000
82 1.0400 1.0400 0.0000
83 1.0125 1.0125 -0.0000
84 1.0115 1.0115 -0.0000
85 1.0284 1.0284 -0.0000
86 1.0281 1.0281 -0.0000
87 1.0115 1.0115 -0.0000
88 1.0115 1.0115 -0.0000
89 1.0216 1.0216 -0.0000
90 1.0215 1.0215 -0.0000
91 1.0151 1.0151 -0.0000
92 1.0141 1.0141 -0.0000
93 1.0216 1.0216 -0.0000
94 1.0215 1.0215 -0.0000
95 1.0154 1.0154 -0.0000
96 1.0149 1.0149 -0.0000
97 1.0219 1.0219 -0.0000
98 1.0213 1.0213 -0.0000
99 1.0260 1.0260 -0.0000
100 1.0258 1.0258 -0.0000
101 1.0260 1.0260 -0.0000
102 1.0173 1.0173 -0.0000
103 1.0167 1.0167 -0.0000
104 1.0209 1.0209 -0.0000
105 1.0196 1.0196 -0.0000
106 1.0226 1.0226 -0.0000
107 1.0224 1.0224 -0.0000
108 1.0189 1.0189 -0.0000
109 1.0182 1.0182 -0.0000
110 1.0319 1.0319 0.0000
111 1.0310 1.0310 0.0000
112 1.0164 1.0163 -0.0000
113 1.0157 1.0157 -0.0000
114 1.0216 1.0216 -0.0000
115 1.0216 1.0216 -0.0000
116 1.0216 1.0216 -0.0000
117 1.0240 1.0240 -0.0000
118 1.0238 1.0238 -0.0000
119 0.9993 0.9993 -0.0000
120 0.9989 0.9989 -0.0000
121 1.0148 1.0148 -0.0000
122 1.0139 1.0139 -0.0000
123 1.0322 1.0322 -0.0000
124 1.0300 1.0300 -0.0000
125 1.0322 1.0322 -0.0000
126 1.0322 1.0322 -0.0000
127 1.0400 1.0400 0.0000
128 1.0400 1.0400 0.0000
129 0.9977 0.9977 -0.0000
130 0.9961 0.9961 -0.0000
131 1.0240 1.0240 -0.0000
132 1.0236 1.0236 -0.0000
133 1.0162 1.0162 -0.0000
134 1.0156 1.0156 -0.0000
135 1.0191 1.0191 -0.0000
136 1.0187 1.0187 -0.0000
137 1.0194 1.0194 -0.0000
138 1.0182 1.0182 -0.0000
139 1.0200 1.0200 -0.0000
140 1.0200 1.0200 -0.0000
141 1.0148 1.0148 -0.0000
142 1.0122 1.0122 -0.0000
143 1.0235 1.0235 0.0000
144 1.0400 1.0400 -0.0000
145 1.0400 1.0400 0.0000
146 1.0299 1.0299 -0.0000
147 1.0201 1.0201 -0.0000
148 1.0191 1.0191 -0.0000
149 1.0262 1.0262 -0.0000
150 1.0262 1.0262 -0.0000
151 1.0319 1.0318 -0.0000
152 1.0312 1.0312 -0.0000
153 1.0300 1.0300 -0.0000
154 1.0143 1.0143 -0.0000
155 1.0134 1.0134 -0.0000
156 0.9937 0.9937 -0.0000
157 0.9921 0.9921 -0.0000
158 1.0124 1.0123 -0.0000
159 1.0120 1.0120 -0.0000
160 1.0217 1.0217 -0.0000
161 1.0213 1.0213 -0.0000
162 1.0302 1.0302 -0.0000
163 1.0263 1.0263 -0.0000
164 1.0246 1.0246 -0.0000
165 1.0323 1.0323 -0.0000
166 1.0272 1.0272 -0.0000
167 1.0400 1.0400 0.0000
168 1.0400 1.0400 0.0000
169 1.0400 1.0400 0.0000
170 1.0149 1.0149 -0.0000
171 1.0147 1.0147 -0.0000
172 1.0158 1.0158 -0.0000
173 1.0157 1.0157 -0.0000
174 1.0157 1.0157 -0.0000
175 1.0365 1.0365 -0.0000
176 1.0278 1.0278 -0.0000
177 1.0276 1.0276 -0.0000
178 1.0277 1.0277 -0.0000
179 1.0060 1.0060 -0.0000
180 1.0058 1.0058 -0.0000
181 1.0281 1.0281 -0.0000
182 1.0268 1.0268 -0.0000
183 1.0291 1.0291 0.0000
184 1.0287 1.0287 0.0000
185 1.0133 1.0133 -0.0000
186 1.0127 1.0127 -0.0000
187 1.0206 1.0206 -0.0000
188 1.0205 1.0205 -0.0000
189 1.0301 1.0301 -0.0000
190 1.0299 1.0299 -0.0000
191 1.0223 1.0223 -0.0000
192 1.0220 1.0220 -0.0000
193 1.0175 1.0175 -0.0000
194 1.0170 1.0170 -0.0000
195 1.0301 1.0301 -0.0000
196 1.0209 1.0209 -0.0000
197 1.0400 1.0400 0.0000
198 1.0400 1.0400 0.0000
199 1.0196 1.0196 -0.0000
200 1.0217 1.0217 -0.0000
201 1.0216 1.0216 -0.0000
202 1.0120 1.0120 -0.0000
203 1.0119 1.0119 -0.0000
204 1.0113 1.0113 -0.0000
205 1.0242 1.0242 -0.0000
206 1.0241 1.0241 -0.0000
207 1.0238 1.0238 -0.0000
208 1.0170 1.0170 -0.0000
209 1.0158 1.0158 -0.0000
210 1.0075 1.0074 -0.0000
211 1.0069 1.0069 -0.0000
212 1.0218 1.0218 -0.0000
213 1.0203 1.0203 -0.0000
214 1.0164 1.0164 -0.0000
215 1.0150 1.0150 -0.0000
216 1.0167 1.0167 -0.0000
217 1.0159 1.0159 -0.0000
218 1.0286 1.0286 -0.0000
219 1.0286 1.0286 -0.0000
220 1.0375 1.0375 -0.0000
221 1.0261 1.0261 0.0000
222 1.0400 1.0400 0.0000
223 1.0400 1.0400 0.0000
224 1.0400 1.0400 0.0000
225 1.0400 1.0400 0.0000
226 1.0154 1.0154 -0.0000
227 1.0153 1.0153 -0.0000
228 1.0053 1.0053 -0.0000
229 1.0045 1.0045 -0.0000
230 1.0248 1.0248 -0.0000
231 1.0400 1.0400 0.0000
232 1.0302 1.0302 -0.0000
233 1.0203 1.0203 -0.0000
234 1.0194 1.0194 -0.0000
235 1.0314 1.0314 -0.0000
236 1.0244 1.0244 -0.0000
237 1.0244 1.0244 -0.0000
238 1.0199 1.0199 -0.0000
239 1.0193 1.0193 -0.0000
240 1.0253 1.0253 -0.0000
241 1.0250 1.0250 -0.0000
242 1.0221 1.0221 -0.0000
243 1.0209 1.0209 -0.0000
244 1.0197 1.0197 -0.0000
245 1.0195 1.0195 -0.0000
246 1.0351 1.0351 -0.0000
247 1.0302 1.0302 0.0000
248 1.0351 1.0351 -0.0000
249 1.0351 1.0351 -0.0000
250 1.0351 1.0351 -0.0000
251 0.9968 0.9967 -0.0000
252 0.9947 0.9947 -0.0000
253 1.0020 1.0020 -0.0000
254 1.0020 1.0020 -0.0000
255 1.0020 1.0020 -0.0000
256 1.0020 1.0020 -0.0000
257 1.0020 1.0020 -0.0000
258 1.0400 1.0400 0.0000
259 1.0019 1.0019 -0.0000
260 1.0280 1.0280 -0.0000
261 1.0277 1.0277 -0.0000
262 1.0296 1.0295 -0.0000
263 1.0240 1.0240 -0.0000
264 1.0296 1.0295 -0.0000
265 1.0296 1.0295 -0.0000
266 1.0231 1.0231 -0.0000
267 1.0182 1.0182 -0.0000
268 1.0165 1.0165 -0.0000
269 1.0207 1.0207 -0.0000
270 1.0206 1.0206 -0.0000
271 1.0308 1.0308 -0.0000
272 1.0257 1.0257 -0.0000
273 1.0256 1.0256 -0.0000
274 1.0307 1.0307 0.0000
275 1.0307 1.0307 0.0000
276 1.0186 1.0185 -0.0000
277 1.0185 1.0185 -0.0000
278 1.0185 1.0185 -0.0000
279 1.0069 1.0069 -0.0000
280 1.0066 1.0066 -0.0000
281 1.0078 1.0078 -0.0000
282 1.0060 1.0060 -0.0000
283 1.0190 1.0190 -0.0000
284 1.0184 1.0184 -0.0000
285 1.0243 1.0243 -0.0000
286 1.0242 1.0241 -0.0000
287 1.0169 1.0169 -0.0000
288 1.0160 1.0160 -0.0000
289 1.0213 1.0213 -0.0000
290 1.0212 1.0212 -0.0000
291 1.0091 1.0091 -0.0000
292 1.0090 1.0090 -0.0000
293 1.0224 1.0224 -0.0000
294 1.0221 1.0221 -0.0000
295 1.0164 1.0164 0.0000
296 1.0163 1.0163 0.0000
297 1.0260 1.0260 -0.0000
298 1.0256 1.0256 -0.0000
299 1.0254 1.0254 -0.0000
300 1.0241 1.0241 -0.0000
301 1.0400 1.0400 0.0000
302 1.0400 1.0400 0.0000
303 1.0226 1.0226 -0.0000
304 1.0037 1.0037 -0.0000
305 1.0400 1.0400 -0.0000
306 1.0400 1.0400 0.0000
307 1.0026 1.0026 -0.0000
308 1.0221 1.0221 -0.0000
309 1.0220 1.0220 -0.0000
310 1.0173 1.0173 -0.0000
311 1.0167 1.0167 -0.0000
312 1.0216 1.0216 -0.0000
313 1.0214 1.0214 -0.0000
314 1.0166 1.0166 0.0000
315 1.0166 1.0166 0.0000
316 1.0241 1.0241 0.0000
317 1.0240 1.0240 0.0000
318 1.0279 1.0279 -0.0000
319 1.0400 1.0400 0.0000
320 1.0117 1.0117 -0.0000
321 1.0093 1.0093 -0.0000
322 1.0069 1.0069 -0.0000
323 1.0062 1.0062 -0.0000
324 0.9995 0.9995 -0.0000
325 0.9979 0.9979 -0.0000
326 1.0156 1.0156 -0.0000
327 1.0139 1.0139 -0.0000
328 1.0086 1.0086 -0.0000
329 1.0083 1.0083 -0.0000
330 1.0189 1.0189 -0.0000
331 1.0180 1.0180 -0.0000
332 1.0329 1.0329 -0.0000
333 1.0243 1.0243 -0.0000
334 1.0243 1.0243 -0.0000
335 1.0243 1.0243 -0.0000
336 1.0236 1.0236 -0.0000
337 1.0235 1.0235 -0.0000
338 1.0170 1.0170 0.0000
339 1.0159 1.0159 0.0000
340 1.0063 1.0062 -0.0000
341 1.0059 1.0059 -0.0000
342 1.0215 1.0215 -0.0000
343 1.0214 1.0214 -0.0000
344 1.0214 1.0214 -0.0000
345 1.0040 1.0039 -0.0000
346 1.0033 1.0033 -0.0000
347 0.9977 0.9977 -0.0000
348 0.9965 0.9965 -0.0000
349 1.0314 1.0314 -0.0000
350 1.0400 1.0400 0.0000
351 1.0400 1.0400 0.0000
352 1.0400 1.0400 0.0000
353 1.0400 1.0400 0.0000
354 1.0229 1.0229 -0.0000
355 1.0226 1.0226 -0.0000
356 1.0140 1.0140 -0.0000
357 1.0132 1.0132 -0.0000
358 1.0198 1.0198 -0.0000
359 1.0196 1.0196 -0.0000
360 1.0073 1.0073 -0.0000
361 1.0073 1.0073 -0.0000
362 1.0073 1.0073 -0.0000
363 1.0073 1.0073 -0.0000
364 1.0069 1.0069 -0.0000
365 1.0096 1.0096 -0.0000
366 1.0087 1.0087 -0.0000
367 1.0176 1.0176 -0.0000
368 1.0175 1.0175 -0.0000
369 1.0210 1.0210 -0.0000
370 1.0205 1.0205 -0.0000
371 1.0292 1.0292 -0.0000
372 1.0208 1.0208 -0.0000
373 1.0204 1.0204 -0.0000
374 1.0068 1.0068 -0.0000
375 1.0065 1.0064 -0.0000
376 1.0118 1.0118 -0.0000
377 1.0106 1.0106 -0.0000
378 1.0266 1.0266 -0.0000
379 1.0264 1.0264 -0.0000
380 1.0165 1.0165 0.0000
381 1.0163 1.0163 0.0000
382 1.0161 1.0161 -0.0000
383 1.0156 1.0156 -0.0000
384 1.0085 1.0085 -0.0000
385 1.0080 1.0079 -0.0000
386 1.0295 1.0295 -0.0000
387 1.0243 1.0243 -0.0000
388 1.0295 1.0295 -0.0000
389 1.0295 1.0295 -0.0000
390 1.0295 1.0295 -0.0000
391 1.0295 1.0295 -0.0000
392 1.0295 1.0295 -0.0000
393 1.0295 1.0295 -0.0000
394 1.0295 1.0295 -0.0000
395 1.0295 1.0295 -0.0000
396 1.0295 1.0295 -0.0000
397 1.0237 1.0237 -0.0000
398 1.0240 1.0240 -0.0000
399 1.0246 1.0246 -0.0000
400 1.0242 1.0242 -0.0000
401 1.0126 1.0126 0.0000
402 1.0113 1.0113 0.0000
403 1.0246 1.0246 -0.0000
404 1.0243 1.0243 -0.0000
405 1.0207 1.0207 -0.0000
406 1.0207 1.0207 -0.0000
407 1.0356 1.0356 -0.0000
408 1.0295 1.0295 -0.0000
409 1.0356 1.0356 -0.0000
410 1.0400 1.0400 0.0000
411 1.0400 1.0400 0.0000
412 1.0400 1.0400 0.0000
413 1.0400 1.0400 0.0000
414 1.0356 1.0356 -0.0000
415 1.0117 1.0117 -0.0000
416 1.0116 1.0116 -0.0000
417 1.0049 1.0049 -0.0000
418 1.0039 1.0039 -0.0000
419 1.0191 1.0191 -0.0000
420 1.0190 1.0190 -0.0000
421 1.0141 1.0141 -0.0000
422 1.0134 1.0134 -0.0000
423 1.0069 1.0069 -0.0000
424 1.0042 1.0042 -0.0000
425 1.0257 1.0257 -0.0000
426 1.0250 1.0250 -0.0000
427 1.0248 1.0248 -0.0000
428 1.0347 1.0347 -0.0000
429 1.0245 1.0245 -0.0000
430 1.0400 1.0400 0.0000
431 1.0400 1.0400 0.0000
432 1.0400 1.0400 0.0000
433 1.0400 1.0400 0.0000
434 1.0400 1.0400 0.0000
435 1.0347 1.0347 -0.0000
436 1.0347 1.0347 -0.0000
437 1.0400 1.0400 0.0000
438 1.0400 1.0400 -0.0000
439 1.0400 1.0400 0.0000
440 1.0219 1.0219 -0.0000
441 1.0131 1.0131 -0.0000
442 1.0219 1.0219 -0.0000
443 1.0219 1.0219 -0.0000
444 1.0219 1.0219 -0.0000
445 1.0110 1.0110 -0.0000
446 1.0096 1.0096 -0.0000
447 1.0165 1.0165 0.0000
448 1.0156 1.0156 0.0000
449 1.0157 1.0157 -0.0000
450 1.0157 1.0157 -0.0000
451 1.0157 1.0157 -0.0000
452 1.0213 1.0213 -0.0000
453 1.0318 1.0318 -0.0000
454 1.0165 1.0165 -0.0000
455 1.0400 1.0400 0.0000
456 1.0400 1.0400 0.0000
457 1.0268 1.0268 -0.0000
458 1.0400 1.0400 0.0000
459 1.0162 1.0162 0.0000
460 1.0085 1.0085 -0.0000
461 1.0062 1.0062 -0.0000
462 1.0234 1.0234 -0.0000
463 1.0234 1.0234 -0.0000
464 1.0320 1.0320 -0.0000
465 1.0275 1.0275 -0.0000
466 1.0226 1.0226 -0.0000
467 1.0183 1.0183 -0.0000
468 1.0088 1.0088 -0.0000
469 0.9972 0.9972 -0.0000
470 1.0195 1.0195 -0.0000
471 1.0327 1.0327 -0.0000
472 1.0256 1.0256 -0.0000
473 1.0196 1.0196 -0.0000
474 0.9908 0.9908 -0.0000
475 1.0205 1.0204 -0.0000
476 1.0186 1.0186 -0.0000
477 1.0171 1.0171 0.0000
478 1.0251 1.0251 0.0000
479 1.0341 1.0341 -0.0000
480 1.0400 1.0400 0.0000
481 1.0400 1.0400 -0.0000
482 1.0400 1.0400 0.0000
483 1.0400 1.0400 0.0000
484 1.0400 1.0400 0.0000
485 1.0137 1.0137 -0.0000
486 1.0166 1.0166 -0.0000
487 1.0254 1.0254 -0.0000
488 0.9966 0.9966 -0.0000
489 1.0253 1.0253 -0.0000
490 1.0238 1.0238 -0.0000
491 1.0182 1.0182 -0.0000
492 1.0169 1.0169 -0.0000
493 1.0209 1.0209 0.0000
494 1.0209 1.0209 0.0000
495 1.0133 1.0133 -0.0000
496 1.0205 1.0205 0.0000
497 1.0400 1.0400 0.0000
498 1.0400 1.0400 0.0000
499 1.0054 1.0054 -0.0000
500 1.0238 1.0238 -0.0000

--- c) Dòng điện nhánh khi có EVs (Hien thi 5 nhanh dau tien) ---
(Xem chi tiet trong bien I_abs_B)
Tu Nut Den Nut Dong Moi (A)
2 1 2001.62
44 1 1670.22
421 1 297.44
4 3 6230.58
62 3 857.92

--- d) Tổn thất khi có EVs ---
Tong ton that P: 91.2463 MW (Tang them 0.0239 MW)
Tong ton that Q: 1601.2970 MVAr

[SharkpJIM]

% --- CHƯƠNG TRÌNH MÔ PHỎNG LƯỚI ĐIỆN LỚN VỚI 500 EVs ---
% Tự động xuất kết quả ra file Excel/CSV để báo cáo
clear; clc;
define_constants; % Nạp thư viện MATPOWER

% ============================================================
% 1. CẤU HÌNH NGƯỜI DÙNG
% ============================================================
CASE_NAME = 'case_ACTIVSg2000'; % Tên lưới điện (2000 bus)
NUM_EVS = 500; % Số lượng trạm sạc muốn lắp
P_EV_KW = 100; % Công suất mỗi trạm (kW)
Q_EV_KVAR = 0; % Công suất phản kháng (kVAr)
V_BASE_KV = 12.66; % Điện áp cơ sở (để tính ra kV thực tế)

fprintf('Đang khởi tạo lưới %s với %d trạm sạc...\n', CASE_NAME, NUM_EVS);

% ============================================================
% 2. PHẦN A: CHẠY LƯỚI CƠ SỞ (BASE CASE)
% ============================================================
try
mpc = eval(CASE_NAME);
catch
error('Không tìm thấy file dữ liệu lưới "%s".', CASE_NAME);
end

fprintf('\n>>> Đang chạy phân bố công suất Phần A...\n');
rA = runpf(mpc, mpoption('verbose', 0, 'out.all', 0));

if ~rA.success
error('Lỗi: Lưới cơ sở không hội tụ!');
end

% --- TÍNH TOÁN KẾT QUẢ PHẦN A ---
% 1. Điện áp
bus_ids = rA.bus(:, BUS_I);
Vm_A = rA.bus(:, VM);
Va_A = rA.bus(:, VA);
Vm_kV_A = Vm_A * V_BASE_KV;

% 2. Dòng điện nhánh
% Lưu ý: Cần lấy đúng index của nút để tra điện áp
S_flow_A = sqrt(rA.branch(:, PF).^2 + rA.branch(:, QF).^2); % MVA
% Tạo map để tra cứu điện áp nút From nhanh chóng
map_V = containers.Map(bus_ids, Vm_A);
V_from_list = zeros(size(rA.branch, 1), 1);
for k = 1:size(rA.branch, 1)
f_node = rA.branch(k, F_BUS);
V_from_list(k) = map_V(f_node);
end
I_pu_A = (S_flow_A / rA.baseMVA) ./ V_from_list;

% 3. Tổn thất
loss_A = get_losses(rA);
Ploss_A = sum(real(loss_A));
Qloss_A = sum(imag(loss_A));

% --- XUẤT FILE PHẦN A ---
fprintf('-> Đang xuất file kết quả Phần A...\n');

% File 1: Điện áp Phần A
Header_Bus = {'Bus_ID', 'V_pu', 'Angle_deg', 'V_kV'};
Data_Bus_A = table(bus_ids, Vm_A, Va_A, Vm_kV_A, 'VariableNames', Header_Bus);
writetable(Data_Bus_A, 'KetQua_PhanA_DienAp.csv');

% File 2: Dòng điện Phần A
Header_Branch = {'From_Bus', 'To_Bus', 'I_pu', 'S_MVA_Flow'};
Data_Branch_A = table(rA.branch(:, F_BUS), rA.branch(:, T_BUS), I_pu_A, S_flow_A, 'VariableNames', Header_Branch);
writetable(Data_Branch_A, 'KetQua_PhanA_DongDien.csv');

fprintf(' [OK] Đã lưu: KetQua_PhanA_DienAp.csv & KetQua_PhanA_DongDien.csv\n');
fprintf(' Tổng tổn thất P (A): %.4f MW\n', Ploss_A);

% ============================================================
% 3. PHẦN B: LẮP ĐẶT 500 EVS & TÍNH TOÁN
% ============================================================
fprintf('\n>>> Đang thiết lập %d trạm sạc (Random)...\n', NUM_EVS);

mpc_EV = mpc;
% Chỉ chọn các nút Tải (PQ Bus) để gắn EV (tránh gắn vào nút nguồn)
candidate_buses = mpc.bus(mpc.bus(:, BUS_TYPE) == 1, BUS_I);

% Random trọn 500 vị trí
rng('shuffle'); % Đảo ngẫu nhiên
shuffled_buses = candidate_buses(randperm(length(candidate_buses)));
if length(shuffled_buses) < NUM_EVS
warning('Số lượng nút tải (%d) ít hơn số EV yêu cầu (%d). Sẽ gắn full nút tải.', length(shuffled_buses), NUM_EVS);
selected_buses = shuffled_buses;
else
selected_buses = shuffled_buses(1:NUM_EVS);
end

% Cập nhật tải EV vào lưới
% Dùng containers.Map để update nhanh mà không cần vòng lặp find (tránh chậm với 2000 bus)
bus_idx_map = containers.Map(mpc.bus(:, BUS_I), 1:size(mpc.bus, 1));

for i = 1:length(selected_buses)
b_id = selected_buses(i);
row = bus_idx_map(b_id); % Tìm dòng trong ma trận
mpc_EV.bus(row, PD) = mpc_EV.bus(row, PD) + P_EV_KW/1000;
mpc_EV.bus(row, QD) = mpc_EV.bus(row, QD) + Q_EV_KVAR/1000;
end

% --- CHẠY PHÂN BỐ CÔNG SUẤT PHẦN B ---
rB = runpf(mpc_EV, mpoption('verbose', 0, 'out.all', 0));

if ~rB.success
fprintf('\n[CẢNH BÁO] Hệ thống KHÔNG HỘI TỤ sau khi thêm 500 EVs!\n');
fprintf('Có thể do sụt áp quá lớn. Vui lòng giảm công suất hoặc giảm số lượng EV.\n');
else
fprintf('-> Chạy Phần B thành công!\n');

% --- TÍNH TOÁN KẾT QUẢ PHẦN B ---
% 1. Điện áp & So sánh
Vm_B = rB.bus(:, VM);
V_drop = Vm_A - Vm_B; % Độ sụt áp

% 2. Dòng điện Phần B
S_flow_B = sqrt(rB.branch(:, PF).^2 + rB.branch(:, QF).^2);

% Update map điện áp mới để tính dòng
map_V_B = containers.Map(rB.bus(:, BUS_I), Vm_B);
V_from_list_B = zeros(size(rB.branch, 1), 1);
for k = 1:size(rB.branch, 1)
    V_from_list_B(k) = map_V_B(rB.branch(k, F_BUS));
end
I_pu_B = (S_flow_B / rB.baseMVA) ./ V_from_list_B;

% 3. Tổn thất
loss_B = get_losses(rB);
Ploss_B = sum(real(loss_B));
Qloss_B = sum(imag(loss_B));

% --- XUẤT FILE PHẦN B ---
fprintf('-> Đang xuất file kết quả Phần B...\n');

% File 3: Điện áp Phần B (Kèm so sánh)
% Đánh dấu nút nào có EV
Is_EV_Node = ismember(bus_ids, selected_buses); 
Header_Bus_B = {'Bus_ID', 'V_Old_pu', 'V_New_pu', 'Voltage_Drop_pu', 'Has_EV_Station'};
Data_Bus_B = table(bus_ids, Vm_A, Vm_B, V_drop, Is_EV_Node, 'VariableNames', Header_Bus_B);
writetable(Data_Bus_B, 'KetQua_PhanB_SoSanhDienAp.csv');

% File 4: Dòng điện Phần B
Header_Branch_B = {'From_Bus', 'To_Bus', 'I_New_pu', 'I_Increase_pu'};
Data_Branch_B = table(rB.branch(:, F_BUS), rB.branch(:, T_BUS), I_pu_B, I_pu_B - I_pu_A, 'VariableNames', Header_Branch_B);
writetable(Data_Branch_B, 'KetQua_PhanB_DongDien.csv');

% File 5: Danh sách vị trí EV (để bạn biết nó đặt ở đâu)
writematrix(selected_buses, 'DanhSach_ViTri_500EV.csv');

fprintf('   [OK] Đã lưu 3 file CSV cho Phần B.\n');

% --- HIỂN THỊ TÓM TẮT TRÊN MÀN HÌNH ---
fprintf('\n================================================\n');
fprintf('          TỔNG HỢP KẾT QUẢ (SUMMARY)\n');
fprintf('================================================\n');
fprintf('1. TỔN THẤT HỆ THỐNG:\n');
fprintf('   - P Loss: %.4f MW  -> %.4f MW (Tăng %.4f MW)\n', Ploss_A, Ploss_B, Ploss_B - Ploss_A);
fprintf('   - Q Loss: %.4f MVAr -> %.4f MVAr\n', Qloss_A, Qloss_B);

fprintf('\n2. ĐIỆN ÁP (Top 5 nút sụt áp nhiều nhất):\n');
[~, idx_sort_drop] = sort(V_drop, 'descend');
fprintf('   Bus ID\tV_Cu\tV_Moi\tSut Ap\n');
for k = 1:5
    idx = idx_sort_drop(k);
    fprintf('   %d\t\t%.4f\t%.4f\t%.4f\n', bus_ids(idx), Vm_A(idx), Vm_B(idx), V_drop(idx));
end

fprintf('\n3. DÒNG ĐIỆN (Top 5 nhánh tải cao nhất):\n');
[~, idx_sort_I] = sort(I_pu_B, 'descend');
fprintf('   From\tTo\tI (p.u.)\n');
for k = 1:5
    idx = idx_sort_I(k);
    fprintf('   %d\t%d\t%.4f\n', rB.branch(idx, F_BUS), rB.branch(idx, T_BUS), I_pu_B(idx));
end

fprintf('\n[THÀNH CÔNG] Bạn hãy mở thư mục hiện tại để xem 5 file CSV chi tiết.\n');

end

[SharkpJIM]

Đã giải quyết được vấn đề cho case_ACTIVSg2000

[rdzman]

I'm sorry, but it wasn't clear to me (using Google Translate) whether you had a specific question.