椭圆方程的有限元法(共6页).doc

精选优质文档-倾情为你奉上椭圆形方程的有限元 上机实习报告专心-专注-专业两点边值问题有限元法（必做）从Galerkin原理出发用线性元解两点边值问题：精确解：。1.1变分形式从Galerkin原理出发推导出两点边值问题的变分形式，将积分区间等分为N份，则步长，记为。写出有限元方程及系数矩阵元素。解： 由题可知p(x)=1,q(x)=1,f(x)= 所以有限元方程为,j=1,2,.,n 其中， ， ， ， 计算有： 1.2利用MATLAB求解问题的过程依次取用MATLAB求解并图形比较数值解与精确解，用表格列出不同剖分时的误差。程序如下：function u= bianzhi(p,q,N)h=1/N;x=0:h:1;A=zeros(N-1);for i=2:N-1 a3=(t)-p./h+h.*q.*t.*(1-t);