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

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，方便大家进行学习！亲民！！！还有我开了一个专栏给女朋友的，很浪漫的喔，代码学累的时候去瞧一瞧，看一看：女朋友的浪漫邂逅。有问题可以私密博主，博主看到会在第一时间回复。
📝目前更新：🌟🌟🌟电力系统相关知识，期刊论文，算法，机器学习和人工智能学习。
🚀支持：🎁🎁🎁如果觉得博主的文章还不错或者您用得到的话，可以关注一下博主，如果三连收藏支持就更好啦！这就是给予我最大的支持！

                           ​

👨‍🎓博主课外兴趣：中西方哲学，送予读者：

👨‍💻做科研，涉及到一个深在的思想系统，需要科研者逻辑缜密，踏实认真，但是不能只是努力，很多时候借力比努力更重要，然后还要有仰望星空的创新点和启发点。当哲学课上老师问你什么是科学，什么是电的时候，不要觉得这些问题搞笑，哲学就是追究终极问题，寻找那些不言自明只有小孩子会问的但是你却回答不出来的问题。在我这个专栏记录我有空时的一些哲学思考和科研笔记：科研和哲思。建议读者按目录次序逐一浏览，免得骤然跌入幽暗的迷宫找不到来时的路，它不足为你揭示全部问题的答案，但若能让人胸中升起一朵朵疑云，也未尝不会酿成晚霞斑斓的别一番景致，万一它居然给你带来了一场精神世界的苦雨，那就借机洗刷一下原来存放在那儿的“真理”上的尘埃吧。

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

​

📋📋📋本文目录如下：⛳️⛳️⛳️

目录

1 概述

2 数学模型 

3 运行结果

4 Matlab代码+数据+文档讲解 

---

1 概述

未来的配电网预计将承载大量分布式微电网，并满足日益增长的需求，例如，受电动汽车更广泛普及的推动。这些网格预计会表现出拥塞现象，仅通过基于最坏情况分析设计超大型网络是无法充分解决的。例如，这些低压和中压网络的电压分布将受到双向有功潮流的影响，预计过压和欠压情况将越来越频繁地发生。

本文为微型发电机提供传​​感和计算能力，并利用电力电子接口的灵活性来注入（或退出）来自电网的无功功率。如果控制得当，这些设备可以充当精细分布的无功功率补偿器网络，为配电网提供有价值的辅助服务，并最终防止可再生能源的弃电，实现广泛的电动交通，以及推迟网格加固。

​

2 数学模型 

见第四部分文件夹，有对文章详细讲解。

3 运行结果

​

​ ​

 ​

 ​

 ​

​ 

​ 

 ​

4 Matlab代码+数据+文档讲解 

部分代码：

%% 可视化 
disp('机组可视化')

close all

colorone = [217,95,2]/255;
colortwo = [117,112,179]/255;
colorgrey = 200*[1,1,1]/255;

figure(1)
    hold on
    line([0 t(end)/60/60], [1 1], 'Color', 'black', 'LineStyle', '-')
    line([0 t(end)/60/60], [params.VMAX params.VMAX], 'Color', 'black', 'LineStyle', '--')
    plot(t/60/60,v_nocontrol, 'Color', colorgrey, 'LineWidth', 0.5)
    plot(t/60/60,v_nocontrol(:,c(1)), 'Color', colorone, 'LineWidth', 1)
    plot(t/60/60,v_nocontrol(:,c(2)), 'Color', colortwo, 'LineWidth', 1)
    xlabel('时间[h]')
    ylabel('电压[p.u.]')
    ylim([0.98 1.1])
    xlim([0 12])
    box on
    title('非受控')
    
    
figure(2)
    hold on
    line([0 t(end)/60/60], [params.VMAX params.VMAX], 'Color', 'black', 'LineStyle', '--')
    line([0 t(end)/60/60], [1 1], 'Color', 'black', 'LineStyle', '-')
    plot(t/60/60,v_static, 'Color', colorgrey, 'LineWidth', 0.5)
    plot(t/60/60,v_static(:,c(1)), 'Color', colorone, 'LineWidth', 1)
    plot(t/60/60,v_static(:,c(2)), 'Color', colortwo, 'LineWidth', 1)
    xlabel('时间[h]')
    ylabel('电压[p.u.]')
    ylim([0.98 1.1])
    xlim([0 12])
    box on
    title('完全分散控制（静态）')
    
    
