编写MATLAB代码，计算二维平面热传导，要求使用绝热壁面条件

以下是一个基本的二维平面热传导的MATLAB代码，使用有限差分法：

% 定义参数
nx = 100; % 网格点数（x方向）
ny = 100; % 网格点数（y方向）
dx = 0.01; % 网格间距（x方向）
dy = 0.01; % 网格间距（y方向）
dt = 0.01; % 时间步长
alpha = 0.01; % 热扩散系数
T_ambient = 25; % 环境温度
T_left = 100; % 左侧温度
T_right = 0; % 右侧温度
T_top = 50; % 顶部温度
T_bottom = 75; % 底部温度

% 初始化温度矩阵
T = ones(nx,ny)*T_ambient;

% 设置边界条件
T(:,1) = T_left; % 左侧
T(:,end) = T_right; % 右侧
T(1,:) = T_top; % 顶部
T(end,:) = T_bottom; % 底部

% 进行时间演化
for n = 1:1000
    T_old = T;
    for i = 2:nx-1
        for j = 2:ny-1
            T(i,j) = T_old(i,j) + alpha*dt*(T_old(i+1,j)-2*T_old(i,j)+T_old(i-1,j))/dx^2 ...
                + alpha*dt*(T_old(i,j+1)-2*T_old(i,j)+T_old(i,j-1))/dy^2;
        end
    end
    % 应用绝热壁面条件
    T(:,1) = T_old(:,1);
    T(:,end) = T_old(:,end);
    T(1,:) = T_old(1,:);
    T(end,:) = T_old(end,:);
end

% 绘制结果
[X,Y] = meshgrid(1:ny,1:nx);
surf(X,Y,T);
xlabel('y');
ylabel('x');
zlabel('Temperature');

这个代码使用有限差分法对热传导方程进行离散化，并通过循环进行时间演化。在每个时间步长中，我们根据边界条件和绝热壁面条件更新温度矩阵。最终，我们使用 `surf` 函数绘制了温度分布。