请给我写一段matlab代码来求解热传导方程

好的，以下是一个简单的Matlab代码来求解热传导方程：

% 定义参数
alpha = 0.1; % 热扩散系数
L = 1; % 空间长度
T = 1; % 时间长度
dx = 0.1; % 空间步长
dt = 0.01; % 时间步长

% 初始化温度分布
x = 0:dx:L; % 空间网格
u0 = sin(pi*x/L); % 初始温度分布
u = u0;

% 迭代求解
for t = 0:dt:T
    % 使用向前差分逼近时间偏导数
    unew(1) = u(1) + alpha*dt/dx^2 * (u(2) - 2*u(1) + u(1));
    for i = 2:length(x)-1
        unew(i) = u(i) + alpha*dt/dx^2 * (u(i+1) - 2*u(i) + u(i-1));
    end
    % 使用向后差分逼近时间偏导数
    unew(end) = u(end) + alpha*dt/dx^2 * (u(end-1) - 2*u(end) + u(end));
    
    u = unew; % 更新温度分布
end

% 绘制结果
plot(x, u)
xlabel('x')
ylabel('u')
title('Temperature Distribution')

这个代码使用向前差分逼近时间偏导数，并使用了简单的欧拉迭代方法来求解热传导方程。您也可以尝试使用其他差分方法或数值求解器来解决此问题。