python 求定积分和不定积分示例
作者:福斯基 时间:2021-06-15 02:17:39
求f(x) = sin(x)/x 的不定积分和负无穷到正无穷的定积分
sin(x)/x 的不定积分是信号函数sig ,负无穷到正无穷的定积分为pi
import math
import numpy as np
import matplotlib.pyplot as plt
from sympy import * #用于求导积分等科学计算
def draw_plot_set():#设置画图格式
ax = plt.gca()
#改变坐标轴位置
ax.spines['right'].set_color('none')#删除原来轴
ax.spines['top'].set_color('none')#删除原来轴
ax.xaxis.set_ticks_position('bottom')#在0点处增加轴
ax.spines['bottom'].set_position(('data',0))
ax.yaxis.set_ticks_position('left')#在0点处增加轴
ax.spines['left'].set_position(('data',0))
#设置坐标名
plt.ylabel('f(x)')
plt.xlabel('x')
plt.grid(True)#打开网格
def dif(left,right,step):#求导 左右区间以及间隔
x,y = symbols('x y')#引入x y变量
#expr = x*pow(E,x)#计算表达式
expr = sin(x)/x
x_value = [] #save x value
y_value = [] #save x f(x) value
y_value_int = [] #save x f(x)_dot value
expr_int = integrate(expr,x)#求函数的不定积分 c=0
print(integrate(expr,(x,-oo,oo)))#对x求定积分 负无穷到正无穷
for i in np.arange(left,right,step):
x_value.append(i)
y_value.append(expr.subs('x',i))#将i值代入表达式
y_value_int.append(expr_int.subs('x',i))#将i值代入积分表达式
draw_plot_set()#设置画图格式
plt.plot(x_value,y_value,"b-",linewidth=1,label='f(x)='+str(expr)) #画图
plt.plot(x_value,y_value_int,"r-",linewidth=1,label='F(x)='+str(expr_int)) #画图
plt.legend()#显示图例
plt.show()#显示图像
if __name__ == '__main__':
draw_plot_set()#设置画图格式
dif(-30,30,0.1)
结果:
来源:https://blog.csdn.net/cool_bot/article/details/91044132
标签:python,定积分,不定积分
![](/images/zang.png)
![](/images/jiucuo.png)
猜你喜欢
python监控文件或目录变化
2023-09-05 16:08:28
讲解SQL Server海量数据导入的最快方法
2008-12-05 16:21:00
一文教你实现Python重试装饰器
2022-04-23 08:01:07
![](https://img.aspxhome.com/file/2023/9/85049_0s.jpg)
ASP万用分页程序
2007-09-21 12:45:00
PHP实现的微信公众号扫码模拟登录功能示例
2023-11-22 18:13:08
banner字体设计与应用
2009-07-06 14:42:00
![](https://img.aspxhome.com/file/UploadPic/20097/6/01-44s.jpg)
python实现简单通讯录管理系统
2021-05-02 10:41:23
Go语言集成开发环境之VS Code安装使用
2023-08-29 13:06:38
![](https://img.aspxhome.com/file/2023/9/101039_0s.png)
Python中使用第三方库xlutils来追加写入Excel文件示例
2022-05-23 10:04:11
PHP使用Swagger生成好看的API文档
2023-05-25 09:56:30
Python手写回归树的实现
2021-08-26 23:03:56
![](https://img.aspxhome.com/file/2023/3/65043_0s.png)
分享13款非常有用的jQuery插件
2011-05-16 19:07:00
PyTorch上实现卷积神经网络CNN的方法
2023-10-30 11:37:05
一百多行python代码实现抢票助手
2022-05-20 00:55:52
提高网页加载显示速度的方法
2007-08-10 13:17:00
原生JS实现左右箭头选择日期实例代码
2023-08-06 04:55:27
![](https://img.aspxhome.com/file/2023/8/55938_0s.png)
python实现单机五子棋
2021-04-26 05:34:35
![](https://img.aspxhome.com/file/2023/9/66459_0s.jpg)
解决Python3.7.0 SSL低版本导致Pip无法使用问题
2023-08-22 01:46:55
python tkinter模块的简单使用
2021-11-20 16:52:33
![](https://img.aspxhome.com/file/2023/5/75325_0s.png)
解决python 读取npy文件太大不能完全显示的问题
2021-08-20 12:07:19
![](https://img.aspxhome.com/file/2023/4/82974_0s.jpg)