不平衡电压跌落场景下分布式并网变流器序电流多目标优化 LVRT 控制方法（Matlab代码、Simulink实现）

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

🎁完整资源、论文复现、期刊合作、论文辅导及科研仿真定制事宜点击：

👉👉👉本文完整资源下载

⛳️座右铭：行百里者，半于九十。

⛳️赠与读者

👨‍💻做科研，涉及到一个深在的思想系统，需要科研者逻辑缜密，踏实认真，但是不能只是努力，很多时候借力比努力更重要，然后还要有仰望星空的创新点和启发点。建议读者按目录次序逐一浏览，免得骤然跌入幽暗的迷宫找不到来时的路，它不足为你揭示全部问题的答案，但若能解答你胸中升起的一朵朵疑云，也未尝不会酿成晚霞斑斓的别一番景致，万一它给你带来了一场精神世界的苦雨，那就借机洗刷一下原来存放在那儿的“躺平”上的尘埃吧。

或许，雨过云收，神驰的天地更清朗.......🔎🔎🔎

💥第一部分——内容介绍

不平衡电网下变流器型分布式电源多目标优化

一、文献基础概况与行业研究背景

（一）文献基础信息

论文聚焦不平衡电网故障穿越场景，面向光伏、风机等变流器接口分布式电源（DG），提出一套完整、可实时在线运算的多目标最优参考电流生成控制方案，全部性能指标与最优控制参数均具备解析闭式解，无需数值迭代寻优。

（二）行业研究背景扩充

全球新能源装机规模持续攀升，分布式光伏、直驱风机全功率变流器大量接入配电网，电网故障下的系统稳定问题日益凸显。各国电网公司陆续更新严苛并网规范，核心要求集中在低电压穿越（LVRT）：当电网发生单相接地、相间短路等不平衡电压跌落故障时，并网变流器不得主动脱网，还需主动向电网注入序电流，完成电压抬升、频率支撑。

在工程实际中，不平衡电压会带来多重负面效应：正负序电压与电流耦合产生二倍频功率振荡，一方面造成直流母线电压剧烈波动，增大直流电容损耗、缩短设备寿命；另一方面诱发交流侧功率振荡，加剧弱电网谐振风险。与此同时，不平衡故障会导致三相输出电流严重不对称，某一相电流极易突破变流器功率器件额定电流阈值，触发过流保护，造成新能源机组大面积脱网，进一步恶化电网功率缺额。

过往研究的控制方案存在明显短板：传统参考电流生成策略仅拥有两个及以内控制自由度，只能单一完成功率振荡抑制或无功电压支撑，无法同步兼顾多重约束；现有最优控制算法高度依赖离线仿真迭代，计算量大、延时高，难以在工业级低算力 DSP 控制器中实时运行；正负序有功、无功电流分配逻辑相互割裂，无法协同实现频率支撑与电压治理双重需求。

针对上述痛点，本文对经典正负序参考电流生成（RCG）架构进行拓展，引入有功分配系数、无功分配系数、有功功率参考、无功功率参考四维独立控制自由度，搭建完整的不平衡工况数学分析体系，推导不同优化目标下的最优控制参数解析计算逻辑，形成五类单目标最优控制策略与一套兼顾电压约束、最大功率输出的复合协同控制方案，并通过完整仿真平台与实物样机完成全部方案有效性验证。

二、不平衡工况下参考电流生成（RCG）基础理论体系

（一）正负序电流四维分解架构

论文基于瞬时无功功率理论，将不平衡电网下并网变流器总输出电流拆分为有功、无功两大基础分量，再进一步拆解为正序、负序两组子分量，由此形成四组独立可控电流通道。四组电流分别对应正序有功、负序有功、正序无功、负序无功，通过两个分配系数调节有功、无功功率在正负序通道之间的分配比例，再配合有功、无功功率基准指令，构成完整四维控制自由度，具备极高控制灵活性。

其中有功分配系数负责调节总有功功率在正序、负序通路的分配权重，无功分配系数同理管控总无功功率的序分量配比，两个系数的取值区间均为 0 到 1：系数取 1 代表全部功率仅通过正序通道输出，取 0 则全部功率经由负序通道注入电网。正交电压矢量作为无功电流生成基准，与原电压矢量保持 90 度相位差，以此实现有功、无功电流完全解耦，互不干扰。该四维架构是全文所有优化策略的底层载体，所有多目标优化逻辑均通过调整四维自由度参数实现。

（二）不平衡工况功率振荡机理分析

