切比雪夫求积分法(附全过程代码)
前言
这是一篇数值计算算法的开胃菜,主要介绍切比雪夫求积分法的使用;如果你遇到编程和数学的问题,强烈建议找到你身边的算法工程师或者计算机、电子、通信、数学专业的朋友,往往能够事半功倍。
求解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
切比雪夫求积分法(附全过程代码)相关推荐
- 这篇Nature子刊文章的蛋白组学数据PCA分析竟花费了我两天时间来重现|附全过程代码...
复现PCA原图之蛋白组学数据 NGS系列文章包括NGS基础.转录组分析 (Nature重磅综述|关于RNA-seq你想知道的全在这).ChIP-seq分析 (ChIP-seq基本分析流程).单细胞测序 ...
- 机器学习——神经网络实验(附全过程代码)
实验名称:实验二.神经网络 一.实验目的 (1)掌握神经网络的 BP 算法原理与实现方法: (2)神经网络的构建.训练和测试方法. 二.实验内容 (1)BP算法 使用 Python 语言编程实现标准 ...
- 翻车实录之Nature Medicine新冠单细胞文献|附全代码
前言 NGS系列文章包括NGS基础.转录组分析 (Nature重磅综述|关于RNA-seq你想知道的全在这).ChIP-seq分析 (ChIP-seq基本分析流程).单细胞测序分析 (重磅综述:三万字 ...
- 利用OpenCV、Python和Ubidots构建行人计数器程序(附完整代码)
作者 | Jose Garcia 译者 | 吴振东 校对 | 张一豪.林亦霖,编辑 | 于腾凯 来源 | 数据派(ID:datapi) 导读:本文将利用OpenCV,Python和Ubidots来编写 ...
- yolov3网络结构图_目标检测——YOLO V3简介及代码注释(附github代码——已跑通)...
GitHub: liuyuemaicha/PyTorch-YOLOv3github.com 注:该代码fork自eriklindernoren/PyTorch-YOLOv3,该代码相比master分 ...
- php 3d animation,css3D+动画的例子(附完整代码)
本篇文章给大家带来的内容是关于css3D+动画的例子(附完整代码),有一定的参考价值,有需要的朋友可以参考一下,希望对你有所帮助. 前言 最近玩了玩用css来构建3D效果,写了几个demo,所以博客总 ...
- Linux创建线程读取双口数据,linux环境下读写一次双口ram尽然要十几个毫秒。(附驱动代码)...
linux环境下读写一次双口ram尽然要十几个毫秒.(附驱动代码) 我用的双口ram是IDT70V28,手册上说的读写时间应该是几个纳秒,我写了个linux驱动,然后做测试,发现读写一次的时间尽然是十 ...
- 从零开始的AI·决策树原来这么好理解(附实例代码)
从零开始的AI系列 从零开始的AI·机器学习の基本概念 从零开始的AI·吃透kNN算法,学完我悟了(附实例代码) 从零开始的AI·朴素贝叶斯?拿来吧你(附实例代码) 文章目录 从零开始的AI系列 前言 ...
- 从零开始的AI·朴素贝叶斯?拿来吧你(附实例代码)
从零开始的AI系列 从零开始的AI·机器学习の基本概念 从零开始的AI·吃透kNN算法,学完我悟了(附实例代码) 从零开始的AI·决策树原来这么好理解(附实例代码) 文章目录 从零开始的AI系列 前言 ...
最新文章
- Python 核心设计理念27个问题及解答
- Office2019完美配置mathtype7.4
- Qt Creator将纹理附加到材质
- python中的魔法属性和方法
- ASP.NET中高级程序员 面试题
- mysql error 1045 的解决方法
- 第一章 概率论的基本概念
- Qt4_Laying Out Widgets Widgets
- 【热门技术】三种SEO方式
- PHP+jQuery实现翻板抽奖
- Archlinux 简明安装指南
- 大家都在发圣诞树,我偏偏要发一个圣诞小游戏给大家玩【内附源码】
- Unit01: Servlet基础 、 HTTP协议
- 【转】Windows XP打印共享四大问题及解决方法
- (转载)Stackoverflow让我们变懒了?
- 尝试Adam代替梯度下降
- Atitit.软件gui按钮and面板---通讯子系统(区)-- github 的使用....
- Footprint:Aave VS Compound之战,谁更能脱颖而出?
- 最网最全python框架--scrapy(体系学习,爬取全站校花图片),学完显著提高爬虫能力(附源代码),突破各种反爬
- URP教务系统自动登录