python应用应用-scipy,numpy,sympy计算微积分计算微积分
python应用应用-scipy,numpy,sympy计算微积分计算微积分
今天来讲一下使用python进行微积分运算，python有很多科学计算库都可以进行微积分运算，当然如果知晓微积分计算的原理
也可以自己编程实现。
下面我们用三种方式进行积分运算圆周率pi
numpy计算pi
import os
import numpy as np
#pi=4(1-1/3+1/5-1/7+1/9-.......)
n = 100000
print(np.sum(4.0 / np.r_[1:n:4, -3:-n:-4]))
#3.141572693
讲解一下上面的代码
首先，这里的编程思路来源于一个公司，代码中也有注解
pi=4((1-1/3+1/5-1/7+1/9-1/11+1/13…)
这个公式怎么来的，实话说我也不知道
上面的代码思路就是，不可能说按照上面的公式计算无数项从1-1/3+1/5…到无穷项，所以我们就尽可能地去计算，去逼近真
实结果，事实上，别的微积分计算函数也不能说就能计算无穷多项，其实也只是一种尽可能地逼近。上面的代码是取了1-
1/3…1/99999。取到这个范围。
还要提到一个方法np.r_[]是指按行连接两个矩阵
所以如果我们呢自己去编程其实也是这样的一个思路
下面我给大家举一个很典型的例子。
计算微函数y=x2y=x^2y=x2从0-5的积分，积分思路：将0-5的y=x2y=x^2y=x2图像进行切分，切成1000000个小矩形再求他们
的面积和，说到这里很多人可能想起了高等数学教材上有关于这一节的积分知识。
下面给出代码：
import os
import numpy as np
#pi=4(1-1/3+1/5-1/7+1/9-.......)
n=1000000
print(sum((5/1000000)*(np.linspace(0,5,1000000)**2)))
如果你用数学积分的知识计算，积分结果应该为125/3，上述结果的精确度已经达到0.00001了，下面给出图形并给出绘图代
码方便大家形象化了解
x=np.linspace(0,5,1000)
y=x**2
plt.figure()
plt.xlabel('x')
plt.ylabel('y')
plt.plot(x,y,label="$y=x^2$")
plt.legend()
x=np.linspace(0,5,100)
for i in x:
x_z=[] y_z=[] x_z.append(i)
y_z.append(0)
x_z.append(i)
y_z.append(i**2)
plt.plot(x_z,y_z)
plt.show()
运行结果：
weixin_38673237
- 粉丝: 2
- 资源: 844
我的内容管理
收起
- 我的资源
快来上传第一个资源
我的收益 登录查看自己的收益
我的积分
登录查看自己的积分
我的C币
登录后查看C币余额
我的收藏 
我的下载 
下载帮助
评论0