电网不平衡的本质是电网同时存在正序、负序两组旋转电压矢量，两组矢量频率相同、旋转方向相反，二者与变流器输出的正负序电流发生交叉耦合，最终在瞬时有功、无功功率上产生 2 倍工频周期性振荡。振荡幅值直接决定直流母线纹波大小与交流系统稳定程度，也是故障工况下需要重点抑制的核心指标。

论文通过矢量点积完整推导得到有功、无功功率最大振荡幅值的量化计算逻辑，仅依靠实时采集的有功功率、无功功率、两组分配系数、正负序电压幅值比值，即可瞬时计算当前工况下功率波动的峰值，无需复杂傅里叶分解、滤波运算。该量化模型为后续最小有功振荡（MOP）、最小无功振荡（MOQ）两类优化策略提供核心优化目标函数，是抑制直流侧损耗、规避电网谐振的理论基础。

（三）单相接地故障三相电流量化模型

论文选取配电网最高发故障 ——A 型单相接地不平衡故障作为核心研究工况，推导静止 αβ 坐标系下正负序标准电压矢量表达式，同步给出用于生成无功电流的正交电压矢量数学形式。借助克拉克坐标变换，将 αβ 静止坐标系合成电流转换为工程可直接观测的 abc 三相实际输出电流，提取各相电流瞬时幅值，以此分析故障下三相电流不对称程度。

为简化三相电流峰值的复杂三角函数表达式，论文定义四组复合中间等效参数，整合正负序电压幅值、有功无功功率、两组分配系数等全部实时变量，大幅简化三相电流峰值的计算复杂度。基于等效参数得到三相各自电流峰值的完整表达式，取三相峰值中的最大值作为故障工况下变流器全局最大相电流，该指标直接绑定功率器件额定电流限值，是最小故障电流（MFC）、最大允许有功（MAP）、最大允许无功（MAQ）三类限流型控制策略的核心约束条件。

论文配套给出多组不同有功、无功、无功分配系数工况下，三相电流峰值随有功分配系数变化的可视化曲线。曲线直观反映：三相电流峰值随有功分配系数呈现非线性变化，全局最大相电流的最小值出现在三条相电流曲线的交点或单条曲线谷底处，为后文多候选点寻优的 MFC 策略提供可视化理论支撑。

三、五类单目标解析最优控制策略设计与工程落地说明

论文依托前文完整的量化数学模型，针对性设计五套独立单目标最优控制方案，所有方案最优参数均通过解析求极值得到，不存在迭代运算，计算延时极低，适配新能源变流器毫秒级实时控制需求，每套策略拥有明确的适用故障场景与工程优化价值。

（一）MOP 最小有功功率振荡策略

该策略以削减瞬时有功功率二倍频振荡幅值为唯一优化目标，通过对振荡幅值表达式求取偏导、令导数等于零，推导出无约束理论最优有功、无功分配系数。从理论解可以看出，不平衡故障下有功分配系数的计算值通常会超出 0~1 的物理可行区间，无法直接投入工程使用。

针对该问题论文给出标准化工程修正方案：将超出上限的有功分配系数强制截断为 1，即故障期间仅向电网注入正序有功电流，仅保留无功分配系数采用理论最优解析值。投入该策略后，直流母线二倍频电压纹波会大幅衰减，能够降低直流侧滤波电容容量选型标准，减少电容发热损耗，提升全功率变流器整体运行效率。仿真与实验波形清晰对比故障前后、投入策略前后的有功功率波形，二倍频脉动几乎完全消除，优化效果直观显著。

（二）MOQ 最小无功功率振荡策略

该策略以最小化瞬时无功功率振荡幅值为优化目标，同样通过求导极值法得到无约束最优分配系数解析解。与 MOP 策略规律一致，故障工况下无功分配系数理论计算值普遍超出 0~1 区间，工程落地时强制限定无功分配系数为 1，仅优化有功分配系数抑制无功波动。

该策略主要适配高阻抗弱电网场景：弱电网中无功功率振荡极易引发公共耦合点（PCC）电压闪变，造成厂区、居民用电设备频闪故障，投入 MOQ 策略能够平稳交流侧电压波形，提升配电网供电质量。实物实验中，故障初期无功功率存在剧烈周期性起伏，激活 MOQ 策略后无功脉动快速收敛，验证方案的振荡抑制能力。

（三）MFC 最小故障相电流策略

