python麦克劳林级数展开

系统 1960 0

我们利用python的sympy模块能够很方便的进行科学计算,可以利用它对变量表达式进行泰勒展开并绘制图像,下面是我对sin(x)进行麦克劳林展开并绘制图像

            
              from matplotlib import pyplot as plt
import numpy as np
from mpl_toolkits.mplot3d import Axes3D
from sympy import * #用于求导积分等科学计算
from sympy.plotting import plot3d
from sympy.functions import exp

x = Symbol('x')#x 变量
print(help(plot))#对x = 0点泰勒展开 至x20次方无穷小 

p0 = plot(sin(x),(x,-2*pi,2*pi),line_color='b',ylim=(-1.5,1.5),title='funion',show=False)
#画出sin(x),范围-2*pi,2*pi,蓝色线,y刻度范围-1.5,1.5,暂不显示
exper1 = series(sin(x),x,n=4)#泰勒展开至第O(x**4)
p1 = plot(exper1.subs(O(x**4),0),(x,-2*pi,2*pi),line_color='br',title='O(x**4)',show=False)
#第O(x**4)项用0替代 颜色为蓝红色
exper2 = series(sin(x),x,n=6)#泰勒展开至第O(x**6)
p2 = plot(exper2.subs(O(x**6),0),(x,-2*pi,2*pi),line_color='g',title='O(x**6)',show=False)
#第O(x**6)项用0替代 颜色为绿
exper3 = series(sin(x),x,n=8)#泰勒展开至第O(x**8)
p3 = plot(exper3.subs(O(x**8),0),(x,-2*pi,2*pi),line_color='r',title='O(x**8)',show=False)
#第O(x**8)项用0替代 颜色为红
exper4 = series(sin(x),x,n=10)#泰勒展开至第O(x**10)
p4 = plot(exper4.subs(O(x**10),0),(x,-2*pi,2*pi),line_color='y',title='O(x**10)',show=False)
#第O(x**10)项用0替代 颜色为黄
p0.extend(p1)
p0.extend(p2)
p0.extend(p3)
p0.extend(p4)
p0.show()#显示图像

            
          
python麦克劳林级数展开_第1张图片 运行结果

可以看到,随着展开的项数越多,展开后的图形就越接近蓝线sin(x) ,黄线是展开至O(x**10)无穷小项的图像

 

 


更多文章、技术交流、商务合作、联系博主

微信扫码或搜索:z360901061

微信扫一扫加我为好友

QQ号联系: 360901061

您的支持是博主写作最大的动力,如果您喜欢我的文章,感觉我的文章对您有帮助,请用微信扫描下面二维码支持博主2元、5元、10元、20元等您想捐的金额吧,狠狠点击下面给点支持吧,站长非常感激您!手机微信长按不能支付解决办法:请将微信支付二维码保存到相册,切换到微信,然后点击微信右上角扫一扫功能,选择支付二维码完成支付。

【本文对您有帮助就好】

您的支持是博主写作最大的动力,如果您喜欢我的文章,感觉我的文章对您有帮助,请用微信扫描上面二维码支持博主2元、5元、10元、自定义金额等您想捐的金额吧,站长会非常 感谢您的哦!!!

发表我的评论
最新评论 总共0条评论