问题标签 [splines]

我正在尝试使用默认数据集（汽车）找到适合的 MSE smooth.spline（R并将其与其他方法进行比较）。但是使用predict函数会减少我的数据点的数量。换句话说，我有 50 对数据点 (x,y) 但predict函数给了我 35 个点 (yhatsp)。如何获得样条曲线的全部 50 分？谢谢

库（数据集）

x=汽车[,2]

y=汽车[,1]

yhatsp=predict(smooth.spline(x,y))$y

MSE=mean((y-yhatsp)^2)

csaps()在 matlab 中，根据平滑参数的特定定义进行三次样条p。这是一些matlab代码及其结果：

我想在 R 中得到相同的结果。我试过base::smooth.spline()了，但是平滑参数spar以不同的方式指定，我似乎无法与 matlab 相关p（你可以吗？）。我能得到的最接近的结果是包的smooth.Pspline()功能。pspline这是一些让事情在 R 中滚动的代码：

csaps() 帮助页面在这里

smooth.spline()可以在这里spar找到帮助（未给出代码，因为我认为和之间的关系可能p很复杂，所以可能不值得走这条路）

pspline::smooth.Pspline()帮助在这里

这个其他人从 2008 年开始的任务似乎没有得到答复，让我觉得自己像这个人。

R充满了样条函数，所以如果你们中的聪明人可以指出我与matlab做同样的事情csaps()（或沿着这些路线的技巧），我会非常感激。

[EDIT 19-8-2013] spar需要指定为(1-p)/p （而不是1/p），然后结果将与数值精度一致。请参阅下面的答案。

我正在拟合自然样条模型，但出现以下错误。您能帮我理解并解决以下错误吗？

Error in double(nx * ncoef) : vector size cannot be NA

In addition: Warning message: In nx * ncoef : NAs produced by integer overflow

不幸的是，我无法上传数据集，不幸的是，您无法复制我的代码。相反，我会尽可能多地解释我的代码：

下面是我运行以适应自然样条拟合的代码：

在我的数据集中，x 是一个从 1 到数据行数 424742 的序列。

我知道这个问题有点模糊，但是在花了很多精力之后，我没有运气弄清楚这个错误是什么以及如何解决它。

非常感谢，

我有兴趣使该图中的所有边都笔直。

据我了解，将样条设置为 false 应该这样做。

在一种情况下，一条边没有样条曲线，并且正好越过另一个节点的顶​​部。

在类似 html 的标签的情况下，会出现样条线，导致边缘看起来像一个大循环。

有没有办法强制“无样条线”？

[在此处查看输出图像。]

提前感谢您的任何帮助或建议！

让我看看我是否可以在没有图片的情况下充分描述这个问题。

假设我有两个变量都线性影响速度，a 和 b。随着 a 的增加，速度线性增加，反之亦然。b- 也是如此，随着它的增加，速度线性增加，反之亦然。现在假设已经确定了随时间变化的速度，并且是一个很好的平滑样条曲线，称为 v(t)。我们还知道给定时间 v = f(a,b)，其中 f 是一些基本的线性函数，可根据 a 和 b 确定 v。最后，我们有一些成本函数 c(a,b,t)，它定义了特定时间的成本以及 a 和 b 的值。

我要做的是绘制一个成本最小化样条曲线，其中每个点在特定时间确定 a 和 b，硬条件是 f(a,b,t) = v(t) 始终存在，并且我们试图最小化 c(a,b,t) 的软条件。如果您将其展平为二维，a 在一个轴上，b 在另一个轴上，在特定时间，您会看到为了满足硬约束，我们必须在 ab 平面中有一条线be on，但是我们应该在哪条线上放置一个点取决于成本函数。

如果成本函数很简单，这将是一个相对容易解决的问题——在每个 t 处，只需确定 a 和 b 以最小化成本，然后就可以完成。但是，成本可能会在保留边界处突然变化（例如，对于 t >= 5，a < 0.6 的成本会急剧增加），我希望我的样条能够预测到这一点并在我们到达 t = 之前开始增加 a 5、这样才能顺利进行。

对我来说，它崩溃的地方是我能找到的所有样条公式都需要 n 空间中的固定点。他们可能不会通过这些点，但他们确实需要它们。我的案例不需要样条通过 [a,b,t] 中的特定点，但确实需要它们通过特定 t 值的线（其余的是最小化）。有没有办法通过查看导数等将这个问题简化为基本样条问题？

本文描述了如何解决类似的问题，但它似乎要求样条曲线通过点而不是线进行最佳拟合。http://www.cs.berkeley.edu/~ravir/dspline.pdf

感谢您的任何帮助，您可以提供。

I'm making a simple 3D racing game for my project. I'm wondering about using the Catmull-Rom splines.

I've set the texture of it. I've also created a functional which takes 4 points and t (the distance between the 2 nearest control points). I've initialized the spline to take 4 CVector3f's, but now i don't know what to do, i.e. how to render it, how to create curves...

我已经计算了散点图的最佳拟合曲线，我想将结果绘制为平滑曲线，类似于 SAS 的样条曲线。

经过一番谷歌搜索后，我发现我应该先在我的数据上使用 interpolate.interp1d，然后再绘制线条。但是，当我尝试根据文档中的教程执行此操作时出现错误。提前感谢您的任何帮助或资源！

我尝试为Petzold Book Block绘制曲线参考

我尝试在 Windows Phone 8 中使用它。这些库不支持StreamGeometry命名空间中的类System.Windows.Media。

如何尝试在 WP8 上绘制它？谢谢

来自 EDIT 的解决方案现在作为答案发布。

老问题

我想将 的样条插值包含Eigen::Spline到一个更大的公式中，并希望借助 来确定该公式的导数Eigen::AutoDiff。

我尝试了以下代码：

遗憾的是，代码的编译给出了以下错误消息：

我怎样才能避免错误？在此先感谢您提供任何有用的提示。

特征版本是 3.2.1。对于 Eigen 的 3.1 和 3.0 版本，编译也会失败。

编译器版本为：

gcc version 4.8.1 (Ubuntu/Linaro 4.8.1-10ubuntu9)

我想在 Java Swing 应用程序中实现 B 样条。从 4 岁开始。我可以尝试几乎所有事情，并度过我空闲时间的每一分钟，但我没有把它们做对 - 从昨天开始我就头疼：/

为了实现 B 样条，我使用了来自wikipedia的规范

我制作了（与MCVE相关的）一个演示版本，可以在这里找到：gist.github.com/soraphis/b-spline

该代码显示了该问题，它可以在摆动 JPanel drawcomponent 方法中调用。但是您可以评论这一行并取消评论 71。

一些附加信息：

• 我的基函数中的 a 和 b 返回值 < 0 或 >1，这不应该是（参见维基百科）
• 我想自己实现 b 样条曲线 - 我不想使用库
• 参考维基百科文章：
  • 基础函数是 B(x)
  • DeBoor 是 S(x)

我真的需要关于基函数的帮助，我想知道如何“正确”建立结向量

感谢您的各种回复
，感谢您阅读本文