浮点计算数值误差及PI的蒙特卡罗近似计算

2023-12-07 14:59

本文主要是介绍浮点计算数值误差及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的蒙特卡罗近似计算的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/466304

相关文章

poj 1113 凸包+简单几何计算

题意: 给N个平面上的点,现在要在离点外L米处建城墙,使得城墙把所有点都包含进去且城墙的长度最短。 解析: 韬哥出的某次训练赛上A出的第一道计算几何,算是大水题吧。 用convexhull算法把凸包求出来,然后加加减减就A了。 计算见下图: 好久没玩画图了啊好开心。 代码: #include <iostream>#include <cstdio>#inclu

uva 1342 欧拉定理(计算几何模板)

题意: 给几个点,把这几个点用直线连起来,求这些直线把平面分成了几个。 解析: 欧拉定理: 顶点数 + 面数 - 边数= 2。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc

uva 11178 计算集合模板题

题意: 求三角形行三个角三等分点射线交出的内三角形坐标。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vector>#include <

XTU 1237 计算几何

题面: Magic Triangle Problem Description: Huangriq is a respectful acmer in ACM team of XTU because he brought the best place in regional contest in history of XTU. Huangriq works in a big compa

音视频入门基础:WAV专题(10)——FFmpeg源码中计算WAV音频文件每个packet的pts、dts的实现

一、引言 从文章《音视频入门基础:WAV专题(6)——通过FFprobe显示WAV音频文件每个数据包的信息》中我们可以知道,通过FFprobe命令可以打印WAV音频文件每个packet(也称为数据包或多媒体包)的信息,这些信息包含该packet的pts、dts: 打印出来的“pts”实际是AVPacket结构体中的成员变量pts,是以AVStream->time_base为单位的显

计算数组的斜率,偏移,R2

模拟Excel中的R2的计算。         public bool fnCheckRear_R2(List<double[]> lRear, int iMinRear, int iMaxRear, ref double dR2)         {             bool bResult = true;             int n = 0;             dou

单位权中误差 详细介绍

单位权中误差(Unit Weight Error, UWE)是用于描述测量数据不确定性的一个统计量,特别是在地理信息系统(GIS)、导航和定位系统中。它主要用于评估和比较不同测量系统或算法的精度。以下是对单位权中误差的详细介绍: 1. 基本概念 单位权中误差(UWE): 定义:单位权中误差表示每个观测值(测量值)在估算中的标准误差。它是误差的一个统计量,主要用于评估测量系统的精度。单位:通常

GPU 计算 CMPS224 2021 学习笔记 02

并行类型 (1)任务并行 (2)数据并行 CPU & GPU CPU和GPU拥有相互独立的内存空间,需要在两者之间相互传输数据。 (1)分配GPU内存 (2)将CPU上的数据复制到GPU上 (3)在GPU上对数据进行计算操作 (4)将计算结果从GPU复制到CPU上 (5)释放GPU内存 CUDA内存管理API (1)分配内存 cudaErro

Java - BigDecimal 计算分位(百分位)

日常开发中,如果使用数据库来直接查询一组数据的分位数,就比较简单,直接使用对应的函数就可以了,例如:         PERCENT_RANK() OVER(PARTITION BY 分组列名 ORDER BY 目标列名) AS 目标列名_分位数         如果是需要在代码逻辑部分进行分位数的计算,就需要我们自己写一个工具类来支持计算了 import static ja

用Python实现时间序列模型实战——Day 14: 向量自回归模型 (VAR) 与向量误差修正模型 (VECM)

一、学习内容 1. 向量自回归模型 (VAR) 的基本概念与应用 向量自回归模型 (VAR) 是多元时间序列分析中的一种模型,用于捕捉多个变量之间的相互依赖关系。与单变量自回归模型不同,VAR 模型将多个时间序列作为向量输入,同时对这些变量进行回归分析。 VAR 模型的一般形式为: 其中: ​ 是时间  的变量向量。 是常数向量。​ 是每个时间滞后的回归系数矩阵。​ 是误差项向量,假