matlab_pde.pptVIP

* 有限元方法——介绍 有限元方法是数值求解偏微分方程边值问题的一种方法，此方法首先于20世纪50年代初由工程师提出，并用于求解简单的结构问题。有限元方法是这一种系统的数值方法，并奠定其数学基础，是在60年代中期以冯康先生为代表的中国学者与西方学者独立并行完成的。 有限元方法不同于差分方法，主要有以下三大特点： (1)从数学物理问题的变分原理出发，而不是从微分方程出发，因此事从问题的整体描述而不是从问题的局部描述出发。 (2)对所考虑问题的区域(以二维情形为例)作三角形（或其他简单多边形）剖分，而不是仅作矩形剖分。 (3)用剖分区域上的简单函数(如分片多项式)去逼近原问题的解，而不是只在剖分节点上的数值逼近。 * 有限元方法——一维边值问题算例 用有限元方法求解如下一维边值问题： 解：一维边值问题线性有限元数值解法的MATLAB程序如下： %线性有限元方法 %参数设置 N=10; X=0:1/(N+1):1; b=zeros(N+1,1); A=zeros(N+1,N+1); for i=2:N+1 F1=@(x)2*(N+1)*(x-X(i-1)).*sin(pi*x/2);%句柄函数 R1=quad(F1,X(i-1),X(i)); F2=@(x)2*(N+1)*(X(i+1)-x).*sin(pi*x/2); R2=quad(