私の問題は、タイトルにあるようvに、第 2 種修正ベッセル関数の導関数 ( に関して)を使用しようとしていますK_v(x)が、成功していません。
派生物として機能するドキュメントの1つを読みましたがbesselDK(v,x)、明らかにこれはRで認識されている関数ではありません。派生物の展開を使用しようとしました。
besselK(v,x)*(1- (1/2v) -log(e*x/2v))
しかし、これは私に正しいプロットを与えるためにも機能しません。これを含む関数をプロットしようとしています。
P <- function(x) (1/2)*log(exp(1)/(2*pi*x^(2)))+(3*exp(1/x^(2))/(sqrt(2*pi*x^(2))))*besselK((1/x^(2)),1/2)*(log(exp(1)/x^(2)))
x <- seq(0.1,2,0.01)
plot(x, P(x), xlim=c(0,2), ylim=c(0,1.2), type="l")
上記のコードから、プロットとして直線を取得します。正しいプロットでは、それは 1 から 1.5 の間で曲がる曲線である必要があります。誰かが正しい方法を教えてくれませんか?