Python日记（6）——数值积分

Python日记（6）——数值积分
每天做一个Python小练习，顺便记录一些小技巧。
在做化学实验中，一般测得的数据都是离散的，所以编写了一些程序来计算离散数据的积分。
首先，粗讲一下已知函数式的数值积分。
以下面的式子举例

这里需要用到numpy
简单地导入以下

import numpy as np

首先，表达出要积分的函数表达式

y = np.cos(2*np.pi*x)*np.exp(-x)+1.2

原始方法
就是把函数积分化成很多且宽很小的矩形，再求所有矩形的面积总和，便可以得到积分的值。这里把x从0.7到4分解成1000分，利用linspace方法

x = np.linspace(0, 6, 1000)

然后dx作宽y作高，计算并总加所有矩形的面积

dx = x[1] - x[0]
fArea = np.sum(y*dx)
print(f'Integral area:{fArea}')

所有代码如下：

import numpy as np

x = np.linspace(0.7, 4.0, 1000)
y = np.cos(2*np.pi*x)