数值分析
From Star.
I'm going on an adventure!
展开
-
【数值分析】python实现四阶亚当姆斯预测矫正方法
python代码import numpy as npimport matplotlib.pyplot as pltDT = np.array(input("份数:").split(' '))y0 = int(input("函数起点:"))def f(y): return np.exp(-y)def Adams(t, y0): n = len(t) y = np.zeros(n) y[0] = y0 for i in range(0, n-1):.原创 2021-04-04 00:48:40 · 1071 阅读 · 0 评论 -
【数值分析】python实现四阶龙格库塔法
python代码import numpy as npimport matplotlib.pyplot as pltDT = np.array(input("份数:").split(' '))y0 = int(input("函数起点:"))def f(y): return np.exp(-y)def Runge_Kutta(t, y0): n = len(t) y = np.zeros(n) y[0] = y0 for i in range(0, .原创 2021-04-03 21:49:33 · 4597 阅读 · 0 评论 -
【数值分析】python实现复化高斯积分
python代码import pandas as pdimport numpy as npimport math# 初始化表格x1, k1 = np.array([[0], [2]])x2, k2 = np.array([[0.5773502692], [1]])x3, k3 = np.array([[0.7745966692, 0], [0.5555555556, 0.8888888889]])x4, k4 = np.array([[0.861.原创 2021-04-03 20:56:12 · 1096 阅读 · 2 评论 -
【数值分析】python实现复化辛普森积分
python代码import numpy as npa, b = input("积分区间:").split(' ')n = int(input("子区间个数:"))a = int(a)b = int(b)x_i = np.linspace(a, b, n+1) # 区间结点h = (b-a)/n# 原函数def f(x): return 4/(1+x**2)def piece(x_i, i, h): return h/6 * (f(x_i[i-1]) + 4.原创 2021-04-03 19:48:02 · 2568 阅读 · 0 评论 -
【数值分析】python实现拉格朗日插值及作图
python代码import numpy as npX = input("x的值:").split(' ')Y = input("y的值:").split(' ')x = input("要预测的值:")X = np.array(X).astype(np.float64)Y = np.array(Y).astype(np.float64)x = np.array(x).astype(np.float64)n = len(X)# 累乘函数def T(x, i, X): T.原创 2021-04-02 19:27:57 · 3228 阅读 · 2 评论