python泰勒公式_泰勒公式和Gamma函数

今天带大家玩下数学中的编程，难度可能有点大，数学不好的人请离开。

泰勒公式

大家知道泰勒公式吗？对它的理解有多深呢？数学中，泰勒公式是一个用函数在某点的信息描述其附近取值的公式。如果函数足够平滑的话，在已知函数在某一点的各阶导数值的情况之下，泰勒公式可以用这些导数值做系数构建一个多项式来近似函数在这一点的邻域中的值。泰勒公式还给出了这个多项式和实际的函数值之间的偏差

泰勒公式，也称泰勒展开式。是用一个函数在某点的信息，描述其附近取值的公式。如果函数足够平滑，在已知函数在某一点的各阶导数值的情况下，泰勒公式可以利用这些导数值来做系数，构建一个多项式近似函数，求得在这一点的邻域中的值。

就是我不知道那函数表达式，我想通过泰勒公式来计算值。

公式

其中，

表示f(x)的n阶导数，等号后的多项式称为函数f(x)在x0处的泰勒展开式，剩余的

是泰勒公式的余项，是

的高阶无穷小。

有名的泰勒级数：

下面咱们来用泰勒公式模拟

函数

代码编写import numpy as np

import math

import matplotlib as mpl

import matplotlib.pyplot as plt

def calc_e_small(x):

'''

计算前10个

:param x:

:return:

'''

n = 10

f = np.arange(1, n+1).cumprod() #阶乘

b = np.array([x]*n).cumprod() # 算的是x的n次方

# 这个就是公式

return np.sum(b / f) + 1

大家可能对np.cumprod()函数不知道，就是计算阶乘的

但是这有个问题，如果x是负数，不就完了吗？可是

有数的，虽然很小。那怎么算？

下面用matplotlib来绘制出来t1 = np.linspace(-2, 0, 10, endpoint=False)

t2 = np.linspace(0, 4, 20)

t = np.concatenate((t1, t2))

print(t) # 横轴数据

y = np.empty_like(t)

for i, x in enumerate(t):

y[i] = calc_e(x)

print('e^', x, ' = ', y[i], '(近似值)\t', math.exp(x), '(真实值)')

plt.figure(facecolor='w')

mpl.rcParams['font.sans-serif'] = ['SimHei']

mpl.rcParams['axes.unicode_minus'] = False

plt.plot(t, y, 'r-', t, y, 'go', linewidth=2, markeredgecolor='k')

plt.title('Taylor展式的应用 - 指数函数', fontsize=18)

plt.xlabel('X', fontsize=15)

plt.ylabel('exp(X)', fontsize=15)

plt.grid(True, ls=':')

plt.show()

图像如下

绘制正正弦函数import numpy as np

import math

import matplotlib as mpl

import matplotlib.pyplot as plt

def calc_sin_small(x):

x2 = -x ** 2

t = x

f = 1

sum = 0

for i in range(10):

sum += t / f

t *= x2

f *= ((2*i+2)*(2*i+3))

return sum

def calc_sin(x):

a = x / (2*np.pi)

k = np.floor(a)

a = x - k*2*np.pi

return calc_sin_small(a)

if __name__ == "__main__":

t = np.linspace(-2*np.pi, 2*np.pi, 100, endpoint=False)

print(t) # 横轴数据

y = np.empty_like(t)

for i, x in enumerate(t):

y[i] = calc_sin(x)

print('sin(', x, ') = ', y[i], '(近似值)\t', math.sin(x), '(真实值)')

print('误差：', y[i] - math.sin(x))

mpl.rcParams['font.sans-serif'] = ['SimHei']

mpl.rcParams['axes.unicode_minus'] = False

plt.figure(facecolor='w')

plt.plot(t, y, 'bo', t, y, 'r-', linewidth=2)

plt.title('Taylor展式的应用 - 正弦函数', fontsize=18)

plt.xlabel('X', fontsize=15)

plt.ylabel('sin(X)', fontsize=15)

plt.xlim((-7, 7))

plt.ylim((-1.1, 1.1))

plt.grid(True)

plt.show()

Gamma函数

伽玛函数（Gamma函数），也叫欧拉第二积分，是阶乘函数在实数与复数上扩展的一类函数。该函数在分析学、概率论、偏微分方程和组合数学中有重要的应用。与之有密切联系的函数是贝塔函数，也叫第一类欧拉积分。可以用来快速计算同伽马函数形式相类似的积分。

