不通风情况下的日光温室温度动态模型，输入参数温室结构、位置、墙体厚度和导热系数、温室内外温湿度、光照强度、土壤温湿度、植物蒸腾、蒸发散热、大气透明度、风速、风向等参数，利用控制模型三大守恒定律、三大流体力学基本方程实现，输出温度随时间变化曲线图和温度分布图，MATLAB完整代码案例实现

以下是一个简单的示例代码，供参考：

% 温室参数
L = 10; % 温室长度
W = 6; % 温室宽度
H = 3; % 温室高度
d = 0.05; % 墙体厚度
k = 0.5; % 墙体导热系数
A = 2*L*H + 2*W*H + W*L; % 温室外表面积
V = L*W*H; % 温室体积

% 初始条件
Tin = 25; % 室内温度
Tout = 20; % 室外温度
q = 100; % 太阳辐射热通量
Q = 0; % 热量损失
dt = 60; % 时间步长
tmax = 86400; % 模拟时间
t = 0:dt:tmax; % 时间数组
N = length(t); % 时间步数

% 空气参数
rho = 1.2; % 空气密度
Cp = 1005; % 空气比热容

% 循环模拟
for i = 2:N
    % 太阳辐射热通量
    Qs = q*(1-exp(-0.0001*t(i)));
    
    % 热量平衡方程
    Q = A*k*(Tin-Tout)/d + Qs;
    
    % 温室内空气温度变化
    Tin = Tin + (Q/(rho*Cp*V))*dt;
end

% 绘制温度随时间变化曲线图
plot(t, Tin*ones(size(t)));
xlabel('时间');
ylabel('温度');
title('温度随时间变化曲线图');

% 绘制温度分布图
[X, Y, Z] = meshgrid(1:L, 1:W, 1:H);
T = Tin*ones(L, W, H);
slice(X, Y, Z, T, L/2, W/2, H/2);
xlabel('X');
ylabel('Y');
zlabel('Z');
title('温度分布图');

请注意，这只是一个示例，并且省略了许多重要的因素和方程。如果您需要更详细的模型和代码实现，请咨询相关领域的专业人士。