R语言高斯消去法解线性方程组

最新推荐文章于 2024-09-02 07:53:48 发布

酒月空鸣

最新推荐文章于 2024-09-02 07:53:48 发布

阅读量1.3k

点赞数

CC 4.0 BY-SA版权

分类专栏： R语言

本文链接：https://blog.youkuaiyun.com/qq_42681040/article/details/103808916

R语言专栏收录该内容

1 篇文章

订阅专栏

本文介绍了一种使用R语言实现的高斯消元法，该方法用于求解线性方程组。通过输入未知量数目和增广矩阵，算法能够输出解向量。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

guass <- function(n,B){
  if(ncol(B) != n+1){
    stop("data input error")
  }
  U <- B
  for (i in 1:(n-1)) {
    for (j in (i+1):n) {
      for (k in i:(n+1)) {
        U[j,k] <- B[j,k] - (B[j,i]*B[i,k])/B[i,i]
      }
    }
    B <- U
  }
  x <- array(,n)
  x[n] <- U[n,n+1]/U[n,n]
  for (p in (n-1):1) {
    for (q in n:(p+1)) {
      U[p,n+1] <- U[p,n+1] - U[p,q]*x[q]
    }
    x[p] <- U[p,n+1]/U[p,p]
  }
  x
}

R语言高斯消去法，输入未知量数目和线性方程组的增广矩阵，输出为解向量