该策略核心目标为降低不平衡故障下变流器输出最大相电流，避免功率器件过流损坏、机组保护脱网。策略运行前提为：无功分配系数、有功无功功率基准由上层电压支撑控制器给定，仅将有功分配系数作为唯一优化变量。

完整寻优流程分为三步：首先分别推导单独使 A、B、C 三相电流最小的有功分配系数解析解；其次联立两相电流峰值方程，求解两条电流曲线交点对应的分配系数，同时补充分配系数区间边界值，共计六组候选优化点；最后将六组候选参数依次代入全局最大相电流计算模型，选取使三相最大电流最小的有功分配系数作为实时控制参数。

实物实验数据充分验证方案价值：未启用 MFC 优化时，故障峰值电流达到 9.5A；激活策略、切换至最优有功分配系数后，故障峰值电流下降至 8A，限流效果明显，能够显著降低故障冲击对变流器功率器件的热应力，提升机组故障穿越可靠性。

（四）MAP 电流约束下最大允许有功输出策略

该策略面向电网故障后频率跌落场景，核心需求为在不突破变流器额定电流限值的前提下，向电网输送尽可能多的有功功率，支撑系统频率稳定。策略固定无功功率指令，以三相电流全部低于额定限流为硬性约束，分别根据 A、B、C 三相电流极限方程求解各相单独约束下的有功功率上限。

为保障三相电流全部满足安全限值，将三组有功上限数值取最小值作为变流器当前工况下安全输出的最大有功基准。在风电、光伏电站大面积故障场景中，该策略可以充分挖掘变流器有功输送潜力，缓解电网有功缺额，抑制频率持续跌落，同时严格遵守设备热保护阈值，兼顾调频需求与设备运行安全。仿真多梯度电压跌落工况下，投入 MAP 策略后有功功率始终维持在电流限值对应的最大值，无任何一相电流超限。

（五）MAQ 电流约束下最大允许无功输出策略

该策略适配深度电压跌落故障场景，匹配各国电网 LVRT 导则强制无功支撑要求：固定有功功率输出，在额定电流约束下输出极限无功功率，最大化抬升 PCC 故障相电压。策略推导逻辑与 MAP 完全对称，通过三相电流极限方程分别求解单通道无功上限，三组数值取最小值作为当前工况最大安全无功基准。

当电网发生严重单相短路、故障相电压大幅跌落时，MAQ 策略可释放变流器全部剩余电流容量增发无功，快速抬高故障相电压，避免电压持续过低触发下级负荷脱网。实物测试中，故障后 PCC 电压出现明显跌落，启用 MAQ 策略增发极限无功后，PCC 电压快速回升至安全区间，同时三相输出电流稳定维持在预设限流值以内。

四、复合多目标协同控制策略 VSS-MAP：电压支撑与最大功率联合优化

前述五类单目标策略仅能实现单一性能指标优化，工程实际故障工况下变流器需要同步满足三重硬性约束：公共耦合点三相相电压维持在电网规范限定的安全区间（0.9~1.1 标幺值）、三相输出电流不超过设备额定限值、尽可能增发有功功率支撑电网频率。针对多重约束耦合需求，论文提出 VSS-MAP 复合分层协同控制架构，融合电压支撑模块（VSS）与最大有功输出模块（MAP），同步完成电压治理、故障限流、有功功率最大化三重目标。

（一）PCC 正负序电压支撑底层模型

建模过程完整计及配网线路阻感阻抗，推导公共耦合点正负序电压偏差与变流器注入正负序电流的耦合关系，分别区分电阻、电感对电压调节的贡献分量。针对工业配电网普遍高感抗特性，线路电感远大于电阻，电阻对电压调节的影响可近似忽略，由此得到简化纯无功调压模型：仅依靠正负序无功电流即可完成 PCC 序电压调节，有功电流不再参与调压，有功通道可完全用于输送有功功率支撑频率。

论文给出简化模型对应的离散控制框图：分别采集当前正负序电压与目标参考序电压的差值，经线路电感增益换算后，搭配 5 毫秒低通滤波环节抑制高频噪声，直接输出正负序无功电流参考值，控制结构简单、动态响应速度快，易于在嵌入式控制器中实现。

（二）三相相电压安全约束求解逻辑

不平衡工况下三相相电压幅值由正负序电压幅值、正负序矢量相位差共同决定，论文推导出 A、B、C 三相电压幅值完整表达式，定义三相电压余弦相位项的最大值、最小值，以此快速求解故障下三相电压全局最高值与全局最低值。

