h=L/N;%空间步长
z=0:h:L;%空间点坐标数组
z=z';%矩阵转置
tao=lambda*h^2/alfa;%时间步长
tm=M*tao;%热传导的总时间tm
t=0:tao:tm;%时间点坐标数组
t=t';%矩阵转置
%计算初值和边值
T=zeros(N+1,M+1);%赋温度矩阵初值
Ti=init_fun(z);%各空间点处的初始条件
To=border_funo(t);%内边界条件
Te=border_fune(t);%外边界条件
T(:,1)=Ti;
T(1,:)=To;
T(N+1,:)=Te;% 赋予温度矩阵边界条件值与初始条件值
%用差分法求出温度T与杆长L、时间t的关系
for k=1:M
    m=2;
    while m<=N
    T(m,k+1)=lambda*(T(m+1,k)+T(m-1,k))+(-2*lambda+1)*T(m,k);
    m=m+1;
    end;
end;
%设置立体网格
for i=1:M+1
    X(:,i)=z;
end;
for j=1:N+1
    Y(j,:)=t;
end
mesh(X,Y,T);
view([1 -1 1]);
xlabel('Z');
ylabel('t');
zlabel('T');

function y=init_fun(z)%初值条件
y=1-z.^2;
return

function y=border_funo(t)%z=0的边界条件
y=1+t.*0;
return

function y=border_fune(t)%z=L的边界条件
y=t*.0;
return

%

2楼: Originally posted by 11jxli1 at 2012-12-24 10:24:46
你把资料截图上传大家方便浏览
以为热传导问题MAT喇叭有现成程序
你查阅MATLAB使用指南

2楼: Originally posted by 11jxli1 at 2012-12-24 10:24:46
你把资料截图上传大家方便浏览
以为热传导问题MAT喇叭有现成程序
你查阅MATLAB使用指南