figure(3)
    hold on
    line([0 t(end)/60/60], [params.QMIN(1) params.QMIN(1)], 'Color', colorone, 'LineStyle', '--')
    line([0 t(end)/60/60], [params.QMIN(2) params.QMIN(2)], 'Color', colortwo, 'LineStyle', '--')
    line([0 t(end)/60/60], [0 0], 'Color', 'black', 'LineStyle', '-')
    plot(t/60/60,q_static(:,1), 'Color', colorone, 'LineWidth', 1)
    plot(t/60/60,q_static(:,2), 'Color', colortwo, 'LineWidth', 1)
    xlabel('时间[h]')
    ylabel('无功功率[MVAR]')
    ylim([min(params.QMIN)-0.1 0.1])
    xlim([0 12])
    box on
    title('完全分散控制（静态）')
    
    
figure(4)
    hold on
    line([0 t(end)/60/60], [params.VMAX params.VMAX], 'Color', 'black', 'LineStyle', '--')
    line([0 t(end)/60/60], [1 1], 'Color', 'black', 'LineStyle', '-')
    plot(t/60/60,v_incremental, 'Color', colorgrey, 'LineWidth', 0.5)
    plot(t/60/60,v_incremental(:,c(1)), 'Color', colorone, 'LineWidth', 1)
    plot(t/60/60,v_incremental(:,c(2)), 'Color', colortwo, 'LineWidth', 1)
    xlabel('时间[h]')
    ylabel('电压[p.u.]')
    ylim([0.98 1.1])
    xlim([0 12])
    box on
    title('完全分散控制（增量）')
    
    
figure(5)
    hold on
    line([0 t(end)/60/60], [params.QMIN(1) params.QMIN(1)], 'Color', colorone, 'LineStyle', '--')
    line([0 t(end)/60/60], [params.QMIN(2) params.QMIN(2)], 'Color', colortwo, 'LineStyle', '--')
    line([0 t(end)/60/60], [0 0], 'Color', 'black', 'LineStyle', '-')
    plot(t/60/60,q_incremental(:,1), 'Color', colorone, 'LineWidth', 1)
    plot(t/60/60,q_incremental(:,2), 'Color', colortwo, 'LineWidth', 1)
    xlabel('时间[h]')
    ylabel('无功功率[MVAR]')
    ylim([min(params.QMIN)-0.1 0.1])
    xlim([0 12])
    box on
    title('完全分散控制(增量)')
    
figure(6)
    hold on
    line([0 t(end)/60/60], [params.VMAX params.VMAX], 'Color', 'black', 'LineStyle', '--')
    line([0 t(end)/60/60], [1 1], 'Color', 'black', 'LineStyle', '-')
    plot(t/60/60,v_networked, 'Color', colorgrey, 'LineWidth', 0.5)
    plot(t/60/60,v_networked(:,c(1)), 'Color', colorone, 'LineWidth', 1)
    plot(t/60/60,v_networked(:,c(2)), 'Color', colortwo, 'LineWidth', 1)
    xlabel('时间[h]')
    ylabel('电压[p.u.]')
    ylim([0.98 1.1])
    xlim([0 12])
    box on
    title('网络控制（DVS）')
    
figure(7)
    hold on
    line([0 t(end)/60/60], [params.QMIN(1) params.QMIN(1)], 'Color', colorone, 'LineStyle', '--')
    line([0 t(end)/60/60], [params.QMIN(2) params.QMIN(2)], 'Color', colortwo, 'LineStyle', '--')
    line([0 t(end)/60/60], [0 0], 'Color', 'black', 'LineStyle', '-')
    plot(t/60/60,q_networked(:,1), 'Color', colorone, 'LineWidth', 1)
    plot(t/60/60,q_networked(:,2), 'Color', colortwo, 'LineWidth', 1)
    xlabel('时间[h]')
    ylabel('无功功率 [MVAR]')
    ylim([min(params.QMIN)-0.1 0.1])
    xlim([0 12])
    box on
    title('网络控制（DVS）')
    
figure(8)
    hold on
    line([0 t(end)/60/60], [0 0], 'Color', 'black', 'LineStyle', '-')
    plot(t/60/60,pvproduction(:,1), 'Color', colorone, 'LineWidth', 1)
    plot(t/60/60,pvproduction(:,2), 'Color', colortwo, 'LineWidth', 1)
    xlabel('时间[h]')
    ylabel('有功功率[MW]')
    xlim([0 12])
    box on
    title('发电')
    
    
figure(9)
    hold on
    line([0 t(end)/60/60], [0 0], 'Color', 'black', 'LineStyle', '-')
    plot(t/60/60,p, 'Color', colorgrey)
    plot(t/60/60,p(:,15), 'Color', 'black', 'LineWidth', 1)
    xlabel('时间[h]')
    ylabel('-负荷[MW]')
    xlim([0 12])
    box on
    title('功率需求')