依据并网导则设定电压安全上下边界，将三相最低相电压目标值固定为 0.9 标幺值，最高相电压目标值兼顾上限 1.1 标幺值与三相电压不平衡差值约束，得到标准化电压参考边界计算规则。反向联立电压极值方程，解析求解能够使三相电压全部落在安全区间内的目标正序、负序参考电压，为电压支撑模块提供精准调节目标。

（三）VSS-MAP 分层协同完整控制流程

整套复合控制架构分为五层递进运算逻辑，层层解耦、实时在线运算：

1. 采样序分离层：实时采集 PCC 三相电压，通过序分量分离算法计算当前正负序电压幅值、矢量相位差；
2. 电压约束求解层：代入三相电压极值模型，计算满足 0.9~1.1 标幺值区间约束的目标正负序参考电压；
3. 无功调压运算层：依托高感抗简化调压模型，计算抬升电压所需的正负序无功电流，换算得到总无功功率与无功分配系数；
4. 极限有功求解层：将无功功率、无功分配系数输入 MAP 策略计算逻辑，求解当前电流限值下可输出的最大有功功率基准；
5. 参考电流生成层：整合最优有功分配系数、无功分配系数、极限有功、调压无功四组参数，送入四维 RCG 电流生成架构，输出四序电流参考指令，送入双 dq 坐标系 PI 电流闭环控制环完成执行。

多梯度深度单相电压跌落仿真测试结果证明：在不同程度电压跌落故障全程，VSS-MAP 策略可稳定将三相相电压钳位在 0.9 至 1.1 标幺值安全区间，有功功率持续保持在额定电流允许的最大值，三相输出电流均未突破限流阈值，三重冲突优化目标同步达成，完美适配复杂电网故障规范要求。

五、仿真平台与实物实验验证体系完整分析

（一）数字仿真平台与测试用例设计

论文搭建 210kVA 工业级三相并网电压源变流器仿真模型，额定线电压 690V、工频 60Hz；直流侧采用理想直流源等效光伏、风机新能源输入，输出端搭配 LC 滤波电路与升压变压器接入大电网，通过故障阻抗模拟 A 型单相接地不平衡跌落故障。电流控制采用正负序独立双 dq 坐标系 PI 调节器，配套完整锁相环模块，电流环、锁相环比例积分参数全部给出标准化设计数值，仿真模型具备完全可复现性。

仿真设计四类差异化测试用例，分别验证不同策略的优化性能：

1. 测试案例 A：验证 MOP、MOQ 功率振荡抑制策略，设置两段不同深度单相电压跌落，对比无优化、投入 MOP、投入 MOQ 三种工况下有功、无功功率波动幅值；
2. 测试案例 B：验证 MFC 最小故障电流策略，对比优化前后三相输出电流峰值，直观体现限流优化效果；
3. 测试案例 C：验证 MAP 最大有功输出策略，分多档无功功率、无功分配系数工况，观测电压跌落过程中有功功率极限输出特性；
4. 测试案例 D：验证 VSS-MAP 复合协同控制策略，四段逐级加深的单相电压跌落工况，同步观测三相电压稳定能力、有功功率输出上限、三相电流限流效果。

全部仿真波形从上层至下层依次展示公共耦合点三相电压、正负序序电流、abc 三相输出电流、瞬时有功无功功率，分层清晰、便于对比优化前后性能差异。

（二）实物实验平台与测试结果说明

为验证解析理论与仿真结论的实际工程可行性，论文搭建基于 dSPACE1104 实时控制器的 1500VA 小功率实物样机，配套 TMS320F240 数字信号协处理器生成开关驱动脉冲，变流器开关频率 10kHz，通过降压变压器接入 110V 相电压 60Hz 工频实验电网，串联故障阻抗模拟单相接地不平衡故障。

实验分为四组独立测试，分别验证 MOP、MOQ、MFC、MAQ、MAP 五类单目标策略：

1. MOP/MOQ 实验：观测故障后功率振荡、投入优化策略后的振荡收敛过程，验证二倍频脉动抑制效果；
2. MFC 实验：对比优化前后三相电流峰值，量化故障限流优化收益；
3. MAQ 实验：固定有功功率，观测增发极限无功后 PCC 电压抬升效果与电流限值约束；
4. MAP 实验：模拟电网导则无功支撑需求，在固定无功基础上输出最大有功，验证电流边界约束有效性。

全部实验波形变化趋势与仿真结果高度匹配，定量数据误差极小，充分证明论文整套解析计算逻辑、控制策略不存在理论与工程偏差，能够直接落地到实际新能源并网变流器控制器中使用。

