Python实现数值积分方式
作者:MachineRandy 发布时间:2022-01-23 14:00:37
标签:Python,数值,积分
原理:
利用复化梯形公式,复化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)) + "个步数")
来源:https://blog.csdn.net/MachineRandy/article/details/82634757
0
投稿
猜你喜欢
- 先来看看基本的定义:channel是Go语言中的一个核心类型,可以把它看成管道。并发核心单元通过它就可以发送或者接收数据进行通讯,这在一定程
- wheel文件Wheel和Egg都是python的打包格式,目的是支持不需要编译或制作的安装过程,实际上也是一种压缩文件,将.whl的后缀改
- pandas模块pandas是一个强大的分析结构化数据的工具集;它的使用基础是Numpy(提供高性能的矩阵运算);用于数据挖掘和数据分析,同
- 一个不错的绿色下划线的简洁CSS导航代码,纯css不用图片,效果图如下:<html> <head> <meta
- 在某些情况下,程序需要在定义函数时为一个或多个形参指定默认值,这样在调用函数时就可以省略为该形参传入参数值,而是直接使用该形参的默认值。为形
- 一、 for 循环根据变量赋值的次数进行循环for item in ["tom","bob",&qu
- 如下所示:plt.subplots_adjust(left=None,bottom=None,right=None,top=None,wsp
- 本文实例为大家分享了python自动发送邮件的具体代码,供大家参考,具体内容如下#coding=utf8 ''&
- python常见的错误有1.NameError变量名错误2.IndentationError代码缩进错误3.AttributeError对象
- 首先说一下需求,我需要将数据以分钟为单位进行分组,然后每一分钟内的数据作为一行输出,因为不同时间的数据量不一样,所以所有数据按照最长的那组数
- 概述OpenCV 是一个跨平台的计算机视觉库, 支持多语言, 功能强大. 今天小白就带大家一起携手走进 OpenCV 的世界. 高
- 在来回切换中英文输入法的时候连按两下shift总是会蹦出来全局搜索框真的很是麻烦,现在是把这个框给禁用掉1.按ctrl+shift+a,弹出
- 一、备份数据库1、打开SQL企业管理器,在控制台根目录中依次点开Microsoft SQL Server2、SQL Server组-->
- 1、Dreamweaver中的复制我在网页中复制的文字,粘贴到Dreamweaver中时,它总是带有原来网页的格式,请问如何只复制其中的文本
- 在默认的情况下,MySQL搜索不区分大小写(但某些字符集始终区分大小写,如czech)。这意味着,如果你使用col_name LIKE
- Paramiko是一个用于执行SSH命令的Python第三方库,使用该库可实现自动化运维的所有任务,如下是一些常用代码的封装方式,多数代码为
- 1. 安装yaml库想要使用python实现yaml与json格式互相转换,需要先下载pip,再通过pip安装yaml库。如何下载以及使用p
- 这可是个综合性的问题,看看下面对文件操作的集大成代码:<% 'Set file i/
- Python中有一个有趣的语法,只要定义类型的时候,实现__call__函数,这个类型就成为可调用的。换句话说,我们可以把这个类型的对象当作
- 现在正在搞三层开发,用ASP和VB6.0,但是现在苦于没有找到合适的方法来调试自己写的DLL文件,效率相当低。 &n