temper.m

来自「关于matlab解热传导问题的解法」· M 代码 · 共 30 行

function T=Temper(L,t)
L=1;t=2.5;
NumX=10;
Numt=1000;
t0=0;
dx=L/NumX;
dt=(t-t0)/Numt;
T(Numt+1,NumX+1)=zeros;
%初始条件
for j = 0:1:NumX
    T(1,j+1)=sin(pi*j*dx);
end
%边界条件
T(:,1)=0;
T(:,NumX+1)=0;
%循环计算
for i =1:1:Numt
    for j = 2:1:NumX
    T(i+1,j)=0.25*dt/dx/dx*(T(i,j+1)-2*T(i,j)+T(i,j-1))+T(i,j);
    end
end
%绘图
x=0:dx:L
plot(x,T(1:10:Numt,:),'-');