Matlab代码书写

% 设置参数
L = 5; % 室内长度（m）
H = 3; % 室内高度（m）
W = 4; % 室内宽度（m）
V = L*H*W; % 室内体积（m^3）
dx = 0.2; % 离散化步长（m）
dy = 0.2;
dz = 0.2;
dt = 1; % 时间步长（s）
total_time = 3600*8; % 总时间（s）
alpha = 1.8e-7; % 热扩散系数（m^2/s）

% 初始化温度场
T = ones(L/dx, H/dy, W/dz) * 293; % 初始温度为293K

% 初始化边界条件
T(:, :, 1) = 273; % 左边界
T(:, :, end) = 293; % 右边界
T(:, 1, :) = 283; % 下边界
T(:, end, :) = 283; % 上边界
T(1, :, :) = 283; % 前边界
T(end, :, :) = 283; % 后边界

% 迭代求解
for t = dt:dt:total_time
    % Block模型部分
    T_mean = mean(T(:));
    T = T + alpha * dt / V * (T_mean - T) * ones(size(T)); % Block模型

    % Gao-Yong湍流方程部分
    [Tx, Ty, Tz] = gradient(T, dx, dy, dz); % 计算梯度
    Txx = diff(Tx, 2, 1) ./ dx^2; % 计算二阶导数
    Tyy = diff(Ty, 2, 2) ./ dy^2;
    Tzz = diff(Tz, 2, 3) ./ dz^2;
    T = T(2:end-1, 2:end-1, 2:end-1) + alpha * dt* (Txx + Tyy + Tzz); % 更新温度场
end

% 可视化结果
[x, y, z] = meshgrid(dx:dx:L, dy:dy:H, dz:dz:W);
slice(x, y, z, T, L/2, H/2, W/2);
xlabel('X (m)');
ylabel('Y (m)');
zlabel('Z (m)');
colorbar;

这个代码Matlab代码有语法错误吗