冷热电联供系统CCHP经济优化运行多能源系统优化MATLAB程序 （1）该程序为冷热电联供系统CCHP经济优化运行

冷热电联供系统CCHP经济优化运行多能源系统优化MATLAB程序 （1）该程序为冷热电联供系统CCHP经济优化运行，多能源系统优化，硕士学位论文源程序，配有该论文。 （2）通过该程序可得到冷热电联供系统的经济运行的相关数据，为能量合理优化调度提供依据和参考。程序与论文包含的内容有冷热电联供系统目标函数与约束条件的建立、粒子群算法求解、算例仿真。 本程序有对应说明，注释全面，可有效帮助大家快速入门该研究领域，减少大量学习成本与时间，让你少走弯路，并且对论文写作与架构的学习都有很大帮助。同时大家也可在该程序基础上进一步拓展程序功能与内容，发表新论文。

该程序实现了一个冷热电联产系统的夏季供需调度优化模型，采用粒子群优化算法对燃气轮机出力与电网购电进行优化，以最小化系统运行成本。以下是各文件的功能说明：

---

1. `calobjvalue.m`

功能：计算粒子群中每个粒子的目标函数值（总运行成本）。

输入：

• pop：粒子群位置矩阵（3维：2×24×80）

输出：

• objvalue：每个粒子的总运行成本（1×80）

计算内容：

• cost1：燃气轮机的燃料成本、运行成本、排污成本
• cost2：从大电网购电的成本（基于分时电价）
• objvalue = cost1 + cost2

---

2. `calvalue.m`

功能：计算给定最优粒子位置 best 对应的总运行成本。

输入：

• best：最优粒子位置（2×24）

输出：

• value：该最优解的总运行成本

计算内容：

• 与 calobjvalue 类似，但针对单个粒子（最优解）

---

3. `initpop.m`

功能：初始化粒子群的位置。

输入：

• popsize：粒子数量（80）
• size：粒子维度（2）
• scope：时间范围（24小时）
• PMTmax、PMTmin：燃气轮机出力上下限

输出：

• pop：初始化后的粒子群位置（2×24×80）

初始化策略：

• 前半部分粒子：按燃气轮机出力均匀分布
• 后半部分粒子：按电网购电功率均匀分布
• 同时满足冷、电功率平衡约束

---

4. `main.m`

功能：主程序，执行粒子群优化算法，求解最优调度策略，并绘制结果图。

流程：

1. 初始化粒子群
2. 计算初始适应度
3. 设置 PSO 参数（惯性因子、加速因子、最大迭代次数）
4. 迭代优化：
   - 更新个体最优 pbest
   - 更新全局最优 gbest
   - 更新粒子速度和位置
   - 检查并修正约束条件（燃气轮机出力上下限、功率平衡）
5. 输出最优解及其运行成本
6. 绘制三张图：
   - 图1：夏季各设备电功率曲线（电负荷、燃机出力、电网交换功率、电空调功率）
   - 图2：夏季各设备冷功率曲线（冷负荷、烟气余热制冷、电空调制冷）
   - 图3：电负荷与冷负荷曲线对比

---

系统模型与约束说明

• 燃气轮机模型：
• 出力 PMT = VMT × 2.91
• 制冷量 QMTco = (19/12) × PMT
• 燃料成本、运行成本、排污成本合计为 3.85 + 0.038×2.91 + 0.2155×2.91 元/kWh

• 电空调模型：
• 制冷效率 COP = 2.3
• 电功率 Pairco = Qairco / 2.3

• 功率平衡约束：
• 电功率平衡：PGrid = PL + P_airco - PMT
• 冷功率平衡：Q_airco = QL - QMTco

• 优化目标：
• 最小化总运行成本 = 燃气轮机成本 + 电网购电成本

---