% 加载33节点电力系统数据文件
mpc = loadcase('case33');

% 修改电力系统数据文件，将18节点的电源改为光伏发电机
mpc.gen(18, 1) = struct('bus', 18, 'pg', 0.5, 'qg', 0, 'vg', 1.0, 'mbase', 100, 'status', 1, 'gen_type', 2, 'r', 0.0, 'x', 0.0, 'b', 0.0, 'mpc', 0, 'contype', 1, 'opf', 0, 'dc', 0);

% 运行最优潮流分析
results = runopf(mpc);

% 分析最优潮流结果，输出各个节点的电压和潮流分布
disp(results.bus(:, [1, 8, 9, 13]));

% 绘制电力系统的潮流分布图
plotpf(mpc, results);