精华内容
下载资源
问答
  • 《复变函数(第五版)/面向21世纪课程教材》内容包括:复数及复平面、复变函数、复变函数的积分、级数、留数、保形映射、解析开拓以及调和函数共八章,其中除单性定理外,均属于复变函数课程的一般内容。...
  • 本书是俄罗斯综合大学和高等技术学校使用的复变函数论教材。它基于前苏联著名数学家、科学院院士拉夫连季耶夫讲稿,由沙巴特补充整理,并经过次修订,使内容更为合理,应用实例更为丰富,已成为该领域一本经典...
  • 第6章 多值函数 第二部分 积分 第7章 积分 第8章 复合闭路积分 第三部分 复级数 第9章 幂函数 第10章 Taylor级数和Laurent级数 第四部分 留数 第11章 留数定理(I) 第12章 留数定理(II) 第13章 留数应用(I...
  • 类题或自我检查题第三章 复变函数的积分 I.重点、要求与例题 §1.复积分的概念及其简单性质(例3.1.1—3.1.11) §2.柯西积分定理(例3.2.1—3.2.9) §3.柯西积分公式及其推论(例3.3.1—3.3.16) §4.解析函数与...
  • 1.1 关于复变函数的若干问答 1.2 函数可导的充分必要条件 1.3 Cauchy定理与Cauchy积分公式 第二章 无穷级数 2.1 无穷级数的收敛性 2.2 幂级数的收敛半径 2.3 无穷级数的Cesaro和与Abel和 2.4 解析函数的幂级数展开 ...
  • 用程序来求积分方法有很,这篇文章主要是有关牛顿-科特斯公式。...这样会将原来求积分问题带到另一个求积分问题:如何求n次多项式积分,而且当次数高时,会出现龙悲歌现象,误差反而可能会增大,并且高次...

    用程序来求积分的方法有很多,这篇文章主要是有关牛顿-科特斯公式。

    学过插值算法的同学最容易想到的就是用插值函数代替被积分函数来求积分,但实际上在大部分场景下这是行不通的。

    插值函数一般是一个不超过n次的多项式,如果用插值函数来求积分的话,就会引进高次多项式求积分的问题。这样会将原来的求积分问题带到另一个求积分问题:如何求n次多项式的积分,而且当次数变高时,会出现龙悲歌现象,误差反而可能会增大,并且高次的插值求积公式有可能会变得不稳定:详细原因不赘述。

    牛顿-科特斯公式解决这一问题的办法是将大的插值区间分为一堆小的插值区间,使得多项式的次数不会太高。然后通过引入参数函数

    104a43c97fd8b20d0dcfc184a7ec33e1.jpg

    将带有幂的项的取值范围固定在一个固定范围内,这样一来就将多项式带有幂的部分的求积变为一个固定的常数,只需手工算出来即可。这个常数可以直接带入多项式求积函数。

    上式中x的求积分区间为[a, b],h = (b - a)/n, 这样一来积分区间变为[0, n],需要注意的是从这个公式可以看出一个大的区间被分为n个等长的小区间。 这一部分具体请参见任意一本有关数值计算的书!

    n是一个事先确定好的值。

    又因为一个大的插值区间需要被分为等长的多个小区间,并在这些小区间上分别进行插值和积分,因此此时的牛顿-科特斯公式被称为:复化牛顿-科特斯公式。

    并且对于n的不同取值牛顿-科特斯有不同的名称: 当n=1时,叫做复化梯形公式,复化梯形公式也就是将每一个小区间都看为一个梯形(高为h,上底为f(t), 下底为f(t+1))。这与积分的本质:无限分隔 相同。

    当n=2时,复化牛顿-科特斯公式被称为复化辛普森公式(非美国法律界著名的那个辛普森)。

    我这篇文章实现的是复化梯形公式:

    4e54624ec7440b2e6c2e7df0aa5342d9.jpg

    首先写一个函数求节点函数值求和那部分:

    """

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

    相当的简单

    试验:

    当把大区间分为两个小区间时:

    f55c8c8c3a71c82ab65b406609ba4988.jpg

    分为20个小区间时:

    0031a4a55acedcc3b238942aef18c601.jpg

    求的积分值就是这些彩色的梯形面积之和。

    测试代码:

    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进行数值计算就是小编分享给大家的全部内容了,希望能给大家一个参考

    您可能感兴趣的文章:

    展开全文
  • 《简明复分析(中国科学技术大学精品教材)》试图用近代数学观点和方法处理复变函数内容,并强调数学统一性。例如,用微分几何初步知识,对Picard大、小定理给出简洁证明;强调变换群概念,利用Pompeiu公式...
  • 9 与多复变函数联系 10 复函数逼近 11 与插值问题的联系 12 Caratheodory和Kobayashi度量及其在复分析中应用 1 序言 2 单化定理 3 源自于Schwarz引理和Schwarz-Pick引理推动 4 关于小林度量基本事实 5...
  • 范例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抖string类来处理文本输入和 输出。最后,...
  • 章讨论其他形式,包括数组(存储个同类型的值)、结构(存储个不同类型的值)、指针(标识内存位置) 。读者还将学习如何创建和存储文本字符串及如何使用C—风格字符数组和C抖string类来处理文本输入和 输出。最后,...
  • 章讨论其他形式,包括数组(存储个同类型的值)、结构(存储个不同类型的值)、指针(标识内存位置) 。读者还将学习如何创建和存储文本字符串及如何使用C—风格字符数组和C抖string类来处理文本输入和 输出。最后,...
  • 用算法程序集(C语言描述)(第五版)+源代码 第1章 多项式计算 1.1 一维多项式求 1.2 一维多项式组求 1.3 二维多项式求 1.4 系数多项式求 1.5 多项式相乘 1.6 系数多项式相乘 1.7 多项式相除 1.8 ...
  • 第14章 特殊函数的计算 14.1 伽马函数 14.2 不完全伽马函数 14.3 误差函数 14.4 第一类整数阶贝塞尔函数 14.5 第二类整数阶贝塞尔函数 14.6 变形第一类整数阶贝塞尔函数 14.7 变形第二类整数阶贝塞尔函数 14.8 不...
  • 第14章 特殊函数的计算 14.1 伽马函数 14.2 不完全伽马函数 14.3 误差函数 14.4 第一类整数阶贝塞尔函数 14.5 第二类整数阶贝塞尔函数 14.6 变形第一类整数阶贝塞尔函数 14.7 变形第二类整数阶贝塞尔函数 14.8 不...
  • 第14章 特殊函数的计算 14.1 伽马函数 14.2 不完全伽马函数 14.3 误差函数 14.4 第一类整数阶贝塞尔函数 14.5 第二类整数阶贝塞尔函数 14.6 变形第一类整数阶贝塞尔函数 14.7 变形第二类整数阶贝塞尔函数 14.8 不...
  • 但是,算法仅仅给出了一个笼统实现策略,对于具体数据类型、运算操作符及优先级情况、怎样支持函数调用、自定义内存变量如何实现等具体问题并无现成解决方案,需要开发人员自己实际设计相应算法解决,而这也...
  • 中文版Excel.2007高级VBA编程宝典 1/2

    热门讨论 2012-04-06 16:00:16
     10.11.4 了解更有关API 函数的信息  第11章 VBA编程示例和技巧  11.1 处理单元格区域  11.1.1 复制单元格区域  11.1.2 移动单元格区域  11.1.3 复制大小可的单元格区域  11.1.4 选中或者识别各种类型的...
  • 第14章 特殊函数的计算 14.1 伽马函数 14.2 不完全伽马函数 14.3 误差函数 14.4 第一类整数阶贝塞尔函数 14.5 第二类整数阶贝塞尔函数 14.6 变形第一类整数阶贝塞尔函数 14.7 变形第二类整数阶贝塞尔函数 14.8 不...
  • <br>【目录信息】 第1章 多项式计算 1.1 一维多项式求 1.2 一维多项式组求 1.3 二维多项式求 1.4 系数多项式求 1.5 多项式相乘 1.6 系数多项式相乘 1.7 多项式相除 ...
  • MATLAB数值分析与应用

    2017-05-13 21:43:13
    第二部分系统地介绍了符号计算在微积分和复变函数两门大学数学基础课程中应用,以及线性方程组、非线性方程与最优化方法、特征与特征向量、插值与函数逼近、估计方法和数据拟合、积分计算、常微分方程等数值方法...
  • 复变函数创立 实变函数论 泛函分析 函数逼近论 傅立叶分析 非标准分析 几何学 欧几里得几何学 非欧几里得几何 解析几何 二次曲线和曲面 三角学 三角函数 反三角函数 仿射几何学 射影几何学 微分...
  • Visual C++ 2008入门经典(中文版).pdf【第一部分】

    千次下载 热门讨论 2011-06-20 09:31:46
    4.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 声明...
  • 4.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 声明...
  • 4.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:15
    IvorHorton还著有关于C、C++和Java的多部入门级好书,如《C语言入门经典(第4版)》和《C++入门经典(第3版)》。 目录 封面 -19 前言 -18 目录 -10 第1章 使用Visual C++ 2008编程 1 1.1 .NET Framework 1 1.2 CLR 2...
  • IvorHorton还著有关于C、C++和Java的多部入门级好书,如《C语言入门经典(第4版)》和《C++入门经典(第3版)》。 目录 封面 -19 前言 -18 目录 -10 第1章 使用Visual C++ 2008编程 1 1.1 .NET Framework 1 1.2 CLR 2...

空空如也

空空如也

1 2 3 4
收藏数 70
精华内容 28
关键字:

复变函数的多值函数问题