matpower怎样调用求解器求解交流潮流

MatPOWER是一个用于电力系统的开源软件包，它主要用于研究电力系统调度、规划和运行等问题。在MatPOWER中，求解交流潮流通常通过调用内置的`runpf()`函数来完成。这个函数基于PQ分解算法，可以计算出电力网络在给定运行条件下的稳态运行状态。

以下是基本的调用流程：

% 导入MatPOWER库
addpath('toolbox/matpower'); % 如果MatPOWER安装在默认路径下

% 加载测试案例数据
case = 'case9'; % 更改为你需要分析的案例，如case9, case14, 等等
[baseM, baseT, bus, gen, branch] = ppci; % 获取基础信息

% 设置运行条件（例如，指定哪些发电机为PQ节点）
p_set = zeros(length(gen), 1); % 初始功率设定
q_set = zeros(length(gen), 1); % 初始无功设定
for i = 1:length(gen)
    if strcmp(gen.type(i), 'PQ') % 只有PQ类型的发电机参与优化
        p_set(i) = gen.Pg(i);
        q_set(i) = gen.Qg(i);
    end
end

% 调用runpf()求解
options = pfOptions(); % 初始化选项结构体
options.dual = true; % 显示灵敏度信息（可选）
options.verbose = false; % 静默模式（可选）
[t, tf, sol] = runpf(baseM, baseT, bus, gen, branch, p_set, q_set, options); % 求解并获取结果

% 查看结果，sol包含了电压、功率等各种变量
disp(sol);