R 에서 미분계수 구하는 법

R 에서 미분계수 구하는 법

개요

R 에서 미분계수를 구하기 위해선 numDeriv 패키지의 grad() 함수를 사용할 수 있다.

코드

예를 들어 $f(x) = x^2 + 4x + 1$ 과 $g(x) = e^{-x}$ 의 미분계수는 다음과 같이 구할 수 있다.

install.packages("numDeriv")
library(numDeriv)
 
f<-function(x) {x^2 + 4*x + 1}
g<-function(x) {exp(-x)}
 
grad(f,2)
grad(g,0)

20190320\_151607.png

실제로 계산해보면 $f'(2) = 2 \cdot 2 + 4 = 8$ 이고 $g'(0) = - e^{0} = -1$ 인 것을 확인할 수 있다.

참고로 스칼라 함수의 경우에도 x 옵션에 벡터를 넣어주면 그래디언트를 잘 계산해준다.

같이보기

댓글