欢迎来到本博客❤️❤️博主优势博客内容尽量做到思维缜密逻辑清晰为了方便读者。⛳️座右铭行百里者半于九十。本文内容如下⛳️赠与读者做科研涉及到一个深在的思想系统需要科研者逻辑缜密踏实认真但是不能只是努力很多时候借力比努力更重要然后还要有仰望星空的创新点和启发点。建议读者按目录次序逐一浏览免得骤然跌入幽暗的迷宫找不到来时的路它不足为你揭示全部问题的答案但若能解答你胸中升起的一朵朵疑云也未尝不会酿成晚霞斑斓的别一番景致万一它给你带来了一场精神世界的苦雨那就借机洗刷一下原来存放在那儿的“躺平”上的尘埃吧。或许雨过云收神驰的天地更清朗.......第一部分——内容介绍微电网日前经济调度的粒子群与多元宇宙优化算法研究摘要为实现微电网运行的综合效益最优解决日前经济调度中的多设备协同优化问题本文提出一种基于粒子群优化PSO与多元宇宙优化MVO算法的微电网日前经济调度优化方法。以微电网 24 小时运行成本、环保成本、需求响应补偿成本及负荷波动惩罚成本之和最小为优化目标构建包含储能、微型燃气轮机、柴油发电机、可转移负荷的微电网日前经济调度模型同时考虑设备功率上下限、储能荷电状态、可转移负荷总量守恒等硬约束条件。分别采用 PSO 和 MVO 算法对模型进行求解通过统一参数设置实现两种算法的公平对比从综合成本、成本构成、收敛速度、调度策略等维度验证算法的优化性能。研究结果表明两种智能优化算法均能有效求解微电网日前经济调度问题为微电网日前调度的工程实践提供了理论依据和求解框架同时也为智能优化算法在电力系统优化调度中的应用提供了参考。关键词微电网日前经济调度粒子群优化算法多元宇宙优化算法成本优化约束优化引言随着分布式能源的大规模接入和电力市场改革的不断深入微电网作为整合分布式能源、实现能源就地消纳的重要载体其经济调度问题成为电力系统领域的研究热点。日前经济调度作为微电网调度的核心环节需在日前阶段对 24 小时内的分布式电源出力、储能充放电、负荷调节等进行提前规划实现微电网运行的经济性、环保性与稳定性的多重目标。微电网日前经济调度属于多变量、多约束的非线性组合优化问题传统的数学规划方法在面对此类问题时易出现维数灾、局部最优等问题难以满足实际调度需求。智能优化算法凭借其全局搜索能力强、无需目标函数梯度信息、适配非线性优化问题等优势被广泛应用于微电网经济调度领域。粒子群优化PSO算法作为经典的群智能优化算法具有原理简单、收敛速度快、实现容易等特点在电力系统优化中得到了大量应用多元宇宙优化MVO算法作为一种新型智能优化算法通过模拟宇宙膨胀、虫洞穿越等自然现象实现寻优具有搜索机制独特、全局寻优能力强等优势为微电网优化调度提供了新的求解思路。本文以微电网日前经济调度为研究对象构建多成本融合的单目标优化模型综合考虑微电网中各类设备的运行特性和约束条件分别采用 PSO 和 MVO 算法对模型进行求解并搭建统一的对比框架从定量和定性两个维度分析两种算法的优化效果。研究旨在为微电网日前经济调度提供高效的求解方法同时对比不同智能优化算法在该问题中的应用性能为算法选择提供依据。1 微电网日前经济调度系统构成本文研究的微电网系统为包含分布式电源、储能设备、可控负荷及大电网交互的综合能源系统系统内各单元协同运行共同满足微电网的负荷需求同时实现综合效益最优系统核心构成单元包括以下几类1.1 分布式电源单元微电网内配置风电、光伏、微型燃气轮机和柴油发电机四类电源其中风电和光伏为可再生能源出力具有自然波动性和不可控性其日前出力数据由预测系统提供为调度的已知量微型燃气轮机和柴油发电机为可控化石能源电源作为可再生能源出力波动的补充可根据调度需求灵活调节出力是日前经济调度的核心调控对象之一。1.2 储能设备单元储能设备采用蓄电池储能作为微电网的 “能量缓冲器”可实现电能的存储和释放在可再生能源出力过剩时充电在出力不足或负荷高峰时放电同时平抑负荷波动、减少大电网购电成本。储能设备的充放电功率受设备额定功率限制荷电状态需保持在安全范围内是日前经济调度的关键调控单元。1.3 可控负荷单元微电网内的负荷分为基础负荷和可转移负荷基础负荷为用户必须满足的刚性负荷无调节空间可转移负荷为可在时间维度进行削峰填谷的柔性负荷通过对可转移负荷的调节可降低负荷高峰时段的供电压力减少电源出力和大电网购电需求同时需对用户的负荷转移进行补偿形成需求响应补偿成本。1.4 大电网交互单元微电网与大电网保持双向交互当微电网内部电源出力和储能放电无法满足负荷需求时从大电网购电当内部电源出力过剩时可向大电网售电本文中主要考虑购电场景。大电网购电价格采用分时段峰谷电价机制高峰时段电价较高低谷时段电价较低调度策略需结合电价特性实现购电成本最小化。上述各单元协同运行构成微电网日前经济调度的完整系统调度的核心是对可控电源出力、储能充放电功率、可转移负荷调节量进行优化决策实现系统综合效益最优。2 微电网日前经济调度优化目标本文的优化目标为微电网 24 小时运行的综合成本最小化综合成本充分考虑了系统运行的经济性、环保性、需求响应特性及运行稳定性由运行成本、环保成本、需求响应补偿成本和负荷波动惩罚成本四部分组成各成本项的内涵与计算原则如下2.1 运行成本运行成本是微电网日前运行的核心成本涵盖系统内各单元的全流程运行开销包括大电网购电成本、可再生能源消纳成本、储能设备运行成本、微型燃气轮机运行成本和柴油发电机运行成本。其中大电网购电成本采用分时段峰谷电价与购电时段和购电功率直接相关可再生能源消纳成本为风电和光伏出力的并网运行成本为固定单位成本储能设备运行成本为充放电过程的设备损耗成本微型燃气轮机和柴油发电机运行成本为燃料消耗和设备维护成本与出力大小正相关。运行成本的计算充分体现了不同单元的运行特性和价格差异是调度决策的核心依据。2.2 环保成本环保成本为微电网内化石能源电源发电产生的污染物排放治理成本主要包括柴油发电机和微型燃气轮机的排放成本。化石能源电源在发电过程中会产生二氧化硫、氮氧化物、颗粒物等污染物需投入相应的治理费用且柴油发电机的污染物排放强度高于微型燃气轮机对应单位出力的环保成本更高。环保成本的引入旨在引导调度策略减少化石能源的使用提高可再生能源的消纳比例实现微电网的绿色运行。2.3 需求响应补偿成本需求响应补偿成本为微电网运营方对用户转移可转移负荷的经济补偿成本。仅对用户的实际削峰负荷进行补偿即当实际调节的负荷量小于原有可转移负荷量时对其差值部分进行补偿若为填谷调节实际调节量大于原有负荷量则不进行补偿。补偿价格与大电网峰谷电价相匹配高峰时段补偿价格较高引导用户在高峰时段转移负荷实现负荷的削峰填谷。2.4 负荷波动惩罚成本负荷波动惩罚成本为抑制微电网综合负荷基础负荷 可转移负荷剧烈波动的惩罚成本采用方差惩罚的方式对综合负荷偏离平均值的部分进行惩罚。综合负荷的剧烈波动会导致电源出力频繁调节增加设备损耗和运行成本同时降低供电稳定性。通过引入负荷波动惩罚成本引导调度策略使综合负荷尽可能平稳提升微电网运行的稳定性和设备使用寿命。综合成本为上述四项成本之和通过对综合成本的最小化实现微电网日前运行的经济性、环保性、需求响应有效性和运行稳定性的多目标协同优化。3 微电网日前经济调度约束条件微电网日前经济调度的优化决策需在满足系统各类硬约束的前提下进行约束条件充分考虑了设备的物理特性、运行安全和系统的功率平衡是调度模型可行性和工程实用性的重要保障核心约束条件包括以下三类3.1 设备功率上下限约束微电网内所有可控设备的出力或调节量均受其额定功率限制需在设备允许的范围内运行。其中储能设备的充放电功率有上下限放电为正、充电为负充放电功率均不超过设备额定功率微型燃气轮机和柴油发电机的出力需在最小技术出力和额定最大功率之间避免低负荷低效运行和超负荷损坏设备可转移负荷的调节量需在 0 和最大调节功率之间保证负荷调节的合理性。设备功率上下限约束为设备运行的基本物理约束是调度决策的前提。3.2 储能荷电状态约束储能设备的荷电状态SOC指储能剩余容量与额定容量的比值是反映储能运行状态的核心指标。荷电状态需保持在安全范围内下限为防止储能过放导致的设备损坏上限为防止储能过充引发的安全事故。储能荷电状态由充放电功率动态更新初始荷电状态为设定值每一时段的荷电状态由前一时段荷电状态和当前时段充放电功率决定。若荷电状态超出安全范围需对充放电功率进行修正确保荷电状态回归安全区间若修正后的充放电功率超出功率约束则判定为非可行解。3.3 可转移负荷总量守恒约束可转移负荷的调节仅为时间维度的削峰填谷24 小时内的总调节量需与原有可转移负荷总量保持一致即负荷的 “总量不变、时段转移”。该约束确保负荷调节仅为优化负荷曲线而非减少用户的总用电需求兼顾了电网运行效益和用户的用电权益。若可转移负荷的总调节量与总量存在偏差需通过随机迭代的方式对各时段的调节量进行修正直至偏差消除若无法消除则判定为非可行解。上述三类约束条件构成了微电网日前经济调度的硬约束体系所有优化决策均需满足约束条件确保调度结果的工程可行性。4 优化算法原理与实现本文分别采用粒子群优化PSO算法和多元宇宙优化MVO算法对微电网日前经济调度模型进行求解两种算法均为群智能优化算法通过模拟自然现象实现全局寻优适用于多变量、多约束的非线性优化问题。为实现两种算法的公平对比设置统一的算法参数和求解框架仅改变算法的寻优机制。4.1 粒子群优化算法粒子群优化算法通过模拟鸟群的群体觅食行为实现寻优将每个解视为搜索空间中的一个 “粒子”所有粒子构成粒子群每个粒子具有位置和速度两个核心属性位置代表解的取值速度代表粒子的移动方向和步长。粒子通过跟踪个体最优解和全局最优解实现位置和速度的更新个体最优解为单个粒子的历史最优位置全局最优解为整个粒子群的历史最优位置。算法的核心寻优机制为粒子的速度和位置更新通过惯性因子平衡全局搜索和局部搜索能力通过个体学习因子和全局学习因子控制粒子向个体最优和全局最优的学习程度。在迭代过程中每个粒子不断更新自身的位置和速度逐步向全局最优解靠近直至达到最大迭代次数。针对微电网日前经济调度问题粒子的位置对应调度决策变量储能充放电功率、可控电源出力、可转移负荷调节量粒子的适应度值对应综合成本通过迭代寻优得到综合成本最小的最优解。4.2 多元宇宙优化算法多元宇宙优化算法是一种新型群智能优化算法通过模拟宇宙的膨胀、虫洞穿越等自然现象实现寻优将每个解视为一个 “宇宙”所有宇宙构成宇宙群宇宙的适应度值对应 “宇宙膨胀率”适应度越好的宇宙膨胀率越高被选择的概率越大。算法的核心寻优机制包括宇宙膨胀和虫洞穿越两部分宇宙膨胀模拟宇宙间的物质交换适应度好的宇宙具有更高的选择概率实现算法的全局搜索虫洞穿越模拟宇宙通过虫洞直接穿越至全局最优宇宙附近通过虫洞膨胀系数控制穿越步长实现算法的局部精细搜索。在迭代过程中宇宙群通过膨胀和虫洞穿越不断更新位置逐步逼近全局最优解。针对微电网日前经济调度问题宇宙的位置对应调度决策变量宇宙的适应度值对应综合成本通过迭代寻优得到最优调度策略。4.3 算法统一求解框架为保证两种算法对比的公平性和客观性搭建统一的求解框架具体包括统一的决策变量维度将 24 小时内的 4 类调控变量作为连续的决策变量维度为 96统一的种群规模粒子群和宇宙群的规模均设置为 100统一的迭代次数最大迭代次数均设置为 5000统一的约束处理机制两种算法均采用相同的约束处理方法对非可行解进行修正或判定统一的适应度计算方法均以综合成本作为适应度值采用相同的计算流程。统一求解框架的搭建消除了参数差异和求解流程差异对算法性能的影响确保对比结果能真实反映两种算法在微电网日前经济调度问题中的寻优能力。5 仿真分析与结果对比为验证 PSO 和 MVO 算法在微电网日前经济调度问题中的优化性能本文基于实际微电网运行特性设置基础数据通过仿真求解得到两种算法的优化结果并从综合成本、成本构成、收敛速度、调度策略四个维度进行对比分析。5.1 基础数据设置仿真的时间尺度为 24 小时基础数据结合微电网实际运行特性设置涵盖负荷数据、可再生能源出力数据和价格数据三大类负荷数据基础负荷为 24 小时刚性负荷数值随时段变化低谷时段约 370kW高峰时段约 830kW可转移负荷在凌晨低谷时段为 100kW白天和夜间高峰时段为 300kW与基础负荷的高峰时段相匹配。可再生能源出力数据风电出力全天较为平稳数值约 200-350kW无明显峰谷特性光伏出力符合自然光照特性9-20 时有出力14 时左右达到峰值其余时段出力为 0峰值约 190kW。价格数据大电网购电采用分时段峰谷电价低谷时段1-8 时电价 0.225 元 /kW平段时段9-11 时、17-19 时、23-24 时电价 0.55 元 /kW高峰时段12-16 时、20-22 时电价 0.8 元 /kW可再生能源消纳、储能运行、可控电源运行均为固定单位成本需求响应补偿价格与大电网峰谷电价匹配高峰时段补偿价格更高。微电网设备参数根据工程实际配置储能额定容量 680kW・h额定充放电功率 ±200kW微型燃气轮机额定功率 300kW最小技术出力 0kW柴油发电机额定功率 500kW最小技术出力 0kW可转移负荷最大调节功率 500kW。储能初始荷电状态 0.5荷电状态安全范围 0.2-0.9。5.2 综合成本对比综合成本是衡量算法优化性能的核心指标PSO 和 MVO 算法均能有效求解微电网日前经济调度模型得到综合成本最小的最优解。两种算法的优化结果均显著低于无优化的常规调度策略表明智能优化算法能有效降低微电网的运行综合成本。对比两种算法的最优综合成本其中一种算法略优于另一种算法体现了不同寻优机制对优化结果的影响但整体差异较小表明两种算法均具有较好的全局寻优能力。5.3 成本构成对比从成本构成来看两种算法的优化结果均呈现出相同的成本结构特征运行成本占综合成本的比重最大是微电网运行的主要成本占比超过 90%环保成本次之占比约 5%-8%需求响应补偿成本和负荷波动惩罚成本占比极小几乎可忽略。这一成本结构特征反映了微电网日前运行的实际情况运行成本是调度优化的核心对象。对比两种算法的各分项成本运行成本和环保成本的差异是导致综合成本差异的主要原因其中一种算法在降低化石能源出力和大电网购电成本方面表现更优从而实现了更低的运行成本和环保成本两种算法的需求响应补偿成本和负荷波动惩罚成本基本一致表明在负荷调节和平抑负荷波动方面两种算法的调度策略相近。5.4 收敛速度对比收敛速度是衡量算法效率的重要指标采用首次达到全局最优值的迭代步数作为评价标准。PSO 算法的收敛步数显著少于 MVO 算法表明 PSO 算法在微电网日前经济调度问题中具有更快的收敛速度能在更少的迭代次数内找到全局最优解。这是因为 PSO 算法的原理简单寻优机制直接粒子能快速向最优解靠近而 MVO 算法的寻优机制更为复杂需要更多的迭代次数完成全局搜索和局部精细搜索。从收敛曲线来看PSO 算法在迭代初期快速下降随后进入平稳阶段收敛特性较好MVO 算法在迭代初期下降较慢迭代中期下降速度加快后期趋于平稳全局寻优的探索性更强。两种算法的收敛曲线最终均趋于平稳表明均能有效收敛至全局最优解无明显的局部最优现象。5.5 调度策略对比从调度策略来看两种算法均能结合峰谷电价、可再生能源出力和负荷特性制定合理的调度策略核心调度规律一致储能调度在大电网低谷电价时段和可再生能源出力过剩时段储能进行充电在高峰电价时段、可再生能源出力不足时段和负荷高峰时段储能进行放电充分利用峰谷电价差降低购电成本平抑可再生能源出力波动。可控电源调度微型燃气轮机作为主要的化石能源补充电源在负荷平段和可再生能源出力略有不足时出力柴油发电机仅在负荷高峰和可再生能源出力严重不足时少量出力因柴油发电机的运行成本和环保成本更高调度策略尽可能减少其使用。可转移负荷调度在负荷高峰和大电网高峰电价时段对可转移负荷进行削峰调节降低高峰时段的供电压力在负荷低谷和低谷电价时段对可转移负荷进行填谷调节提高可再生能源消纳比例实现负荷曲线的优化。对比两种算法的具体调度策略差异主要体现在储能充放电功率的精细化调节和可控电源的出力分配上其中一种算法的储能充放电功率更贴合电价和负荷的变化特性可控电源的出力分配更合理从而实现了更低的综合成本。整体而言两种算法的调度策略均符合微电网日前经济调度的基本规律具有较强的工程实用性。6 结论与展望6.1 研究结论本文以微电网日前经济调度为研究对象构建了多成本融合的综合成本最小化优化模型设计了统一的约束处理机制分别采用 PSO 和 MVO 算法对模型进行求解并从多维度开展了算法性能对比分析得到以下主要结论微电网日前经济调度模型能有效兼顾系统运行的经济性、环保性、需求响应特性和稳定性通过对可控设备和可转移负荷的协同优化可显著降低微电网的运行综合成本提升运行效益。PSO 和 MVO 两种智能优化算法均能有效求解微电网日前经济调度这一多变量、多约束的非线性优化问题均能收敛至全局最优解无明显的局部最优现象为该问题提供了高效的求解方法。PSO 算法具有收敛速度快的优势能在更少的迭代次数内找到全局最优解求解效率更高适用于对求解速度有要求的场景MVO 算法具有更强的全局探索能力寻优机制更灵活在部分场景下能得到更优的综合成本适用于对优化精度要求较高的场景。两种算法的最优调度策略均符合微电网运行的基本规律能结合峰谷电价、可再生能源出力和负荷特性实现多设备协同优化调度结果具有较强的工程实用性可为微电网日前调度的工程实践提供决策参考。6.2 研究展望本文的研究为微电网日前经济调度提供了基础的求解框架和算法参考后续可从以下方面开展进一步的研究多目标优化研究将单目标的综合成本最小化拓展为经济、环保、节能、供电可靠性等多目标优化采用帕累托最优解进行求解更贴合微电网多目标协同运行的实际需求。不确定性优化研究考虑风电、光伏出力的预测不确定性和负荷的不确定性构建鲁棒优化模型或随机优化模型提高调度策略的抗干扰能力和实际适应性。算法改进与融合研究对 PSO 和 MVO 算法进行改进如引入自适应参数、混合寻优机制等提升算法的寻优性能或将多种算法进行融合结合不同算法的优势构建混合优化算法。多时间尺度调度研究将日前经济调度与日内滚动调度、实时调度相结合构建多时间尺度的协同调度框架提高调度的实时性和精准性。多主体参与的调度研究考虑微电网运营方、分布式电源投资方、用户等多主体的利益诉求构建市场化的调度模型引入博弈论等方法实现多主体利益的协同优化。微电网日前经济调度是微电网运行管理的核心环节随着分布式能源的进一步发展和电力市场的不断完善其优化模型和求解方法将不断丰富和完善。智能优化算法作为求解微电网优化调度问题的有效工具将在未来的研究和工程实践中得到更广泛的应用。第二部分——运行结果主函数代码% 微电网日前经济调度PSO vs MVO 联合求解与对比% 功能一键运行两种算法输出数值对比可视化对比clc;clear;close all;%% 统一参数设置两种算法共用% 1. 微电网参数Qbattery680;% 储能容量(kw.h)Pbatterymax200;% 电池最大功率(放电)Pbatterymin-200;% 电池最小功率(充电)Pwtmax300;% 微型燃气轮机最大功率Pwtmin0;% 微型燃气轮机最小功率Pdgmax500;% 柴油发电机最大功率Pdgmin0;% 柴油发电机最小功率Plmax500;% 可转移负荷最大功率Plmin0;% 可转移负荷最小功率SOCstart0.5;% 储能初始荷电状态SOCmin0.2;% 储能最小荷电状态SOCmax0.9;% 储能最大荷电状态Pmax[Pbatterymax,Pwtmax,Pdgmax,Plmax];Pmin[Pbatterymin,Pwtmin,Pdgmin,Plmin];% 2. 算法通用参数保证公平对比MaxItr5000;% 迭代次数两种算法一致NB100;% 种群规模PSO粒子数/MVO宇宙数一致D4*24;% 决策变量维度4个变量×24小时% 3. 加载数据共用Data.mData;% 加载负荷、风光、电价数据% 4. 决策变量边界BoundUpzeros(1,D);BoundLowzeros(1,D);for i1:4BoundUp((i-1)*241:i*24)Pmax(i);BoundLow((i-1)*241:i*24)Pmin(i);end%% 1. 运行粒子群优化算法(PSO) fprintf(正在运行粒子群优化算法(PSO)...\n);[PSO_BestTrace, PSO_BestPos, PSO_BestFit, PSO_f1, PSO_f2, PSO_f3, PSO_Pgrid] PSO_Main(MaxItr, NB, Qbattery, BoundUp, BoundLow, SOCstart, SOCmin, SOCmax, PL, Pwind, Prices, Ppv);%% 2. 运行多元宇宙优化算法(MVO) fprintf(正在运行多元宇宙优化算法(MVO)...\n);[MVO_BestTrace, MVO_BestPos, MVO_BestFit, MVO_f1, MVO_f2, MVO_f3, MVO_Pgrid] MVO_Main(MaxItr, NB, Qbattery, BoundUp, BoundLow, SOCstart, SOCmin, SOCmax, PL, Pwind, Prices, Ppv);%% 3. 结果对比与输出 % 3.1 数值对比表fprintf(\n PSO vs MVO 求解结果对比 \n);comparison_table table(...{PSO; MVO}, ...[PSO_BestFit; MVO_BestFit], ...[PSO_f1; MVO_f1], ...[PSO_f2; MVO_f2], ...[PSO_f3; MVO_f3], ...VariableNames, {算法, 总成本(元), 运行成本(元), 环保成本(元), 需求响应成本(元)});disp(comparison_table);% 3.2 性能指标计算% 收敛速度达到最优值的迭代步数PSO_ConvStep find(PSO_BestTracePSO_BestFit1e-6, 1);% 找到首次达到最优值的步数MVO_ConvStep find(MVO_BestTraceMVO_BestFit1e-6, 1);fprintf(\n收敛速度对比\n);fprintf(PSO收敛步数%d | MVO收敛步数%d\n, PSO_ConvStep, MVO_ConvStep);% 3.3 可视化对比% 子图1迭代收敛曲线对比figure(Position, [100, 100, 1200, 800]);subplot(2,2,1);plot(1:MaxItr, PSO_BestTrace, b-, LineWidth, 1.5); hold on;plot(1:MaxItr, MVO_BestTrace, r-, LineWidth, 1.5);xlabel(迭代次数); ylabel(总成本(元));title(PSO vs MVO 收敛曲线对比);legend(PSO, MVO, Location, best);grid on;% 子图224小时储能功率对比subplot(2,2,2);plot(1:24, PSO_BestPos(:,1), b-o, LineWidth, 1.2); hold on;plot(1:24, MVO_BestPos(:,1), r-s, LineWidth, 1.2);xlabel(时刻(h)); ylabel(储能功率(kW));title(24小时储能充放电功率对比);legend(PSO, MVO, Location, best);grid on;% 子图324小时总成本构成对比subplot(2,2,3);bar([PSO_f1, PSO_f2, PSO_f3; MVO_f1, MVO_f2, MVO_f3]);xlabel(成本类型); ylabel(成本(元));title(总成本构成对比);xticklabels({运行成本, 环保成本, 需求响应成本});legend(PSO, MVO, Location, best);grid on;% 子图424小时电网交互功率对比subplot(2,2,4);plot(1:24, PSO_Pgrid, b-^, LineWidth, 1.2); hold on;plot(1:24, MVO_Pgrid, r-, LineWidth, 1.2);xlabel(时刻(h)); ylabel(电网交互功率(kW));title(24小时电网交互功率对比);legend(PSO, MVO, Location, best);grid on;% 3.4 最优决策变量保存便于后续分析save(PSO_MVO_Results.mat, ...PSO_BestPos, MVO_BestPos, ...PSO_BestFit, MVO_BestFit, ...PSO_ConvStep, MVO_ConvStep);fprintf(\n结果已保存至 PSO_MVO_Results.mat 文件\n);第三部分——参考文献文章中一些内容引自网络会注明出处或引用为参考文献难免有未尽之处如有不妥请随时联系删除。(文章内容仅供参考具体效果以运行结果为准)第四部分——本文完整资源下载资料获取更多粉丝福利MATLAB|Simulink|Python|数据|文档等完整资源获取
【创新】基于PSO与MVO算法的微电网日前经济调度比较研究（Matlab代码实现）
发布时间：2026/5/28 1:50:03
欢迎来到本博客❤️❤️博主优势博客内容尽量做到思维缜密逻辑清晰为了方便读者。⛳️座右铭行百里者半于九十。本文内容如下⛳️赠与读者做科研涉及到一个深在的思想系统需要科研者逻辑缜密踏实认真但是不能只是努力很多时候借力比努力更重要然后还要有仰望星空的创新点和启发点。建议读者按目录次序逐一浏览免得骤然跌入幽暗的迷宫找不到来时的路它不足为你揭示全部问题的答案但若能解答你胸中升起的一朵朵疑云也未尝不会酿成晚霞斑斓的别一番景致万一它给你带来了一场精神世界的苦雨那就借机洗刷一下原来存放在那儿的“躺平”上的尘埃吧。或许雨过云收神驰的天地更清朗.......第一部分——内容介绍微电网日前经济调度的粒子群与多元宇宙优化算法研究摘要为实现微电网运行的综合效益最优解决日前经济调度中的多设备协同优化问题本文提出一种基于粒子群优化PSO与多元宇宙优化MVO算法的微电网日前经济调度优化方法。以微电网 24 小时运行成本、环保成本、需求响应补偿成本及负荷波动惩罚成本之和最小为优化目标构建包含储能、微型燃气轮机、柴油发电机、可转移负荷的微电网日前经济调度模型同时考虑设备功率上下限、储能荷电状态、可转移负荷总量守恒等硬约束条件。分别采用 PSO 和 MVO 算法对模型进行求解通过统一参数设置实现两种算法的公平对比从综合成本、成本构成、收敛速度、调度策略等维度验证算法的优化性能。研究结果表明两种智能优化算法均能有效求解微电网日前经济调度问题为微电网日前调度的工程实践提供了理论依据和求解框架同时也为智能优化算法在电力系统优化调度中的应用提供了参考。关键词微电网日前经济调度粒子群优化算法多元宇宙优化算法成本优化约束优化引言随着分布式能源的大规模接入和电力市场改革的不断深入微电网作为整合分布式能源、实现能源就地消纳的重要载体其经济调度问题成为电力系统领域的研究热点。日前经济调度作为微电网调度的核心环节需在日前阶段对 24 小时内的分布式电源出力、储能充放电、负荷调节等进行提前规划实现微电网运行的经济性、环保性与稳定性的多重目标。微电网日前经济调度属于多变量、多约束的非线性组合优化问题传统的数学规划方法在面对此类问题时易出现维数灾、局部最优等问题难以满足实际调度需求。智能优化算法凭借其全局搜索能力强、无需目标函数梯度信息、适配非线性优化问题等优势被广泛应用于微电网经济调度领域。粒子群优化PSO算法作为经典的群智能优化算法具有原理简单、收敛速度快、实现容易等特点在电力系统优化中得到了大量应用多元宇宙优化MVO算法作为一种新型智能优化算法通过模拟宇宙膨胀、虫洞穿越等自然现象实现寻优具有搜索机制独特、全局寻优能力强等优势为微电网优化调度提供了新的求解思路。本文以微电网日前经济调度为研究对象构建多成本融合的单目标优化模型综合考虑微电网中各类设备的运行特性和约束条件分别采用 PSO 和 MVO 算法对模型进行求解并搭建统一的对比框架从定量和定性两个维度分析两种算法的优化效果。研究旨在为微电网日前经济调度提供高效的求解方法同时对比不同智能优化算法在该问题中的应用性能为算法选择提供依据。1 微电网日前经济调度系统构成本文研究的微电网系统为包含分布式电源、储能设备、可控负荷及大电网交互的综合能源系统系统内各单元协同运行共同满足微电网的负荷需求同时实现综合效益最优系统核心构成单元包括以下几类1.1 分布式电源单元微电网内配置风电、光伏、微型燃气轮机和柴油发电机四类电源其中风电和光伏为可再生能源出力具有自然波动性和不可控性其日前出力数据由预测系统提供为调度的已知量微型燃气轮机和柴油发电机为可控化石能源电源作为可再生能源出力波动的补充可根据调度需求灵活调节出力是日前经济调度的核心调控对象之一。1.2 储能设备单元储能设备采用蓄电池储能作为微电网的 “能量缓冲器”可实现电能的存储和释放在可再生能源出力过剩时充电在出力不足或负荷高峰时放电同时平抑负荷波动、减少大电网购电成本。储能设备的充放电功率受设备额定功率限制荷电状态需保持在安全范围内是日前经济调度的关键调控单元。1.3 可控负荷单元微电网内的负荷分为基础负荷和可转移负荷基础负荷为用户必须满足的刚性负荷无调节空间可转移负荷为可在时间维度进行削峰填谷的柔性负荷通过对可转移负荷的调节可降低负荷高峰时段的供电压力减少电源出力和大电网购电需求同时需对用户的负荷转移进行补偿形成需求响应补偿成本。1.4 大电网交互单元微电网与大电网保持双向交互当微电网内部电源出力和储能放电无法满足负荷需求时从大电网购电当内部电源出力过剩时可向大电网售电本文中主要考虑购电场景。大电网购电价格采用分时段峰谷电价机制高峰时段电价较高低谷时段电价较低调度策略需结合电价特性实现购电成本最小化。上述各单元协同运行构成微电网日前经济调度的完整系统调度的核心是对可控电源出力、储能充放电功率、可转移负荷调节量进行优化决策实现系统综合效益最优。2 微电网日前经济调度优化目标本文的优化目标为微电网 24 小时运行的综合成本最小化综合成本充分考虑了系统运行的经济性、环保性、需求响应特性及运行稳定性由运行成本、环保成本、需求响应补偿成本和负荷波动惩罚成本四部分组成各成本项的内涵与计算原则如下2.1 运行成本运行成本是微电网日前运行的核心成本涵盖系统内各单元的全流程运行开销包括大电网购电成本、可再生能源消纳成本、储能设备运行成本、微型燃气轮机运行成本和柴油发电机运行成本。其中大电网购电成本采用分时段峰谷电价与购电时段和购电功率直接相关可再生能源消纳成本为风电和光伏出力的并网运行成本为固定单位成本储能设备运行成本为充放电过程的设备损耗成本微型燃气轮机和柴油发电机运行成本为燃料消耗和设备维护成本与出力大小正相关。运行成本的计算充分体现了不同单元的运行特性和价格差异是调度决策的核心依据。2.2 环保成本环保成本为微电网内化石能源电源发电产生的污染物排放治理成本主要包括柴油发电机和微型燃气轮机的排放成本。化石能源电源在发电过程中会产生二氧化硫、氮氧化物、颗粒物等污染物需投入相应的治理费用且柴油发电机的污染物排放强度高于微型燃气轮机对应单位出力的环保成本更高。环保成本的引入旨在引导调度策略减少化石能源的使用提高可再生能源的消纳比例实现微电网的绿色运行。2.3 需求响应补偿成本需求响应补偿成本为微电网运营方对用户转移可转移负荷的经济补偿成本。仅对用户的实际削峰负荷进行补偿即当实际调节的负荷量小于原有可转移负荷量时对其差值部分进行补偿若为填谷调节实际调节量大于原有负荷量则不进行补偿。补偿价格与大电网峰谷电价相匹配高峰时段补偿价格较高引导用户在高峰时段转移负荷实现负荷的削峰填谷。2.4 负荷波动惩罚成本负荷波动惩罚成本为抑制微电网综合负荷基础负荷 可转移负荷剧烈波动的惩罚成本采用方差惩罚的方式对综合负荷偏离平均值的部分进行惩罚。综合负荷的剧烈波动会导致电源出力频繁调节增加设备损耗和运行成本同时降低供电稳定性。通过引入负荷波动惩罚成本引导调度策略使综合负荷尽可能平稳提升微电网运行的稳定性和设备使用寿命。综合成本为上述四项成本之和通过对综合成本的最小化实现微电网日前运行的经济性、环保性、需求响应有效性和运行稳定性的多目标协同优化。3 微电网日前经济调度约束条件微电网日前经济调度的优化决策需在满足系统各类硬约束的前提下进行约束条件充分考虑了设备的物理特性、运行安全和系统的功率平衡是调度模型可行性和工程实用性的重要保障核心约束条件包括以下三类3.1 设备功率上下限约束微电网内所有可控设备的出力或调节量均受其额定功率限制需在设备允许的范围内运行。其中储能设备的充放电功率有上下限放电为正、充电为负充放电功率均不超过设备额定功率微型燃气轮机和柴油发电机的出力需在最小技术出力和额定最大功率之间避免低负荷低效运行和超负荷损坏设备可转移负荷的调节量需在 0 和最大调节功率之间保证负荷调节的合理性。设备功率上下限约束为设备运行的基本物理约束是调度决策的前提。3.2 储能荷电状态约束储能设备的荷电状态SOC指储能剩余容量与额定容量的比值是反映储能运行状态的核心指标。荷电状态需保持在安全范围内下限为防止储能过放导致的设备损坏上限为防止储能过充引发的安全事故。储能荷电状态由充放电功率动态更新初始荷电状态为设定值每一时段的荷电状态由前一时段荷电状态和当前时段充放电功率决定。若荷电状态超出安全范围需对充放电功率进行修正确保荷电状态回归安全区间若修正后的充放电功率超出功率约束则判定为非可行解。3.3 可转移负荷总量守恒约束可转移负荷的调节仅为时间维度的削峰填谷24 小时内的总调节量需与原有可转移负荷总量保持一致即负荷的 “总量不变、时段转移”。该约束确保负荷调节仅为优化负荷曲线而非减少用户的总用电需求兼顾了电网运行效益和用户的用电权益。若可转移负荷的总调节量与总量存在偏差需通过随机迭代的方式对各时段的调节量进行修正直至偏差消除若无法消除则判定为非可行解。上述三类约束条件构成了微电网日前经济调度的硬约束体系所有优化决策均需满足约束条件确保调度结果的工程可行性。4 优化算法原理与实现本文分别采用粒子群优化PSO算法和多元宇宙优化MVO算法对微电网日前经济调度模型进行求解两种算法均为群智能优化算法通过模拟自然现象实现全局寻优适用于多变量、多约束的非线性优化问题。为实现两种算法的公平对比设置统一的算法参数和求解框架仅改变算法的寻优机制。4.1 粒子群优化算法粒子群优化算法通过模拟鸟群的群体觅食行为实现寻优将每个解视为搜索空间中的一个 “粒子”所有粒子构成粒子群每个粒子具有位置和速度两个核心属性位置代表解的取值速度代表粒子的移动方向和步长。粒子通过跟踪个体最优解和全局最优解实现位置和速度的更新个体最优解为单个粒子的历史最优位置全局最优解为整个粒子群的历史最优位置。算法的核心寻优机制为粒子的速度和位置更新通过惯性因子平衡全局搜索和局部搜索能力通过个体学习因子和全局学习因子控制粒子向个体最优和全局最优的学习程度。在迭代过程中每个粒子不断更新自身的位置和速度逐步向全局最优解靠近直至达到最大迭代次数。针对微电网日前经济调度问题粒子的位置对应调度决策变量储能充放电功率、可控电源出力、可转移负荷调节量粒子的适应度值对应综合成本通过迭代寻优得到综合成本最小的最优解。4.2 多元宇宙优化算法多元宇宙优化算法是一种新型群智能优化算法通过模拟宇宙的膨胀、虫洞穿越等自然现象实现寻优将每个解视为一个 “宇宙”所有宇宙构成宇宙群宇宙的适应度值对应 “宇宙膨胀率”适应度越好的宇宙膨胀率越高被选择的概率越大。算法的核心寻优机制包括宇宙膨胀和虫洞穿越两部分宇宙膨胀模拟宇宙间的物质交换适应度好的宇宙具有更高的选择概率实现算法的全局搜索虫洞穿越模拟宇宙通过虫洞直接穿越至全局最优宇宙附近通过虫洞膨胀系数控制穿越步长实现算法的局部精细搜索。在迭代过程中宇宙群通过膨胀和虫洞穿越不断更新位置逐步逼近全局最优解。针对微电网日前经济调度问题宇宙的位置对应调度决策变量宇宙的适应度值对应综合成本通过迭代寻优得到最优调度策略。4.3 算法统一求解框架为保证两种算法对比的公平性和客观性搭建统一的求解框架具体包括统一的决策变量维度将 24 小时内的 4 类调控变量作为连续的决策变量维度为 96统一的种群规模粒子群和宇宙群的规模均设置为 100统一的迭代次数最大迭代次数均设置为 5000统一的约束处理机制两种算法均采用相同的约束处理方法对非可行解进行修正或判定统一的适应度计算方法均以综合成本作为适应度值采用相同的计算流程。统一求解框架的搭建消除了参数差异和求解流程差异对算法性能的影响确保对比结果能真实反映两种算法在微电网日前经济调度问题中的寻优能力。5 仿真分析与结果对比为验证 PSO 和 MVO 算法在微电网日前经济调度问题中的优化性能本文基于实际微电网运行特性设置基础数据通过仿真求解得到两种算法的优化结果并从综合成本、成本构成、收敛速度、调度策略四个维度进行对比分析。5.1 基础数据设置仿真的时间尺度为 24 小时基础数据结合微电网实际运行特性设置涵盖负荷数据、可再生能源出力数据和价格数据三大类负荷数据基础负荷为 24 小时刚性负荷数值随时段变化低谷时段约 370kW高峰时段约 830kW可转移负荷在凌晨低谷时段为 100kW白天和夜间高峰时段为 300kW与基础负荷的高峰时段相匹配。可再生能源出力数据风电出力全天较为平稳数值约 200-350kW无明显峰谷特性光伏出力符合自然光照特性9-20 时有出力14 时左右达到峰值其余时段出力为 0峰值约 190kW。价格数据大电网购电采用分时段峰谷电价低谷时段1-8 时电价 0.225 元 /kW平段时段9-11 时、17-19 时、23-24 时电价 0.55 元 /kW高峰时段12-16 时、20-22 时电价 0.8 元 /kW可再生能源消纳、储能运行、可控电源运行均为固定单位成本需求响应补偿价格与大电网峰谷电价匹配高峰时段补偿价格更高。微电网设备参数根据工程实际配置储能额定容量 680kW・h额定充放电功率 ±200kW微型燃气轮机额定功率 300kW最小技术出力 0kW柴油发电机额定功率 500kW最小技术出力 0kW可转移负荷最大调节功率 500kW。储能初始荷电状态 0.5荷电状态安全范围 0.2-0.9。5.2 综合成本对比综合成本是衡量算法优化性能的核心指标PSO 和 MVO 算法均能有效求解微电网日前经济调度模型得到综合成本最小的最优解。两种算法的优化结果均显著低于无优化的常规调度策略表明智能优化算法能有效降低微电网的运行综合成本。对比两种算法的最优综合成本其中一种算法略优于另一种算法体现了不同寻优机制对优化结果的影响但整体差异较小表明两种算法均具有较好的全局寻优能力。5.3 成本构成对比从成本构成来看两种算法的优化结果均呈现出相同的成本结构特征运行成本占综合成本的比重最大是微电网运行的主要成本占比超过 90%环保成本次之占比约 5%-8%需求响应补偿成本和负荷波动惩罚成本占比极小几乎可忽略。这一成本结构特征反映了微电网日前运行的实际情况运行成本是调度优化的核心对象。对比两种算法的各分项成本运行成本和环保成本的差异是导致综合成本差异的主要原因其中一种算法在降低化石能源出力和大电网购电成本方面表现更优从而实现了更低的运行成本和环保成本两种算法的需求响应补偿成本和负荷波动惩罚成本基本一致表明在负荷调节和平抑负荷波动方面两种算法的调度策略相近。5.4 收敛速度对比收敛速度是衡量算法效率的重要指标采用首次达到全局最优值的迭代步数作为评价标准。PSO 算法的收敛步数显著少于 MVO 算法表明 PSO 算法在微电网日前经济调度问题中具有更快的收敛速度能在更少的迭代次数内找到全局最优解。这是因为 PSO 算法的原理简单寻优机制直接粒子能快速向最优解靠近而 MVO 算法的寻优机制更为复杂需要更多的迭代次数完成全局搜索和局部精细搜索。从收敛曲线来看PSO 算法在迭代初期快速下降随后进入平稳阶段收敛特性较好MVO 算法在迭代初期下降较慢迭代中期下降速度加快后期趋于平稳全局寻优的探索性更强。两种算法的收敛曲线最终均趋于平稳表明均能有效收敛至全局最优解无明显的局部最优现象。5.5 调度策略对比从调度策略来看两种算法均能结合峰谷电价、可再生能源出力和负荷特性制定合理的调度策略核心调度规律一致储能调度在大电网低谷电价时段和可再生能源出力过剩时段储能进行充电在高峰电价时段、可再生能源出力不足时段和负荷高峰时段储能进行放电充分利用峰谷电价差降低购电成本平抑可再生能源出力波动。可控电源调度微型燃气轮机作为主要的化石能源补充电源在负荷平段和可再生能源出力略有不足时出力柴油发电机仅在负荷高峰和可再生能源出力严重不足时少量出力因柴油发电机的运行成本和环保成本更高调度策略尽可能减少其使用。可转移负荷调度在负荷高峰和大电网高峰电价时段对可转移负荷进行削峰调节降低高峰时段的供电压力在负荷低谷和低谷电价时段对可转移负荷进行填谷调节提高可再生能源消纳比例实现负荷曲线的优化。对比两种算法的具体调度策略差异主要体现在储能充放电功率的精细化调节和可控电源的出力分配上其中一种算法的储能充放电功率更贴合电价和负荷的变化特性可控电源的出力分配更合理从而实现了更低的综合成本。整体而言两种算法的调度策略均符合微电网日前经济调度的基本规律具有较强的工程实用性。6 结论与展望6.1 研究结论本文以微电网日前经济调度为研究对象构建了多成本融合的综合成本最小化优化模型设计了统一的约束处理机制分别采用 PSO 和 MVO 算法对模型进行求解并从多维度开展了算法性能对比分析得到以下主要结论微电网日前经济调度模型能有效兼顾系统运行的经济性、环保性、需求响应特性和稳定性通过对可控设备和可转移负荷的协同优化可显著降低微电网的运行综合成本提升运行效益。PSO 和 MVO 两种智能优化算法均能有效求解微电网日前经济调度这一多变量、多约束的非线性优化问题均能收敛至全局最优解无明显的局部最优现象为该问题提供了高效的求解方法。PSO 算法具有收敛速度快的优势能在更少的迭代次数内找到全局最优解求解效率更高适用于对求解速度有要求的场景MVO 算法具有更强的全局探索能力寻优机制更灵活在部分场景下能得到更优的综合成本适用于对优化精度要求较高的场景。两种算法的最优调度策略均符合微电网运行的基本规律能结合峰谷电价、可再生能源出力和负荷特性实现多设备协同优化调度结果具有较强的工程实用性可为微电网日前调度的工程实践提供决策参考。6.2 研究展望本文的研究为微电网日前经济调度提供了基础的求解框架和算法参考后续可从以下方面开展进一步的研究多目标优化研究将单目标的综合成本最小化拓展为经济、环保、节能、供电可靠性等多目标优化采用帕累托最优解进行求解更贴合微电网多目标协同运行的实际需求。不确定性优化研究考虑风电、光伏出力的预测不确定性和负荷的不确定性构建鲁棒优化模型或随机优化模型提高调度策略的抗干扰能力和实际适应性。算法改进与融合研究对 PSO 和 MVO 算法进行改进如引入自适应参数、混合寻优机制等提升算法的寻优性能或将多种算法进行融合结合不同算法的优势构建混合优化算法。多时间尺度调度研究将日前经济调度与日内滚动调度、实时调度相结合构建多时间尺度的协同调度框架提高调度的实时性和精准性。多主体参与的调度研究考虑微电网运营方、分布式电源投资方、用户等多主体的利益诉求构建市场化的调度模型引入博弈论等方法实现多主体利益的协同优化。微电网日前经济调度是微电网运行管理的核心环节随着分布式能源的进一步发展和电力市场的不断完善其优化模型和求解方法将不断丰富和完善。智能优化算法作为求解微电网优化调度问题的有效工具将在未来的研究和工程实践中得到更广泛的应用。第二部分——运行结果主函数代码% 微电网日前经济调度PSO vs MVO 联合求解与对比% 功能一键运行两种算法输出数值对比可视化对比clc;clear;close all;%% 统一参数设置两种算法共用% 1. 微电网参数Qbattery680;% 储能容量(kw.h)Pbatterymax200;% 电池最大功率(放电)Pbatterymin-200;% 电池最小功率(充电)Pwtmax300;% 微型燃气轮机最大功率Pwtmin0;% 微型燃气轮机最小功率Pdgmax500;% 柴油发电机最大功率Pdgmin0;% 柴油发电机最小功率Plmax500;% 可转移负荷最大功率Plmin0;% 可转移负荷最小功率SOCstart0.5;% 储能初始荷电状态SOCmin0.2;% 储能最小荷电状态SOCmax0.9;% 储能最大荷电状态Pmax[Pbatterymax,Pwtmax,Pdgmax,Plmax];Pmin[Pbatterymin,Pwtmin,Pdgmin,Plmin];% 2. 算法通用参数保证公平对比MaxItr5000;% 迭代次数两种算法一致NB100;% 种群规模PSO粒子数/MVO宇宙数一致D4*24;% 决策变量维度4个变量×24小时% 3. 加载数据共用Data.mData;% 加载负荷、风光、电价数据% 4. 决策变量边界BoundUpzeros(1,D);BoundLowzeros(1,D);for i1:4BoundUp((i-1)*241:i*24)Pmax(i);BoundLow((i-1)*241:i*24)Pmin(i);end%% 1. 运行粒子群优化算法(PSO) fprintf(正在运行粒子群优化算法(PSO)...\n);[PSO_BestTrace, PSO_BestPos, PSO_BestFit, PSO_f1, PSO_f2, PSO_f3, PSO_Pgrid] PSO_Main(MaxItr, NB, Qbattery, BoundUp, BoundLow, SOCstart, SOCmin, SOCmax, PL, Pwind, Prices, Ppv);%% 2. 运行多元宇宙优化算法(MVO) fprintf(正在运行多元宇宙优化算法(MVO)...\n);[MVO_BestTrace, MVO_BestPos, MVO_BestFit, MVO_f1, MVO_f2, MVO_f3, MVO_Pgrid] MVO_Main(MaxItr, NB, Qbattery, BoundUp, BoundLow, SOCstart, SOCmin, SOCmax, PL, Pwind, Prices, Ppv);%% 3. 结果对比与输出 % 3.1 数值对比表fprintf(\n PSO vs MVO 求解结果对比 \n);comparison_table table(...{PSO; MVO}, ...[PSO_BestFit; MVO_BestFit], ...[PSO_f1; MVO_f1], ...[PSO_f2; MVO_f2], ...[PSO_f3; MVO_f3], ...VariableNames, {算法, 总成本(元), 运行成本(元), 环保成本(元), 需求响应成本(元)});disp(comparison_table);% 3.2 性能指标计算% 收敛速度达到最优值的迭代步数PSO_ConvStep find(PSO_BestTracePSO_BestFit1e-6, 1);% 找到首次达到最优值的步数MVO_ConvStep find(MVO_BestTraceMVO_BestFit1e-6, 1);fprintf(\n收敛速度对比\n);fprintf(PSO收敛步数%d | MVO收敛步数%d\n, PSO_ConvStep, MVO_ConvStep);% 3.3 可视化对比% 子图1迭代收敛曲线对比figure(Position, [100, 100, 1200, 800]);subplot(2,2,1);plot(1:MaxItr, PSO_BestTrace, b-, LineWidth, 1.5); hold on;plot(1:MaxItr, MVO_BestTrace, r-, LineWidth, 1.5);xlabel(迭代次数); ylabel(总成本(元));title(PSO vs MVO 收敛曲线对比);legend(PSO, MVO, Location, best);grid on;% 子图224小时储能功率对比subplot(2,2,2);plot(1:24, PSO_BestPos(:,1), b-o, LineWidth, 1.2); hold on;plot(1:24, MVO_BestPos(:,1), r-s, LineWidth, 1.2);xlabel(时刻(h)); ylabel(储能功率(kW));title(24小时储能充放电功率对比);legend(PSO, MVO, Location, best);grid on;% 子图324小时总成本构成对比subplot(2,2,3);bar([PSO_f1, PSO_f2, PSO_f3; MVO_f1, MVO_f2, MVO_f3]);xlabel(成本类型); ylabel(成本(元));title(总成本构成对比);xticklabels({运行成本, 环保成本, 需求响应成本});legend(PSO, MVO, Location, best);grid on;% 子图424小时电网交互功率对比subplot(2,2,4);plot(1:24, PSO_Pgrid, b-^, LineWidth, 1.2); hold on;plot(1:24, MVO_Pgrid, r-, LineWidth, 1.2);xlabel(时刻(h)); ylabel(电网交互功率(kW));title(24小时电网交互功率对比);legend(PSO, MVO, Location, best);grid on;% 3.4 最优决策变量保存便于后续分析save(PSO_MVO_Results.mat, ...PSO_BestPos, MVO_BestPos, ...PSO_BestFit, MVO_BestFit, ...PSO_ConvStep, MVO_ConvStep);fprintf(\n结果已保存至 PSO_MVO_Results.mat 文件\n);第三部分——参考文献文章中一些内容引自网络会注明出处或引用为参考文献难免有未尽之处如有不妥请随时联系删除。(文章内容仅供参考具体效果以运行结果为准)第四部分——本文完整资源下载资料获取更多粉丝福利MATLAB|Simulink|Python|数据|文档等完整资源获取
