1.SciPy是什么
SciPy是python的一个著名的开源科学库,SciPy一般都是操纵NumPy数组来进行科学计算,统计分析,可以说是基于NumPy之上。SciPy提供了许多科学计算的库函数,如线性代数,微分方程,信号处理,图像处理,系数矩阵计算等,
2.Scipy安装
官网:http://www.scipy.org/SciPy

安装:假定python安装目录在C:\Python27,则在C:\Python27\Scripts下打开cmd执行

执行命令:pip install scipy

python之SciPy简介

运用sicpy模块中integrate函数求定积分,不定积分,多重积分 ,以及含有断点如何求积分 并运用matplotlib画出含有断点函数的图形
代码如下

py import integrate
import numpy as np
import matplotlib.pyplot as plt
'''求定积分'''
#建立定积分函数含有 一重积分
def f(x):return x**2+1
v,err=integrate.quad(f,1,2)#求1到2的定积分 err误差
print("函数f的一重积分积分为{:0.4}".format(v))
#建立含有参数a、b的积分函数  一重积分
def F(x,a,b):return a*x**2+b
#可以用参数args传入a=2、b=1的值
v,err=integrate.quad(F,1,2,args=(2,1))
print("函数F的一重积分为{:0.4}".format(v))
#建立有断点的函数求积分 如 1/根号下x的绝对值
def w(x):return 1/np.sqrt(abs(x))
#w函数在x=0处存在断点 用quad()函数中的points可跳过求解
v,err=integrate.quad(w,-1,1,points=[0])
print("函数w的一重积分为{:0.4}".format(v))
#绘制图像
fig,ax=plt.subplots(figsize=(8,3))
x=np.linspace(-1,1,10000)
ax.plot(x,w(x),lw=2)
ax.fill_between(x,w(x),color='green',alpha=0.5)
ax.set_xlabel("$x$",fontsize=18)
ax.set_ylabel("$w(x)=1/sqrt(abs(x))$",fontsize=18)
ax.set_ylim(0,25)
# plt.show()
plt.draw()
plt.pause(1)# 间隔的秒数:1s
plt.close(fig)
#建立函数求二重定积分
def W(x,y):return x*y
def h(x):return x
v,err1=integrate.dblquad(W,1,2,1,h) #外层积分1到2 内层积分1到x
print("函数W(x,y)=x*y的二重积分为{:0.4}".format(v))
#建立函数求三重积分
f=lambda x,y,z:x
g=lambda x:(1-x)/2
h=lambda x,y:(1-x-2*y)
v,err=integrate.tplquad(f, 0, 1, 0, g, 0, h)
print("函数f=x的三重积分为{:0.4}".format(v))

运行结果

