上个月对Math.NET的基本使用进行了介绍，主要内容有矩阵，向量的相关操作，解析数据格式，数值积分，数据统计，相关函数，求解线性方程组以及随机数发生器的相关内容。这个月接着深入发掘Math.NET的各种功能，并对源代码进行分析，使得大家可以尽可能的使用Math.NET在.NET平台下轻易的开发数学...
这两个函数都是MATLAB中的内置函数，是基于梯形法则的数值积分公式 例如我们有函数y=x^3-2x-3，为了计算在[0,1]上的积分，可以这么做： 其中x和y分别是自变量和对应的值，trapz其实就是Trapezoidal（梯形的简写），cumtrapz函数和trapz函数使用...
分类：
其他好文 时间：
2015-04-03 14:56:22
阅读次数：
430
原文:【原创】开源Math.NET基础数学类库使用(08)C#进行数值积分 本博客所有文章分类的总目录：http://www.cnblogs.com/asxinyu/p/4288836.html 开源Math.NET基础数学类库使用总目录：http://www.cnblogs.com/asxinyu...
hdu 1724
题意：
求椭圆面积, x^2/a^2+y^2/b^2=1, 给出a,b,l,r, 求椭圆被直线x=l, x=r所夹面积。思路：
1. 自适应simpson，精度不是很高，所以eps要设小一点。
2. romberg，精度很好，但效率不是很高。自适应simpson：
{CSDN:CODE:617152}romberg：
{CSDN:CODE:617154}...
分类：
其他好文 时间：
2015-03-11 19:42:58
阅读次数：
138
hdu 4498 自适应simpson
题意：
给出k1,k2,…,kn, a1,a2,…,an 和 b1,b2,…,bn
求函数：
F(x)=min{100,min{ki*(x-ai)^2+bi | 0 < i <= n}}
在坐标上画出的曲线的长度。限制：
1 <= n <= 50; 0 <= ai,bi < 100; 0 < ki < 100思路：
先求出所有交点，然后排序，把函...
分类：
其他好文 时间：
2015-03-11 19:41:44
阅读次数：
116
在数值计算的需求中，数值积分也是比较常见的一个。我们也知道像Matlab，Mathematics等软件的积分求解功能非常高大上，不仅能求解定积分，还能求解不定积分，甚至多重积分等等。而Math.NET这个组件没有如此高级的功能，目前也只提供了比较件的闭区间上的定积分求解功能。今天就一起来看看，因为不...
（一）数值积分 一、数值积分的MATLAB实现方法： 1、变步长辛普生法(quad)法： （1）调用格式： [I,n]=quad('fname',a,b,tol,trace); fname是被积函数； a,b是积分上下限； tol来控制积分精度，默认为0.001； trace控制是否展现积分过程，默...
分类：
其他好文 时间：
2015-02-20 00:07:54
阅读次数：
275
对于积分：
只要找到被积公式的原函数F(x)，利用牛顿莱普利兹公式有：
但是，实际使用这种求积分的方法往往是有困难的，因为大量的被积函数的原函数是不能用初等函数表示的；另外，当f(x)是由测量或数值计算给出的一张数据表时，牛顿莱普利兹公式也无法直接运用，因此有必要研究积分的数值计算问题。
对于一些理论的推导，大家可以看看维基百科，下面我主要给出牛顿-科特斯公式在n...
分类：
其他好文 时间：
2015-02-02 12:38:58
阅读次数：
243
SciPy以NumPy为基础，提供了应用更加广泛的科学计算工具。
其在以下方面有着优秀的函数库：
1.线性代数
2.数值积分
3.插值
4.优化
5.随机数生成
6.信号处理
7.其他
与NumPy一样，SciPy有着稳定，成熟，且应用广泛的数值运算库。
许多SciPy函数仅仅是给诸如LAPACK，BLAS这样的Fortran数值计算工业标准库提供了接口。
在本文中，我们仅仅...
分类：
其他好文 时间：
2015-01-14 23:01:17
阅读次数：
1055
