求解热传导方程matlab

这是非稳态一维热传导的方法，也叫古典显格式。

如果是做数学建模，就别用了，这种方法计算量比较大，算的很慢，而且收敛不好。

但是如果实在没办法也能凑合用。

该改的地方我都用？？？代替了。

给个详细解释https://wenku.baidu.com/view/78a359d43b3567ec112d8a77.html?qq-pf-to=pcqq.group

 function rechuandao() %

 Llist = ？？？
 N = ; % 空间点数
 M = ;
 alfa = ？？？ % 导热 / (密度*比热)
 lambda = 0.4; % 稳定条件
 %************
 zone=Llist()/N;%空间步长
 z=:zone:Llist();%空间点坐标数组
 z=z';%矩阵转置

 n = lambda * (Llist()/N)^ / alfa();

 TM = M * n;   % 总时间

 t=:n:TM;    %时间点坐标数组
 t=t';        %矩阵转置

 %计算初值和边值
 T=zeros(N+,M+);  %赋温度矩阵初值
 Ti=init_fun(z);    %各空间点处的初始条件
 To=border_funo(t); %内边界条件
 Te=border_fune(t); %外边界条件
 T(:,)=Ti;         %(初始条件 温度都为37)
 T(,:)=To;         %（边界条件x=0处温度为）
 T(N+,:)=Te;       %（边界条件x=L 处温度为）

 %用差分法求出温度T与杆长L、时间t的关系
 for k=:M  %时间点数
     n=;
     while n<=N %空间点数
     T(n,k+)=lambda*(T(n+,k)+T(n-,k))+(-*lambda+)*T(n,k);
     n=n+;
     end
 end
 %设置立体网格
 for i=:M+
     X(:,i)=z;
 end
 for j=:N+
     Y(j,:)=t;
 end
 mesh(X,Y,T);
 view([ - ]);
 xlabel('长度');
 ylabel('时间');
 zlabel('温度');
 csvwrite('csv.csv',T)

 function y=init_fun(z)%初值条件
 y=？？？
 return

 function y=border_funo(t)%z=0的边界条件
 y=？？？
 return

 function y=border_fune(t)%z=L的边界条件
 y= ？？？
 return
 %