含分布式电源的配电网可靠性评估(matlab代码)

目录

[1 主要内容](#1 主要内容)

模型节点图

[2 部分代码](#2 部分代码)

[3 程序结果](#3 程序结果)

[4 下载链接](#4 下载链接)


1 主要内容

该程序参考《基于仿射最小路法的含分布式电源配电网可靠性分析》文献方法,通过概率模型和时序模型分别进行建模,实现基于概率模型+最小路法的含分布式电源配电网可靠性评估以及时序模型+序贯蒙特卡洛模拟法的含分布式电源配电网可靠性评估。程序采用matlab编写,注释清楚,方便学习!

模型节点图

采用IEEERBTSBUS6_F4结构图:

2 部分代码

复制代码
% 读取数据
[line_num,line,load_num,load] = IEEE_RBTS_BUS6_F4;%可以自己修改参数
% 生成图
g=graph(line(:,1),line(:,2),line(:,3));
% plot(g)%生成节点图
% DG接入位置,可修改
DG_position=[13,27];
% 初始化负荷点可靠性指标
Lambda_load=zeros(1,load_num);
Gamma_load=zeros(1,load_num);
U_load=zeros(1,load_num);
P_load=zeros(1,load_num);
load_users=zeros(1,load_num);
% DG四状态概率模型,DG出力大小可以自己修改,如果四个状态均为零就相当于不含DG的系统
[P1,S1]=DG_probability_model(0,240,640,800);
[P2,S2]=DG_probability_model(0,240,640,800);
% 各个状态下可形成孤岛供电的负荷及恢复比例
node_island1=cell(2,4);
node_island2=cell(2,4);
% 可恢复负荷的搜索顺序
search_order1=[9 10 8 11 12 7 13 6];
search_order2=[19 18 20 17 16 15 14];
% 确定可以形成孤岛的负荷
for k=1:4
    medalS1=S1(k);
    medalS2=S2(k);
    medal_index1=1;
    while medalS1~=0
        if medalS1>=load(search_order1(medal_index1),4)
            medalS1=medalS1-load(search_order1(medal_index1),4);
            node_island1{1,k}=[node_island1{1,k},search_order1(medal_index1)];
            node_island1{2,k}=[node_island1{2,k},load(search_order1(medal_index1),4)];
        else
            node_island1{1,k}=[node_island1{1,k},search_order1(medal_index1)];
            node_island1{2,k}=[node_island1{2,k},medalS1];
            medalS1=0;
        end
        medal_index1=medal_index1+1;
    end
    medal_index2=1;
    while medalS2~=0
        if medalS2>=load(search_order2(medal_index2),4)
            medalS2=medalS2-load(search_order2(medal_index2),4);
            node_island2{1,k}=[node_island2{1,k},search_order2(medal_index2)];
            node_island2{2,k}=[node_island2{2,k},load(search_order2(medal_index2),4)];
        else
            node_island2{1,k}=[node_island2{1,k},search_order2(medal_index2)];
            node_island2{2,k}=[node_island2{2,k},medalS2];
            medalS2=0;
        end
        medal_index2=medal_index2+1;
    end
end

3 程序结果

概率模型+最小路法得到的结果如下:

时序模型+序贯蒙特卡洛模拟法得到结果如下:

4 下载链接

相关推荐
斯普信专业组2 分钟前
fluent-bit使用kafka作为数据源采集问题
分布式·kafka·fluent-bit
熊猫_豆豆1 小时前
用MATLAB画一只可爱的小熊
前端·matlab·画图
工程师小星星1 小时前
消息队列Apache Kafka教程
分布式·kafka·apache
熊猫_豆豆3 小时前
MATLAB画出湖面波纹相遇所形成的现象
开发语言·matlab·仿真
风跟我说过她3 小时前
Hadoop HA (高可用) 配置与操作指南
大数据·hadoop·分布式·zookeeper·centos
还是大剑师兰特4 小时前
Kafka 面试题及详细答案100道(66-80)-- 运维与部署
分布式·kafka·大剑师·kafka面试题
机器学习之心6 小时前
基于RNN循环神经网络的锂电池剩余寿命预测Matlab实现
rnn·matlab·锂电池剩余寿命预测·rnn循环神经网络
机器学习之心9 小时前
多目标鲸鱼优化算法(NSWOA),含46种测试函数和9个评价指标,MATLAB实现
算法·matlab·多目标鲸鱼优化算法·46种测试函数·9个评价指标
Insist75314 小时前
基于OpenEuler部署kafka消息队列
分布式·docker·kafka
在未来等你15 小时前
Elasticsearch面试精讲 Day 20:集群监控与性能评估
大数据·分布式·elasticsearch·搜索引擎·面试