[求助]一维热传导方程Matlab求助

小弟初学FDTD   现在想用Matlab 编写一个一维热传导方程的程序，哪位高手能帮忙一下。。。谢谢 )wNP( @$L  
9V5-
%Iv  
du/dt=d^2/dx^2; ooQQ-?"m  
u(0,t)=0;  u(1,t)=0; .?l\g-;=  
u(x,0)=  x  for   0<=x<=0.5;  u(x,0)=1-x,    for  0.5=<x<=1;
xt*u4%  
4R\Hpt  
对众多高手可能很简单， .. YY!Rz[/  
x7$}8LZ"B  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

我也是新手，好好的学习。。。

简单点的吧。 S;286[oq@  
d(u)/dx = u(n) - u(n-1) / dx .E8_Oz  
d^2(u) / dx^2 = (u(n) + u(n-2) - 2 * u(n-1)) / dx^2 4?q
<e*W  
所以du/dt = (u(n)|t+1 - u(n)|t ) / dt I!Z_[M  
IvtJ0  
所以递推公式： u(n)  = u(n) + dt / dx^2 * (u(n)  + u(n-2) - 2*u(n-1)) -^hWM}F  
!`Hd-&}bYz  
如果是编程的话，大概就是这样： 7%|~>   
vkEiOFU!u  
int numGrid = 100; P 'od`  
int xStart = 0; }%{LJ}\Px  
int xStop = 1; ^Xq 6:
 
x[1:numGrid] =[0:numGrid-1] / numGrid * (xStop - xStart) ; &keR~~/  
/*** 初始化U***/ hRD=Y<>A  
***** j~jV'f.:H  
int dt = 0.01; [R^iF  
int timeCout = 10000; `3WFjU5a  
int curtime = 0; (Hb:?(  
EWkLXU6t  
for curtime = 0 to timeCount \=ux atw  
... SPY|K  
end sl`s_$J  
3@$,s
~+ 3  
肚子饿了，写不动了，很简单的。

版主太强悍了， ' V^6XI  
谢谢了