-
复变函数 第五版 余家荣 高清 带详细目录
2018-04-14 07:26:48《复变函数(第五版)/面向21世纪课程教材》内容包括:复数及复平面、复变函数、复变函数的积分、级数、留数、保形映射、解析开拓以及调和函数共八章,其中除单值性定理外,均属于复变函数课程的一般内容。... -
复变函数论方法 第6版-拉夫连季耶夫 沙巴特pdf
2019-01-03 14:21:06本书是俄罗斯综合大学和高等技术学校使用的复变函数论教材。它基于前苏联著名数学家、科学院院士拉夫连季耶夫的讲稿,由沙巴特补充整理,并经过多次修订,使内容更为合理,应用实例更为丰富,已成为该领域一本经典... -
《复变函数札记》作者: 梁昌洪 出版年: 2011年
2019-05-27 22:30:45第6章 多值函数 第二部分 复积分 第7章 复积分 第8章 复合闭路积分 第三部分 复级数 第9章 复幂函数 第10章 Taylor级数和Laurent级数 第四部分 留数 第11章 留数定理(I) 第12章 留数定理(II) 第13章 留数应用(I... -
《复变函数学习指导书》作者: 钟玉泉 编 出版年: 1996年
2019-05-15 18:06:47类题或自我检查题第三章 复变函数的积分 I.重点、要求与例题 §1.复积分的概念及其简单性质(例3.1.1—3.1.11) §2.柯西积分定理(例3.2.1—3.2.9) §3.柯西积分公式及其推论(例3.3.1—3.3.16) §4.解析函数与... -
数学物理方法专题复变函数与积分变换 [吴崇试 编] 2013年版
2019-05-07 22:25:561.1 关于复变函数的若干问答 1.2 函数可导的充分必要条件 1.3 Cauchy定理与Cauchy积分公式 第二章 无穷级数 2.1 无穷级数的收敛性 2.2 幂级数的收敛半径 2.3 无穷级数的Cesaro和与Abel和 2.4 解析函数的幂级数展开 ... -
python梯形法计算积分值_复化梯形求积分实例――用Python进行数值计算
2021-02-10 16:52:42用程序来求积分的方法有很多,这篇文章主要是有关牛顿-科特斯公式。...这样会将原来的求积分问题带到另一个求积分问题:如何求n次多项式的积分,而且当次数变高时,会出现龙悲歌现象,误差反而可能会增大,并且高次...用程序来求积分的方法有很多,这篇文章主要是有关牛顿-科特斯公式。
学过插值算法的同学最容易想到的就是用插值函数代替被积分函数来求积分,但实际上在大部分场景下这是行不通的。
插值函数一般是一个不超过n次的多项式,如果用插值函数来求积分的话,就会引进高次多项式求积分的问题。这样会将原来的求积分问题带到另一个求积分问题:如何求n次多项式的积分,而且当次数变高时,会出现龙悲歌现象,误差反而可能会增大,并且高次的插值求积公式有可能会变得不稳定:详细原因不赘述。
牛顿-科特斯公式解决这一问题的办法是将大的插值区间分为一堆小的插值区间,使得多项式的次数不会太高。然后通过引入参数函数
将带有幂的项的取值范围固定在一个固定范围内,这样一来就将多项式带有幂的部分的求积变为一个固定的常数,只需手工算出来即可。这个常数可以直接带入多项式求积函数。
上式中x的求积分区间为[a, b],h = (b - a)/n, 这样一来积分区间变为[0, n],需要注意的是从这个公式可以看出一个大的区间被分为n个等长的小区间。 这一部分具体请参见任意一本有关数值计算的书!
n是一个事先确定好的值。
又因为一个大的插值区间需要被分为等长的多个小区间,并在这些小区间上分别进行插值和积分,因此此时的牛顿-科特斯公式被称为:复化牛顿-科特斯公式。
并且对于n的不同取值牛顿-科特斯有不同的名称: 当n=1时,叫做复化梯形公式,复化梯形公式也就是将每一个小区间都看为一个梯形(高为h,上底为f(t), 下底为f(t+1))。这与积分的本质:无限分隔 相同。
当n=2时,复化牛顿-科特斯公式被称为复化辛普森公式(非美国法律界著名的那个辛普森)。
我这篇文章实现的是复化梯形公式:
首先写一个函数求节点函数值求和那部分:
"""
@brief: 求和 ∑f(xk) : xk表示等距节点的第k个节点,不包括端点
xk = a + kh (k = 0, 1, 2, ...)
积分区间为[a, b]
@param: xk 积分区间的等分点x坐标集合(不包括端点)
@param: func 求积函数
@return: 返回值为集合的和
"""
def sum_fun_xk(xk, func):
return sum([func(each) for each in xk])
然后就可以写整个求积分函数了:
"""
@brief: 求func积分 :
@param: a 积分区间左端点
@param: b 积分区间右端点
@param: n 积分分为n等份(复化梯形求积分要求)
@param: func 求积函数
@return: 积分值
"""
def integral(a, b, n, func):
h = (b - a)/float(n)
xk = [a + i*h for i in range(1, n)]
return h/2 * (func(a) + 2 * sum_fun_xk(xk, func) + func(b))
相当的简单
试验:
当把大区间分为两个小区间时:
分为20个小区间时:
求的积分值就是这些彩色的梯形面积之和。
测试代码:
if __name__ == "__main__":
func = lambda x: x**2
a, b = 2, 8
n = 20
print integral(a, b, n, func)
''' 画图 '''
import matplotlib.pyplot as plt
plt.figure("play")
ax1 = plt.subplot(111)
plt.sca(ax1)
tmpx = [2 + float(8-2) /50 * each for each in range(50+1)]
plt.plot(tmpx, [func(each) for each in tmpx], linestyle = '-', color='black')
for rang in range(n):
tmpx = [a + float(8-2)/n * rang, a + float(8-2)/n * rang, a + float(8-2)/n * (rang+1), a + float(8-2)/n * (rang+1)]
tmpy = [0, func(tmpx[1]), func(tmpx[2]), 0]
c = ['r', 'y', 'b', 'g']
plt.fill(tmpx, tmpy, color=c[rang%4])
plt.grid(True)
plt.show()
注意上面代码中的n并不是上文开篇提到的公式中的n,开篇提到的n是指将每一个具体的插值区间(也就是小区间)等距插n个节点,复化梯形公式的n是固定的为1.
而代码中的n指将大区间分为n个小区间。
以上这篇复化梯形求积分实例――用Python进行数值计算就是小编分享给大家的全部内容了,希望能给大家一个参考
您可能感兴趣的文章:
-
《简明复分析(第2版)》作者: 龚昇 出版年: 2009年
2019-05-17 19:02:27《简明复分析(中国科学技术大学精品教材)》试图用近代数学的观点和方法处理复变函数内容,并强调数学的统一性。例如,用微分几何的初步知识,对Picard大、小定理给出简洁的证明;强调变换群的概念,利用Pompeiu公式... -
《数学中的小问题大定理(第2辑):许瓦兹引理 从一道加利福尼亚大学伯克利分校数学系博士生试题谈起》作者: ...
2019-05-26 20:32:269 与多复变函数论的联系 10 复函数的逼近 11 与插值问题的联系 12 Caratheodory和Kobayashi度量及其在复分析中的应用 1 序言 2 单值化定理 3 源自于Schwarz引理和Schwarz-Pick引理的推动 4 关于小林度量的基本事实 5... -
C开发金典随书源码:含数据结构 数值计算分析 图形图像处理 目录和文件操作 系统调用方面的范例
2013-10-25 13:12:12范例1-52 双向链表元素值的查询 129 ∷相关函数:GetElemP函数 1.3.22 稀疏矩阵的建立 136 范例1-53 稀疏矩阵的建立 136 ∷相关函数:Create函数 1.3.23 稀疏矩阵的删除 138 范例1-54 稀疏矩阵的删除 138 ∷... -
C语言通用范例开发金典.part2.rar
2012-08-31 14:18:18范例1-52 双向链表元素值的查询 129 ∷相关函数:GetElemP函数 1.3.22 稀疏矩阵的建立 136 范例1-53 稀疏矩阵的建立 136 ∷相关函数:Create函数 1.3.23 稀疏矩阵的删除 138 范例1-54 稀疏矩阵的删除 138 ∷... -
C 开发金典
2013-06-20 16:20:03范例1-52 双向链表元素值的查询 129 ∷相关函数:GetElemP函数 1.3.22 稀疏矩阵的建立 136 范例1-53 稀疏矩阵的建立 136 ∷相关函数:Create函数 1.3.23 稀疏矩阵的删除 138 范例1-54 稀疏矩阵的删除 138 ∷... -
最权威的C++教程_C++_Primer_Plus中文第五版+C++_Primer中文第四版(都含源码+习题)(共4分卷)分卷1
2010-06-23 17:33:55章讨论其他形式,包括数组(存储多个同类型的值)、结构(存储多个不同类型的值)、指针(标识内存位置) 。读者还将学习如何创建和存储文本字符串及如何使用C—风格字符数组和C抖string类来处理文本输入和 输出。最后,... -
最权威的C++教程_C++_Primer_Plus中文第五版+C++_Primer中文第四版(都含源码+习题)(共4分卷)分卷2
2010-06-23 17:47:19章讨论其他形式,包括数组(存储多个同类型的值)、结构(存储多个不同类型的值)、指针(标识内存位置) 。读者还将学习如何创建和存储文本字符串及如何使用C—风格字符数组和C抖string类来处理文本输入和 输出。最后,... -
最权威的C++教程_C++_Primer_Plus中文第五版+C++_Primer中文第四版(都含源码+习题)(共4分卷)分卷3
2010-06-23 18:03:39章讨论其他形式,包括数组(存储多个同类型的值)、结构(存储多个不同类型的值)、指针(标识内存位置) 。读者还将学习如何创建和存储文本字符串及如何使用C—风格字符数组和C抖string类来处理文本输入和 输出。最后,... -
常用算法程序集(C语言及C++描述)第五版(徐士良) 高清版 pdf 文档 第5版
2019-05-25 15:13:00用算法程序集(C语言描述)(第五版)+源代码 第1章 多项式的计算 1.1 一维多项式求值 1.2 一维多项式多组求值 1.3 二维多项式求值 1.4 复系数多项式求值 1.5 多项式相乘 1.6 复系数多项式相乘 1.7 多项式相除 1.8 ... -
徐士良常用算法程序集(C语言描述)第三版(带书签)及程序源码
2017-12-01 15:03:41第14章 特殊函数的计算 14.1 伽马函数 14.2 不完全伽马函数 14.3 误差函数 14.4 第一类整数阶贝塞尔函数 14.5 第二类整数阶贝塞尔函数 14.6 变形第一类整数阶贝塞尔函数 14.7 变形第二类整数阶贝塞尔函数 14.8 不... -
徐士良C常用算法程序集第三版高清电子书+源代码
2018-01-11 17:27:42第14章 特殊函数的计算 14.1 伽马函数 14.2 不完全伽马函数 14.3 误差函数 14.4 第一类整数阶贝塞尔函数 14.5 第二类整数阶贝塞尔函数 14.6 变形第一类整数阶贝塞尔函数 14.7 变形第二类整数阶贝塞尔函数 14.8 不... -
常用算法程序集(C语言描述)(第三版)+完整源代码
2010-12-05 13:50:20第14章 特殊函数的计算 14.1 伽马函数 14.2 不完全伽马函数 14.3 误差函数 14.4 第一类整数阶贝塞尔函数 14.5 第二类整数阶贝塞尔函数 14.6 变形第一类整数阶贝塞尔函数 14.7 变形第二类整数阶贝塞尔函数 14.8 不... -
VC++科学计算器的设计V1_4Beta.rar
2010-03-23 04:54:59但是,算法仅仅给出了一个笼统的实现策略,对于具体的数据类型、运算操作符及优先级情况、怎样支持函数调用、自定义内存变量如何实现等具体问题并无现成的解决方案,需要开发人员自己实际设计相应的算法解决,而这也... -
中文版Excel.2007高级VBA编程宝典 1/2
2012-04-06 16:00:1610.11.4 了解更多有关API 函数的信息 第11章 VBA编程示例和技巧 11.1 处理单元格区域 11.1.1 复制单元格区域 11.1.2 移动单元格区域 11.1.3 复制大小可变的单元格区域 11.1.4 选中或者识别各种类型的... -
常用算法程序集(C语言描述)(第三版)+源代码
2008-11-04 22:14:05第14章 特殊函数的计算 14.1 伽马函数 14.2 不完全伽马函数 14.3 误差函数 14.4 第一类整数阶贝塞尔函数 14.5 第二类整数阶贝塞尔函数 14.6 变形第一类整数阶贝塞尔函数 14.7 变形第二类整数阶贝塞尔函数 14.8 不... -
常用算法程序集(C语言第三版)只要2个积分
2008-01-30 13:55:41<br>【目录信息】 第1章 多项式的计算 1.1 一维多项式求值 1.2 一维多项式多组求值 1.3 二维多项式求值 1.4 复系数多项式求值 1.5 多项式相乘 1.6 复系数多项式相乘 1.7 多项式相除 ... -
MATLAB数值分析与应用
2017-05-13 21:43:13第二部分系统地介绍了符号计算在微积分和复变函数两门大学数学基础课程中的应用,以及线性方程组、非线性方程与最优化方法、特征值与特征向量、插值与函数逼近、估计方法和数据拟合、积分计算、常微分方程等数值方法... -
数与形 [崔玉泉,包芳勋 主编] 2013年版
2019-05-29 23:21:24复变函数论的创立 实变函数论 泛函分析 函数逼近论 傅立叶分析 非标准分析 几何学 欧几里得几何学 非欧几里得几何 解析几何 二次曲线和曲面 三角学 三角函数 反三角函数 仿射几何学 射影几何学 微分... -
Visual C++ 2008入门经典(中文版).pdf【第一部分】
2011-06-20 09:31:464.1 处理多个相同类型的数据值 134 4.1.1 数组 135 4.1.2 声明数组 135 4.1.3 初始化数组 138 4.1.4 字符数组和字符串处理 139 4.1.5 多维数组 142 4.2 间接数据存取 145 4.2.1 指针的概念 145 4.2.2 声明... -
Visual C++ 2008入门经典(中文版).pdf【第三部分】
2011-06-20 09:25:514.1 处理多个相同类型的数据值 134 4.1.1 数组 135 4.1.2 声明数组 135 4.1.3 初始化数组 138 4.1.4 字符数组和字符串处理 139 4.1.5 多维数组 142 4.2 间接数据存取 145 4.2.1 指针的概念 145 4.2.2 声明... -
Visual C++ 2008入门经典(中文版).pdf【第二部分】
2011-06-20 09:29:154.1 处理多个相同类型的数据值 134 4.1.1 数组 135 4.1.2 声明数组 135 4.1.3 初始化数组 138 4.1.4 字符数组和字符串处理 139 4.1.5 多维数组 142 4.2 间接数据存取 145 4.2.1 指针的概念 145 4.2.2 声明... -
Excel百宝箱8.0
2011-06-07 21:32:17【筛选唯一值】【建工作表目录】【批量新建复制工作表】【批量加解密】【工作表批量命名】【破解工作表/簿密码】 【批量命名文件】【可还原的合并】【合并列中相同值】【取消合并还原数据】【合并区域自动换行】... -
Visual C++ 2008入门经典--详细书签版
2013-02-02 16:07:15IvorHorton还著有关于C、C++和Java的多部入门级好书,如《C语言入门经典(第4版)》和《C++入门经典(第3版)》。 目录 封面 -19 前言 -18 目录 -10 第1章 使用Visual C++ 2008编程 1 1.1 .NET Framework 1 1.2 CLR 2... -
Visual C++ 2008入门经典--源代码及课后练习答案
2013-02-02 16:13:25IvorHorton还著有关于C、C++和Java的多部入门级好书,如《C语言入门经典(第4版)》和《C++入门经典(第3版)》。 目录 封面 -19 前言 -18 目录 -10 第1章 使用Visual C++ 2008编程 1 1.1 .NET Framework 1 1.2 CLR 2...
-
xxljob源码分析
-
Reply-tools-analysis-2020:计算参与度和采用度指标以评估答复工具的影响-源码
-
毛概课件 第七章改革和对外开放.ppt
-
linux 查看docker 容器创建时命令——筑梦之路
-
laravel-源码
-
kds:多平台Kotlin 1.3的数据结构库-源码
-
Exercise.ppt
-
MacOS基础教程——最全的Mac键盘快捷键分享
-
朱老师鸿蒙系列课程第1期-2鸿蒙系统Harmonyos源码架构分析
-
MySQL 查询与高级查询(多表、嵌套和正则表达式)
-
Oracle_11g_Linux到Linux_DataGuard部署
-
MySQL 高可用工具 DRBD 实战部署详解
-
Java中的ThreadLocal详解
-
MPAP_c579ebd1135a4f42bf1140fe0a1c1b77_001.provxml
-
1902_20191611321_韩玉满.zip
-
linux常用基本命令
-
Chap3 Theories of FPI.ppt
-
MySQL 数据类型和运算符
-
it_blog-源码
-
华为1+X——网络系统建设与运维(高级)