Python 求高斯误差函数 erf 和 erfc

文章目录

  • Part.I Introduction
  • Part.II 概念定义
    • Chap.I 误差函数 erf
    • Chap.II 误差函数补 erfc
  • Part.II 求值与绘图
    • Chap.I 求取方式
    • Chap.II 绘图
  • Reference

Part.I Introduction

本文将对误差函数(ERror Function, ERF)进行简单的介绍,并探索如何用 Python 来求误差函数。

在这里插入图片描述

更多与之相关的理论知识可以参看 关于高斯函数 / 高斯分布 / 高斯误差函数的探讨

Part.II 概念定义

Chap.I 误差函数 erf

误差函数(ERror Function, ERF)的定义式如下
erf ⁡ ( x ) = 2 π ∫ 0 x e − t 2 d t \operatorname{erf}(x)=\frac{2}{\sqrt{\pi}}\int_{0}^{x} e^{-t^2}dt erf(x)=π 20xet2dt

Chap.II 误差函数补 erfc

有的地方叫做『误差补函数』或『互补误差函数』,这里我随意取成了误差函数补(ERror Function Complement, ERFC)。一是参考了『舒尔补』这一名字,二是和字母顺序保持一致。

误差函数补(ERror Function Complement, ERFC)的定义式如下

erfc ⁡ ( x ) = 1 − erf ⁡ ( x ) = 2 π ∫ x ∞ e − t 2 d t \operatorname{erfc}(x)=1-\operatorname{erf}(x)=\frac{2}{\sqrt{\pi}}\int_{x}^{\infty} e^{-t^2}dt erfc(x)=1erf(x)=π 2xet2dt

PS: 2 π ∫ 0 ∞ e − t 2 d t = 1 \frac{2}{\sqrt{\pi}}\int_{0}^{\infty} e^{-t^2}dt=1 π 20et2dt=1

Part.II 求值与绘图

Chap.I 求取方式

其实很简单

import mathmath.erf(x)
math.erfc(x)

比如

>>> import math
>>> math.erf(1)
0.8427007929497149
>>> math.erfc(1)
0.1572992070502851
>>> 1-math.erf(1)
0.1572992070502851

GeoGebra 验证如下:

在这里插入图片描述

Chap.II 绘图

使用下面的代码可以绘制处误差函数和误差函数补的图像

import math
import matplotlib.pyplot as pltdef plot_erf():x=list(range(0,200))x=[i*0.02 for i in x]y1=[math.erf(i) for i in x]y2=[math.erfc(i) for i in x]plt.plot(x,y1)plt.plot(x,y2)plt.legend(["erf","erfc"])plt.show()plot_erf()

绘图结果为:
在这里插入图片描述

Reference

  • 浅谈误差函数 erf(x)
  • 关于高斯函数 / 高斯分布 / 高斯误差函数的探讨

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/bicheng/11510.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【Flutter】极光推送配置流程(VIVO/OPPO/荣耀厂商通道) 章三

前言 很高兴大家来看小编写的文章~~ 继【Flutter】极光推送配置流程(极光通道/华为厂商/IOS) 章一 继【Flutter】极光推送配置流程(小米厂商通道) 章二 接下配置VIVO/OPPO/华为荣耀的厂商通道 所有截图来源于公司项目,所以会有大量马赛克&am…

【C/C++】详解不同场景下的文件读写操作

🔗 运行环境:C/C 🚩 撰写作者:左手の明天 🥇 精选专栏:《python》 🔥 推荐专栏:《算法研究》 🔐#### 防伪水印——左手の明天 ####🔐 💗 大家好&…

虚拟数字人及AI相关应用分享

一、虚拟数字人 1、简介 虚拟数字人可分为基础类和仿真智能类。可用于直播的,一般是仿真智能类;基础类动作缓慢,体验差,很容易被直播平台封号。 目前各大短视频平台上介绍的数字人,出于营销目的,有夸大宣传…

代码随想录35期Day38-Java(Day37休息)

