numpy
scipy
matplotlib
seaborn
没有合适的资源?快使用搜索试试~ 我知道了~
基于 Python 实现的有限元方程求解程序源码+程序说明.zip
共25个文件
png:20个
py:2个
txt:1个
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
5星 · 超过95%的资源 4 下载量 31 浏览量
2023-08-10
15:07:34
上传
评论 1
收藏 1.89MB ZIP 举报
温馨提示
基于 Python 实现的有限元方程求解程序源码+程序说明.zip 旨在通过有限元方法解决非齐次两点边值问题。它用于近似求解在给定域内的微分方程解。该项目主要利用numpy和scipy库,将域离散化为较小的单元,并使用适当的基函数构建 Ritz-Galerkin 方程。最后分别通过应用高斯消去法和雅可比迭代等数值技术,求解生成的线性方程组。 结果分析 从高斯消去法求解的结果来看,$u_h$函数近似估计精确解的效果很好,节点处的数值解与精确解的值几乎是重合的,而且随着h的减小误差也不断减少,当h=0.01时,误差的尺度为1e-5至1e-4,基本可以忽略不计。 从雅可比迭代法求解的结果来看,$u_h$函数近似估计精确解的效果不太好,节点处的数值解与精确解的值之间误差较大,而且随着h的减少,误差下降到一定程度(1e-2至1e-1)后不再下降。经过程序检验发现造成雅可比迭代不收敛的原因在于对有限元方程构建的总刚度矩阵是一个非对角占优矩阵,即不满足雅可比迭代的收敛要求,所以通过雅可比迭代法求解线性方程组 无法得到收敛的数值解。
资源推荐
资源详情
资源评论
收起资源包目录
基于 Python 实现的有限元方程求解程序源码+程序说明.zip (25个子文件)
Uh-jacobi-0.02.png 92KB
Uh-gauss-0.2.png 115KB
Error-gauss-0.02.png 155KB
Error-jacobi-0.1.png 102KB
Uh-gauss-0.01.png 115KB
Error-gauss-0.2.png 145KB
main.py 356B
Error-jacobi-0.05.png 126KB
LICENSE 1KB
Uh-gauss-0.05.png 114KB
Error-jacobi-0.01.png 141KB
Error-jacobi-0.2.png 122KB
Uh-jacobi-0.2.png 113KB
Uh-jacobi-0.01.png 90KB
Uh-gauss-0.1.png 114KB
Uh-jacobi-0.1.png 103KB
requirements.txt 30B
程序说明.md 5KB
Error-jacobi-0.02.png 138KB
Error-gauss-0.1.png 149KB
Uh-jacobi-0.05.png 101KB
Error-gauss-0.05.png 140KB
Error-gauss-0.01.png 158KB
femsolver.py 5KB
Uh-gauss-0.02.png 114KB
共 25 条
- 1
资源评论
- ywy20202024-03-17非常有用的资源,可以直接使用,对我很有用,果断支持!
- 一只鸟2432024-02-07这个资源值得下载,资源内容详细全面,与描述一致,受益匪浅。
- weixin_498637652024-01-23资源使用价值高,内容详实,给了我很多新想法,感谢大佬分享~
- m0_628933382023-12-13资源很实用,对我启发很大,有很好的参考价值,内容详细。Make程序设计2024-01-16感谢支持和认可,加油啊
Make程序设计
- 粉丝: 5629
- 资源: 3568
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功