3.5!——阶乘的插值(python)
本文的全部理论基础见于 神奇的Gamma函数(scipy)。
k!\Rightarrow x!,\quad x \geq 0
阶乘插值,在更大的范围里说这其实是一个 数列插值的问题,通俗地说,就是把数列的通项公式定义从 整数集合延拓到实数集合。
n\Rightarrow \mathbb{R}
例如数列 1,4,9,16,…1, 4, 9, 16, \ldots 可用通项公式 n2n^2自然表达,即便 nn为实数的时候,这个通项公式也有良好定义。直观地说,也就是可以找到一条平滑的曲线通过 y=x2y=x^2 通过所有的整数点( (n,n2)(n, n^2))这些点,从而可以把定义在整数集上的公式延拓到整个实数集合。可是对于阶乘序列 1,2,6,24,120,720,…,1, 2,6,24,120,720,\ldots,我们是否像计算 2.522.5^2一样计算 3.5!3.5!的呢。
答案是显然地:
k!=\int_0^\infty x^ke^{-x}dx
k∈Rk\in \mathbb{R} 对等式右边有意义,自然可通过等式右边计算正实数的阶乘,下面我们就用python(scipy)简单仿真实践插值:
import numpy as np
import matplotlib.pyplot as plt
from scipy import integratedef fact_int(x, k):return x**k*np.exp(-x)if __name__ == '__main__':x = np.arange(0, 4, .01)plt.plot(x, [integrate.quad(fact_int, 0, np.inf, args=k)[0] for k in x], 'g', lw=2)k = np.arange(0, 5)plt.plot(k, [integrate.quad(fact_int, 0, np.inf, args=k)[0] for k in k], 'r--', marker='o', markersize=10)plt.show()
3.5!——阶乘的插值(python)相关推荐
- python求n的阶乘并输出身份信息_python编程求n的阶乘_使用Python编程的阶乘
python编程求n的阶乘 Before we start implementing factorial using Python, let us first discuss what factori ...
- python编程求n的阶乘_使用Python编程的阶乘
python编程求n的阶乘 Before we start implementing factorial using Python, let us first discuss what factori ...
- python求5的阶乘_用python分享5的阶乘 (咋办)用python计算阶乘
(急)用python计算阶乘 小编需要用python列出一个阶乘的表. 小编不太清楚python的阶乘公式 以阶乘reduce为逐次操作list里的每项,接收的参数为 2个,最后返回的为一个结果 &g ...
- 拉格朗日插值python代码_Python实现的拉格朗日插值法示例
本文实例讲述了Python实现的拉格朗日插值法.分享给大家供大家参考,具体如下: 拉格朗日插值简单介绍 拉格朗日插值法是以法国十八世纪数学家约瑟夫·拉格朗日命名的一种多项式插值方法. 许多实际问题中都 ...
- python缺失值拉格朗日插值_拉格朗日插值-python
在数据库中,有些数据是异常值或者空值,这些值在分析的时候应该特殊处理,比如最简单的忽略掉或者通过算法推测它的值.其中拉格朗日插值就是通过其他已经知道的值,对x位置缺失的值插入的算法. 假定我们已经知道 ...
- 三次样条插值Python实现
函数 y=11+x2y = \frac{1}{1 + x^2}y=1+x21 算法分析 三次样条插值.就是在分段插值的一种情况. 要求: 在每个分段区间上是三次多项式(这就是三次样条中的三次的来源) ...
- C语言大数阶乘取余,python中math模块常用函数介绍 取模(取余)取绝对值 求阶乘 求最大公约数最小公倍数 取对数 取根号 取幂(取次方) 取整函数 三角函数与反三角函数...
前提:import math 两个常用常量 e = 2.718281828459045 pi = 3.141592653589793 >>> import math >> ...
- 十一届蓝桥杯 100阶乘的约数-python
100阶乘的约数 问题描述 定义阶乘 n! = 1 × 2 × 3 × ··· × n. 请问 100! (100 的阶乘)有多少个约数. 答案提交 这是一道结果填空的题,你只需要算出结果后提交即可. ...
- python判断阶乘和数_[C语言训练]阶乘和数*-题解(Python代码)
就是按部就班的遍历 最后字典方式排序输出,就是将每一个数字都转换为字符串,然后排序输出 answer_list = [] # 需要一个求阶乘的函数 def func(number): if numbe ...
最新文章
- 联机分析处理系统与联机事务处理系统的区别和联系_混合事务分析处理“HTAP”的技术要点分析...
- Sql server中如何将表A和表B的数据合并(乘积方式)
- 正则中除了空格还有什么标识间距
- 如何把Linux工具里的“军刀”BusyBox移植到RT-Thread Smart?
- (已解决) centos6.5 yum源 失效 The whole CentOS 6 is dead and shouldn’t be used anywhere at all
- CentOS7.0系统安全加固实施方案
- Rayman的绝顶之路——Leetcode每日一题打卡11
- Excel数据去重(删除重复项)
- 平江南江计算机培训,平江县南江学区: 假期培训促提升,收心归位再出发
- 阿里巴巴-数据平台事业部-数据产品部:寻找真爱
- geogebra动态数学软件,实用工具
- 赛迪数据丨我国新型显示产业累积总投资已超过1.2万亿元
- 全球及中国造纸行业十四五发展决策及运营能力预测报告2021-2027年
- app架构师,HTTPS面试常问全解析,吊打面试官系列!
- 仙剑5前传大吐糟 2013-2-5记录
- 项目日常管理制度大全,都列在这里了(共36大项)
- 微信小程序录音功能实现、录音转文字(接口与插件的实现)
- 【无限互联】iOS开发视频教程—1.iPhone开发概述——必看
- 带有汉字的字符串截断出现半个“汉字”的解决方法-C语言源码
- 将RGB图像转换为二值图像--MATLAB
热门文章
- Redis五大数据类型
- linux删除更新的系统内核,Linux系统如何删除旧内核
- oracle如何获取日期月份差,Oracle获取日期和月份
- mysql数据库设置索引 以下_如何设置MySQL数据库表的索引总长度
- android人脸识别demo_C#开发实录:基于免费SDK实现人脸识别应用开发
- 第八届蓝桥杯第二题等差素数列
- 剑指offer面试题33. 二叉搜索树的后序遍历序列(二叉树)(递归)
- 高性能加速深度学习算法
- HiveQL学习笔记(二):Hive基础语法与常用函数
- js 事件模型 + ( 事件类型 )