Day38题目 LeetCode509.斐波那契数列 核心思想:很简单dp[i]dp[i-1]dp[i-2].这里用了数组存储的形式,也可以递归 class Solution {public int fib(int n) {int[] dp new int[n2];dp[0] 0;dp[1] 1;for(int i 2 ; i < n ; i ){dp[i] dp[i-1] dp[i-2];}return dp[n];} …

go语言冒泡排序算法之泛型实现 支持整型,浮点型和字符串排序 golang bubble sort general

冒泡排序是一个比较经典的基础算法&#xff0c; 是学习编程的同学们所必备的算法之一&#xff0c;今天就介绍一个go语言里面使用泛型实现的冒泡排序算法&#xff0c; 支持字符串 string&#xff0c;数字 int float 等的排序。 在go语言中通常我们在进行数据比较时都是使用的c…

prophet时间序列模型水质预测应用

前言 此前已经分析了&#xff0c;ARIMA 模型在水质预测中的应用&#xff0c;今天用 prophet 模型测试下在水质预测中的效果。 Prophet 简介 Prophet 是 Facebook 于2017年开源的一个时间序列预测框架&#xff0c;特别适合于处理具有明显趋势性和季节性的数据。该模型设计初衷…

领略Java内部类的“内部”

内部类有两种情况&#xff1a; (1) 在类中定义一个类(私有内部类&#xff0c;静态内部类) (2) 在方法中定义一个类(局部内部类&#xff0c;匿名内部类) 1、私有内部类 —— 在方法之间定义的内部类&#xff0c;非静态 我们首先看看类中内部类的两个特点&#xff1a; (1) 在外部…

springboot3项目练习详细步骤(第三部分:文章管理模块)

目录 发布文章 接口文档 业务实现 自定义参数校验 项目参数要求 实现思路 实现步骤 文章列表(条件分页) 接口文档 业务实现 mapper映射 更新文章 接口文档 业务实现 获取文章详情 接口文档 业务实现 删除文章 接口文档 业务实现 文章管理业务表结构…

数据可视化训练第6天(美国人口调查获得关于收入与教育背景的数据,并且可视化)

数据来源 https://archive.ics.uci.edu/dataset/2/adult 过程 首先&#xff1b;关于教育背景的部分翻译有问题。 本次使用字典嵌套记录数据&#xff0c;并且通过lambda在sorted内部进行对某个字典的排序&#xff0c;最后用plotly进行绘图 本次提取数据的时候&#xff0c;用到…

ubuntu server 22.04.4 系统安装详细教程

本教程使用vmware workstation 17创建虚拟机进行安装演示&#xff0c;安装方式和真机安装没有区别。 1、下载镜像 下载ubuntu server版本系统镜像&#xff0c;官网下载地址&#xff1a;https://cn.ubuntu.com/download/server/step1 注意&#xff1a;自己下载时需要确认是否是…

MySQL索引的存储结构

2.3.2. 索引存储结构 索引存储结构有堆表和索引组织表两种方式。 堆表和索引组织表有什么区别&#xff1f; 堆表&#xff1a; 数据和索引是分开存放的&#xff0c;索引是排序后的&#xff0c;但堆数据是无序的进行的都是随机访问&#xff0c;索引的叶子节点中存放的是数据在…

基于springboot + vue 实现的简易博客系统

项目效果图 登陆页面 文章列表 发表文章 用户管理 栏目管理 数据统计 后端技术栈后端主要采用了&#xff1a; 1.SpringBoot 2.SpringSecurity 3.MyBatis 4.部分接口遵循Restful风格 5.MySQL 前端技术栈前端主要采用了&#xff1a; 1.Vue 2.axios 3.Elemen…

每周一算法:传递闭包

题目描述 不等式排序 给定 n n n个变量和 m m m个不等式。其中 n n n小于等于 26 26 26&#xff0c;变量分别用前 n n n 的大写英文字母表示。 不等式之间具有传递性&#xff0c;即若 A > B A>B A>B 且 B > C B>C B>C&#xff0c;则 A > C A>C …

GPU prompt

提问&#xff1a; GPU是如何与CPU协调工作的&#xff1f; GPU也有缓存机制吗&#xff1f;有几层&#xff1f;速度差异是多少&#xff1f; GPU渲染流程有哪些阶段&#xff1f;他们的功能分别是什么&#xff1f; Early-Z技术是什么&#xff1f;发生在哪个阶段&#xff1f;这个…

Cocos 2048从创建到发布上线

二、审核通过之后上线流程 代码通过审核之后&#xff0c;会通过站内信和微信消息发送通知&#xff0c;在管理后台&#xff0c;点击提交发布&#xff0c;去备案

ensp简单ac+ap+sta无线配置和脚本

接入交换机与ap连线配置: interface E0/0/5 port link-type trunk port trunk pvid vlan 10 port trunk allow-pass vlan all配置AC的IP配置&#xff1a; [AC]Vlan 2 创建vlan 2 [AC]interface Vlanif 2 进入vlan 2 [AC-Vlanif2]ip address 192.168.80.3 24 设置一个IP地址…

File类~路径、创建文件对象

路径分为相对路径&#xff08;不带盘符&#xff09;&#xff0c;绝对路径&#xff08;带盘符&#xff09; 路径是可以存在的&#xff0c;也可以是不存在的 创建文件对象的三个方法&#xff1a;

IDEA-SpringBoot项目启动类位置错误导致启动报错

文章目录 前言&#xff1a;本文为SpringBoot项目启动类位置错误导致启动报错解决方案 笔者在初学SpringBoot项目时遇到过因启动类位置错误导致项目启动报错现象&#xff0c;总结如下&#xff1a; SpringBoot启动发送请求后遇到如下错误&#xff1a; Postman界面信息&#xff1…

速了解及使用布隆过滤器

布隆过滤器 介绍 概念&#xff1a;是一种高效查询的数据结构 作用&#xff1a;判断某个元素是否在一个集合中。&#xff08;但是会出现误判的情况&#xff09; 实现原理 加入元素&#xff1a; 当一个元素需要加入到布隆过滤器中时&#xff0c;会使用一组哈希函数对该元素进…

个人模拟面试java2

文章目录 SQL如何进行优化数据库性能调优1. 性能评估2. 确定调优目标 3. 监控和分析4. 硬件和配置调优5. 数据库设计和模式优化6. 查询优化7. 并发控制8. 缓存策略9. 定期维护10. 监控和持续调优 Redis支持哪些数据类型&#xff1f;Redis如何实现持久化&#xff1f;如何解决Red…