MATPOWER中如何计算各线路上与节点上的电流？

MATPOWER是一个用于研究电力系统运行情况的开源软件包，它基于MATLAB环境编写。要计算线路的电流和节点的功率注入，MATPOWER通常涉及到两个关键步骤：首先设置网络模型，然后解决直流或交流潮流方程。

1. 设置网络模型：通过`pfs`函数创建一个系统描述数据结构，包括网络的连接、发电机和负荷等信息。例如：

baseMVA = 100; % 基准容量
 bus = ...; % 节点数据
 gen = ...; % 发电机数据
branch = ...; % 线路数据
[PF, busf, branchf] = pfs(baseMVA, bus, gen, branch);

2. 计算电流和功率：`PF`矩阵包含了线路上的复功率分布，可以直接读取其中的元素来获取电流。每个元素`PF(i,j)`代表从第i个节点流入第j个节点的有功功率（实部），无功功率（虚部）。若需要电流，则可以利用欧姆定律和基尔霍夫定律转换得到。对于三相系统，还需要考虑相间功率分量。

% 电压幅值和角度
Vbus = bus(:,2); % 电压幅值
theta = bus(:,4); % 电压相角

% 计算每条线路的复电流 I = P / (V1*conj(V2) - V2*conj(V1)) * (cos(theta2 - theta1) + j*sin(theta2 - theta1))
Ibranch = PF ./ conj((Vbus(branch(:,1)).^2 - Vbus(branch(:,2)).^2) .* exp(1i*(theta(bus(branch(:,1),4)) - theta(bus(branch(:,2),4)))));