濟(jì)源建設(shè)工程管理處網(wǎng)站網(wǎng)絡(luò)推廣員要怎么做
目錄
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á)!