R中无yacas方程的算术运算

Question

我定义了likef<-"x^2"and g<-"y^2"__方程。我想得到像z<-(x^2)*(y^2)__这样的z方程。

'Class(f)‘、'class(g)’和'class(z)‘值对我来说并不重要。

我试过这个：

> f<-"x^2"
> g<-"y^2"
> z<-f*g

我得到了：

Error in f * g : non-numeric argument to binary operator<br/>

我尝试将f<-expression(f)和g<-expression(g)相乘，但没有结果。

我也试过：

> f<-function(x) x^2
> g<-function(y) y^2
> z<-function(x,y) f*g
> z<-parse(text=z)

我得到了：

Error in as.character(x) : 
  cannot coerce type 'closure' to vector of type 'character'

用paste(z)代替parse(z)

> paste(z)
Error in paste(z) : 
  cannot coerce type 'closure' to vector of type 'character'

有没有一种方法可以在R中不使用像yacas这样的重软件来进行符号运算？

Answer 1

您可以尝试以下几种方法：

f <- expression(x^2)
g <- expression(y^2)
z <- expression(eval(f) * eval(g))
#> eval(z,list(x = 1, y = 2))
#[1] 4
#...and here's another test, just to be sure:
#> identical(eval(z, list(x = 17, y = 21)), 17^2 * 21^2)
#[1] TRUE

或者您可以使用rSymPy包：

library(rSymPy)
x <- Var("x")
y <- Var("y")
sympy("f = x**2")
sympy("g = y**2")
sympy("z = f * g")
#> sympy("z")
#[1] "x**2*y**2"
#> sympy("z.subs([(x,3),(y,2)])")
#[1] "36" 

第二个建议可能不是很有吸引力，因为它可能更容易直接使用Python。

Answer 2

你对符号和功能的使用非常令人困惑。如果您想使用语言对象的“片段”，您可以使用更类似的

f <- quote(x^2)
g <- quote(y^2)
z <- quote(f*g)
do.call("substitute", list(z, list(f=f, g=g)))
# x^2 * y^2

如果您只想要某种形式的函数组合，只需确保您将变量视为函数。

f <- function(x) x^2
g <- function(y) y^2
z <- function(x,y) f(x)*g(y)

但这些只是构建R语言元素的方法。把这些看作代数运算是不明智的。R不会试图简化或类似的东西。