问题标签 [fixed-point-iteration]

我需要在 Julia 中找到多变量函数的不动点。

考虑以下最小示例：

理想情况下，我会使用像 Roots.jl 和 call 这样的包find_zeros(p -> p - example(p))，但我找不到用于多变量函数的类似包。我找到了一个叫IntervalRootFinding，但奇怪的是它需要 unicode 字符并且文档很少，所以我不知道如何使用它。

我正在阅读fix-pointSICP：

从上面的案例中，我了解到，一个本质while是抽象概念“尝试”

所以我可以用有效的函数抽象思维在 python 中编写迭代。

反思时try，不仅仅是“尝试是迭代中的一段时间”，它教会了我什么？

它可以在没有的情况下重新构建try，但return fixed_point(f, nex)直接返回。

那么为什么try在这里介绍SICP，我猜效率可能不是作者的主要考虑因素。

用 elisp 测试

它按预期工作。

似乎 returnfixed-point f next比带有 try 的内部迭代更干净一些。

SICP在这里的考虑是什么，打算教什么？

我newton-method在elisp中写了从Scheme示例中查找root

它有效，但观察扭曲的代码：

三个funcalls，我无法想象如果闭包的深度会更糟。

elisp 的问题是否有替代解决方案？（我猜它贬低了闭包）

我正在编写一个程序，该程序涉及RWS跟踪可变状态并生成一些日志。我的目的是定义一个计算来评估一些动作，收集随后的状态，并根据它在日志的开头附加一些东西。Writer最小的例子：

在这里，我通过在事件发生MonadFix之前写入日志来改变过去。act它完美地返回"1ab"。但是，如果我使用M来遍历状态，那么它会挂起：

这种行为对我来说很奇怪，我不明白为什么这个计算会发散。更难证明其合理性，因为prepend第二个变体中的状态不会在任何程度上改变状态。为什么这个程序不收敛？我可以做些什么来修复（inb4“hehe fix”）吗？

我知道我可以使用Statepart of解决它RWS，但出于某种原因，我想避免它。

我有一组集合，这意味着集合中的项目必须在实际开始之前完成。例如：

我希望每一行都包含递归之前的行。所以在这种情况下，它应该像这样结束：

我已经尝试生成一个变量并从空白创建集合，但我没有设法做到这一点。任何想法我怎么能做到这一点？

我希望能帮助我解决这个问题。我喜欢在这个系统中找到修复点。

我在 Matlab 中编写了如下代码：

当我使用此代码制作图形时，上述系统的解决方案：

这是图形

根据上图，系统的解收敛到不动点。有没有什么办法可以用Matlab数值计算上一个系统的不动点？

我是 Python 的初学者，我对这个任务有疑问：

• 编写一个函数，使用定点迭代找到用户数学函数的根。
• 使用此函数查找以下的根：x^3 + x - 1。
• 用迭代算法的步数画出根近似的依赖关系图。

这是我第一次使用 Python，所以我真的需要帮助。这是我的代码，但它不起作用：

我目前正在从事我的博士论文编码工作。本文基于 Arellano (2008) 处理选择性违约。从@quantecon notebooks 获得灵感。

我正在尝试使用 Numba 在 Python 中创建一个迭代一堆变量的函数，但是其中一个是有条件的。也就是说，根据自然状态的概率，迭代必须不同。我曾尝试使用 if/else 语句，但 Python 指出语法错误。该代码是我博士论文的价值函数迭代。这是我写的：

错误在最后一部分，我收到以下消息：

如果有人提供帮助，我将不胜感激！

谢谢 ：）

我是 Julia 初学者。我想使用 nlsolve 求解以下非线性方程。

然后，我收到以下错误。

我假设发生错误是因为w在解决问题的过程中某些组件变为负数。所以我将上面的方程修改如下，w在解题过程中使 的分量取正值。

但是，迭代没有收敛，所有结果都变成了 NaN。

在问题求解过程中，如何在不取幂函数的负基的情况下求解上述非线性方程？

我有大量的一维非线性不动点问题要解决，最有效的数值求解器是什么？我目前正在使用 scipy.optimize.fixed_point，运行 1000 个任务大约需要 17 秒。感谢您的任何建议。