给出的答案

-4

我试图解决逆问题以下功能R.给出的答案

x + 2 (C1 * y) + C1 * C1 * z = d2

我目前可以进入C1并获得d2但需要输入d2并获得当查找函数的系数C1。变量x,y和z都是已知的并且永远不会改变。

我已经有一些已知的使用C1和d2值。

是否有R功能，这将让我进入功能，以前的结果和d2值和它返回C1系数？

来源

2012-07-16 TrueWheel

代数是1000多年前发明的，用于解决这类问题。 – Andrie 2012-07-16 15:39:30

@Andrie我确定这是。我也认为有人会创建一个R包来帮助解决这个问题。 – TrueWheel 2012-07-16 15:44:19

嗯，我想你可以尝试[聊天[（http://chat.stackoverflow.com/）这种问题。 – 2012-07-16 15:49:40

你有一个二次方程的f ORM：

(x - d2)*C1^0 + (2*y)*C1^1 + (z)*C1^2 = 0

在R中可以解决与功能polyroot()二次方程式（实际上任何多项式方程）：

x <- 1 
y <- 2 
z <- 3 

d <- 0 

polyroot(c(x-d, 2*y, z)) 
[1] -0.3333333+0i -1.0000000+0i

（可提供两种解决方案，正如你所期望）

为了解决一系列的输入值，你需要把它放到你最喜欢的apply函数中，在这种情况下为sapply()：

d <- seq(0, 1, 0.2) 

sapply(d, function(dd)polyroot(c(x-dd, 2*y, z))) 

       [,1]   [,2]   [,3]   [,4]   [,5]   [,6] 
[1,] -0.3333333+0i -0.2450296+0i -0.1722534-0i -0.1088933-0i -0.05203037+0i 0.000000+0i 
[2,] -1.0000000+0i -1.0883037+0i -1.1610799+0i -1.2244400+0i -1.28130296+0i -1.333333+0i

来源

2012-07-16 16:03:43 Andrie

+1代数技能。 – 2012-07-16 16:12:37

谢谢。 :)我想我需要刷上我的代数。 – TrueWheel 2012-07-16 16:16:00

你有

d2 = x + 2 C1 y + C1^2 z

，你可以重新排列，或者只是将其插入Wolfram Alpha得到

z C1^2 + 2 y C1 + x - d2 = 0

这是C1，它可以解决或者使用二次公式quadratic equation得到

C1 = (-sqrt(d2 * z - x * z + y^2) - y)/z

来源

2012-07-16 15:50:44

回答

相关问题