前言

  这是一篇数值计算算法的开胃菜,主要介绍切比雪夫求积分法的使用;如果你遇到编程和数学的问题,强烈建议找到你身边的算法工程师或者计算机、电子、通信、数学专业的朋友,往往能够事半功倍。

求解Gamma函数

Gamma函数对特殊值

  以下资料来自于维基百科

  切比雪夫多项式展开

通过求解方程组的方式解出c[i]的值,即可通过c[I]求解出Gamma函数的值

#写于2019年3月13日
import numpy as npGam = np.array([2.363271801207,-3.544907701811,1.772453850906,1,\0.886226925453,1,1.329340388179,2,3.323350970448,6])
a = [-1.5,-0.5,0.5,1,1.5,2,2.5,3,3.5,4]
n = len(a)
T = np.ones((n,n))
for i in range(1,n):for j in range(n):if 1 == i:T[j][i] = a[j]else:T[j][i] = 2*a[j]*T[j][i-1]-T[j][i-2]
c = np.linalg.solve(T,Gam)
print(c)
x= 1
T_n = np.zeros((n))
T_n[0] = 1
T_n[1] = x
for i in range(2,n):T_n[i] = 2*x*T_n[i-1]-T_n[i-2]
p = 0
for i in range(n):s += c[i]*T_n[i]
print(p)

求解得到的相应c[i]的值
[ 0.24759993 6.92093386 -33.03116357 88.41245909 -132.01201401
125.18957712 -77.98534069 31.11558493 -7.24600984 0.75164499]

关于作者

  我是cclplus,一名算法萌新,请多多指教。
  如果你有任何疑问,可以联系作者cclplus
  如果你想要更深入的学习Python,你可以加入QQ群:916372346

切比雪夫求积分法(附全过程代码)相关推荐

  1. 这篇Nature子刊文章的蛋白组学数据PCA分析竟花费了我两天时间来重现|附全过程代码...

    复现PCA原图之蛋白组学数据 NGS系列文章包括NGS基础.转录组分析 (Nature重磅综述|关于RNA-seq你想知道的全在这).ChIP-seq分析 (ChIP-seq基本分析流程).单细胞测序 ...

  2. 机器学习——神经网络实验(附全过程代码)

    实验名称:实验二.神经网络 一.实验目的 (1)掌握神经网络的 BP 算法原理与实现方法: (2)神经网络的构建.训练和测试方法. 二.实验内容 (1)BP算法 使用 Python 语言编程实现标准 ...

  3. 翻车实录之Nature Medicine新冠单细胞文献|附全代码

    前言 NGS系列文章包括NGS基础.转录组分析 (Nature重磅综述|关于RNA-seq你想知道的全在这).ChIP-seq分析 (ChIP-seq基本分析流程).单细胞测序分析 (重磅综述:三万字 ...

  4. 利用OpenCV、Python和Ubidots构建行人计数器程序(附完整代码)

    作者 | Jose Garcia 译者 | 吴振东 校对 | 张一豪.林亦霖,编辑 | 于腾凯 来源 | 数据派(ID:datapi) 导读:本文将利用OpenCV,Python和Ubidots来编写 ...

  5. yolov3网络结构图_目标检测——YOLO V3简介及代码注释(附github代码——已跑通)...

    GitHub: liuyuemaicha/PyTorch-YOLOv3​github.com 注:该代码fork自eriklindernoren/PyTorch-YOLOv3,该代码相比master分 ...

  6. php 3d animation,css3D+动画的例子(附完整代码)

    本篇文章给大家带来的内容是关于css3D+动画的例子(附完整代码),有一定的参考价值,有需要的朋友可以参考一下,希望对你有所帮助. 前言 最近玩了玩用css来构建3D效果,写了几个demo,所以博客总 ...

  7. Linux创建线程读取双口数据,linux环境下读写一次双口ram尽然要十几个毫秒。(附驱动代码)...

    linux环境下读写一次双口ram尽然要十几个毫秒.(附驱动代码) 我用的双口ram是IDT70V28,手册上说的读写时间应该是几个纳秒,我写了个linux驱动,然后做测试,发现读写一次的时间尽然是十 ...

  8. 从零开始的AI·决策树原来这么好理解(附实例代码)

    从零开始的AI系列 从零开始的AI·机器学习の基本概念 从零开始的AI·吃透kNN算法,学完我悟了(附实例代码) 从零开始的AI·朴素贝叶斯?拿来吧你(附实例代码) 文章目录 从零开始的AI系列 前言 ...

  9. 从零开始的AI·朴素贝叶斯?拿来吧你(附实例代码)

    从零开始的AI系列 从零开始的AI·机器学习の基本概念 从零开始的AI·吃透kNN算法,学完我悟了(附实例代码) 从零开始的AI·决策树原来这么好理解(附实例代码) 文章目录 从零开始的AI系列 前言 ...

最新文章

  1. Python 核心设计理念27个问题及解答
  2. Office2019完美配置mathtype7.4
  3. Qt Creator将纹理附加到材质
  4. python中的魔法属性和方法
  5. ASP.NET中高级程序员 面试题
  6. mysql error 1045 的解决方法
  7. 第一章 概率论的基本概念
  8. Qt4_Laying Out Widgets Widgets
  9. 【热门技术】三种SEO方式
  10. PHP+jQuery实现翻板抽奖
  11. Archlinux 简明安装指南
  12. 大家都在发圣诞树,我偏偏要发一个圣诞小游戏给大家玩【内附源码】
  13. Unit01: Servlet基础 、 HTTP协议
  14. 【转】Windows XP打印共享四大问题及解决方法
  15. (转载)Stackoverflow让我们变懒了?
  16. 尝试Adam代替梯度下降
  17. Atitit.软件gui按钮and面板---通讯子系统(区)-- github 的使用....
  18. Footprint:Aave VS Compound之战,谁更能脱颖而出?
  19. 最网最全python框架--scrapy(体系学习,爬取全站校花图片),学完显著提高爬虫能力(附源代码),突破各种反爬
  20. URP教务系统自动登录

热门文章

  1. python实现二分查找(对新手友好,内容通俗易懂)
  2. ActiveMQ原理分析
  3. 服务网格 ASM 年终总结:最终用户如何使用服务网格?
  4. 解决Microsoft SQL Server 2005 安装失败方法
  5. Android和wp7,Android 4.0与WP7孰强孰弱?
  6. 基于Android Studio的Android内存泄漏检测方法
  7. No such file or directory错误
  8. Excel VBA 高级编程-来自直男的Excel表白
  9. [原创]双硬盘四系统安装全功略
  10. IE浏览器设置兼容模式不生效或者只是当次生效的解决办法