matlab求解三维热传导问题.zip

U0=2; a2=2; N=10; u0=0; theta=linspace(0,2*pi,100); r=eps:0.05:1; t=0.1:0.001:0.2; [RR,TT]=meshgrid(r,t); figure(1) [R,TH]=meshgrid(theta,r); [X,Y]=pol2cart(R,TH); for tt=1:100 un=0; for k=1:N unn =2*(U0-u0)*(-1)^k.*sin(k.*pi.*(X.^2+Y.^2).^0.5).*exp(-k^2*pi^2*a2*t(tt))./(pi *k*((X.^2+Y.^2).^0.5)); un=unn+un; end mesh(X,Y,un); axis([-1 1 -1 1 -0.4 0]) pause(0.1) end figure(2) wn=0; for k=1:N wnn =2*(U0-u0)*(-1)^k.*sin(k.*pi.*RR).*exp(-k^2*pi^2*a2*TT)./(pi *k.*RR); n=wnn+wn; end waterfall(RR,TT,wn) xlabel('r') ylabel('t')