问题标签 [differentiation]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
1108 浏览

maple - 如何强制 maple 进行链式微分?

在微分函数时,我通常不清楚,在哪些情况下 maple 执行链微分,何时不执行。

让我们看一个例子:

然后A(a,b)给出只是

问题:为什么 Maple 不通过应用链式规则的定义来执行微分?我怎样才能让枫树这样做?

更令人费解的是,在这个更简单的示例中,maple 的行为如我所愿:

然后A(a,b)返回

也许这有助于解决问题......

0 投票
4 回答
4120 浏览

python - 可在 Python 模块 Sympy 中以矩阵形式使用的微分运算符

我们需要两个微分算子矩阵,[B]例如[C]

这也可以应用于计算矢量场的旋度,例如:

为了使用 Sympy 解决这个问题,必须创建以下 Python 类:

当微分算子矩阵在左边相乘时,这个类单独解决。diff仅当要区分的功能已知时才执行此处。

当微分运算符的矩阵在右侧相乘时,要解决此问题,必须按以下方式更改__mul__核心类中的方法:Expr

它工作得很好,但是 Sympy 中应该有更好的本地解决方案来处理这个问题。有人知道它可能是什么吗?

0 投票
1 回答
152 浏览

matlab - Matlab不处理第二个for循环

我正在尝试通过将收集到的一些加速器数据与也记录的 LVDT 位移进行比较来验证这些数据。为此,我试图区分 LVDT 数据两次以获得加速。但是,当我运行我的代码时,它有错误

对此进行调查,我发现 Matlab 没有处理第二个 for 循环,因此从未生成变量 DiffDiffLVDT。

为什么它会跳过第二个 for 循环?

0 投票
1 回答
914 浏览

python - 避免在 Python 模块 Sympy 中对 args 进行排序

我目前正在开发一个可以以矩阵形式放置的微分算子。sympy在这种情况下,args创建Mul对象时列表的顺序对于保证仅在需要的地方执行区分非常重要。

问题是,当完成以下操作时:

它返回(t, x, z)是因为args发生了一些重新排列。如何避免args被排序?

0 投票
2 回答
17329 浏览

java - Java - 使用 Apache Commons 数学库计算推导

我在使用 apache commons 数学库时遇到问题。
我只想创建像 f(x) = 4x^2 + 2x 这样的函数,我想计算这个函数的导数
--> f'(x) = 8x + 2

我阅读了有关差异化的文章(http://commons.apache.org/proper/commons-math/userguide/analysis.html,第 4.7 节)。
有一个我不明白的例子:

在第 5 行,当然会发生编译错误。该函数f(x)被调用但未定义。我做错了什么?
有没有人对apache commons数学库的微分/推导有任何经验,或者有人知道另一个可以帮助我的库/框架吗?

谢谢

0 投票
3 回答
461 浏览

matlab - matlab自动微分

如果你们中的某个人可以帮助我或指出我正确的方向,那就太好了。

我有以下公式,其中每个主题有 8 个不同的参数。该等式描述了一种增长模式,即t时间和y增长。m1tom8是每个主题的不同参数

我想做的是插入这个公式并计算(数字)速度曲线。从这些信息中,我想获得以下信息

  • 当时的最大速度和相应的年龄和身高
  • 达到最大速度之前的最小速度。

这样我就可以列出我所有的科目及其参数,然后就可以收到我的结果了?

我完全没有编程和数学方面的背景,但我希望能够做到这一点。我有 matlab 供我使用。我试图弄清楚事情是如何运作的,但我无法抓住它。作为一名生物学家,我明年将学习编程课程;-)

你们中的任何人都可以帮助我吗?

谢谢

更新但仍然不存在(感谢你们:-))

-我制作了一个从 excel 派生的制表符分隔文件,并具有以下结构列 A:主题 B 到 I 列的 ID:不同参数 m1 ... m8 的值(如等式中)每行是一个不同的主题有不同的参数

  • 我选择了导入选择和下一个生成功能(或者我应该选择生成脚本)。

  • 如果我打开它,我会看到一个新的编号选项卡出现。第一行是 afunction[ID,m1,...,m8]=importfile。然后按照 112 行。

  • 我已经复制了 dan 从第 113 行开始给出的文本(感谢 dan)。

    tmin=0;tmax=20;dt=1/12;t=timn:dt:tmax;y = m1。( 1-1./(1+(m2.(t+m8)).^m5+(m3.(t+m8)).^m6+(m4.(t+m8)).^m7 ) );dy=diff(y)./dt; 最大(dy);分钟(dy);imax=find(dy==max(dy))+1; imin=find(dy==min(dy))+1; t(最大);t(最小);y(imax); y(imin);

    • 接下来我点击运行并被要求保存我已经完成的文件。
    • 现在我再次单击运行,但出现错误

谁能指出我正确的方向?非常感谢

问候

0 投票
1 回答
893 浏览

python - sympy 中共轭的导数

当我尝试用 SymPy 区分符号时,我得到以下信息

当我区分符号与其共轭时,结果是

但是,当我尝试区分符号 SymPy 的共轭时,它并没有做到这一点

这仍然是正确的,但结果应该为零。如何让 SimPy 执行共轭的导数?

0 投票
1 回答
526 浏览

maxima - 绝对值函数在 x=0 时是否可微

Wolfram Alpha ( diff abs(x)) 和 Maxima ( diff(abs(x),x,1);) 都说绝对值函数在 x=0 处是可微的,并且导数是 x/abs(x)。

在制定查询时,您将如何从域中排除点 x=0?谢谢!

PS 这发生在 Maxima (build_info("5.27.0","2012-05-09 21:24:21","x86_64-unknown-linux-gnu","GNU Common Lisp (GCL)","GCL 2.6. 7"))。

0 投票
1 回答
276 浏览

matlab - 如何在 Matlab 中区分和评估?

我需要区分exp((s^2*sigma^2)/2 + mu*s)和评估它s=0

任何人都可以就语法或我应该如何进行这种区分提供建议?

如果有帮助,上述函数是标准法线的 mgf。

我想区分和评估s=0,以便获得均值、方差、偏斜和峰度。

谢谢!

0 投票
1 回答
576 浏览

python - Sympy 的导数离散化

有没有办法在 sympy 中离散化未知函数的导数?我正在努力实现以下目标: