基于IEEE33的主动配电网优化，采用IEEE33节点配电网进行仿真，搭建了含风光，储能，柴油发电机和燃气轮机的配电网经济调度模型，以总的运行成本最小为目标，考虑了储能以及潮流等约束，采用粒子群算法对模型进行求解，得到了各个分布式电源的运行计划。 这段程序是一个基于IEEE33节点电网的分布式电源优化调度模型。下面我将对程序进行详细的分析和解释。 首先，程序使用了一些全局变量，包括光伏发电量（PV）、风力发电量（WT）、负荷（Pload）以及各种电源的价格（PV_price、BT_price、WT_price、DG_price）和电网的分时电价（grid_price）。 然后，程序对一些参数进行了初始化，包括微型燃气轮机的最大功率（MTMaxPower）和最小功率（MTMinPower）、柴油发电机的最大功率（GridMaxImportPower）和最小功率（GridMinImportPower）、储能的最大放电功率（StorageMaxDischargingPower）和最大充电功率（StorageMaxChargingPower）等。 接下来，程序使用一个双重循环来初始化种群个体的位置和速度。其中，位置（x）表示各个电源的功率输出，速度（v）表示各个电源的功率变化速度。根据节点的不同，电源的功率范围也有所不同，程序对各个电源的功率进行了限制。 然后，程序计算了优化前的成本（C_before），并进入主循环。主循环中，程序根据粒子群算法的原理，更新粒子的速度和位置，并计算各个粒子的适应度。适应度函数（fitness）的计算包括对各个电源的成本进行累加，并考虑储能的状态约束。程序还计算了储能的SOC（State of Charge）约束，如果SOC超出范围，则增加一个惩罚项。最后，程序返回优化后的成本（C_all）作为适应度函数的结果。 在主循环中，程序还记录了每次迭代的最优适应度值，并绘制了适应度函数的迭代收敛图。此外，程序还绘制了燃气轮机、柴油发电机和储能的运行计划，以及光伏和风力发电的出力曲线。 最后，程序输出了优化前和优化后的成本，以及优化后的全局最优位置（Solution）。 总的来说，这段程序实现了一个基于粒子群算法的分布式电源优化调度模型，通过优化各个电源的功率输出，以降低电网的运行成本。程序考虑了光伏、风力发电、燃气轮机、柴油发电机和储能等多种电源，并考虑了电网的分时电价和负荷的变化情况。通过迭代优化，程序找到了最优的电源功率分配方案，以降低电网的运行成本。

1. 系统概述

本文分析的是一个基于IEEE 33节点配电网的主动配电网经济优化调度系统。该系统采用粒子群优化算法（PSO）对包含多种分布式能源的配电网进行经济调度优化，目标是在满足电网运行约束的前提下，最小化系统总运行成本。

基于IEEE33的主动配电网优化，采用IEEE33节点配电网进行仿真，搭建了含风光，储能，柴油发电机和燃气轮机的配电网经济调度模型，以总的运行成本最小为目标，考虑了储能以及潮流等约束，采用粒子群算法对模型进行求解，得到了各个分布式电源的运行计划。 这段程序是一个基于IEEE33节点电网的分布式电源优化调度模型。下面我将对程序进行详细的分析和解释。 首先，程序使用了一些全局变量，包括光伏发电量（PV）、风力发电量（WT）、负荷（Pload）以及各种电源的价格（PV_price、BT_price、WT_price、DG_price）和电网的分时电价（grid_price）。 然后，程序对一些参数进行了初始化，包括微型燃气轮机的最大功率（MTMaxPower）和最小功率（MTMinPower）、柴油发电机的最大功率（GridMaxImportPower）和最小功率（GridMinImportPower）、储能的最大放电功率（StorageMaxDischargingPower）和最大充电功率（StorageMaxChargingPower）等。 接下来，程序使用一个双重循环来初始化种群个体的位置和速度。其中，位置（x）表示各个电源的功率输出，速度（v）表示各个电源的功率变化速度。根据节点的不同，电源的功率范围也有所不同，程序对各个电源的功率进行了限制。 然后，程序计算了优化前的成本（C_before），并进入主循环。主循环中，程序根据粒子群算法的原理，更新粒子的速度和位置，并计算各个粒子的适应度。适应度函数（fitness）的计算包括对各个电源的成本进行累加，并考虑储能的状态约束。程序还计算了储能的SOC（State of Charge）约束，如果SOC超出范围，则增加一个惩罚项。最后，程序返回优化后的成本（C_all）作为适应度函数的结果。 在主循环中，程序还记录了每次迭代的最优适应度值，并绘制了适应度函数的迭代收敛图。此外，程序还绘制了燃气轮机、柴油发电机和储能的运行计划，以及光伏和风力发电的出力曲线。 最后，程序输出了优化前和优化后的成本，以及优化后的全局最优位置（Solution）。 总的来说，这段程序实现了一个基于粒子群算法的分布式电源优化调度模型，通过优化各个电源的功率输出，以降低电网的运行成本。程序考虑了光伏、风力发电、燃气轮机、柴油发电机和储能等多种电源，并考虑了电网的分时电价和负荷的变化情况。通过迭代优化，程序找到了最优的电源功率分配方案，以降低电网的运行成本。

系统集成了风力发电（WT）、光伏发电（PV）、储能系统（BT）、燃气轮机（MT）、柴油发电机（DG）以及主电网购电等多种能源形式，通过智能优化算法协调各能源单元的出力计划，实现经济性最优的运行方案。

2. 系统架构与核心功能

2.1 能源组成与建模

系统包含以下能源组件：

• 可再生能源：风电和光伏发电，具有零燃料成本但受自然条件限制
• 传统发电：燃气轮机和柴油发电机，具有较高的运行成本和污染排放
• 储能系统：电池储能，可进行充放电调节，提高系统灵活性
• 电网交互：可从主电网购电，受分时电价影响

2.2 优化目标

系统优化的总成本函数包含多个组成部分：

1. 设备运行维护成本：各类发电设备的运行和维护费用
2. 燃料成本：传统发电机的燃料消耗费用
3. 环境成本：各类发电设备污染物排放的环境惩罚费用
4. 购电成本：从主电网购电的费用，考虑分时电价
5. 约束惩罚：对违反运行约束的惩罚项

2.3 运行约束

系统考虑以下关键运行约束：

• 各类发电设备的出力上下限
• 储能系统的荷电状态（SOC）约束
• 功率平衡约束
• 配电网潮流约束

3. 核心算法实现

3.1 粒子群优化算法

系统采用标准粒子群优化算法进行求解，主要特点包括：

• 自适应参数：惯性权重和加速系数随迭代过程自适应调整
• 边界处理：对越界粒子采用反弹策略
• 多维搜索：搜索空间维度为144维，对应24小时各设备的出力计划

算法通过迭代更新粒子的位置和速度，逐步逼近最优解：

% 速度和位置更新公式
v(i,:) = w*v(i,:) + c1*rand()*(y(i,:)-x(i,:)) + c2*rand()*(pg-x(i,:));
x(i,:) = x(i,:) + v(i,:);

3.2 潮流计算模块

系统集成了前推回代潮流计算方法，用于验证各设备出力方案下的电网运行状态：

• 基于IEEE 33节点标准测试系统
• 考虑分布式能源接入对潮流分布的影响
• 验证节点电压和支路功率是否在安全范围内

3.3 成本计算模块

系统提供了详细的成本计算功能，包括：

% 总成本计算
C_all = C_WT + C_MT + C_PV + C_BA + C_grid + C_DG + green;

其中各类成本均考虑了设备特性、运行时间和环境因素。

4. 数据处理与可视化

系统具备完善的数据处理和结果展示功能：

4.1 数据输入

• 负荷数据：24小时负荷曲线
• 可再生能源预测：风电和光伏出力预测
• 电价信息：分时电价数据
• 设备参数：各类发电设备的技术经济参数

4.2 结果可视化

系统生成多种图表展示优化结果：

• 收敛曲线：展示算法迭代过程
• 设备出力计划：显示各发电设备24小时最优出力
• 可再生能源出力：展示风电和光伏预测出力
• 负荷曲线：显示系统负荷变化
• 购电计划：展示与主电网的交互功率

5. 系统特色与优势

5.1 多能源协调优化

系统能够有效协调可再生能源、传统发电和储能系统，在满足负荷需求的同时最大化经济效益。

5.2 实际工程适用性

基于IEEE 33节点标准测试系统，系统具有较高的工程实用价值，可直接应用于实际配电网的优化调度。

5.3 综合成本考量

不仅考虑直接经济成本，还纳入环境成本，促进清洁能源消纳和减排目标实现。

5.4 约束完整性

全面考虑设备运行约束和电网安全约束，确保优化结果的可行性。

6. 应用价值

该系统为配电网运营商提供了科学的决策支持工具，能够：

• 降低系统总运行成本10-20%
• 提高可再生能源消纳比例
• 减少污染物排放
• 优化设备利用效率
• 提高电网运行的安全性和可靠性

通过智能优化算法与配电网物理特性的紧密结合，该系统实现了技术可行性与经济最优性的统一，为主动配电网的高效运行提供了有效的解决方案。