目錄

1?主要內(nèi)容

2?程序難點(diǎn)及問(wèn)題說(shuō)明

3?部分程序

4 下載鏈接

---

1?主要內(nèi)容

該程序分為兩部分，日前優(yōu)化部分——該程序首先根據(jù)《電力系統(tǒng)云儲(chǔ)能研究框架與基礎(chǔ)模型》上面方法，根據(jù)每個(gè)居民的實(shí)際需要得到響應(yīng)儲(chǔ)能充放電功率，優(yōu)化得到整體的儲(chǔ)能充放電功率情況。日內(nèi)滾動(dòng)mpc跟蹤部分——采用《基于MPC的微電網(wǎng)并網(wǎng)優(yōu)化調(diào)度》P31-36頁(yè)相關(guān)內(nèi)容，通過(guò)預(yù)測(cè)模型、滾動(dòng)優(yōu)化、反饋校正得到soc跟蹤情況。

2?程序難點(diǎn)及問(wèn)題說(shuō)明

• 代碼修正

for i=1:96    C=[C;P_C1(i) == P_G1(i)+P_PV1(i);%儲(chǔ)能功率=電網(wǎng)供電+光伏儲(chǔ)能P_PV1(i)<=P_C1(i)<=P_cap;%儲(chǔ)能值約束0<=P_G1(i)<=P_cap;%電網(wǎng)供電約束E1(1)+sum(P_C1)*nc*t==E_cap;%soc約束，到儲(chǔ)能結(jié)束soc=1];
end

該程序含有這段代碼，最后一句代碼和i不相關(guān)，卻重復(fù)循環(huán)96次，這句代碼和下面的代碼功能重復(fù)，在學(xué)習(xí)和修改代碼的時(shí)候增加了復(fù)雜度。此外，還有幾處類似問(wèn)題已經(jīng)修正。

原文獻(xiàn)有儲(chǔ)能充放電功率模型如下，該部分約束能夠保證充放電功率不能同時(shí)不為0，但是在代碼里面卻存在這方面問(wèn)題，已經(jīng)進(jìn)行了修正。

• 部分問(wèn)題

代碼原來(lái)的出圖效果如下：

該部分是第一階段優(yōu)化結(jié)果，對(duì)于云儲(chǔ)能用戶，考慮最簡(jiǎn)單的調(diào)度策略，即根據(jù)電價(jià)詳細(xì)分清楚每個(gè)時(shí)間段的充放電狀態(tài)，所以soc曲線也較為規(guī)律?（該圖應(yīng)該為5段，即充-放-充-放-充，但是考慮到前四段已經(jīng)能夠說(shuō)明問(wèn)題，就只分析前四段）。

得到兩個(gè)階段的soc跟蹤效果圖，但是為什么只是充電階段的跟蹤圖呢?？為什么第一個(gè)圖效果那么好呢?？第一個(gè)問(wèn)題留給大家思考，第二個(gè)問(wèn)題一會(huì)揭曉?。

• 修正后的代碼運(yùn)行結(jié)果

修正代碼實(shí)現(xiàn)了全過(guò)程的soc跟蹤控制。

本模型的跟蹤控制主要是由于光伏波動(dòng)引起的，在不存在光伏波動(dòng)的時(shí)間段就會(huì)出現(xiàn)soc完美的跟蹤效果，也就回答了上面第二個(gè)問(wèn)題?。

3?部分程序

clc;clear all;
%% 數(shù)據(jù)準(zhǔn)備
data = xlsread('預(yù)測(cè)數(shù)據(jù).xlsx');
?
p_pv = data;%正值代表負(fù)荷。負(fù)值代表用戶光伏發(fā)電量高于負(fù)荷的部分
p_pv(data<0) = 0;    %用戶的剩余光伏數(shù)據(jù)
prepv=sum(p_pv,2);
data_d = data;
data_d(data>0) = 0;   %用戶放電負(fù)荷數(shù)據(jù)
data_d = -1.*data_d;
preload=sum(data_d,2);
%% 參數(shù)設(shè)置
e_cap = xlsread('各用戶租賃容量.xlsx');  %各用戶租賃的能量容量
p_cap = 0.5.*e_cap;                                %各用戶租賃的功率容量
E_cap = 175;                                       %CES儲(chǔ)能容量
P_cap = 90;                                        %CES功率容量
?
e0 = 0.2.*e_cap;                                   %各用戶初始電量
E0 = 0.2*E_cap;                                    %CES初始電量
emin = 0.1.*e_cap; 
Emin = 0.1*E_cap;
?
nc = 0.96;  %充電效率
nd = 0.96;  %放電效率
?
r1 = 0.32;   %低谷電價(jià)
r2 = 0.66;   %平谷電價(jià)
r3 = 1.1;    %高峰電價(jià)
r = zeros(288,1);  %各時(shí)段電價(jià)
r(1:96) = r1;
r(97:144) = r3;
r(145:204) = r2;
r(205:252) = r3;
r(253:288) = r2;
?
t = 5/60;  %時(shí)間間隔
%% 定義變量
E1 = sdpvar(96,1);%儲(chǔ)能容量變量
E1(1) = E0;
P_PV1 = sum(p_pv(1:96,:),2);%充電時(shí)段0-8：00整體光伏出力
P_C1 = sdpvar(96,1);%儲(chǔ)能出力變量
P_G1 = sdpvar(96,1);%電網(wǎng)供電變量
C=[];%初始化約束
for i=1:96    C=[C;P_C1(i) == P_G1(i)+P_PV1(i);%儲(chǔ)能功率=電網(wǎng)供電+光伏儲(chǔ)能P_PV1(i)<=P_C1(i)<=P_cap;%儲(chǔ)能值約束0<=P_G1(i)<=P_cap;%電網(wǎng)供電約束
%        E1(1)+sum(P_C1)*nc*t==E_cap;%soc約束，到儲(chǔ)能結(jié)束soc=1];
end

4 下載鏈接

點(diǎn)擊直達(dá)！