matlab有限元解二维抛物方程.docxVIP

%%%%% 真解 u=sin(pi*x)*sin(pi*y)*sin(t) %%%%% 方程 diff(u,t)-Laplace(u)=f %%%%% f=si n(pi*x)*si n(pi*y)*cos(t)+2*p^2*si n(pi*x)*si n( pi*y)*si n(t) %clear all % clc %%%%fi nite eleme nt code for parabolic equatio n with con sta nt coefficie nt %%%mesh%% no de=[0,0;1,0;1,1;0,1]; elem=[2,3,1;4,1,3]; T=1; bdEdge=setboundary(node,elem, Dirichlet ; n=in put( Please in put in itial mesh:); M=input( M=; for i=1: n [no de,elem,bdEdge]=u niformrefi ne( no de,elem,bdEdge); end N=size( no de,1); NT=size(elem,1); S=1/NT; r=1/M; A=zeros(N,N); u=zeros(N,M+1); u1=zeros(N,1); f=inline( sin(pi*xx(1,1))*sin(pi*xx(