Matpower 39 节点频率变化下的有功-频率下垂特性控制措施及曲线绘制

以下是基于 matpower 39 节点的有功-频率下垂特性控制措施的 MATLAB 代码示例：

% 导入 matpower 39 节点数据
mpc = loadcase('case39');
baseMVA = mpc.baseMVA;

% 设置频率范围及其对应的有功功率下垂值
frequencies = [50.1 50.2 50.3 50.4 50.5]; % 频率范围
deltaP = [-2 -1 0 1 2]; % 有功功率下垂值

% 初始化结果变量
results = zeros(length(frequencies), length(deltaP));

% 循环计算不同频率和有功功率下垂值下的电力系统结果
for i = 1:length(frequencies)
    for j = 1:length(deltaP)
        % 设置频率
        mpc.bus(:, 3) = mpc.bus(:, 3) * (frequencies(i) / 50);
        
        % 设置有功功率下垂值
        mpc.gen(:, 2) = mpc.gen(:, 2) - deltaP(j) / baseMVA;
        
        % 运行潮流计算
        results(i, j) = runpf(mpc).success;
    end
end

% 绘制有功-频率下垂曲线
figure;
hold on;
for i = 1:length(deltaP)
    plot(frequencies, results(:, i), 'o-');
end
hold off;
xlabel('频率 (Hz)');
ylabel('有功功率下垂');
legend(cellstr(num2str(deltaP', 'deltaP = %d MW')));
grid on;

该代码首先导入 matpower 39 节点数据，并设置了频率范围和对应的有功功率下垂值。然后，通过循环计算不同频率和有功功率下垂值下的电力系统结果，并将结果存储在 results 变量中。最后，绘制出有功-频率下垂曲线。

请注意，该代码假设 case39 文件是位于当前工作目录中的 MATPOWER 数据文件，您可以根据需要更改数据文件的路径和名称。另外，该代码使用 MATPOWER 的 runpf 函数运行潮流计算，并使用 MATLAB 的 plot 函数绘制曲线。