- 博客(5)
- 收藏
- 关注
原创 共轭梯度法求解线性方程组python代码
从任意给定的初始点出发,沿一组关于矩阵A的共轭方向进行线性搜索,在无舍入误差的假定下,最多迭代n次(其中n为矩阵A的阶数),就可求得二次函数的极小点,也就求得线性方程组Ax= b的解。r = b-np.dot(A,x) #r=b-Ax r也是是梯度方向。f=cg(A,b,x) #调用共轭斜量法。alpha = np.dot(r,r)/pap #直接套用公式。def cg(A,b,x): #共轭梯度法。n=10 #n 个解。
2024-06-10 10:56:28 415
原创 python高斯消元法代码
当用于一个矩阵时,高斯消元法会产生出一个"行梯阵式"。高斯消元法可以用在电脑中来解决数千条。数学上,高斯消元法(或译:高斯消去法),是线性代数中的一个。,可用来为线性方程组求解,求出矩阵的秩,以及求出可逆方阵的。
2024-06-08 17:53:15 724
原创 LU分解法计算线性方程组
for i in range(k + 1, n): # i表示第二层循环,(k+1,n)行,计算该行消元的系数。x[n - 1] = b[n - 1] / a[n - 1][n - 1] #先算最后一位的x解。for k in range(n - 1): # k表示第一层循环,(0,n-1)行。for i in range(1,n-1):# 求U的第i行 L的第i行。for i in range(n-2,-1,-1):#求x[i]sum = 0.0 #求和。#求U[n-1][n-1]
2024-06-07 17:51:21 328
原创 用python实现雅可比矩阵
在向量微积分中,雅可比矩阵是一阶偏导数以一定方式排列成的矩阵,其行列式称为雅可比行列式。还有,在代数几何中,代数曲线的雅可比量表示雅可比簇:伴随该曲线的一个群簇,曲线可以嵌入其中。雅可比矩阵的重要性在于它体现了一个可微方程与给出点的最优线性逼近。因此,雅可比矩阵类似于多元函数的导数。雅可比矩阵定义为向量对向量的微分矩阵。print("精确度等于", eps, "时,需要迭代", times, "次收敛")print("在最大迭代次数内不收敛", "最大迭代次数后的结果为", x)A:为方程组的系数矩阵。
2024-06-07 17:33:41 421
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人