问题标签 [derivative]

我正在尝试将参数拟合到使用optim()R 中的数据集。目标函数需要对方程进行迭代根求解，G以便预测值p将G（嵌套在目标函数中）的值带到 0（或接近 0 到可能；为了稳定性，我使用了 50 次二分法迭代）。

问题是：我真的更愿意为 包含一个分析梯度optim()，但我怀疑迭代函数不可能。但是，在我放弃分析梯度之前，我想在这里的每个人都运行这个问题，看看是否有我忽略的解决方案。有什么想法吗？

注意：在确定二分法之前，我尝试了其他求解根的方法，但所有非括号法（牛顿等）似乎都不稳定。

以下是该问题的可重现示例。使用提供的数据集和 的起始值optim()，该算法在没有分析梯度的情况下收敛得很好，但对于其他数据集和起始值，它的表现并不好。

谢谢！！

我正在研究python上的心电图信号，并希望从参数中捕捉心电图信号的振荡行为——自然频率、阻尼因子和输入信号，这些参数又是从语义挖掘中提取的。该算法在单点使用四次时间导数，表示为 x'、x''、x''' 和 x''''，其中 x 是我的傅立叶变换心电信号。我的问题是如何在频域中获取信号的时间导数？为此，我应该使用 python 中的哪个函数，如 scipy.misc.derivative 或 numpy.diff？

我对 Mathematica 相当陌生，目前正在使用它来研究四元数。在一个特定情况下，我想获得一个涉及四元数范数的导数，如下所示：

Mathematica 实际上能够计算四元数的范数，但有人知道它是否可以扩展范数函数，以便上述输入产生扩展的导数（//FunctionExpand 没有改变输出的任何内容......）。

任何提示表示赞赏，谢谢！

在 R 中，我想要一种方法来获取公式右侧的符号导数，其中可能包括交互项、平方项等。

例如，我希望能够对以下两个 [edit:three] 公式中的每一个求得右手边关于 x 的导数：

y~x+I(x^2)

y~x:z

编辑：y~x*z

我想要一个函数，当输入上述三个公式中的每一个时，分别返回 1+2x、z 和 1+z。

我尝试了以下方法：

有什么方法可以强制 R 将 I(x^2)（或类似地，I(x*z) 等）识别为 x^2（分别为 x*z），x:z 识别为 x*z（在数学意义上）和 x*z（在公式意义上）作为 x+z+x*z（在数学意义上）用于计算导数？

其次，有没有办法从 deriv() 获取输出并将其重塑为看起来像公式的右侧？特别是，我知道 D() 将缓解这个问题并以我想要的形式生成输出（尽管 D() 不能将公式作为输入处理），但是如果我想对多个变量进行导数怎么办？我可以通过对每个我想对其进行导数的变量反复应用 D() 来解决这个问题，但是最好简单地输入所有这些变量的字符串并接收适合放置的输出在公式的右边。

谢谢！

我正在寻找 R 中生存函数的导数并将其存储在一个新函数中。

到目前为止，这是我的代码：

我想在 x=1.0 点计算 Matlab 中以下函数的导数：

我试过了，但它返回 py = []

I have a long 3-columns table of numerical data. For example, a small piece of it looks like this:

To be clear, the first column shows the time(t) and the other two indicate some other variables both as functions of time, let say f(t) and g(t). I want to calculate the time-derivative of each function and save it in a new data file. I'm not sure how should I handle this numerically. Any idea?

您好，我在 C++ 中创建了一个小型运动模拟。我想向学生展示欧拉法、龙格-库塔法和中点法之间的差异。

但是当我切换到 Rungy-Kutta 模式时，我在某个地方犯了一个错误，所有材料点都消失了。我犯的错误是在solveRK4 方法中。附件中的 MinGW Developer Studio 项目，还有库文件夹放入 mingw 编译器目录。： http ://speedy.sh/CvDHj/LABO3.zip

当您按下“R”按钮时，当“E”到欧拉时，当“M”到中点时，它会切换到 RK4。

问题是我在solveRK4函数中哪里出错了？

结果（发布）如下所示：http://speedy.sh/h28VP/zad3.exe按 R 键会从屏幕上消失。

Punkt - Point Wektor - 矢量图

主要的：

结构：

自从我处理一些数学问题以来已经有一段时间了，我有点生疏了，如果我问一个愚蠢的问题，请客气。

问题：我有 n 对线，它们作为二维点数组保存在内存中，因此没有显式函数。我必须检查情侣身上的线条是否平行，这是一项非常容易的任务，因为检查它们的导数是否相同就足够了。

要在算法中做到这一点，我必须检查函数两点之间的直线斜率（我有），因为我不需要极高的精度，我可以使用简单的公式：

m = (y2-y1)/(x2-x1)

但很明显，这导致我遇到了 x2 = x1 的大问题。我无法为这种情况提供默认值......我该如何解决？

目前我有两个 numpy 数组：x并且y大小相同。

我想编写一个函数（如果存在，可能调用 numpy/scipy... 函数）：

其中是result一个x与_ _ _nyxy