六、论文核心创新价值总结

1. 底层理论建模创新：首次构建不平衡单相接地故障下四维自由度参考电流生成完整解析体系，推导出功率振荡幅值、三相故障电流峰值全套量化计算逻辑，摆脱传统控制方案依赖离线数值迭代的缺陷，全部优化参数可微秒级在线实时运算，适配低算力工业控制器；
2. 分层单目标控制体系创新：搭建五套标准化单目标最优控制策略，分别对应功率振荡抑制、故障电流限流、极限有功 / 无功馈送四大电网核心需求，每套策略配套完整理论最优解与工程边界修正方案，模块化程度高，可根据电网故障类型、并网需求自由切换；
3. 多约束协同控制创新：原创 VSS-MAP 复合分层控制架构，首次实现 PCC 三相电压区间约束、变流器额定电流热限值、最大化有功频率支撑三重相互制约目标同步优化，完美匹配当下严苛、多维度的电网 LVRT 并网规范；
4. 工程落地适配性优势：整套控制逻辑仅依靠基础代数运算完成参数求解，无复杂矩阵、高阶迭代运算，不需要高端浮点运算芯片，中小功率光伏、储能、风机变流器均可低成本搭载；
5. 电网工程应用价值：能够大幅降低新能源机组故障脱网概率，减小直流滤波电容容量与发热损耗，提升弱电网电压稳定水平，有效提高分布式新能源并网消纳渗透率，方案适配全类型变流器接口分布式电源。

七、现有研究局限性与后续拓展研究方向

（一）论文现有研究局限

1. 故障覆盖范围有限：全文完整推导、验证仅针对单相接地不平衡故障，未针对两相短路、两相接地等其余常见不平衡故障类型推导对应量化模型与最优控制参数求解逻辑；
2. 约束优化方式简化处理：当有功、无功分配系数理论最优解超出 0~1 可行区间时，仅采用简单边界截断处理，未引入带约束二次规划算法求解区间内全局折中最优分配系数，优化潜力未完全释放；
3. 系统动态耦合因素忽略：建模与优化过程中未考虑锁相环动态响应、电网谐波畸变、变流器开关死区非线性等实际扰动对最优分配参数的耦合干扰，复杂畸变电网下控制精度存在提升空间；
4. 优化目标单一割裂：单目标控制策略仅单独优化一项性能指标，未构建多目标加权代价函数，无法根据电网工况动态权衡功率振荡、故障电流、电压偏差多项指标。

（二）可延伸拓展研究方向

1. 多故障通用解析模型拓展：推导两相短路、两相接地故障下功率振荡、三相故障电流统一量化表达式，形成覆盖全部配网不平衡故障类型的通用最优控制框架；
2. 带约束多目标加权优化：构建包含功率振荡幅值、故障电流峰值、PCC 电压偏差的加权综合代价函数，引入约束优化算法求解区间内全局最优分配系数，实现多指标动态折中；
3. 孤岛微网工况适配改造：将整套四维 RCG 多目标优化策略拓展至孤岛微网带不平衡负载场景，实现离网状态下三相电压平衡、负载功率稳定控制；
4. 模型预测控制融合优化：将本文解析最优解作为模型预测控制（MPC）滚动优化初始值，大幅削减 MPC 在线迭代计算量，兼顾控制精度与运算速度；
5. 风光储协同故障支撑：结合储能变流器快速有功吞吐特性，实现风机、光伏、储能多机协同，故障期间同步完成电压抬升与频率稳定双重支撑。

📚第二部分——运行结果

所有模型和讲解见第4部分。

set(0,'DefaultAxesFontSize',24,'DefaultAxesFontName','Helvetica','DefaultAxesFontWeight','bold','DefaultLineLineWidth',2,'DefaultAxesLineWidth',1)
plot(k1,Imax,'--','color',[0 .45 .74],'LineWidth',6)
grid on
hold on
plot(k1,sqrt(Ia),'color',[0.75 0 .75])
plot(k1,sqrt(Ib),'color',[.85 .33 .1])
plot(k1,sqrt(Ic),'color',[0 .5 0])

legend('Imax','Ia','Ib','Ic')
hold on
plot(k1_opt(O),min(Imax_opt),'s','color',[0 .45 .74],'MarkerSize',25,'MarkerFaceColor',[0 .45 .74])

plot(k1(J),min(Imax),'o','color',[.85 .33 .1],'LineWidth',2,'MarkerSize',15)

