问题标签 [curve-fitting]

我需要在矩阵的每一行都使用傅立叶函数，避免使用for循环。我尝试使用单元格数组，但不知道如何配合使用它。我的程序是

它给出的??? Error using ==> cellfun cfit type is not currently implemented.可能单元格数组不是这个问题的解决方案。如何解决这个问题呢？谢谢

更新

工作没有错误，但 wyn 只有空单元格
我测试了一行

并z包含模型的好值
我在 cellfun 中的语法有问题，但不知道是什么

我有一个 3D 数据立方体，有两个空间维度，第三个是 2D 图像每个点的多波段光谱。

给定一个波长（或波段号），我想提取与该波长对应的二维图像。这将是一个简单的数组切片，如H[:,:,bnd]. 类似地，给定空间位置 (i,j)，该位置的频谱为H[i,j]。

我还想在光谱上“平滑”图像，以抵消光谱中的低光噪声。对于 band bnd，我选择一个大小的窗口wind并将 n 次多项式拟合到该窗口中的频谱。使用 polyfit 和 polyval 我可以在该点找到适合的光谱值 band bnd。

现在，如果我想要拟合值的整个图像bnd，那么我必须在每个(i,j)图像上执行这个窗口拟合。我还想要 的二阶导数图像bnd，即每个点处拟合光谱的二阶导数的值。

在这些点上运行，我可以对每个光谱进行 polyfit-polyval-polyder x*y。虽然这可行，但这是一个逐点操作。有没有一些 pytho-numponic 方法可以更快地做到这一点？

我有一组 XY 值（即散点图），我想要一个 Pascal 例程来生成适合这些点的 N 阶多项式的系数，就像 Excel 所做的那样。

我在 Matlab r2010b 中使用曲面拟合工具箱来曲线拟合 3 维数据。该工具是通过 GUI 进行 3d 曲线拟合的非常有用且简单的工具。到目前为止，我处理 3d 数据，例如；（X 输入：x（nx 1 向量），Y 输入：y（nx 1 向量），Z 输出：z（nxn 矩阵））

但是现在，我正在尝试处理 4 维数据，例如；（X 输入：x（nx 1 向量），Y 输入：y（nx 1 向量），Z 输入：z（nx 1 向量）和 T 输出：t（nxnxn 矩阵））。

我无法通过曲面拟合工具箱来实现这一点。如何使用 4 维数据进行曲线拟合？

将不胜感激任何帮助。

我正在尝试将下面的 matlab/octave 函数转换为 C（传统方式 - 了解 matlab 函数并从头开始用 C 编写代码）。它使用多项式拟合将数据拟合到高斯曲线。

但是当我检查函数 polyfit 时，这似乎需要做很多工作，因为它涉及到对更多 octave 库函数的大量调用。它首先计算 Vandermonde 矩阵，然后对其执行一些 QR 分解，并计算向量的范数等......

1. 我可以利用哪些其他选项/处理来获得类似的功能（上面发生的实际操作的近似值），但使用一些更简单的曲线拟合或插值方法。

任何指针都会很有用。

我想知道是否有一个用于将数据拟合到正弦曲线的类。我在这里找到了一个算法，但首先理解然后编码需要一些时间......在我走这条路之前，我想知道论坛是否知道已经为这种数学建模创建的类。

我一直在寻找一种在 Python 中编写代码片段的方法，该代码片段可以计算任何 n 阶傅里叶级数曲线拟合。要计算傅里叶级数曲线拟合的某个阶数，比如 3 阶非常简单，但是在阶数 n 可变的情况下进行计算，仍然不可行。也许有人已经做到了，但我的搜索还没有找到。我想知道是否有人可以提供帮助。谢谢。

我正在尝试使用 GSL 进行最小二乘拟合，但我什至无法使示例正常运行。它可以编译，但是在运行程序时会出现错误

dyld: lazy symbol binding failed: Symbol not found: _cblas_dnrm2

Referenced from: /usr/local/lib/libgsl.0.dylib

Expected in: dynamic lookup

我相信我已经将问题缩小到了这一行

gsl_multifit_fdfsolver_set (s, &f, &x.vector);

但我一点也不知道为什么。

我在 OSX 上使用 GSL 1.15 运行它（尽管我最初安装了错误的版本，1.9）。我使用以下命令在终端中编译它：gcc unchanged_example.c -o examplefitter -lgsl -lm

编辑：在这里找到了解决方案。编译时添加' -lgslcblas'整理出来

我正在尝试在我生成的散点图上生成线性回归，但是我的数据是列表格式，并且我可以找到使用的所有示例polyfitrequire using arange。arange虽然不接受列表。我已经搜索了有关如何将列表转换为数组的高低搜索，但似乎没有什么清楚的。我错过了什么吗？

接下来，我怎样才能最好地使用我的整数列表作为输入polyfit？

这是我正在关注的 polyfit 示例：

我想在 R 中对正常和双对数图中的数据进行线性回归。

对于普通数据，数据集可能如下：

在那里，我想计算仅为数据点 2、3 和 4 的线性回归绘制一条线。

对于双对数数据，数据集可能如下：

在这里，我想为数据集 1:7 和 8:15 绘制回归线。

我可以计算斜率和y 偏移量以及拟合参数（R^2，p 值）吗？

如何处理正常数据和对数数据？

谢谢你的帮助，

斯文