python正态分布函数计算误差极大_标准正态分布函数的快速计算方法

本文介绍了两种提高Python中正态分布函数计算速度的方法,包括A&S 26.2.17公式和查表法。这两种方法在保证误差在10^-7级别的同时,速度分别约为R中`pnorm()`的3倍和15倍。通过实例展示了这两种方法的实现和性能对比,强调在计算资源有限时,这些传统方法的智慧仍然有价值。
摘要由CSDN通过智能技术生成

标准正态分布的分布函数 Φ(x)Φ(x) 可以说是统计计算中非常重要的一个函数,基本上有正态分布的地方都或多或少会用上它。在一些特定的问题中,我们需要大量多次地计算这个函数的取值,比如我经常需要算正态分布与另一个随机变量之和的分布,这时候就需要用到数值积分,而被积函数就包含 Φ(x)Φ(x)。如果 Z?N(0,1),X?f(x)Z?N(0,1),X?f(x),ff 是 XX 的密度函数,那么 Z+XZ+X 的分布函数就是

P(Z+X≤t)=∫+∞?∞Φ(t?x)f(x)dxP(Z+X≤t)=∫?∞+∞Φ(t?x)f(x)dx

我们知道,Φ(x)Φ(x) 没有简单的显式表达式,所以它需要用一定的数值方法进行计算。在大部分的科学计算软件中,计算的精度往往是第一位的,因此其算法一般会比较复杂。当这个函数需要被计算成千上万次的时候,速度可能就成为了一个瓶颈。

当然有问题就会有对策,一种常见的做法是略微放弃一些精度,以换取更简单的计算。在大部分实际应用中,一个合理的误差大小,例如 10?710?7,一般就足够了。在这篇文章中,给大家介绍两种简单的方法,它们都比R中自带的 pnorm() 更快,且误差都控制在 10?710?7 的级别。

第一种办法来自于经典参考书 Abramowitz and Stegun: Handbook of Mathematical Functions

的公式 26.2.17。其基本思想是把 Φ(x)Φ(x) 表达成正态密度函数 ?(x)?(x) 和一个有理函数的乘积。这种办法可以保证误差小于 7.5×10?87.5×10?8,一段C++

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
对于Vue 3和Vite的大屏自适应,你可以使用CSS的媒体查询和Vue的响应式特性来实现。 首先,为了让你的大屏页面能够自适应不同的屏幕尺寸,你可以使用CSS的媒体查询来设置不同的样式。在你的Vue组件中,通过在<style>标签内添加媒体查询,根据不同的屏幕尺寸应用不同的样式。例如: ```html <template> <div class="container"> <!-- 页面内容 --> </div> </template> <style> .container { /* 普通屏幕样式 */ } @media screen and (min-width: 1200px) { .container { /* 大屏幕样式 */ } } </style> ``` 在上述示例中,`.container` 类选择器设置了普通屏幕下的样式,而`@media` 媒体查询选择器用于设置大屏幕下的样式,当屏幕宽度大于等于1200px时生效。 此外,Vue 3还提供了响应式特性,你可以使用`ref`或`reactive`来定义响应式数据,并在模板中根据不同的屏幕尺寸动态渲染内容。例如: ```html <template> <div> <h1 v-if="isLargeScreen">大屏幕内容</h1> <h1 v-else>普通屏幕内容</h1> </div> </template> <script> import { ref, onMounted, onBeforeUnmount } from 'vue' export default { setup() { const isLargeScreen = ref(false) const handleResize = () => { isLargeScreen.value = window.innerWidth >= 1200 } onMounted(() => { window.addEventListener('resize', handleResize) handleResize() }) onBeforeUnmount(() => { window.removeEventListener('resize', handleResize) }) return { isLargeScreen } } } </script> ``` 上述示例中,我们使用`ref`定义了一个名为`isLargeScreen`的响应式变量,并通过`window`对象的`resize`事件来监听窗口大小变化。在模板中,根据`isLargeScreen`的值来动态渲染不同的内容。 希望以上信息能对你有所帮助!如有更多问题,请随时提问。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值