C++ 计算正态分布积分

python,C#等语言计算正态分布积分的非常的方便。

提到C++,很多人就会觉得这个语言非常的笨, 如果不是为了追求效率,很少人会用C++写程序。 其实最近发现,C++有很多封装好的库, 比如矩阵计算有eigen库等。

今天想用C++计算一个标准正态分布的积分。 查了很多东西,包括boost等。

最后发现C++11可以直接算erf(x),稍微变换一下就可以计算高斯函数的积分了。erf的介绍在  https://en.wikipedia.org/wiki/Error_function

具体代码如下:

const double M_SQRT1_2 = sqrt(0.5);
double normalCFD(double value)
{
  return 0.5 * erfc(-value * M_SQRT1_2);
}

主要参考资料:https://stackoverflow.com/questions/2328258/cumulative-normal-distribution-function-in-c-c

时间: 2024-12-11 02:46:23

C++ 计算正态分布积分的相关文章

[复变函数]第22堂课 6.2 用留数定理计算实积分

3. 函数在 \infty 的留数 (1) 定义: 设 \infty 为 f 的孤立奇点, 则称 \bex \cfrac{1}{2\pi i}\int_{\vGa^-}f(z)\rd z\quad (\vGa:\ |z|=\rho) \eex 为 f 在 \infty 的留数, 记作 \dps{\underset{z=\infty}{\Res}f(z)} . (2) 若 f 在 r<|z|<\infty 内有 Laurent 展式 \bex f(z)=\cdots+\cfrac{c_{-n}}

hihocoder #1083 : 积分

时间限制:1000ms 单点时限:1000ms 内存限制:256MB 描述 在平面上有一个顶点数为N的多边形P,区域 你需要写一个程序计算这个积分 输入 输入包含T (T<=500)个测试用例.数字T在输入文件的第一行给出.每个测试用例的第一行是一个整数N代表多边形的点数.其后跟随N行,每行包含两个点Xi和Yi,表示第i个点的坐标,当我们以给定的顺序连接这些点,我们得到了一个多边形.题目保证多边形不自交,这个多边形的面积也不会是零,也不会出现3个相邻的点共线的情况. [参数说明] 3 <= N

matlab学习笔记第八章——积分

1.使用下面的格式就可以得到f的一个积分表达式: int(f) 2.当然,我们也可以使用int(f, v) 语法来调用int,其中f就是要积分的函数,而v是积分变量. 3.通过调用subs命令,使用给定的条件我们可以计算表达式的数值.要在一个命令中用数代替符号变量,我们要把被代替的符号列表和用来代替的数值列表都用花括号括起来.例如: >> subs(F,{b,x},{2,4}) ans = 23.0831 4.int还可以用来计算定积分,此时要把积分区间传递给它.如果我们输入int(f, a,

时域积分与频域积分 实现及对比

玩陀螺仪的都会遇到一个问题就是,陀螺仪输出的是角速度和线加速度.怎么把加速度转化成位移就值得研究一下. 首先我们讲一下傅立叶变换,傅立叶本身就是一个线性积分变换.主要是将信号在时域和频域中进行变换.因为我们相信任何一个信号都可以分解成sin函数.sin函数的频率,振幅可以组合成很多的信号形式. 傅立叶变换的数学公式是这样的. 简单的有一个示意动画就可以说明问题. 经过傅立叶变换,我们所谓的频域积分也都是基于sin函数的积分. 傅立叶函数有个积分性质,当积分函数进行傅立叶变换的时候有下面这个特性

IMU 预积分推导

给 StereoDSO 加 IMU,想直接用 OKVIS 的代码,但是有点看不懂.知乎上郑帆写的文章<四元数矩阵与 so(3) 左右雅可比 >提到 OKVIS 的预积分是使用四元数,而预积分论文中使用 so(3) 的右雅克比.才疏学浅,先整理好 so(3) 的预积分,写好 StereoDSO 加上 IMU,再考虑其他的东西. 以下的内容参考预积分的的论文,还有它的 Supplementary Material.预积分的论文中有一些 typo 所以看上去还是比较迷的,参考网络上多份预积分论文的

38、生鲜电商平台-会员积分系统的设计与架构

说明:互联网平台积分体系主要用于激励和回馈用户在平台的消费行为和活动行为,一个良好的积分体系可以很好的提升用户的粘性及活跃度. 一.互联网平台积分体系设计必要性 互联网平台积分体系是一个独立.完整的系统模块,主要用于激励和回馈用户在平台的消费行为和活动行为,通过积分体系可以激发与引导用户在平台的活跃行为,逐步形成用户对平台的依赖性和习惯性,提升用户对平台的黏度和重复下单率. 积分体系在保持系统独立性的同时,又与平台会员系统.商品系统.订单系统等具有紧密的关联性,积分体系的规划设计需与平台其他系统

小小知识点(二十)利用MATLAB计算定积分

一重定积分 1. Z = trapz(X,Y,dim) 梯形数值积分,通过已知参数x,y按dim维使用梯形公式进行积分 %举例说明1 clc clear all % int(sin(x),0,pi) x=0:pi/100:pi; %积分区间 y=sin(x); %被积函数 z = trapz(x,y) %计算方式一 z = pi/100*trapz(y) %计算方式二  运行结果 被积函数曲线 2.[q,fcnt]= quad(fun,a,b,tol,trace,p1,p2...) 自适应sim

CS考研_统考大纲

序号 政治 外语 业务课一 业务课二 1 (101)思想政治理论 (201)英语一 (301)数学一 (408)计算机学科专业基础综合 以上是计算机全国统考考试科目,三门公共课非统考基本也都是这三个,大家如果看到非统考的科目如果是三个1,就可以直接来参考我这里列出的大纲了!所以在此,我就直接列出最近的2015年考研这四个的考试大纲: 政治101: Ⅰ.考试性质 思想政治理论考试是为高等院校和科研院所招收硕士研究生而设置的具有选拔性质的全国招生考试科目,其目的是科学.公平.有效地测试考生掌握大学本

复旦考研信息整理

复旦大学考研网:http://www.kaofudan.com/kaoyanjingyan/ 复旦大学考研资料:http://yz.kaoyan.com/fudan/ziliao/ 复旦大学考研真题:http://download.kaoyan.com/list-12-t-1 复旦大学考研网:http://www.kaofudan.com/kaoyanziliao/ 计算机科学与技术专业:http://www.kaofudan.com/kaoyanziliao/jisuanjikexuejish