matlab求解热传导实例

?利 片 di produces the solution to the FEM hrmnuiaiiori oTt^e scalar PDE problemР- V (^Vwfc + ow - f on 口Р別！Рor rhe system PDE problemРQ甘一 f (c ® ^iO+ fiff - f " QР类比可得系数c=1,a=0,f=0,d=1.计算参考程序如下：Р[p,e,t]=i nitmesh('squareg');Р[p,e,t]=refi nemesh('squareg',p,e,t);Рu0=zeros(size(p,2),1);Рix=fi nd(sqrt(p(1,:).A2+p(2,:).A2)<0.8);Рu0(ix)=on es(size(ix));Рtlist=li nspace(0,0.1,20);Рu1=parabolic(u0,tlist,'squareb1',p,e,t,1,0,0,1);Рpdeplot(p,e,t,'xydata',u1(:,8),'mesh','off,'colormap','hot');Рx=li nspace(_1,1,31);y=x;Р[unu sed,t n, a2,a3]=tri2grid(p,t,u0,x,y);Р%制作动画Рn ewplot;Рumax=max(max(u1));Рumi n=min( mi n(u1));Рfor j=1:8Рu=tri2grid(p,t,u1(:,j),t n, a2,a3);Рi=fi nd(is nan( u));Рu(i)=zeros(size(i));Рsurf(x,y,u);caxis([u min umax]);colormap(hot)Рaxis([-1 1 -1 1 0 1]);Рm(j)= getframe;РendРmovie(m);