import numpy as np
import matplotlib.pyplot as pltdef inter_x(a, b, n):""":param a: 插值区间左端点:param b: 插值区间右端点:param n: 插值点个数:return: 插值节点向量"""x np.zeros(n1)h (b - a) / nfor i in range(n1):x[i] a i *…
Axb
系数矩阵A按行严格对角占优
若系数矩阵A按行严格对角占有>Jacobi迭代法和Gauss-Seidel迭代法均收敛
迭代矩阵B谱半径<1
Jacobi迭代法和Gauss-Seidel迭代法收敛<>迭代矩阵B的谱半径<1 其中,Jacobi迭代法的迭代矩阵 B − D − 1 ( L U ) B-D…
反幂法20231226
一种求实矩阵 A {A} A 的按模最小的特征值,及其对应的特征向量 x i {x_i} xi 的方法,只能求一个。 要保证矩阵最小特征值只有一个,有 n {n} n 个线性无关的特征向量,矩阵可逆。 可以通过求矩阵 A − 1 {A^{…
区间折半法
从梯形公式出发,上一步步长为 h {h} h ,则有步长折半后的积分 T 2 n 1 2 T n h 2 ∑ i 0 n − 1 f ( x i 0.5 ) T_{2n} \frac{1}{2}T_n \frac{h}{2} \sum_{i0}^{ n-1}f(x_{i0.5}) T2n21Tn2hi0∑n−1f(xi0.5) matlab实现
…
4. 牛顿法
收敛时牛顿法的收敛速度是二阶的,不低于二阶。如果函数有重根,牛顿法一般不是二阶收敛的。 x k 1 x k − f ( x k ) f ′ ( x k ) x_{k1}x_k- \frac{f(x_k)}{f(x_k)} xk1xk−f′(xk)f(xk) matlab实现
%% 牛顿迭代例子
f (x) x…
例题 代码
import numpy as np
np.set_printoptions(suppressTrue)def h_u_lamb_A(x):h []u []lamb []nx.shape[0]A np.mat(2*np.eye(n-2, n-2))A[0, 1] 1A[n-3, n-4] 1for i in range(0, n - 3):h.append(x[i 2] - x[i1])for i in range(0, n - 4):u.append(h[i] / (h…