MATLAB代码："引号内的标题是：“计及碳排放交易与多元需求响应的微网/虚拟电厂
日前优化调度模型”。
最近在搞虚拟电厂优化调度发现个有意思的现象——把空调遥控器揣兜里的大爷大妈和满街跑的电
驴子，竟然能跟碳中和扯上关系。今儿咱们就扒开这个微网调度代码，看看怎么用MATLAB让空调、电动车和
碳指标在电表上打架。
先说这个空调负荷建模，搞过的人都知道温度延迟效应有多头大。咱这代码里用了个状态空间模型
，比常规的等效热参数模型准多了。看这段约束设置：
```matlab
% 室内温度状态方程
for t=2:24
T_in(t) = T_in(t-1) + dt*(alpha*(T_out(t-1)-T_in(t-1)) + beta*P_ac(t-1)/C_room);
end
% 温度舒适度约束
T_set_low <= T_in <= T_set_high;
```
这里的alpha、beta是建筑热惯性参数，dt取0.5小时步长。有意思的是把空调功率P_ac作为控制变
量，既保证室温在舒适区，又能让空调集群当虚拟储能使。
说到电动车调度，代码里把用户出行规律玩明白了。每个电动车对象都带着出发时间、里程焦虑：
```matlab
classdef EV_model
properties
arrive_time; % 到家时间
depart_time; % 出门时间
soc_now; % 当前电量
soc_max; % 电池容量
mileage; % 日行驶里程
end
methods
function obj = calc_charge_curve(obj)
% 根据出行需求反推充电曲线
