能源调度

%负荷的置信区间求解 sigma=xlsread('D:/gesktop/P4&P5/00融合/论文放图/风光相关性/风光等不确定性结果.xlsx','负荷标准差','B2:Y5'); % bata=[0.05,0.15,0.25,0.95,0.85,0.75]; % a=zeros(4,24); % for i=1:4 % for j=1:24 % a(i,j)=icdf('norm',bata(1,6),0,sigma(i,j)); % end % end a=[6.579414508 -3.289707254 6.579414508 4.934560881 -3.289707254 -3.289707254 8.224268135 3.289707254 6.579414508 3.289707254 1.644853627 8.224268135 3.289707254 4.934560881 6.579414508 1.644853627 1.644853627 1.644853627 4.934560881 8.224268135 1.164735875 8.224268135 -0.822659486 8.224268135 3.289707254 -8.224268135 4.934560881 3.289707254 -1.644853627 8.224268135 8.224268135 4.934560881 -6.579414508 -3.289707254 3.289707254 6.579414508 -4.934560881 -0.800723545 -3.289707254 -8.224268135 -3.289707254 -6.579414508 6.579414508 -3.289707254 6.579414508 -6.579414508 -4.934560881 -1.644853627 0 8.224268135 0 0 -3.289707254 -4.934560881 -4.934560881 -4.934560881 -3.289707254 4.934560881 8.224268135 -4.934560881 -8.224268135 4.934560881 3.289707254 -8 -3.289707254 -1.644853627 -6.579414508 -6.579414508 -4.934560881 5.934532649 -1.644853627 4.934560881 0 0 0 0 0 0 0 0 0 0 0 -4.934560881 0 0 0 0 0 -6.579414508 -8.224268135 0 0 0 0 0 ]; b=zeros(4,24); for i=1:4 for j=1:24 %b(i,j)=abs(normpdf(a(i,j),0,sigma(i,j))-normpdf(0,0,sigma(i,j))); b(i,j)=abs(normcdf(a(i,j),0,sigma(i,j))-normcdf(0,0,sigma(i,j))); end end %bata % xlswrite('D:/gesktop/P4&P5/00融合/论文放图/风光相关性/风光等不确定性结果.xlsx',a,'负荷0.25','B2:Y5'); % xlswrite('D:/gesktop/P4&P5/00融合/论文放图/风光相关性/风光等不确定性结果.xlsx',b,'0.25','B4:Y7'); %1-bata % xlswrite('D:/gesktop/P4&P5/00融合/论文放图/风光相关性/风光等不确定性结果.xlsx',a,'负荷0.25','B12:Y15'); % xlswrite('D:/gesktop/P4&P5/00融合/论文放图/风光相关性/风光等不确定性结果.xlsx',b,'0.25','B14:Y17'); gfj=0 %% % x=-6:0.2:6; % y1=normpdf(x,0,1); % y2=normpdf(x,0,3); % y3=normpdf(x,0,5); % y4=normpdf(x,0,7); % y5=normpdf(x,0,9); % plot(x,y1,'g'); % hold on % plot(x,y2,'b'); % hold on % plot(x,y3,'r'); % hold on % plot(x,y4,'y'); % hold on % plot(x,y5,'m');