CHP可行域图形复现matlab代码

以下是使用Matlab绘制CHP可行域图形的代码：

% CHP可行域图形
clear;clc;

% 参数设置
P_e = 0:0.1:100; % 电功率范围
P_th = 0:0.1:100; % 热功率范围
eta_th = 0.9; % 热功率转换效率
eta_e = 0.3; % 电功率转换效率
P_th_min = 10; % 最小热功率
P_th_max = 90; % 最大热功率
P_e_min = 10; % 最小电功率
P_e_max = 80; % 最大电功率

% 计算CHP可行域
P_th_CHP = P_e * eta_e / (eta_e + eta_th); % CHP热功率
P_e_CHP = P_th_CHP * (eta_e + eta_th) / eta_e; % CHP电功率
index1 = find(P_th_CHP >= P_th_min & P_th_CHP <= P_th_max); % 筛选符合最小和最大热功率的点
index2 = find(P_e_CHP >= P_e_min & P_e_CHP <= P_e_max); % 筛选符合最小和最大电功率的点
index = intersect(index1, index2); % 取交集
P_th_CHP_feasible = P_th_CHP(index); % CHP可行热功率
P_e_CHP_feasible = P_e_CHP(index); % CHP可行电功率

% 绘制CHP可行域图形
figure;
plot(P_e, P_th, 'k--'); % 绘制边界线
hold on;
plot(P_e_CHP_feasible, P_th_CHP_feasible, 'r', 'LineWidth', 2); % 绘制CHP可行域
xlabel('电功率/kW');
ylabel('热功率/kW');
title('CHP可行域');
legend('边界线', '可行域');

运行上述代码可以得到如下的CHP可行域图形：