阶乘对于有数学基础的人来说都不陌生，简单理解就是数的累乘。10的阶乘10！=1098765432*1。但是我们有没有思考过，如分数的阶乘是如何运算的？有没有方法估算一个数的阶乘？

Γ函数是阶乘在实数上的推广import numpy as np

import matplotlib as mpl

import matplotlib.pyplot as plt

from scipy.special import gamma

from scipy.special import factorial

mpl.rcParams['axes.unicode_minus'] = False

mpl.rcParams['font.sans-serif'] = 'SimHei'

if __name__ == '__main__':

N = 5

x = np.linspace(0, N, 50)

y = gamma(x+1)

plt.figure(facecolor='w')

plt.plot(x, y, 'r-', x, y, 'mo', lw=2, ms=7)

z = np.arange(0, N+1)

f = factorial(z, exact=True) # 阶乘

print(f)

plt.plot(z, f, 'go', markersize=9)

plt.grid(b=True)

plt.xlim(-0.1,N+0.1)

plt.ylim(0.5, np.max(y)*1.05)

plt.xlabel('X', fontsize=15)

plt.ylabel('Gamma(X) - 阶乘', fontsize=15)

plt.title('阶乘和Gamma函数', fontsize=16)

plt.show()

最后，祝有所学习，有所成长

python泰勒公式_泰勒公式和Gamma函数相关推荐

python怎么用gamma函数_不完全gamma函数 python有gamma函数吗
matlab中gammainc(不完全伽马函数)是代表什么运算? 如题,matlab中gammainc函数是什么运算,比如A==gammainc(B/2,C/2) 麻烦问一下,图中的红框的下一步是怎么 ...
复平面上gamma函数_神奇的Gamma函数 (中)
Gamma 函数欣赏 Each generation has found something of interest to say about the gamma function. Perhaps ...
Python带_的变量或函数命名，带下划线的方法
python中的标识符可以包含数字.字母和_,但必须以字母或者_开头,其中以_开头的命名一般具有特殊的意义目录单前置下划线的名字: 1.作为类名或函数名时, 2.作为类的属性名或方法名时前后均带 ...
python function gamma_Python中使用Gamma函数的累积分布函数
我正在处理Schechter Luminosity功能,如下所示: phi(L)dL = norm. Factor * (L/Lstar)^(a) * exp (L/Lstar) d(L/Lstar) ...
按键控制led灯python程序_树莓派使用threading函数实现多按键控制LED灯
讲两个知识点.一个是关于上拉和下拉,另一个是关于threading函数. 1.上拉电阻和下拉电阻上拉就是把一个不确定的信号通过一个电阻连接到高电位,这样在开关断开是信号为高电位,开关合上时信号为低电 ...
python怎么用gamma函数_如何通俗的理解伽马(gamma)函数
如何通俗的理解伽马(gamma)函数 - 直觉,求导和示例我为什么要在乎garmma分布? 使用伽马函数定义了许多概率分布,例如伽马分布,Beta分布,狄利克雷分布,卡方分布和学生t分布等. 对于数 ...
python怎么用gamma函数_「gamma函数」Gamma函数 - seo实验室
gamma函数伽玛函数(Gamma函数),也叫欧拉第二积分,是阶乘函数在实数与复数上扩展的一类函数.该函数在分析学.概率论.偏微分方程和组合数学中有重要的应用.与之有密切联系的函数是贝塔函数,也叫第 ...
[转载] gamma函数stiriling公式_数学笔记|特殊函数(1):Gamma函数
参考链接: Python数学库math | gamma()函数使用教材: 王竹溪,郭敦仁.特殊函数概论设 ,积分称为 Gamma函数(第二类Euler积分). 性质1:Gamma函数有递推关系 ...
已知自然常数e的泰勒展开式是_泰勒公式：微分学的顶峰 (数学分析 · 导数的应用 (2))...
微分中值定理是构建微分学等理论的基本工具,但是它只能初步地解决问题.要问仅用导数可以将函数刻画到什么程度,还是要看泰勒公式. 所谓的刻画函数,就是我们试图用足够简单的函数作为尽量多的函数的近似.多项式 ...

python泰勒公式_泰勒公式和Gamma函数

python泰勒公式_泰勒公式和Gamma函数相关推荐

最新文章

热门文章