0

我在 R 中使用严重依赖 fda 包的 CollocInfer 包来估计常微分方程的参数。我能够成功实施广义分析方法来获得这些参数估计值，但我无法弄清楚如何使用 Profile.LS() 函数的输出来评估 t0 时 ODE 的估计解，这将给出估计值基于数据的 ODE 的初始条件。

编辑我之前的例子太简单了，没有捕捉到问题。下面是修改后的例子。

简单示例假设我有数据 x，y，我使用 B 样条基础对其进行平滑处理

require(fda)
x <- seq(0, 25)
y <- 3 + sin(x) + rnorm(1, 0, 1)
b.basis <- create.bspline.basis(range(x), nbasis=8, norder=4)
smooth.y <- smooth.basis(x, y, b.basis)
yfd <- smooth.y$fd
yfdPar <- fdPar(yfd, 1, 1)
yfd1 <- smooth.fd(yfd, yfdPar)

在 x=0 处计算 yfd1 的值的（最简单的）方法是什么？

4

1 回答 1

2

您可以使用指定的参数值评估功能数据对象eval.fd

eval.fd(0,yfd1)
        rep1
[1,] 3.27759
于 2013-01-31T21:51:32.613 回答