2013-10-10 44 views
0

我有amatrix m其包括三次多项式的系数的不同值,并且我想找到这些三次函数的最大实根。我怎样才能做到这一点R中?应用一个矩阵的每一行用于polyroot中的R

 m = matrix(c(1, 3, 4, 6, 7, 8, 9, 2, -3, -5, 6, 2), nrow = 3) 
     colnames(m) <- c("d", "c", "b", "a") 

我想用polyroot为m矩阵中的每一行,我怎么能做到这一点的R'

回答

1

这是不是仅仅是:

apply(m , 1, polyroot) 

...记忆,当然,该结果将在列的顺序。

apply(apply(m, 1,polyroot) , 2, function(x) ifelse(abs(Im(x)) < 0.00001, x,NA)) 
      [,1]   [,2]   [,3] 
[1,]   NA -0.416256-0i -0.4167669+0i 
[2,]   NA   NA   NA 
[3,] 2.34748+0i   NA   NA 
+0

谢谢,但我怎么才能得到最大的真正的根? – rose

+0

嗯,这是不是原来的问题了,是吗? –

+0

不,soory,我忘了提到它。 – rose

相关问题