- 博客(10)
- 收藏
- 关注
原创 python实现函数的最佳平方逼近
先放代码,有时间补内容import sympy as syimport numpy as nyfrom scipy import integratefrom guass_swap import gaussindef InterProduct(f,g,do,up): val=sy.integrate(f*g,(x, #函数 do,#x下界0 up))#x上界pi return valdef SquaresA
2020-05-24 15:35:08 2994
原创 python实现最小二乘方法
先放代码,有时间补内容import numpy as nyimport sympy as syfrom guass_swap import gaussinny.set_printoptions(suppress=True)def LeastSquares(X,Y,m): assert len(X)==len(Y) n=len(X) A = ny.zeros([m,m]) B = ny.transpose(ny.zeros([m])) for i in r
2020-05-24 15:34:11 469
原创 python实现用改进的乘幂法求方阵的按模最大特征值和特征向量
先放代码,有时间补内容import numpy as npdef MaxFabs(v): m,n = v.shape assert n == 1 ans = 0 for i in range(m): if ans < np.fabs(v[i]): ans = v[i] maxn = i return ansdef PowerMethod(a,v,g): m,n = a.shape
2020-05-24 15:32:45 1560
原创 拉格朗日插值法(Lagrange)求数值解的python实现
先放代码,有时间补内容def LagrangeInterpolation(x, y, x0): m = len(x) n = len(y) assert m==n L = 0 l = [1]*n for j in range(n): for i in range(n): if i != j: ...
2020-04-04 16:10:09 940
原创 牛顿迭代法(Newton)解方程数值解的python实现
事情过多,先放代码,有时间再补内容。import sympy as sydef FunValue(x0, f): result = f.subs(x, x0) return resultdef Newton(x0, f, g = 1e-6): times = 0 while True: tempx = x0 x0 = x0 ...
2020-04-04 16:06:08 5977
原创 Gauss-Seidel迭代法的python实现
Gauss-Seidel迭代法求线性方程组的原理和python代码Gauss-Seidel迭代法求线性方程组的介绍举例讲解算法描述Gauss-Seidel迭代法的python代码Gauss-Seidel迭代法求线性方程组的介绍举例讲解算法描述Gauss-Seidel迭代法的python代码import numpy as npdef G_S(a, b, x, g): # a...
2020-03-23 10:56:56 6287 5
原创 矩阵LU分解中,Crout分解法的python实现
数值分析中,LU分解的Crout分解法的原理和python代码LU分解法分为Doolittle分解法和Crout分解法,其中这里只介绍Crout分解。LU分解法的介绍如下图这里有一个具体的实例算法描述Crout分解法的python代码import numpy as npdef Crout(a, b): cout = 0 m, n = a.shape i...
2020-03-23 10:45:08 4084 1
原创 高斯列主元消去法的python实现
数值分析中,高斯列主元消去法的python实现首先要了解高斯消去法的原理和代码,高斯消去法的原理和python实现,高斯列主元消去法是高斯消去法的改进。高斯列主元消去法的python代码import numpy as npdef swap(a, b, k, n): # 找到主元并交换,这仅是一个仅用来交换的函数 ans = 0 for i in ra...
2020-03-23 10:37:36 8909 3
原创 高斯消去法的python实现
数值分析中,高斯消去法的原理和python代码高斯消去法(Gauss)的原理如下:下面是python代码import numpy as np #首先确保python安装了numpy包def gauss(a, b): #自己定义一个函数,自变量是a,b cout = 0 #定义计算次数 m, n = a.shape #矩阵a的行数和列数 ...
2020-03-23 10:10:34 6330 1
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人