python用复化梯形公式积分_复化梯形求积分——用Python进行数值计算

最新推荐文章于 2024-04-25 11:20:08 发布
最新推荐文章于 2024-04-25 11:20:08 发布
阅读量2.6k 收藏 1
点赞数
文章标签： python用复化梯形公式积分
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_32770541/article/details/112827203
版权
本文介绍了使用Python实现复化梯形公式进行数值积分的方法，通过将大区间分解为多个小区间，避免了高次多项式求积分带来的问题。文章提供了详细的代码示例，展示了如何利用函数计算节点值的和以及整个积分过程。此外，还提到了将大区间分为两个或更多小区间时，积分值相当于彩色梯形面积之和，并提供了相应的绘图演示。
摘要由CSDN通过智能技术生成

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

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

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

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

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

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

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

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

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

最低0.47元/天 解锁文章
素心连城
关注
python用复化梯形公式积分_荐数值积分 python代码实现
weixin_39693662的博客
12-18 1242
老规矩，数学原理什么的就不写了。直接贴代码和实例演示，以下代码基于python和numpy。在这里，我将用代码实现复化梯形算法、复化 Simpson 算法、Romberg 积分算法和三点 Gauss-Legendre求积算法。往期博客：复化梯形算法复化梯形公式如下图所示：首先import numpy as np定义函数以下便是我定义的函数：def tx_fh(x0,f,n):a=x0[0]b=x0...
复化梯形积分和龙格积分（python）
weixin_56273009的博客
10-10 1292
通过设计、编制、调试2~3个数值积分与微分算法的程序，加深对其数值计算方法及有关的基础理论知识的理解。用编程语言实现复化梯形积分、Romberg积分的程序。积分区间为0到1 结果为1.718281828459045。1.718281828459045 # 真值。R为： 1.7182818287945303。Romberg积分算法求积分。的自动控制误差算法求积分。
python 实现复合梯度公式及复合辛普森公式
(｡◕‿◕｡)
11-02 3350
1.复合梯形公式 2.复合辛普森公式 python实现 import math def fun(x): return math.sin(x)/(x+1e-16) #加上1e-16避免除零错误 # 复合梯度 def tx(a,b,n): h=(b-a)/n fxi=0 xi=a for i in range(1,n): xi=xi+h...
复化梯形求积分实例——用Python进行数值计算
09-18
今天小编就为大家分享一篇复化梯形求积分实例——用Python进行数值计算，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
复化梯形公式和复化Simpson公式的python实现
☆下山☆的博客
10-27 4497
复化梯形公式和复化Simpson公式的python实现
基于复合梯形公式和复合辛普森求积公式计算积分在python中的实现.txt
05-05
基于复合梯形公式和复合辛普森求积公式计算积分在python中的实现.txt
复化梯形求积分——用Python进行数值计算
anyou7947的博客
12-10 4325
　　用程序来求积分的方法有很多，这篇文章主要是有关牛顿-科特斯公式。 　　学过插值算法的同学最容易想到的就是用插值函数代替被积分函数来求积分，但实际上在大部分场景下这是行不通的。 　　插值函数一般是一个不超过n次的多项式，如果用插值函数来求积分的话，就会引进高次多项式求积分的问题。这样会将原来的求积分问题带到另一个求积分问题：如何求n次多项式的积分，而且当次数变高时，会出现龙悲歌现象...
数值分析实验 实验2-1 复化梯形公式 python3实现
qq_35857125的博客
09-28 1647
题目： 设计复化梯形公式求积算法 代码： # 要积的函数，字符串形式 fx_str = str(input("请输入要积的函数，格式为 x*x+1 （即变量用x表示，乘号用*表示）：")) # 求积的区间[a, b] a, b = list(map(float, input("请输入求积区间[a, b]，格式为 a b ：").split())) # 区间划分的份数 n = int(input("...
复化梯形公式、复化Simpon公式、Romberg算法（python）
weixin_46039719的博客
11-16 2410
import numpy as np PI25DT=3.141592653589793238462643 def funEval(x): fx = 4.0/(1+x**2) return fx def fhtan(x0,f): #复化梯形公式 a=x0[0] b=x0[-1] fx = f(x0) y=2*np.sum(fx)-fx[0]-fx[-1] tn=((b-a)/(fx.shape[0]-1))*y/2 return tn.
C++程序用复化梯形公式计算积分
12-14
C++程序用复化梯形公式计算积分sinx从0到1的值
用梯形公式计算定积分
05-14
分别输入所求定积分的上下限，算法实现梯形公式计算定积分的近似值。
复化梯形公式和复化辛普森公式
11-24
（① 实验目的 1）观察复化梯形公式和复化辛普森公式随区间数n增加各自误差的减少规律；研究广义积分的数值计算如何将其转化为普通积分，再由已有数值积分方法进行计算； 2）利用复化梯形公式和复化辛普森公式计算定积分，编程实现。
复合梯形公式与复合辛普森公式MATLAB
03-18
采用复合梯形公式与复合辛普森公式，计算 sin(x)/x 在[0, 1]范围内的积分。采样点数 目为 5、9、17、33。
定步长复化梯形求解积分方程
DragonLord的学习园地
06-14 2001
/*****定步长复化梯形求解积分方程***                                                     n-1                                   n-1  定步长复化梯形公式： Tn=   ∑  h/2〔 f(a) + f(b) + 2 ∑ f(x(k)) 〕                              
变步长复化梯形法的Python实现
Jaryn62的博客
06-09 1282
变步长复化梯形法的Python
[计算机数值分析]复化梯形的递推公式的变步长算法求积分
Spring-_-Bear 的 CSDN 博客
12-14 4153
变步长求积法的主要思想是利用若干小梯形的面积代替原方程的积分，当精度达不到要求时，可以通过增加点数对已有的区间再次划分，达到所需精度时即可；其中由于新的式子中有原来n点中的部分项，故可以省略一些计算，符合了计算机计算存储的思想。
复化积分公式
huhuandk的博客
04-10 1354
对于积分： 只要找到被积公式的原函数F(x)，利用牛顿莱普利兹公式有： 但是，实际使用这种求积分的方法往往是有困难的，因为大量的被积函数的原函数是不能用初等函数表示的；另外，当f(x)是由测量或数值计算给出的一张数据表时，牛顿莱普利兹公式也无法直接运用，因此有必要研究积分的数值计算问题。 对于一些理论的推导，大家可以看看维基百科，下面我主要给出牛顿-科特斯公式在n=...
数值积分 复化梯形求积公式 | 北太天元 or matlab
最新发布
以算法实现为主
04-25 1468
数值积分——复化梯形求积公式 | 北太天元 主要以算法实现为主，使用软件为北太天元
Python采用复化梯形公式、复化Simpson公式求积分的程序
05-23
以下是使用复化梯形公式和复化Simpson公式求积分的Python程序： 复化梯形公式： ```python def trapezoidal(f, a, b, n): h = (b - a) / n x = a s = f(x) for i in range(1, n): x += h s += 2 * f(x) s += f(b) return h * s / 2 # 示例 def f(x): return x ** 2 print(trapezoidal(f, 0, 1, 100)) # 输出 0.3333333333333333 ``` 复化Simpson公式： ```python def simpson(f, a, b, n): h = (b - a) / n x = a s = f(x) for i in range(1, n): x += h if i % 2 == 0: s += 2 * f(x) else: s += 4 * f(x) s += f(b) return h * s / 3 # 示例 def f(x): return x ** 2 print(simpson(f, 0, 1, 100)) # 输出 0.3333333333333333 ``` 其中，`f`为要求积分的函数，`a`和`b`为积分区间的上下限，`n`为分割的区间数。在示例中，分割区间数为100，可以根据需要进行调整。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值