导热方程求解matlab(共3页).doc

使用差分方法求解下面的热传导方程初值条件：边值条件：使用差分公式上面两式带入原热传导方程令，化简上式的编程MATLAB程序，运行结果如下function mypdesolutionc=1;xspan=0 1;tspan=0 0.2;ngrid=100 10;f=(x)4*x-4*x.2;g1=(t)0;g2=(t)0;T,x,t=rechuandao(c,f,g1,g2,xspan,tspan,ngrid);x,t=meshgrid(x,t);mesh(x,t,T);xlabel(x)ylabel(t)zlabel(T)function U,x,t=rechuandao(c,f,g1,g2,xspan,tspan,ngrid)% 热传导方程：% Ut(x,t)=c2*Uxx(x,t) axb tsttf% 初值条件：% u(x,0)=f(x)% 边