概率算法是以前一直没有接触过的算法,第一眼看见的时候真的觉得很厉害,这个算法的思想很简单但是运用好很难,大概就是将问题转化为几何图形,然后通过点的分布来求解我们需要的信息,这里的例子是求解π,
画的图确实有点丑,大概就是以正方形的变长为1,圆形的半径为1,那么我知道知道第一象限圆的面积为π/4,然后正方形的面积为1,那么我们可以通过落入圆形中的点数占的比例来估计π的值,顺便提一句,这个算法也叫蒙特卡洛概率算法好了,可以贴代码了。
package www.jk.monte;
import java.io.InputStreamReader;
import java.util.Random;
import java.util.Scanner;
/**
*
* @author jk 这个程序演示的是通过概率的方式来求取圆周率,反正我第一看的时候觉得很厉害,首先原点为中心,半径为1画一个圆,
* 然后在第一象限以原点为左下角画一个边长为1的正方形,然后1/4圆和正方形的面积比就是π/4,接下来我们通过函数随机的产生小于
* 小于1的x,y所构成的n个点,然后我们通过计算出其中落在圆形中的点数和一共的点数来计算出π。
*
*/
public class Test {
public static void main(String[] args) {
double PI = 0.0;
int n;
System.out.println("输入需要采样的点的个数");
Scanner input = new Scanner(System.in);
n = input.nextInt();
// 调用方法来得到PI
// invoke the method to get PI
PI = montePI(n);
System.out.println("PI:" + PI);
}
public static double montePI(int n) {
double PI = 0.0;
double x, y;
int res = 0;
for (int i = 0; i < n; i++) {
// 通过函数获取随机的x,y
// get random x,y
x = Math.random();
y = Math.random();
// 如果距离小于1,这说明在圆圈内,结果加1
// if the distance is small than 1,it turns the point is in the
// cycle,res+1
if ((x * x + y * y) < 1) {
res++;
}
}
// 通过比例算出PI
// calculate PI
PI = 4.0 * res / n;
return PI;
}
}