- 博客(3)
- 收藏
- 关注
RSS订阅原创 单步线性定常迭代法
现在，在计算x(k)用x(k-1)的各个分量计算，但当计算x(k)的第二分量时，由于x(k)的第一个分量已经算出，则用它来代替x(k-1)的第一个分量，其余分量不变。同Jacobi迭代法的计量相同，G-S迭代法可以写成 x(k) = (D-L)(逆)*U*x(k-1) + (D-L)(逆)*b。则方程组可以写成 x = B*x + g，其中B = D(逆)*(L+U)，g = D(逆)*b。print("请输入非奇异矩阵(A)（一行一行的输入）：")print("请输入非奇异矩阵的行标(n)：")
2023-05-07 19:27:44
223
1
原创 最小二乘问题的正则化方法求解（QR分解的实验报告）
解决方法：通过人为的变化，例如矩阵出现可忽略不计的值时，可计为零。print("注:由于计算机的精度问题，存在一定误差。1.问题：矩阵的行标列标在进行运算时可能存在不匹配的情况。print("输出的上三角+零矩阵(QR分解的R):")2.问题：由于计算机的精度，计算值和理论值存在一定误差。print("输出的正交矩阵(QR分解的Q):")解决方法：用逐步输出的方法找到问题所在并解决。print("请输入需要求的矩阵：")print("请输入矩阵的行标：")print("请输入矩阵的列标：")
2023-04-05 16:51:23
481
1
原创 全主元三角分解实验报告（数值线性代数）（未使用numpy）
如在找最大主元时，进行到第k步时，只需要定义的函数作用于第k行和第k列后的元素。注: 在全主元三交分解中，可能遇到的问题，如在定义找最大主元的函数、换行函数、换列函数时，会因为作用矩阵的范围的变化而导致错误。jz[i][j]=float(1) #定义一个单位矩阵。a[i][j]=float(input()) #输入一个矩阵。a[i][l] = c #创建列变换函数。a[k] = m #创建行变换函数。print("请输入一个矩阵:(一行一行的输入)")
2023-03-31 21:42:57
121
1
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人