孙开心

摘要： 继续上一篇一维gauss积分的讨论，本文讨论二维gauss数值积分公式的使用，并给出数值实验。 [TOC] 一. 标准区间 按照一维情况类似方式，首先给出如下二重积分公式： $$ I = \int_{ 1}^1dx\int_{ 1}^1 f(x,y)dy=\Sigma_{k=1}^m\Sigma_{ 阅读全文

摘要： 在数值分析中，尤其是 有限元 刚度矩阵、质量矩阵等的计算中，必然要求如下定积分: $$ I=\int_a^b f(x)dx $$学好 gauss 积分也是学好 有限元 的重要基础，学过高等数学的都知道，手动积分能把人搞死（微笑脸），而且有些函数还不存在原函数，使用原始的手动算出原函数几乎是不现实的。 阅读全文

摘要： gallery 函数是 matlab 收集的一些重要的典型矩阵集合，主要用作算法的测试，在我的上一篇博客，我详细介绍了 "矩阵市场" 的使用。矩阵市场相当于是一个矩阵的大润发超市，各种矩阵应有尽有，而 gallery 相当于是矩阵的一个袖珍版杂货铺收集了一些典型的矩阵，供matlab自身算法测试使用 阅读全文

摘要： 矩阵市场收集了很多科学与工程实际问题离散后得到的大型稀疏矩阵，但是为了统一存储方式和提高效率，矩阵市场存储的矩阵格式并不是我们 matlab 可以直接使用的，本文介绍如何将下载的矩阵转化为 matlab 可以直接使用的稀疏矩阵。同时矩阵市场也提供转化为其他语言数据的方式，如Python，C语言、Fo 阅读全文