原理:
利用复化梯形公式,复化Simpson公式,计算积分。
步骤:
import math """测试函数""" def f(x,i): if i == 1: return (4 - (math.sin(x)) ** 2) ** 0.5 if i == 2: if x == 0: return 1 else: return math.sin(x) / x if i == 3: return (math.exp(x)) / (4 + x ** 2) if i == 4: return math.log(1+x,math.e) / (1 + x ** 2) """打印显示函数""" def p(i,n): return "第" + str(i) + "题,n=" + str(n) + "时的积分值为:" """复化Simpson函数""" def Simpson(a, b, n, i): h = (b - a) / (2 * n) F0 = f(a,i) + f(b,i) F1 = 0 F2 = 0 for j in range(1,2 * n): x = a + (j * h) if j % 2 == 0: F2 = F2 + f(x,i) else: F1 = F1 + f(x,i) SN = (h * (F0 + 2 * F2 + 4 * F1)) / 3 print("复化Simpson函数" + p(i,n) + str("%-10.7f"%(SN))) return SN def T(a, b, n, i): h = (b - a) / n F0 = f(a,i) + f(b,i) F = 0 for j in range(1,n): x = a + (j * h) F = F + f(x,i) SN = (h * (F0 + 2 * F)) / 2 print("复化梯形函数" + p(i,n) + str("%-10.7f"%(SN))) return SN def SimpsonTimes(x): n = 1 y = Simpson(0, math.pi/4, n, 1) while(abs(y - 1.5343916) > x): n = n + 1 y = Simpson(0, math.pi/4, n, 1) else: return n def Times(x): n = 1 y = T(0, math.pi/4, n, 1) while(abs(y - 1.5343916) > x): n = n + 1 y = T(0, math.pi/4, n, 1) else: return n """ 测试部分 """ Simpson(0, math.pi/4, 10, 1) Simpson(0, 1, 10, 2) Simpson(0, 1, 10, 3) Simpson(0, 1, 10, 4) Simpson(0, math.pi/4, 20, 1) Simpson(0, 1, 20, 2) Simpson(0, 1, 20, 3) Simpson(0, 1, 20, 4) T(0, math.pi/4, 10, 1) T(0, 1, 10, 2) T(0, 1, 10, 3) T(0, 1, 10, 4) T(0, math.pi/4, 20, 1) T(0, 1, 20, 2) T(0, 1, 20, 3) T(0, 1, 20, 4) print("复化梯形函数求解第一问,精度为0.00001时需要" + str(Times(0.00001)) + "个步数") print("复化Simpson函数求解第一问,精度为0.00001时需要" + str(SimpsonTimes(0.00001)) + "个步数") print("复化梯形函数求解第一问,精度为0.000001时需要" + str(Times(0.000001)) + "个步数") print("复化Simpson函数求解第一问,精度为0.000001时需要" + str(SimpsonTimes(0.000001)) + "个步数")
以上这篇Python实现数值积分方式就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持。
标签:
Python,数值,积分
免责声明:本站文章均来自网站采集或用户投稿,网站不提供任何软件下载或自行开发的软件!
如有用户或公司发现本站内容信息存在侵权行为,请邮件告知! 858582#qq.com
白云城资源网 Copyright www.dyhadc.com
暂无“Python实现数值积分方式”评论...
《魔兽世界》大逃杀!60人新游玩模式《强袭风暴》3月21日上线
暴雪近日发布了《魔兽世界》10.2.6 更新内容,新游玩模式《强袭风暴》即将于3月21 日在亚服上线,届时玩家将前往阿拉希高地展开一场 60 人大逃杀对战。
艾泽拉斯的冒险者已经征服了艾泽拉斯的大地及遥远的彼岸。他们在对抗世界上最致命的敌人时展现出过人的手腕,并且成功阻止终结宇宙等级的威胁。当他们在为即将于《魔兽世界》资料片《地心之战》中来袭的萨拉塔斯势力做战斗准备时,他们还需要在熟悉的阿拉希高地面对一个全新的敌人──那就是彼此。在《巨龙崛起》10.2.6 更新的《强袭风暴》中,玩家将会进入一个全新的海盗主题大逃杀式限时活动,其中包含极高的风险和史诗级的奖励。
《强袭风暴》不是普通的战场,作为一个独立于主游戏之外的活动,玩家可以用大逃杀的风格来体验《魔兽世界》,不分职业、不分装备(除了你在赛局中捡到的),光是技巧和战略的强弱之分就能决定出谁才是能坚持到最后的赢家。本次活动将会开放单人和双人模式,玩家在加入海盗主题的预赛大厅区域前,可以从强袭风暴角色画面新增好友。游玩游戏将可以累计名望轨迹,《巨龙崛起》和《魔兽世界:巫妖王之怒 经典版》的玩家都可以获得奖励。
更新日志
2024年10月06日
2024年10月06日
- 群星《前途海量 电影原声专辑》[FLAC/分轨][227.78MB]
- 张信哲.1992-知道新曲与精丫巨石】【WAV+CUE】
- 王翠玲.1995-ANGEL【新艺宝】【WAV+CUE】
- 景冈山.1996-我的眼里只有你【大地唱片】【WAV+CUE】
- 群星《八戒 电影原声带》[320K/MP3][188.97MB]
- 群星《我的阿勒泰 影视原声带》[320K/MP3][139.47MB]
- 纪钧瀚《胎教古典音乐 钢琴与大提琴的沉浸时光》[320K/MP3][148.91MB]
- 刘雅丽.2001-丽花皇后·EMI精选王【EMI百代】【FLAC分轨】
- 齐秦.1994-黄金十年1981-1990CHINA.TOUR.LIVE精丫上华】【WAV+CUE】
- 群星.2008-本色·百代音乐人创作专辑【EMI百代】【WAV+CUE】
- 群星.2001-同步过冬AVCD【环球】【WAV+CUE】
- 群星.2020-同步过冬2020冀待晴空【环球】【WAV+CUE】
- 沈雁.1986-四季(2012梦田复刻版)【白云唱片】【WAV+CUE】
- 纪钧瀚《胎教古典音乐 钢琴与大提琴的沉浸时光》[FLAC/分轨][257.88MB]
- 《国语老歌 怀旧篇 3CD》[WAV/分轨][1.6GB]