小编Jus*_*n K的帖子

我试图通过在R中使用傅立叶级数来实现2D拉普拉斯反演.我的代码在单个反演情况下工作但不用于双反演.任何帮助将非常感谢.谢谢!

下面是R中的代码,用于f(x,t) = exp(-3(x+t))具有双拉普拉斯变换的简单函数F(p,s) = 1/((p+3)*(s+3)).

# General implementation of 2D Laplace Inversion

# How to use this?
# 1. Change the Laplace Transform function (F)
# 2. Change the Original function(f_1)

# Double Laplace Transform function, parameters p and s
F <- function(p,s) {1/((p+3)*(s+3))}

# Adjust parameters here
a <- 0.1
n <- 1000
# Note that the approximation works accurately up to x=T/2 only
T <- 10

# First inversion - p
    F_1 <- …