描述 |
在三个灵活而简短的MATLAB程序中,用最低阶Raviart-Thomas混合有限元实现了二维非均匀混合边界条件下拉普拉斯方程的数值逼近。第一种,混合实现(LMmfem)基于拉格朗日乘法器技术。第二种直接方法(EBmfem)对最低阶Raviart-Tomas有限元使用edge-basis函数。第三个分析(CRmfem)利用了Crouzeix和Raviart提出的P1非协调有限元方法,然后通过Marini提出的技术对离散通量进行后处理。本文的目的是推导、记录、说明和验证三种MATLAB实现EBmfem、LMmfem和CRmfem,以便在教育和研究中进一步使用和修改。为了监测离散化误差,采用了一种具有可靠和有效平均技术的后验误差控制。其中,重点是对mexed边界条件的正确处理。数值例子说明了所提供软件的一些应用以及误差估计的质量。 |