问题标签 [chaos]

我正在尝试使用三阶辛积分器 Ruth 绘制 Henon-Heiles 哈密顿量的庞加莱曲面部分，但我遇到了一些问题，因为经过几个步骤后，变量开始取巨大值导致错误。我不明白为什么会发生这种情况，因为我使用的积分器应该是辛的。

我曾多次尝试编写此代码，但我总是遇到某种错误，所以如果有人能解释我做错了什么，我将不胜感激。

我用python写的代码是这样的：

我明白了：

“此语句不在任何函数内。（它遵循终止函数“plot_2_components”定义的END。）”

我试过在“结束”之前输入情节，但没有出现情节。假设变量具有分配给它们的值。

谢谢！

是否可以在实验中对单个动作进行多次迭代？

假设我想在每 600 秒后运行几次 terminate-db-pod。

下面是我的迭代程序的代码。我希望能够使用海龟图形来获取每个参数 (k)，并将方程输出绘制在其对应的 k 值上。如果我没记错的话，这应该创建一个 feigenbaum 图？我的问题是，如何让乌龟为每个 k 值绘制这些点，然后将它们连接到相邻 k 值的点等等？

用于对部署在 Openshift 上的 Vertx 应用程序进行混沌测试的任何工具的指针。Chaos Monkey 或任何其他工具会起作用吗？

我正在尝试创建混乱游戏。我有一个 if 语句，它调用不同的混乱事物并指出绘制它们的点。但我认为我的 if 语句有问题。

我不知道该尝试什么...

我希望这能调用我的课程并绘制我的观点。

我正在尝试使用 DynamicalSystems julia 包来绘制由以下方程控制的离散动力系统的轨道图：

y_{n+1} = 1- by_n^2

我可以让它绘制但是我似乎无法让初始条件起作用，我正在尝试使用y_0 = -0.68.

这是我的代码：

这是它生成的数字：

如果可能的话，有人能告诉我是否有办法使用 DynamicalSystems 获得分叉点。

我想通过使用本文第 81 页中描述的标准方法来数值计算洛伦兹系统的 Lyapunov 谱。

一个基本上需要整合洛伦兹系统和切向向量（我为此使用了龙格-库塔方法）。切向向量的演化方程由洛伦兹系统的雅可比矩阵给出。在每次迭代之后，需要对向量应用 Gram-Schmidt 方案并存储其长度。然后由存储长度的平均值给出三个 Lyapunov 指数。

我在 python 中实现了上面解释的方案（使用版本 3.7.4），但我没有得到正确的结果。

我认为错误在于 der 向量的 Rk4 方法，但我找不到任何错误......轨迹 x、y、z 的 RK4 方法工作正常（由图指示）和实施的 Gram-Schmidt 方案也正确执行。

我希望有人可以查看我的短代码并可能找到我的错误

编辑：更新的代码

上面的代码是根据 Lutz 的建议更新的。结果看起来好多了，但它们仍然不是 100% 准确的。

• 正确解 (2.16, 0.0, -32.4)

• 我的解决方案（1.3540301507934012，-0.0021967491623752448，-16.351653561383387）

正确的解决方案来自Wolf's Paper, p.289。在第 290-291 页，他描述了他的方法，该方法看起来与我在本文开头提到的论文中完全相同（论文，第 81 页）。

所以我的代码中一定还有另一个错误......

研究了混沌原理，找了一些开源项目，比如阿里巴巴开源的chaosblade ，vmware开源的mangle。

这些工具都是故障注入工具，对被测系统不做任何分析。

根据混沌原理，我们应该

1.首先将“稳态”定义为系统的一些可测量输出，表明正常行为。

2.假设这种稳定状态将在对照组和实验组中继续存在。

3. 引入反映现实世界事件的变量，例如服务器崩溃、硬盘驱动器故障、网络连接中断等。

4.尝试通过寻找对照组和实验组之间的稳态差异来反驳假设。

那么我们如何进行第 4 步呢？我们是否应该使用监控系统来监控一些主要的指标，以检查故障注入后系统的状态。

有什么好的建议或最佳实践吗？

所以在某个时刻，他从这个愚蠢的方程式开始

他开始谈论旅鼠以及增长率如何波动并可能导致非常不可预测的混乱

Pnext = R * P - R * P * P

简化为

Pnext = R * P(1-P)

威奇是一样的。

其中 R 是繁殖率，他说大多数时候小于 2 P 是当前人口 Pnext 是幸存的后代，我认为

所以下一次 P = P + Pnext

当我看到这个时，我想；等一下

没有办法为下一代提供一个正数。除非增长或人口为负数。

无论我尝试什么，一切都消失了。

我究竟做错了什么？

谢谢