配电网电压调节及通信联系(Matlab代码实现)

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

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

                           ​

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

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

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

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

目录

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('功率需求')
   
  • 1
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

荔枝科研社

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值