返回函数的系数值

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请问：如果告诉一个关于y的函数，在R软件中如何返回y的j次方的系数？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：r软件 函数

函数形式是多项式吗？给个例子。

我要做的那个函数很麻烦，里面还有其他的变量，不过对于具体的实验来说它的值是可以求出来的。可以举个简单的例子：f(y)=1+ay+by^2+3y^3  ,a=2,b=3,怎样求f的各项系数的返回值：1,2,3,3,

qoiqpwqr 发表于 2012-2-13 21:34
函数形式是多项式吗？给个例子。

我要做的那个函数很麻烦，里面还有其他的变量，不过对于具体的实验来说它的值是可以求出来的。可以举个简单的例子：f(y)=1+ay+by^2+3y^3  ,a=2,b=3,怎样求f的各项系数的返回值：1,2,3,3,

I guess you can evaluate your functions at different values and get the coefficients by solving those linear equations.

1. > # In your example
   
2. > a <- 2
   
3. > b <- 3
   
4. > aFun <- function(y) {1 + a*y + b*y^2 + 3*y^3}
   
5. > # Highest order of your polynomial
   
6. > n <- 3
   
7. > # Here, we are going to evaluate the function
   
8. > # at 1,2, ..., n + 1
   
9. > # Independent variable at different orders
   
10. > ivs <- outer(c(1:(n + 1)), c(0:n), "^")
    
11. > # Values of functions at correcponding y values
    
12. > # e.g., aFun(1), aFun(2), aFun(3), aFun(4)
    
13. > funVal <- sapply(1:(n + 1), aFun)
    
14. > # Find the coefficients of the function by
    
15. > # solving the equations
    
16. > coeff <- solve(ivs, funVal)
    
17. > coeff
    
18. [1] 1 2 3 3

复制代码

The only thing you need to know before hand is the highest oder of your polynomial. I can not think of any good way to find it out in case you do not know it.

zhangyangsmith 发表于 2012-2-15 15:06
I guess you can evaluate your functions at different values and get the coefficients by solving thos ...

谢谢！  对于变量的不同取值最高次项的次数是不同的，不过我在matlab里找到了一个命令可以返回其系数值，函数中也可以有其他的变量（该变量的值必须事先求出）。