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

小弟初学FDTD   现在想用Matlab 编写一个一维热传导方程的程序，哪位高手能帮忙一下。。。谢谢 *-+C<2"  
Y{`3`Pg&N  
du/dt=d^2/dx^2; s-fKh`  
u(0,t)=0;  u(1,t)=0; wbo{JQ  
u(x,0)=  x  for   0<=x<=0.5;  u(x,0)=1-x,    for  0.5=<x<=1; |YJ$c@  
0,+EV,  
对众多高手可能很简单， .. 1fmSk$ y.9  
uT#Acg  

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

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

简单点的吧。 1A)wbH)  
d(u)/dx = u(n) - u(n-1) / dx 8`rAE_n`%  
d^2(u) / dx^2 = (u(n) + u(n-2) - 2 * u(n-1)) / dx^2 lyOrM7Gs  
所以du/dt = (u(n)|t+1 - u(n)|t ) / dt {
3F}Slb  
n41\y:CAo  
所以递推公式： u(n)  = u(n) + dt / dx^2 * (u(n)  + u(n-2) - 2*u(n-1)) K\Y6 cj  
^)%wq@Hi  
如果是编程的话，大概就是这样： !)51v {  
Ue>A  
int numGrid = 100; )P(d66yq'u  
int xStart = 0; RJ`/qXL  
int xStop = 1; LNyL>VHkK  
x[1:numGrid] =[0:numGrid-1] / numGrid * (xStop - xStart) ; 4/mj"PBKL  
/*** 初始化U***/ q)z1</B-  
***** am$-sh72  
int dt = 0.01; 7Da^Jv k  
int timeCout = 10000; gl(6m`a>  
int curtime = 0; U(Hq4D  
.x-Z+Rs{
g  
for curtime = 0 to timeCount 4%c7#AX[T  
... RL]lt0O{  
end ?SsRN jeL  
oN1wrf}Sh  
肚子饿了，写不动了，很简单的。

版主太强悍了， )+{omQ7v  
谢谢了