您好,登錄后才能下訂單哦!
原理:
利用復(fù)化梯形公式,復(fù)化Simpson公式,計(jì)算積分。
步驟:
import math """測(cè)試函數(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 """ 測(cè)試部分 """ 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ù)求解第一問(wèn),精度為0.00001時(shí)需要" + str(Times(0.00001)) + "個(gè)步數(shù)") print("復(fù)化Simpson函數(shù)求解第一問(wèn),精度為0.00001時(shí)需要" + str(SimpsonTimes(0.00001)) + "個(gè)步數(shù)") print("復(fù)化梯形函數(shù)求解第一問(wèn),精度為0.000001時(shí)需要" + str(Times(0.000001)) + "個(gè)步數(shù)") print("復(fù)化Simpson函數(shù)求解第一問(wèn),精度為0.000001時(shí)需要" + str(SimpsonTimes(0.000001)) + "個(gè)步數(shù)")
以上這篇Python實(shí)現(xiàn)數(shù)值積分方式就是小編分享給大家的全部?jī)?nèi)容了,希望能給大家一個(gè)參考,也希望大家多多支持億速云。
免責(zé)聲明:本站發(fā)布的內(nèi)容(圖片、視頻和文字)以原創(chuàng)、轉(zhuǎn)載和分享為主,文章觀點(diǎn)不代表本網(wǎng)站立場(chǎng),如果涉及侵權(quán)請(qǐng)聯(lián)系站長(zhǎng)郵箱:is@yisu.com進(jìn)行舉報(bào),并提供相關(guān)證據(jù),一經(jīng)查實(shí),將立刻刪除涉嫌侵權(quán)內(nèi)容。