Matlab|【需求响应】空调负荷需求响应模型

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，即空调停止工作elses(1,j)=1;                        %忽视这一行即可endTr(1,j)=Tint(j)+(-Ke(j)*Pr*s(1,j)./lint-Hr(j)*(Tint(j)-Tw(j)))/Cr;        %室内温度计算，即论文中的式(2.4)，Hr为室内-外部环境热能交互的等效热导，Cr为室内的等效热质量
end

3程序结果