偏微分方程的有限元法求解(共22页).doc

精选优质文档-倾情为你奉上 %对于d2u/dx2=f的FEM解算器，其中f=x*(1-x)%边界条件u(0)=0, u(1)=0.%精确解用以比对xx=linspace(0,1,101);%产生0-1之间的均分指令，101为元素个数uex=(1/6).*xx.3-(1/12).*xx.4-(1/12).*xx;%对力项设置高斯点的数目NGf=2;if (NGf=2)xiGf=-1/sqrt(3);1/sqrt(3);%1、2的值aGf=1 1;else,NGf=1;xiGf=0.0;aGf=2.0;end%单元数目Ne=5;%建立网格节点x=linspace(0,1,Ne+1);%零刚性矩阵K=zeros(Ne+1,Ne+1);b=zeros(Ne+1,1);%对所有单元循环计算刚性和残差