快捷導(dǎo)航

Python實(shí)現(xiàn)數(shù)值積分方式

更新時(shí)間：2019年11月20日 10:09:23 作者：MachineRandy

今天小編就為大家分享一篇Python實(shí)現(xiàn)數(shù)值積分方式，具有很好的參考價(jià)值，希望對大家有所幫助。一起跟隨小編過來看看吧

原理：

利用復(fù)化梯形公式，復(fù)化Simpson公式，計(jì)算積分。

步驟：

import math

"""測試函數(shù)"""
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)

"""打印顯示函數(shù)"""   
def p(i,n):
  return "第" + str(i) + "題,n=" + str(n) + "時(shí)的積分值為："

"""復(fù)化Simpson函數(shù)"""
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("復(fù)化Simpson函數(shù)" + 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("復(fù)化梯形函數(shù)" + 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("復(fù)化梯形函數(shù)求解第一問，精度為0.00001時(shí)需要" + str(Times(0.00001)) + "個步數(shù)")
print("復(fù)化Simpson函數(shù)求解第一問，精度為0.00001時(shí)需要" + str(SimpsonTimes(0.00001)) + "個步數(shù)")
print("復(fù)化梯形函數(shù)求解第一問，精度為0.000001時(shí)需要" + str(Times(0.000001)) + "個步數(shù)")
print("復(fù)化Simpson函數(shù)求解第一問，精度為0.000001時(shí)需要" + str(SimpsonTimes(0.000001)) + "個步數(shù)")

以上這篇Python實(shí)現(xiàn)數(shù)值積分方式就是小編分享給大家的全部內(nèi)容了，希望能給大家一個參考，也希望大家多多支持腳本之家。

您可能感興趣的文章: