码迷,mamicode.com
首页 > 编程语言 > 详细

概率算法 蒙特卡罗算法

时间:2014-12-12 00:05:55      阅读:205      评论:0      收藏:0      [点我收藏+]

标签:概率算法   蒙特卡罗算法   

/***Date:2014.12.11***/
//概率算法思想:统计学思路。
//基本过程:1)将问题转化为应的容易计算面积的几何图形S,问题结果对应几何图形中的某一部分S1;
////////////2)向几何图形中随机撒点;
////////////3)统计几何图形S、S1中的点数,根据二者面积关系以及二者中的点数来计算得到结果;
////////////4)判断结果是否达到需要精度结果,若符合则输出,否则继续(2)步骤。
//蒙特卡罗(Monte Carlo)算法
//均匀撒点:利用随机函数实现,产生(0,1)之间随机坐标值[x,y];
//区域判断:通过计算 x^2 + y^2 <=1 实现。

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>

double MontePI(int n)
{
double PI,x,y;
int i,sum;
sum = 0;
srand((unsigned)time(NULL));//64位要强转一下,要不然虽可执行,但不安全,PS,我电脑VS都64位
for (i = 1;i < n;i ++)
{
x=(double)rand()/RAND_MAX;
y=(double)rand()/RAND_MAX;
if (x*x + y*y <= 1)
{
sum ++;
}
}
PI = 4.0 * sum /n;
return PI;
}

int main()
{
int n;
double PI;

printf("蒙特卡罗概率算法计算 π:\n");
printf("输入点的数量(数量越多,结果越接近真实值):");
scanf_s("%d",&n);
PI = MontePI(n);
printf("π =  %f\n",PI);

system("pause");
return 0;
}

概率算法 蒙特卡罗算法

标签:概率算法   蒙特卡罗算法   

原文地址:http://blog.csdn.net/sdgtliuwei/article/details/41878265

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!