numpy 定积分案例

定积分： 直观的说，对于一个给定的正实值函数，在一个实数区间上的定积分可以理解为坐标平面上由曲线、直线以及轴围城的去边梯形的面积值(一个确定的实数值)

案例

求二次函数y = 2x^2 + 3x + 4在[-5,5]区间内的积分

import numpy as np
import matplotlib.pyplot as mp  #画图用
import matplotlib.patches as mc #图像打补丁用
import scipy.integrate as si    #numpy求定积分用#声明曲线函数
def f(x):return 2 * x ** 2 + 3 * 4 + 4#1. 在区间[-5,5]间拆出1000个满足f(x)的点,画出该函数曲线
a, b = -5, 5
x1 = np.linspace(a, b, 1001)
y1 = f(x1)
mp.figure('Integral',facecolor='lightgray')
mp.title('Integral',fontsize=20)
mp.xlabel('x',fontsize=14)
mp.ylabel('y',fontsize=14)
#设置刻度参数，labelsize设置刻度字体大小
mp.tick_params(labelsize=10）
#显示网格
mp.grid(linestyle=':')
mp.plot(x1,y1,c='orangered',linewidth=6,label=r'$y=2x^2+3x+4$',zorder=0)#2. 微元法求积分
#把区间拆成50份(份数拆的越多，积分求得越精细)
n = 50
x2 = np.linespace(a,b,n+1)
#将x2传入f(x)获取y2
y2 = f(x2)
#通过计算拆出来的50个梯形的面积和，来求积分
area = 0
for i in range(n):#每个梯形的面积 = (上底+下底)*高/2area += (y2[i + 1] + y2[i] * (x2[1] - x2[0]) / 2
print(area) #用微元法求得的积分#将微元法分出的梯形在图中呈现出来
for i in range(n):mp.gca().add_patch(mc.Polygon([[x2[i],0],[x2[i],y2[i]],   #梯形的左下点和左上点[x2[i+1],y2[i+1]],[x2[i+1],0]],   #梯形的右上点和右下点fc='deepskyblue', #填充色ec='dodgerblue', #边缘色alpha=0.5))#3. numpy直接求定积分的API
#利用quad求定积分，给出函数f,积分下限和积分上限[a,b]，返回值为(积分值,最大误差)
area = si.quad(f,a,b)[0]
print(area)mp.legend()
mp.show()

numpy 定积分案例相关推荐

Matlab在线IDE：MATLAB Online介绍与计算定积分案例
目录 1.MATLAB Online介绍功能与特点命令行窗口和编辑器窗口适用场景计费方式使用方法 2.注册登录 3.计算定积分 1.MATLAB Online介绍 MATLAB Online ...
【译】Understanding Batch Normalization with Examples in Numpy and Tensorflow with Interactive Code
因此,今天,我将探讨批量标准化( 批量标准化: 通过 Sergey Ioffe 和 Christian Szegedy 减少内部协变量来加速深度网络培训 ). 但是,为了加强我对数据预处理的理解,我将 ...
人工智能数学基础---定积分1：定积分的概念以及近似计算
一.引言在日常计算中,需要进行一些非线性的计算,如曲边型的面积和变速直线运动的总里程等,由于非线性,导致这些计算不能使用常规的方法来进行.但如果将这些计算涉及的函数在其定义区间上细分成n(n-> ...
python图形界面实践_数据可视化之实践篇——python
一.10个可视化例子 import matplotlib.pyplot as plt import seaborn as sns import numpy as np import pandas as ...
jupyter notebook第三章numpy库的一些案例分析加相关函数的解析
NumPy是在1995年诞生的Python库Numeric的基础上建立起来的,但真正促使NumPy的发行的是Python的SciPy库.但SciPy中并没有合适的类似于Numeric中的对于基础数据对 ...
【ML】numpy ravel、c_、r_函数使用说明（结合实际案例）
[ML]numpy ravel.c_.r_函数使用说明(结合实际案例) 基本使用 ravel()函数用法 c_函数用法 r_函数用法实践基本使用 ravel()函数用法 ravel函数将一个多维数 ...
波士顿房价预测python代码_使用Python和Numpy构建神经网络模型——波士顿房价预测案例...
原标题:使用Python和Numpy构建神经网络模型--波士顿房价预测案例
numpy 获取特定数值的索引_入门Numpy、Pandas及案例实践
用Python进行数据分析,需要用到数据分析相关的包,如本文中用到的pandas和numpy包. 本文目录: 用python分析一维数据用python分析二维数据案例:销售数据分析一.一维数据分 ...
科学计算库NumPy案例：酒鬼漫步
通过前面对NumPy的学习,相信大家一定对NumPy这个科学计算包有了一定的了解,接下来,本节将通过酒鬼漫步的案例来介绍如何运用NumPy随机数模块与数据处理. 下面先为大家描述一下场景,在一片空旷的 ...

numpy 定积分案例

案例

numpy 定积分案例相关推荐

最新文章

热门文章