R-借助fft函数计算2个多维连续密度的卷积
我需要计算两个独立的多维连续密度的卷积,只由网格上的数值给出 为了理解如何计算卷积,我尝试使用R中的fft函数计算一些简单的例子,但是,它不能正常工作。例如,我想计算两个独立的标准标准rv的卷积,代码如下R-借助fft函数计算2个多维连续密度的卷积,r,fft,convolution,R,Fft,Convolution,我需要计算两个独立的多维连续密度的卷积,只由网格上的数值给出 为了理解如何计算卷积,我尝试使用R中的fft函数计算一些简单的例子,但是,它不能正常工作。例如,我想计算两个独立的标准标准rv的卷积,代码如下 #values from a standard normal density eps <- 1e-1 points <- seq(-5,5, by = eps) x <- dnorm(points) y <- dnorm(points) #padded
#values from a standard normal density
eps <- 1e-1
points <- seq(-5,5, by = eps)
x <- dnorm(points)
y <- dnorm(points)
#padded values x, y with 0 to use fft function
xPad <- c(rep(0, length(y)), x)
yPad <- c(y, rep(0, length(x)))
zPad <- fft(fft(xPad)*fft(yPad), inverse = 'TRUE')/length(fft(xPad)*fft(yPad))
plot(Re(zPad))
max(Re(zPad))
dnorm(0, sd = sqrt(2))
标准正常密度的值
EPS