matpower修改自导纳重新计算潮流的方法

修改自导纳并重新计算潮流
在电力系统分析工具 MatPower 中，可以通过修改节点的自导纳矩阵 (self-admittance matrix) 来调整系统的参数，并进而重新计算潮流分布。具体操作涉及以下几个方面：
准备工作
为了确保能够顺利进行自导纳的修改以及后续的潮流计算，在开始之前需加载所需的案例数据。
% 加载测试用例
mpc = loadcase('case9');
登录后复制
修改自导纳
对于特定母线上的自导纳值可以基于实际需求进行更改。假设要增加编号为 bus_id 的母线上电抗器的影响，则可通过如下方式实现：
bus_id = 3; % 假设目标母线ID为3
additional_shunt_admittance = 0.02 + j*0.05; % 新增的串联阻抗值(实部G,虚部B)

for i = 1:length(mpc.bus)
if mpc.bus(i, BUS_I) == bus_id
mpc.bus(i, GS) = real(additional_shunt_admittance); % 更新实数部分(G)
mpc.bus(i, BS) = imag(additional_shunt_admittance)*1e6; % 更新虚数部分(B)，单位转换成uS
break;
end
end
登录后复制
这里需要注意的是，Matpower中的BS是以微西门子(uS)作为默认单位存储的，因此当输入新的虚部时需要乘以(1 \times 10^{6})[^1]。
执行潮流计算
完成上述配置之后就可以调用内置函数执行潮流求解:
options = mpoption('out.all', 0);
results = runpf(mpc, options);
disp(results.gen(:, PG)); % 显示发电机有功功率输出结果
登录后复制
此段代码会根据最新的电网模型再次运行直流或交流潮流算法，并返回相应的计算结果。