python乘积函数_龙贝格求积公式（Python实现）

#Author：glm233

#这个程序就是一个黑箱接口，把要进行龙贝格求积的函数放在func函数里，可以自行修改，然后范围就是在Romberg里改

import math

'''给定一个函数，如：f(x)= x^(3/2），和积分上下限a,b，用机械求积Romberg公式求积分。'''

import numpy as np

def func(x):

return x**2+math.sin(x)/x

class Romberg:

def __init__(self, integ_dowlimit, integ_uplimit):

'''初始化积分上限integ_uplimit和积分下限integ_dowlimit输入一个函数，输出函数在积分上下限的积分'''

self.integ_uplimit = integ_uplimit

self.integ_dowlimit = integ_dowlimit

def calc(self):

'''计算Richardson外推算法的四个序列'''

t_seq1 = np.zeros(5, 'f')

s_seq2 = np.zeros(4, 'f')

c_seq3 = np.zeros(3, 'f')

r_seq4 = np.zeros(2, 'f')

# 循环生成hm间距序列

hm = [(self.integ_uplimit - self.integ_dowlimit) / (2 ** i) for i in range(0,5)]

print(hm)

# 循环生成t_seq1

fa = func(self.integ_dowlimit)

fb = func(self.integ_uplimit)

t0 = (1 / 2) * (self.integ_uplimit - self.integ_dowlimit) * (fa+fb)

t_seq1[0] = t0

for i in range(1, 5):

sum = 0

# 多出来的点的累加和

for each in range(1, 2**i,2):

sum =sum + hm[i]*func( self.integ_dowlimit+each * hm[i])#计算两项值

temp1 = 1 / 2 * t_seq1[i - 1]

temp2 =sum

temp = temp1 + temp2

# 求t_seql的1-4位

t_seq1[i] = temp

print('T序列：'+ str(list(t_seq1)))

# 循环生成s_seq2

s_seq2 = [round((4 * t_seq1[i + 1] - t_seq1[i]) / 3,6) for i in range(0, 4)]

print('S序列：' + str(list(s_seq2)))

# 循环生成c_seq3

c_seq3 = [round((4 ** 2 * s_seq2[i + 1] - s_seq2[i]) / (4 ** 2 - 1),6) for i in range(0, 3)]

print('C序列：' + str(list(c_seq3)))

# 循环生成r_seq4

r_seq4 = [round((4 ** 3 * c_seq3[i + 1] - c_seq3[i]) / (4 ** 3 - 1),6) for i in range(0, 2)]

print('R序列：' + str(list(r_seq4)))

r_seq5 = [round((4 ** 4 * r_seq4[i + 1] - r_seq4[i]) / (4 ** 4 - 1),6) for i in range(0, 1)]

print('A序列：' + str(list(r_seq5)))

return 'end'

rom = Romberg(0.3, 0.8)

print(rom.calc())

python乘积函数_龙贝格求积公式（Python实现）相关推荐

数值积分：梯形规则--复合梯形规则--辛普森规则--复合辛普森规则--龙贝格求积公式
数值积分:梯形规则--复合梯形规则--辛普森规则--复合辛普森规则--龙贝格求积公式 1.问题描述微积分方法求积有很大的局限性,当碰到被积函数很复杂时,找不到相应的原函数.积分值在几何上可解释为由 ...
Romberg(龙贝格)求积公式求解数值积分时的注意事项
<数值分析>第5版(李庆扬编著)的第四章课后习题第8-(2)题中,要求使用Romberg(龙贝格)求积公式求解f(x)=xsinx在区间[0,2pi]上的积分,要求误差小于10^(-5). ...
python fact函数_第5天：Python 函数
函数是组织好的,可重复使用的,用来实现单一,或相关联功能的代码段,所以我经常说函数是程序员规模化使用的基础. 函数能提高应用的模块性,和代码的重复利用率.在程序设计中,常将一些常用的功能模块编写成函数 ...
Romberg算法（龙贝格求积公式C语言实现）
目录 1.原理 2.案例 3.代码 4.结果 5.误差分析 1.原理 2.案例 3.代码 #include<stdio.h> #include<math.h> float f( ...
python事件处理函数_事件驱动的简明讲解(python实现)
关键词:编程范式,事件驱动,回调函数,观察者模式作者:码匠信龙举个简单的例子: 有些人喜欢的某个公众号,然后去关注这个公众号,哪天这个公众号发布了篇新的文章,没多久订阅者就会在微信里收到这个公众号 ...
python纯函数_理想国真恵玩Python从入门到精通第006天_纯函数写游戏管理系统
原标题:理想国真恵玩Python从入门到精通第006天_纯函数写游戏管理系统前面已经带大家学习了函数,高级数据类型,比如说字典,今天带大家用函数加字典做一个游戏管理系统,希望大家喜欢.废话不多说,直 ...
python数学函数_「分享」关于Python整理的常用数学函数整理
原标题:「分享」关于Python整理的常用数学函数整理 1.函数说明 abs(number)返回数字的绝对值,如abs(-10)返回10 pow(x,y[,z]) 返回x的y次幂(所得结果对z取模), ...
python初始化函数_当你学会了Python爬虫，网上的图片素材就免费了
前言本文的文字及图片来源于网络,仅供学习.交流使用,不具有任何商业用途,版权归原作者所有,如有问题请及时联系我们以作处理. 加入作者的python学习圈子:1156465813 即可免费获取,资料全 ...
python初始化函数_当你学会了Python爬虫，网上的图片素材就开始免费了
前言本文的文字及图片来源于网络,仅供学习.交流使用,不具有任何商业用途,版权归原作者所有,如有问题请及时联系我们以作处理. 加企鹅群695185429即可免费获取,资料全在群文件里.资料可以领取包括 ...

python乘积函数_龙贝格求积公式（Python实现）

python乘积函数_龙贝格求积公式（Python实现）相关推荐

最新文章

热门文章