![](https://img-blog.csdnimg.cn/20201014180756928.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数值分析实验
SArthas
这个作者很懒,什么都没留下…
展开
-
【数值分析实验】矩阵特征值及特征向量计算:幂法、反幂法、幂法加速(python)
矩阵特征值及特征向量计算调包import numpy as np幂法#幂法算法 def M(A,u0,epsilon = 1e-6,M=100): uk = u0 mk_1 = 1 k = 0 while True: k += 1 vk = A*uk mk = vk[np.argmax(np.abs(vk))] uk = vk / mk if np.abs(mk - mk_原创 2022-01-03 23:11:34 · 1540 阅读 · 0 评论 -
【数值分析实验】线性方程组的迭代解法:雅可比迭代、高斯赛德尔迭代、SOR迭代算法(c#)
线性方程组的迭代解法class Operation { double[][] a; //系数矩阵 double[] b; //常数向量 double[] x00,x0, x; //初始迭代向量和迭代向量 int n; //矩阵大小 public double epsilon = 1e-5; //最大精度 int M, m; //最大迭代向量次数和迭代次数原创 2022-01-03 23:07:36 · 1016 阅读 · 0 评论 -
【数值分析实验】线性代数方程组的直接解法:列主元高斯消去法、LU分解法、LU分解法求矩阵的逆(python)
线性代数方程组的直接解法调包消元函数和矩阵求解高斯消去法列主元高斯消去法LU分解法LU分解求矩阵的逆下三角矩阵求逆LU分解求逆演示调包import numpy as np消元函数和矩阵求解#消元def KillYuan(A,i): try: A[i] /= A[i][i] except: print("矩阵奇异") A[i] /= A[i][i] for j in range(i + 1,len(A)): A[j]原创 2022-01-03 23:00:01 · 1520 阅读 · 0 评论 -
【数值分析实验】常微分方程初值问题:显示欧拉法、隐式欧拉法、欧拉改进法、四阶龙格库塔(python)
常微分方程初值问题的数值解法调包import mathimport numpy as npimport matplotlib.pyplot as plt显示欧拉法fStr为函数str名#显式欧拉法def EulerExplicit(x0,y0,h,fStr): xn = x0 yn = y0 n = 0 ns = [n] xs = [xn] ys = ['%.8f'%yn] while n < 50: n +=原创 2022-01-03 22:53:22 · 2605 阅读 · 0 评论 -
【数值分析实验】数值积分:复化梯形、复化辛普森、龙贝格求积公式、自适应求积公式(python)
数值积分调包import mathimport numpy as npimport matplotlib.pyplot as pltfStr为被积函数函数名(str)梯形求积#梯形求积def Trapezium(a,b,fStr): return (b - a) / 2 * (globals()[fStr](a) + globals()[fStr](b))复化梯形#复化梯形def TrapeziumComplex(a,b,epsilon,fStr): def F(原创 2021-12-29 22:38:47 · 2891 阅读 · 0 评论 -
【数值分析实验】插值与拟合:拉格朗日插值、牛顿插值、分段插值;线性拟合、最小二乘拟合(python)
插值与拟合插值拉格朗日插值牛顿插值分段Hermite插值拟合最小二乘法拟合例题最小二乘拟合变式例题调包import mathimport numpy as npimport sympy as spimport matplotlib.pyplot as plt插值拉格朗日插值#拉格朗日插值def Lagrange(x,xs,ys): n = len(xs) ls = [1 for i in range(n)] for i in range(n): f原创 2021-12-29 22:14:28 · 1834 阅读 · 0 评论 -
【数值分析实验】非线性方程求根:二分法、牛顿迭代法、弦截法(python)
非线性方程求根二分法#二分法def Dichotomy(a,b,fStr,epsilon = 1e-6,eta = 1e-9): def f(x): return globals()[fStr](x) n = 0 fa = f(a) fb = f(b) if fa * fb >= 0: print("算法失败") return while abs(b - a) > epsilon:原创 2021-12-29 21:53:51 · 2005 阅读 · 0 评论