本文主要是介绍浮点计算数值误差及PI的蒙特卡罗近似计算,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
看JAVA语言程序设计基础,因为以前学习的都是C++, 总体基本上相似。不过这本书里有几个比较有意思的小地方,比如书中写的最小化数值误差(4.7),书中说在大数之前先增加小数是减少误差的方法。如:
//第一种方法求0.01到1之间递增的数列之和
for(float i = 0.01f; i<= 1.0f; i = i + 0.01f)
{
sum +=i;
}//第二种
for(double i = 0.01; i<= 1.0; i = i + 0.01)
{
sum +=i;
}//第三种
currentValue = 0.01;
for(int i = 0; i< 100; i++)
{
sum += currentValue;
currentValue += 0.01;
}//第四种
currentValue = 1.0;
for(int i = 0; i< 100; i++)
{
sum += currentValue;
currentValue -= 0.01;
}
第三种方法优于第四种,因为加大数前先加小数会提高精度。
近似的计算PI的值(4.8.3),画一个长度为2的矩形,中间内接一个圆,则理论上落在圆内的概率为PI/4。用蒙特卡罗模拟随机数落在圆内的点数。判断是否落在圆内,用坐标x^2+y^2 <=1 来判断。
public class ComputeArea {public static void main(String[] args) {final int numPoint = 10000000;int rHitNum = 0;for (int i = 0; i < numPoint; i++) {double x = Math.random() * 2 - 1;double y = Math.random() * 2 - 1;if ((Math.pow(x, 2) + y * y) <= 1) {rHitNum++;}}double pi = 4 * ((double)rHitNum / (double)numPoint); //记得强制类型转化为double再除,书中好像有问题System.out.printf("the PI is : %f", pi);}}
这篇关于浮点计算数值误差及PI的蒙特卡罗近似计算的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!