蒙特卡罗论坛：基于蒙特卡洛的电动车有序充放电（Matlab代码实现）

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

📝目前更新：🌟🌟🌟电力系统相关知识，期刊论文，算法，机器学习和人工智能学习。
🚀支持：🎁🎁🎁如果觉得博主的文章还不错或者您用得到的话，可以关注一下博主，如果三连收藏支持就更好啦！这就是给予我最大的支持！ 

 

 

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

​

目录

1 概述

2 蒙特卡洛模拟方法介绍

3 基于蒙特卡洛法的规模化电动汽车充电负荷预测（Python&Matlab实现）

4 （Matlab实现）基于蒙特卡诺和拉格朗日乘子法的电动车调度【有序、无序充放电】

5 规模化电动汽车充电负荷预测计算方法

6 仿真结果

7 兴趣

8 Matlab代码实现

---

​
式中SOC1为完成充电时的电池荷电状态,SOC2为上一次完成充电时的电池荷电状态。

计算电动汽车充电所需时间T,可以通过将电池容量C、起始荷电状态SOC和充电功率Р求
得，具体公式如下:

​
(2)计算某一种用途电动汽车在第i个充电负荷计算点时的总充电负荷。本文将每天计算为1440分钟，每15分钟计算一次充电负荷，共计96个充电负荷计算点，计算预测北京市各用途的电动汽车充电负荷，得到各用途电动汽车充电负荷。某一种用途电动汽车的充电负荷的预测计算方法如下:

​

其中

​
 

​
 

式中，Nc，Nt，Nb分别表示在i时刻充电的电动私家车、电动出租车、电动公交车的数量;

部分代码：

Load=xlsread('数据.xls')for i=1:1:5000 S=20+10*rand(1);%容量 P=rand(1)+2;%充电功率 finish(i,1)=17.6+3.4*randn(1);%最后出行结束时刻 if finish(i,1)>24 finish(i,1)=finish(i,1)-24; end while(finish(i,1)<0) finish(i,1)=17.6+3.4*randn(1); end finish(i,1)=ceil(finish(i,1)); distance(i,1)=lognrnd(3.2,0.88); while(distance(i,1)>160/1.61) distance(i,1)=lognrnd(3.2,0.88);%日行驶里程 end %% 放电 temp=9*rand(1); if (temp<=6) disch(i,1)=10+temp; else disch(i,1)=20-6+temp; end

Load=xlsread('数据.xls')
for i=1:1:5000
    S=20+10*rand(1);%容量
    P=rand(1)+2;%充电功率
    finish(i,1)=17.6+3.4*randn(1);%最后出行结束时刻
    if finish(i,1)>24
        finish(i,1)=finish(i,1)-24;
    end
    while(finish(i,1)<0)
        finish(i,1)=17.6+3.4*randn(1);
    end
    finish(i,1)=ceil(finish(i,1));
    distance(i,1)=lognrnd(3.2,0.88);
    while(distance(i,1)>160/1.61)
        distance(i,1)=lognrnd(3.2,0.88);%日行驶里程
    end
    %% 放电
    temp=9*rand(1);
    if (temp<=6)
        disch(i,1)=10+temp;
    else
        disch(i,1)=20-6+temp;
    end

Load=xlsread('数据.xls')
for i=1:1:5000
    S=20+10*rand(1);%容量
    P=rand(1)+2;%充电功率
    finish(i,1)=17.6+3.4*randn(1);%最后出行结束时刻
    if finish(i,1)>24
        finish(i,1)=finish(i,1)-24;
    end
    while(finish(i,1)<0)
        finish(i,1)=17.6+3.4*randn(1);
    end
    finish(i,1)=ceil(finish(i,1));
    distance(i,1)=lognrnd(3.2,0.88);
    while(distance(i,1)>160/1.61)
        distance(i,1)=lognrnd(3.2,0.88);%日行驶里程
    end
    %% 放电
    temp=9*rand(1);
    if (temp<=6)
        disch(i,1)=10+temp;
    else
        disch(i,1)=20-6+temp;
    end 

7 兴趣

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

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

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

8 Matlab代码实现