python 求定积分和不定积分示例


Posted in Python onNovember 20, 2019

求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)

结果:

python 求定积分和不定积分示例

以上这篇python 求定积分和不定积分示例就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
windows下python模拟鼠标点击和键盘输示例
Feb 28 Python
pyqt5自定义信号实例解析
Jan 31 Python
PyQt5主窗口动态加载Widget实例代码
Feb 07 Python
python实现傅里叶级数展开的实现
Jul 21 Python
利用pandas进行大文件计数处理的方法
Jul 25 Python
python中pip的使用和修改下载源的方法
Jul 08 Python
Django ORM 查询管理器源码解析
Aug 05 Python
python实现文件批量编码转换及注意事项
Oct 14 Python
基于python调用jenkins-cli实现快速发布
Aug 14 Python
python 下划线的不同用法
Oct 24 Python
从np.random.normal()到正态分布的拟合操作
Jun 02 Python
Python使用DFA算法过滤内容敏感词
Apr 22 Python
python 计算积分图和haar特征的实例代码
Nov 20 #Python
Python 函数绘图及函数图像微分与积分
Nov 20 #Python
python抓取多种类型的页面方法实例
Nov 20 #Python
Python超越函数积分运算以及绘图实现代码
Nov 20 #Python
python自动化实现登录获取图片验证码功能
Nov 20 #Python
python通过链接抓取网站详解
Nov 20 #Python
python爬虫之遍历单个域名
Nov 20 #Python
You might like
用文本文件实现的动态实时发布新闻的程序
2006/10/09 PHP
php 正则表达式小结
2009/08/31 PHP
php中的PHP_EOL换行符详细解析
2013/10/26 PHP
PHP的fsockopen、pfsockopen函数被主机商禁用的解决办法
2014/07/08 PHP
CentOS下与Apache连接的PHP多版本共存方案实现详解
2015/12/19 PHP
php mysql 封装类实例代码
2016/09/18 PHP
Yii2框架实现登录、退出及自动登录功能的方法详解
2017/10/24 PHP
php判断IP地址是否在多个IP段内
2020/08/18 PHP
IE6中使用position导致页面变形的解决方案(js代码)
2011/01/09 Javascript
JS cookie中文乱码解决方法
2014/01/28 Javascript
详解JavaScript数组的操作大全
2015/10/19 Javascript
EXT中单击button按钮grid添加一行(光标位置可设置)的实例代码
2016/06/02 Javascript
AngularJS指令中的绑定策略实例分析
2016/12/14 Javascript
jQuery模拟淘宝购物车功能
2017/02/27 Javascript
jQuery实现在HTML文档加载完毕后自动执行某个事件的方法
2017/05/08 jQuery
微信小程序之前台循环数据绑定
2017/08/18 Javascript
vue router demo详解
2017/10/13 Javascript
JavaScript数据结构之双向链表和双向循环链表的实现
2017/11/28 Javascript
jQuery实现的下雪动画效果示例【附源码下载】
2018/02/02 jQuery
vue技术分享之你可能不知道的7个秘密
2018/04/09 Javascript
npm scripts 使用指南详解
2018/10/08 Javascript
详解VSCode配置启动Vue项目
2019/05/14 Javascript
解决layui批量传值到后台操作时出现传值为空的问题
2019/09/28 Javascript
js new Date()实例测试
2019/10/31 Javascript
[49:54]Ti4 循环赛第三日 LGD vs Titan
2014/07/12 DOTA
[02:31]2014DOTA2国际邀请赛2009专访:干爹表现出乎意料 看好DK杀回决赛
2014/07/20 DOTA
python基于phantomjs实现导入图片
2016/05/13 Python
python 生成器协程运算实例
2017/09/04 Python
python3 selenium自动化 下拉框定位的例子
2019/08/23 Python
tensorflow -gpu安装方法(不用自己装cuda,cdnn)
2020/01/20 Python
tensorflow ckpt模型和pb模型获取节点名称,及ckpt转pb模型实例
2020/01/21 Python
Hotels.com南非:酒店预订
2017/11/02 全球购物
Linux如何修改文件和文件夹的权限
2013/09/05 面试题
驾驶员管理制度范本
2015/08/06 职场文书
Java基于字符界面的简易收银台
2021/06/26 Java/Android
java中为什么说子类的构造方法默认访问的是父类的无参构造方法
2022/04/13 Java/Android