2013-08-28 5 views
8

मैं एक बिवाएरेट फ़ंक्शन के संख्यात्मक डेरिवेटिव ढूंढना चाहता हूं।मनमाने ढंग से परिभाषित फ़ंक्शन के संख्यात्मक डेरिवेटिव

  • समारोह अपने आप से परिभाषित किया गया है
  • मैं प्रत्येक तर्क के संबंध में पहले डेरिवेटिव की जरूरत है और दूसरा व्युत्पन्न

वहाँ एक पैकेज या निर्मित समारोह यह करने के लिए है पार?

उत्तर

11

numDeriv पैकेज को स्थापित और लोड करें।

library(numDeriv) 
f <- function(x) { 
    a <- x[1]; b <- x[2]; c <- x[3] 
    sin(a^2*(abs(cos(b))^c)) 
} 
grad(f,x=1:3) 
## [1] 0.14376097 0.47118519 -0.06301885 
hessian(f,x=1:3) 
##   [,1]  [,2]  [,3] 
## [1,] 0.1422651 0.9374675 -0.12538196 
## [2,] 0.9374675 1.8274058 -0.25388515 
## [3,] -0.1253820 -0.2538852 0.05496226 

(मेरे उदाहरण बल्कि द्विचर से trivariate है, लेकिन यह स्पष्ट रूप से और साथ ही एक द्विचर समारोह के लिए काम करेंगे।) कैसे ढाल और विशेष रूप से हेस्सियन संगणना किया जाता है के बारे में अधिक जानकारी के लिए सहायता पृष्ठ देखें।

संबंधित मुद्दे