Python 基于sicpy求解定积分 ,不定积分以及多重积分相关推荐

  1. Matlab求解定积分/不定积分/微分

    使用Matlab求解定积分/不定积分 https://blog.csdn.net/qq_34374664/article/details/79186465 用MATLAB求定积分 https://bl ...

  2. matlab积分求解(定积分/不定积分)总结

    Matlab求积分(定积分/不定积分)@TOC 一些最近学习的总结,希望可以帮助到有需要的童鞋,不对的地方也欢迎批评指正 首先是求定积分: 方法一: 利用intergral 数值积分(用于求定积分): ...

  3. Python 基于积分原理计算定积分并可视化数值积分计算的动画过程

    文章目录 一.问题描述 二.代码实现 1.  并基于积分原理计算  ∫ 0 1 x 3 + 1  的值  1. \text { 并基于积分原理计算 } \int_{0}^{1} x^{3}+1 \te ...

  4. 使用Matlab求解定积分/不定积分

    https://blog.csdn.net/qq_34374664/article/details/79186465 https://blog.csdn.net/qq_27008079/article ...

  5. python辗转相除法求最大公约数的递归函数_Python基于辗转相除法求解最大公约数的方法示例...

    本文实例讲述了Python基于辗转相除法求解最大公约数的方法.分享给大家供大家参考,具体如下: 之前总结过一次高德纳TAOCP中的最大公约数求解,其实课后题中的算法修改要求实现的是辗转相除法求解最大公 ...

  6. python辗转相除_Python基于辗转相除法求解最大公约数的方法示例

    Python基于辗转相除法求解最大公约数的方法示例 本文实例讲述了Python基于辗转相除法求解最大公约数的方法.分享给大家供大家参考,具体如下: 之前总结过一次高德纳TAOCP中的最大公约数求解,其 ...

  7. python辗转相除法求最大公约数最小公倍数_Python基于辗转相除法求解最大公约数的方法示例...

    本文实例讲述了Python基于辗转相除法求解最大公约数的方法.分享给大家供大家参考,具体如下: 之前总结过一次高德纳TAOCP中的最大公约数求解,其实课后题中的算法修改要求实现的是辗转相除法求解最大公 ...

  8. 基于蒙特卡洛法的定积分求解

    基于蒙特卡洛法的定积分求解 0.摘要 本文首先介绍了蒙特卡洛法的起源和应用,接着详细推导了基于蒙特卡洛法求解定积分的数学理论公式.然后进行相关的实验研究.本文的实验共有三个,分别为蒙特卡洛法可行性研究 ...

  9. 基于遗传算法求解TSP问题(旅游路径规划,Python实现,超详细,可视化,结果分析)

    ps:作者是很用心写的,如果觉得不错,请给作者一点鼓励噢!(点赞收藏评论噢) 基于遗传算法求解TSP问题 摘要 巡回旅行商问题(TSP)是组合优化中的经典问题.常见的TSP问题求解算法例如穷举法.贪心 ...

最新文章

  1. React——Flow代码静态检查
  2. 学业水平测试计算机都考什么,高中学业水平考试与高中会考有什么区别
  3. 深度学习100例-卷积神经网络(CNN)花朵识别 | 第4天
  4. NYOJ 485 A*B Problem
  5. audio标签不展示_设计标签式PPT:排版简洁明了,强烈信息秩序感,你想试试吗?...
  6. 使用SAP UI5 Web Components开发React应用
  7. 如何查看某个employee被assign到了哪个sales organization上
  8. java中的多线程的示例
  9. jQuery框架学习第六天:jQuery中的Ajax应用
  10. Python3安装cx_Oracle连接oracle数据库实操总结
  11. ajax在success方法中打点不能调用
  12. 软件项目需求调研报告模板下载_软件项目需求分析报告模板
  13. winform直接控制云台_这款达到广播级机器 功能需求的 轻量级云台!真香!
  14. 淘客APP开发定制系统推荐
  15. 目前微型计算机硬件主要采用,目前使用的微型计算机硬件主要采用的电子器件是()。 A. 真空管 B. 晶体管 C. 大规模和超大规模集成电路...
  16. 【毕业设计】深度学习卫星遥感图像检测与识别 -opencv python 目标检测
  17. 高代第2章第1节n元排列:奇排列、偶排列的判断
  18. 点云配准论文复现:Robust generalized point cloud registration with orientational data based on expectation ma
  19. iOS地图定位导航与大头针的简单使用
  20. 【结构型模式】《大话设计模式》——读后感 (9)牛市股票还会亏钱?——外观模式

热门文章

  1. 智能停车场[简易版]
  2. C++程序设计课程设计报告——自助点餐系统
  3. \t\t门萨智商测试 有些面试的公司喜欢玩这个
  4. TensorFlow RNN MNIST字符识别Demo快速了解TF RNN核心框架
  5. 数学黑洞(二)任何数都逃不出的西西弗斯黑洞
  6. 享有盛誉的PHP高级教程
  7. GUI编程自学,编写自己的计算器!
  8. Struts2-052 漏洞复现
  9. 是的你没看错,js生成word文档
  10. 张朝阳开课手推E=mc²,李永乐现场狂做笔记!CEO当太久都忘了他是MIT物理博士