偏微分方程程序.doc

成 绩 2009级数学与应用数学和信息与计算科学专业 偏微分方程数值解上机实验 实验题目　 利用有限元方法和有限差分方法求解偏微分方程 完成日期 2012年12月17日 学生姓名　 张灵刚 所在班级　 1102090 任课教师　 王晓东 西北工业大学理学院应用数学系 目录 一．实验目的………………………………………..…….（2） 二．实验要求…………………………………….…...….（2） 三．实验题目……………………………………….………（3） 四．实验二……………………………………………….…（4） 1.实验内容…………………………………….…..（4） 2.实验原理……………………………………….…（4）. 3算法流程………………………………….………（5） 4结果分析…………………………………………..（5） 5总结讨论…………………………………………..（6） 6源程序………………………………………..……（6） 五. 实验三 1.实验内容…………………………………………（17） 2.实验原理…………………………………….….（17）. 3算法流程………………………………………….（18） 4结果分析……………………………………….….（18）. 5总结讨论………………………………………….（21） 6源程序……………………………………………（21） 偏微分方程数值解上机实验报告 实验地点：数学系机房 实验时间：第13—15周，周一、四下午5、6节 实验分数：占期末考试成绩的30% 一、实验目的及意义 掌握有限元方法和有限差分方法的程序实现；学会选择合适的有限差分格式求解一维非线性对流占优的非定常对流扩散问题；学会使用三角线性元和四边形线性元的有限元方法求解二维椭圆方程边值问题，并对计算结果进行收敛性分析；尝试采用有限元方法或有限差分方法实现二维初边值抛物型方程的大规模数值求解。通过实验可以提高学生的动手能力，加深学生对算法的理解。 二、实验要求 在下列给出的三个问题中，最少选择两个问题进行编程实现。要求给出格式的推导过程、算法流程、实现程序、选取的网格参数、以表格或图形的方式给出计算结果、对计算结果进行分析、最后对实验进行总结和讨论。  问题2：用三角线性元和四边形线性元的有限元方法求解方程： 取比较两种方法的计算精度，并给出数值收敛率. 问题3：选用合适的数值方法求解方程： 求时，点、、、、、、、、处的数值解。 上机实验（二） 一、实验内容 用三角线性元和四边形线性元的有限元方法求解方程： 取比较两种方法的计算精度，并给出数值收敛率 二、实验原理 由于计算机只能存储有限个数据和做有限次运算，所以必须把连续问题离散化，有限元法通过把求解偏微分问题转化为变分问题，实质上就是Ritz-Galerkin法，利用有穷维空间近似代替无穷维空间，从而转化为求多元二次函数的极值问题。然后选定单元形状，对求解域进行剖分。构造基函数或单元形状函数，形成有限元空间了，再形成有限元方程，并提供求解有限元方程的有效方法。 三、算法流程 四、计算结果及分析 空间步长h 三角元实际误差 四边形元实际误差 0.129816 0.027308 0.016676 0.003659 0.002081 0.000457 0.000260 0.000057 0.000032 0.000007 结果分析：通过给出的上述结果可以发现三种方法相比四边形元的实际求解误差最小，这是由于四边形型元中的基函数中含有项，使得四边形元的误差远比三角形元小，且当空间步长不断加细时，向前差分的误差逐渐和四边形元接近。而对于收敛率的求解，如果舍去误差，三种方法的收敛率都约为3. 五、总结和讨论 有限元计算的有关问题有：把初值问题化为变分形式，对求解域作网络分割，构造基函数（或单元形状函数），形成有限元方程。通过本次实验，我懂得了用有限元方法求解初值问题的基本数学思想,理解了线性有限元法的基本原理和方法。另外，我也懂得了按Galerkin方法推导有限元方程的优点，它比Ritz法更加方便直接。我也对虚