微分方程数值方法实验报告.docxVIP

一、实验题目 (1)用Ritz-Galerkin方法求解边值问题u + w (1)用Ritz-Galerkin方法求解边值问题 w(0) = ?(1) = 0 精确解仏=空兰一 ％ sinl (2)用有限元方法求解L+之+小血仏)以5(加) I w(0)=弘⑴=0 实验目的 运用MATLAB数学软件编写Ritz-Galerkin方法和有限元方法程序，进一步熟悉MATLAB 的应用及掌握偏微分方程数值方法中Ritz-Galerkin方法和有限元方法，対各个方法求解精 度进一步明确。 实验原理 (1) Ritz-Galerkin 方法 Ritz方法 设给出二次泛函 J(v) = — A (v,v)-F(v), 2 (1) 其屮人(叫°是Hilbert空间V上的双线性泛函，而且满足对称性、有界性、正定性；F(v) 是V上的有界线性泛函。考虑一下变分问题：V求满足 J(w) = min J(v). ? (2) 设Hilbert空间V是可分的，即V中有可数多个元素构成一个稠密集。在V屮取N个线 性无关的元素5?，,他们张成一个N维子空间* ,即 N S/v ={Wu =》CS；C],???,Cn为实数} 台 ，或记成SN=spanW詡2,…。 ⑶ 上述元素％ 輕，…,0V称为空间的基 用h代替V,在S上求泛函丿(弘)的极小，即求UN G SN满足 J(uN) = min J(v). 氏Sn (4) 显然，原先的变分问题(2)与后面的变分问题(4)是不同的，他们的解弘与也是不同 的。如果V的子空间必充分接近V,那么，用?代替V而得到的解柿也就充分接近u, 从而把作为原变分问题的近似解，亦即把无穷维空间V中的极值问题近似为有限维空间 S中的极值问题，这就是Ritz方法得基本思想。 问题（4）很容易化为一个线性代数方稈组问题。实际上，每一个中的元素都可 N 二 Y Cjipi 以用°,02，???，0n的线性组合表示，故 /=!,这时，由A （^v）及FW）的性质 TOC \o 1-5 \h \z I N N N 丿他）=牙A（工g,X5?） - F（工训） L /=1 ;=1 /=1 ]N N N = A0（P）CiCj -工 F?）Cj 2 Z=1 J=1 /=1 , 因此，丿（?）是一个以55…,Cn为变量的二次多项式Q（C1，C2，???,CQ。由于双线性泛函 ° （仏Q是对称、正定的，则二次项 tLZ A?（p）CiCj = - A（?，） 厶 /=1 ;=1 / 是一个以55…工N为变量的正定二次型，即矩阵[A（%9j）“N是对称正定的。因此， /⑷v）= °（C[,C2，???,Cn）在点（cf,c；,???,J）取得极小的充分必要条件是 器 l（d = O,心 1,2,…,N / O N 工 A@?）c； =F@）,d = l,2,???,N 可以算出，上式即戶 ’ ，即N维向量 C° =时,￡,???,氏卩是线性方程组Kc = f的解，其中 K = [A（?,伟）]爲,—[q,C2，...,c舁,/ =[尸（％）用輕）,???,尸（孤）]二求出（冲,￡,???,氏） 之后，即得近似变分问题（4）的解 N UN = W? i=\ 总之，Ritz法吧泛函的极小值问题化为多元二次函数的极小值问题，最后通过求解一 个线性代数方程组而得出变分问题的近似解。 Galerkin 方法 当双线性泛函A (^v)对称时，变分问题(2)等价于变分方程：求眈wV使得 A (w,v) = F(v),Vvg V ⑸ TOC \o 1-5 \h \z Galerkin方法是一种求解变分方程(5)的近似方法。假设V是可分的Hi lbert空间，仍然取有限维子空间为SN=span{(p[y(p2^(pN] 考虑近似变分方程：求n w 满足A (w?v,v) = F(v),Vvg Sn ⑹ 需要注意的是，近似变分方程(6)的解“N在Sn中，且检验函数V只要求对中的元素成 立而不是对V中的一切元素成立。 这时，求解n仍然化为求解一个线性方程组。因为中的元素仍是通过{0}的线性 N N 你二工弓%」二工 组合来表示，可以设 冃 i=l ,代入式(6)得 N N N 工 Z A((Pj,(P疋5 -工 )q = o Z=1 ；=1 /=1 N N 乞送 A?，?)c；-F@Jki=O 或者 日月 由 的任意性，即 “2,…,Cn得任意性，可得出线性代数方程组 N 工 = F(%)(1 iN) ;=1 由于双线性泛函A (^v)是对称的，所以该线性方程组与Ritz得到的线性方程组是相同 的。 (1)有限元方法 考虑两点边值问题(P) r / 1 p(x)— 4-q{x)u = /(x),a xhdx\ dx) v u(a)= 0 p(b)『(b)