函數將線性方程轉化爲R中的矩陣形式？

我想知道是否存在能夠將線性方程組轉換爲矩陣形式（例如，通過Gauss Seidel algorithm解決的問題）的R的任何包或其他預構建解決方案，類似於equationsToMatrix(eqns,vars)function in Matlab？函數將線性方程轉化爲R中的矩陣形式？

從MATLAB一個例子：

[A, b] = equationsToMatrix([x - y == 0, x + 2*y == 3, [x, y]) 

A = 
[ 1, -1] 
[ 1, 2] 

b = 
    0 
    3

建議積木將是非常有益的，太。

來源

2014-11-22 dreamon

你會用這個做什麼？只需直接輸入矩陣即可輕鬆打字。 – Roland 2014-11-22 12:10:59

1）這是你問不正是對，但也許它會幫助反正：

library(Ryacas) 
x <- Sym("x") 
y <- Sym("y") 
Simplify(Solve(List(x - y == 0, x + 2*y == 3), List(x, y)))

，並提供：

expression(list(list(x - y == 0, y - 1 == 0)))

2）如果我們知道這些都是線性那麼問題中所示形式的方程就會嘗試這一點。這兩個strapply調用執行正則表達式與args組件匹配，捕獲括號內正則表達式部分匹配的字符串，並將那些捕獲的字符串作爲參數調用指定爲第三個參數的函數。我們將strapply輸出結合使用rbind.fill，並將其生成的所有NAs替換爲零。

library(gsubfn) # strapply 
library(plyr) # rbind.fill 

eqn <- function(...) { 
    args <- c(...) 
    x2num <- function(x, y) { # determine coefficient value as a numeric 
     z <- gsub(" ", "", x) 
     setNames(if (z == "-") -1 else if (z == "") 1 else as.numeric(z), y) 
    } 
    lhs <- strapply(args, "(-? *\\d*)[ *]*([a-z])", x2num) 
    lhs <- do.call(rbind.fill, lapply(lhs, function(x) as.data.frame(t(x)))) 
    lhs <- as.matrix(lhs) 
    lhs[] <- ifelse(is.na(lhs), 0, lhs) 
    list(lhs = lhs, rhs = strapply(args, "== *(\\d)", as.numeric, simplify = TRUE)) 
} 

# test it out 
eqn("x - y == 0", "2*y == 3")

，並提供：

$lhs 
    x y 
[1,] 1 -1 
[2,] 0 2 

$rhs 
[1] 0 3

更新：廣義所以現在不是所有的變量需要在每個方程中，也變量可以在不同的公式不同的訂單。

來源

2014-11-22 13:30:52

這與我所尋找的非常接近，謝謝！ – dreamon 2014-11-30 13:25:26

無法從您的示例中看出您是僅僅想要簡單的線性方程式求解還是更普遍的系統求解器。如果是後者，請看看BB和nleqslv包。

您可能還對包裝ktsolve中的一些嚴重扭曲的頭腦:-)編寫的「包裝器」工具感興趣。這最後一個工具可以讓您設置任意一組方程，併爲任何所需的一組變量進行反解。

來源

2014-11-22 14:25:32

我不知道這存在。棒極了！非常感謝！ – dreamon 2014-11-30 13:24:51

函數將線性方程轉化爲R中的矩陣形式？

回答

相關問題