title(['P=', num2str(P(p)), ' pu Q=', num2str(Q(q)), ' pu Vp=', num2str(Vp(vp)), ' pu Vn=', num2str(Vn(vn)), ' pu k_q=', num2str(k2(k)),' pu k_p_,_o_p_t= ', num2str(k1_opt(O(1))), ' I_m_a_x_,_o_p_t= ', num2str(min(Imax_opt)), ' pu'])

xlabel('k_p');
ylabel('I_m_a_x (pu)');

%%
A=3*P(p)^2*n; B=-3*P(p)^2*n+sqrt(3)*P(p)*Q(q)*n*(2*k2(k)-1); C=3*n*k2(k)*Q(q)^2*(1-k2(k))-sqrt(3)*P(p)*Q(q)*n*k2(k);
delta=B^2-4*A*C;
if delta>=0
k1_intsec_IaIb=(-B+sqrt(delta))/2/A;
o_num=4;
k1_opt(o_num)=k1_intsec_IaIb;


if k1_intsec_IaIb<1 && k1_intsec_IaIb>0

K1_intsec=P(p)/Vn(vn)*((n+1)*k1_intsec_IaIb-1);
K2_intsec=Q(q)/Vn(vn)*((n-1)*k2(k)+1);
K3_intsec=P(p)/Vn(vn)*((n-1)*k1_intsec_IaIb+1);
K4_intsec=Q(q)/Vn(vn)*((n+1)*k2(k)-1);

Ia_intsec=K1_intsec^2+K2_intsec^2;
Ic_intsec=(-.5*K1_intsec-sqrt(3)/2*K4_intsec)^2+(.5*K2_intsec-sqrt(3)/2*K3_intsec)^2;

plot(k1_intsec_IaIb,sqrt(Ia_intsec),'s','color',[0 .45 .74],'MarkerSize',25)

Imax_intsec=sqrt(max(Ia_intsec,Ic_intsec));
end
end



%%
A=3*P(p)^2*n; B=-3*P(p)^2*n-sqrt(3)*P(p)*Q(q)*n*(2*k2(k)-1); C=3*n*k2(k)*Q(q)^2*(1-k2(k))+sqrt(3)*P(p)*Q(q)*n*k2(k);
delta_IaIc=B^2-4*A*C;
if delta_IaIc>=0
k1_intsec_IaIc=(-B+sqrt(delta_IaIc))/2/A;

if k1_intsec_IaIc<1 && k1_intsec_IaIc>0

K1_intsec=P(p)/Vn(vn)*((n+1)*k1_intsec_IaIc-1);
K2_intsec=Q(q)/Vn(vn)*((n-1)*k2(k)+1);
K3_intsec=P(p)/Vn(vn)*((n-1)*k1_intsec_IaIc+1);
K4_intsec=Q(q)/Vn(vn)*((n+1)*k2(k)-1);

Ia_intsec=K1_intsec^2+K2_intsec^2;
Ib_intsec=(-.5*K1_intsec+sqrt(3)/2*K4_intsec)^2+(.5*K2_intsec+sqrt(3)/2*K3_intsec)^2;

plot(k1_intsec_IaIc,sqrt(Ia_intsec),'s','color',[0 .45 .74],'MarkerSize',25)

end

end
end

i=i+1;

end
end
end
end
% end

🎉第三部分——参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。(文章内容仅供参考，具体效果以运行结果为准)

[1]M. M. Shabestary and Y. A. -R. I. Mohamed, "Analytical Expressions for Multiobjective Optimization of Converter-Based DG Operation Under Unbalanced Grid Conditions," in IEEE Transactions on Power Electronics, vol. 32, no. 9, pp. 7284-7296, Sept. 2017, doi: 10.1109/TPEL.2016.2628405.

[2]谭畅舒,李艳,田杰,明威宇,李妍,王少荣.计及分布式能源的主动配电网恢复力综合评估[J].电力科学与技术学报,2023,38(01):108-113.DOI:10.19781/j.issn.1673-9140.2023.01.012.

[3]葛钦. 电网故障下电力电子变压器故障穿越机理与控制策略研究[D].湖南大学,2020.DOI:10.27135/d.cnki.ghudu.2020.000151.

[4]庞永恒. 基于故障失配度的含新能源配电网接地故障定位方法研究[D].东北大学,2017.

​​​​​​🌈第四部分——本文完整资源下载

资料获取，更多粉丝福利，MATLAB|Simulink|Python|数据|文档等完整资源获取

本文完整资源下载