蒙特卡洛法求多重积分-Cpp

0~1随机分布的点列

C++类的写法：

class RND
{
private:int R;
public:RND(int r = 1){R = r;}double rnd1(){int m;double s, u, v, p;s = 65536.0;u = 2053.0;v = 13849.0;m = (int)(R / s);m = R / s;R = R - m * s;R = u * R + v;m = (int)(R / s);R = (int)(R - m * s);p = R / s;return(p);}
};

这样，调用rnd1就可以产生0~1随机分布的随机数。

蒙特卡洛算法的核心函数，

C++写法：

double mtml(int n, double a[], double b[], double (*f) (int, double[]))
{int m, i;double s, d, *x;RND r(1.0);x = new double[n];d = 65536.0;s = 0.0;for (m = 0; m <= 65535; m++){for (i = 0; i <=n - 1; i++)x[i] = a[i] + (b[i] - a[i])*r.rnd1();s = s + (*f) (n, x) / d;}for (i = 0; i <= n - 1; i++) s=s*(b[i] - a[i]);delete []x;return (s);}

接下来在主函数中写被积函数的写法的形式为：

double f(int n,double x[])
{double z;z=f(x0,x1,x2,...,xn-1)表达式；return z;
}

积分结果为：

mtml（n,a,b,f)

本文参考了清华大学出版社的《常用函数算法集（C++描述）》

完整代码请参考我的GitHubhttps://github.com/TYduoduo

建模：蒙特卡洛法求多重积分相关推荐

蒙特卡洛法求圆周率 c语言,c++蒙特卡洛法求圆周率
在linux下用 g++ pi.cpp -o pi -lm 编译. //蒙特卡洛法求圆周率 #include #include #include #include using namespace st ...
用Matlab实现蒙特卡洛法求心形线面积
用Matlab实现蒙特卡洛法求心形线面积思路: 1.生成一定数量的随机数,其大小在[0,1)之间: 2.用随机数确定一个长宽已知的矩形,并保证心形线在矩形中: 3.计算随机点落在心形线中的概率,然后 ...
蒙特卡洛法求非线性方程组
蒙特卡罗法的概念及应用蒙特卡洛法 (又称统计试验法)是描述装备运用过程中各种随机现象的基本方法,而且它特别适用于一些解析法难以求解甚至不可能求解的问题,因而在装备效能评估中具有重要地位. 用蒙特卡洛 ...
蒙特卡洛求圆周率c语言代码,蒙特卡洛法求圆周率
利用蒙特卡洛算法求圆周率是一个概率的方法,关于这方面的内容很多,而且也很容易理解,更多具体分析过程可以参考如下文章: http://blog.csdn.net/zzhflyqh/article/det ...
蒙特卡洛法求概率密度函数matlab,matlab的蒙特卡洛算法-----用概率算法求pi值
问题补充: 能提供一个例子看看吗?我就不懂MATLAB,想知道具体代码.因为现在写论文急,也没时间仔细看书了蒙特卡罗模拟就是随机数相关的东西,你只要知道随机数是怎么得到.其它的事就要好办了. ra ...
利用蒙特卡洛法求π的近似值：
蒙特卡洛法是一种非常重要的数值计算方法. 对于边长为2的正方形.半径为1的圆,在正方形内随机投掷大小可忽略不计的点,设落在圆内的概率为P,则有P=π/4. 对于投掷点的个数n,当n趋近于无穷时,可以 ...
用蒙特卡洛法求圆周率pi
思路正方形的面积为: 圆的面积为: 圆的面积比上正方形的面积为:π / 4 使用蒙特卡洛法在正方形内随机撒点,落在圆内的点 / 落在正方形内的点(全部的点),就约等于圆的面积 / 正方形的面积 = ...
蒙特卡洛法求圆周率100亿数据
代码 import time import random hits=0 pi=0 DARTS=100000*100000 start=time.perf_counter() for i in rang ...
用随机投点法计算pi值matlab,用蒙特卡洛法求圆周率PI的近似值，其算法如下：使用随机函数Rnd,随机生成两个[0,1）之间的随机数，由此构成一个投点坐标（x,y），然后根据x...
(2) 实现上述功能的VB程序如下,请在划线处填入合适代码. Private Sub Command1_Click() Dim m1 As String, m2 As String, sum As S ...

建模：蒙特卡洛法求多重积分

蒙特卡洛法求多重积分-Cpp

建模：蒙特卡洛法求多重积分相关推荐

最新文章

热门文章