基于蒙特卡洛生成电动汽车充电负荷曲线程序

clc; clear all; %本程序基于蒙特卡洛思想生成10000台充电汽车充电负荷曲线 %具体题目要求及数据来源为第十一届电工杯 %具体疑问qq:1427167261(有偿） rate=[0.1 0.4 0.5]; %三种充电汽车比例 P_h=[1.4 1.9; 7.7 25.6;40 100];%三种充电功率取值范围 M=10000; %10000辆车 r=0; Bh=zeros(M,96); for i=1:M Ph=zeros(96,1); r=r+1; % for j=1:96 T_star=get_star(); %充电开始时间 qua=get_qua(); %充电电量 P_hi=get_i(rate,P_h);%充电功率 T_char=round(4*qua/P_hi); %一小时有4个15分钟 T_sum=T_star+T_char; %充电结束时间 if T_sum>96 %时间不能大于96 new_T_sum=T_sum-96; Ph(T_star:96)=P_hi; Ph(1:new_T_sum)=P_hi; else Ph(T_star:T_sum)=P_hi; end % end for i=1:96 Bh(r,i)=Ph(i); end end B=sum(Bh,1); %%%%%寻找最大功率值，最大同时充电数量%%%% max_P=max(B); %最大功率 Max_moun=zeros(96,1); %每时刻最大充电数量 for i=1:96 for j=1:M if Bh(j,i)~=0 Max_moun(i)=Max_moun(i)+1;%某时刻最大同时充电数量 end end end max_moun=max(Max_moun); %最大同时充电数量 plot(B,'k'); % legend('柱状图','拟合曲线','1','1'); axis([0 96 0 2*10^4]); xlabel('时间/(15min)'); ylabel('充电功率/kW'); title('充电负荷曲线'); set(gcf,'color',[1 1 1]); set(gca,'xtick',0:4:96);