1 主要内容
程序主要复现《溫控负荷的需求响应潜力评估及其协同优化管理研究_谢敦见》2.5部分章节的内容,建立空调负荷的聚合模型,考虑调节空调温度对空调响应潜力的影响,程序结果充分说明随着上调温度增大,响应程度逐渐增大。
具体模型如下:
2 部分代码
代码句句注释,非常方便学习!
%% 设置基础参数
N=1000; %设置空调数量
lint=60; %1小时=60分钟
Hr=normrnd(75.52,50,[1,N]); %设置正态分布的室内-外部环境热能交互的等效热导,生成服从正态分布(均值75.52,标准差50)的1行N列随机数矩阵,单位为m^2
Hr=abs(Hr); %abs()数字取绝对值
Cr=288; %房间热容参数设置为288,单位为KJ/m^2℃
Ke=3+0.6.*rand(1,N); %能效比,设置为(3-3.6)之间的均匀分布随机数,rand(1,N)为随机1行N列的数发生器
Pr=2000; %每台空调额定功率设为2,单位为kW
Te=21+7.*rand(1,N); %用户设定温度的参数为(21-28)之间的均匀分布随机数,单位为℃
tt=0.2; %温度调整间隔为0.2,单位为℃
Te=Te+tt; %用户设定的温度+调节的温度
Tint=24.*ones(1,N); %初始的室内温度为24,单位为℃
Tw=30.*ones(1,N); %外界温度的变化率
deta=3; %设置data为温控负荷的温度迟滞区间为3
td=60; %
for j=1:N %循环N次,因为有N个空调
if Tint(j)>=Te(j)+deta %当室内温度>=用户设定的温度+温控负荷的温度迟滞区间
s(1,j)=1; %s(1,j)为温控负荷的实时状态为1,即空调工作
elseif Tint(j)<=Te(j)-deta %当室内温度<=用户设定的温度-温控负荷的温度迟滞区间
s(1,j)=0; %s(1,j)为温控负荷的实时状态为0,即空调停止工作
else
s(1,j)=1; %忽视这一行即可
end
Tr(1,j)=Tint(j)+(-Ke(j)*Pr*s(1,j)./lint-Hr(j)*(Tint(j)-Tw(j)))/Cr; %室内温度计算,即论文中的式(2.4),Hr为室内-外部环境热能交互的等效热导,Cr为室内的等效热质量
end
3 程序结果
4 下载链接