用R语言求函数的导函数
看中文版《R导论》的时候,惊奇地发明R语言中有一个D()函数可以用来对函数求导函数,巧合的是前两天建宁兄也问到这个问题,不妨在这里写一写吧。
D()函数的用法很是简朴。好比我们要求函数 的导函数的话,我们可以先界说一个表达式(expression):
> fun=expression(x^2+2*x+3) #界说函数为表达式
> fun #查察表达式
expression(x^2 + 2 * x + 3)
接下来用D()函数对fun举办求导就可以了:
> D(fun,”x”)
2 * x + 2
也许,你会问,能求偏导数吗?那就试试看:
> fun=expression(x^2+y^2+x*y)
> D(fun,”x”)
2 * x + y
> D(fun,”y”)
2 * y + x
> D(D(fun,”x”),”x”)
[1] 2
> D(D(fun,”x”),”y”)
[1] 1
假如要求函数 在(2,3)处的关于 的一阶偏导数的值呢?
> x=2;y=3
> eval(D(expression(x^2+y^2+x*y),”x”))
[1] 7