我正在使用 GNU plot (gnuplot) 来绘制来自一个奇怪吸引子的集成的一些输出数据（参见 Wikipedia） - 因为这是一个奇怪的吸引子，所以线“在 3D 空间中围绕自身滚动”。

这使得在绘制大量数据时很难看到会发生什么。看这张图片：

有什么方法可以绘制从开头的红色到结尾的蓝色的过渡或褪色颜色？ 我知道 gnuplot 是一个强大的工具，但也许我问的太多了？

我有3个方程如下

我尝试使用 scipy.integrate.odeint 从 t=0 解决到 t=10 我的初始条件是 u(0)=0, v(0)=1.0 和 w(0)=2.0 我没有太多有用的注释scipy.integrate.odeint。所以可以使用一些建议来解决我的问题。

我正在尝试在 MATLAB 中模拟 Lorenz 系统中的轨迹，目前使用以下代码 -

函数 g、g2 的定义方式相同：

这就是洛伦兹系统。整个代码的目的是制作两个初始状态的轨迹变化非常微小的电影，以展示该系统的混沌行为。代码本身确实有效，但占用了我计算机的所有内存，并且在尝试制作轨迹的 .avi 文件时，它抱怨超过 7.5 GB——这对于这个模拟来说当然太多了。

我的问题包括两部分：

(1) 我该如何管理这段代码以使其运行更顺畅？

(2) 如何制作轨迹的.avi 文件？我在网上找了很久，但是每次都是MATLAB或者我的电脑都放弃了。

提前致谢！

我正在使用 MATLAB 绘制 Lorenz 吸引子，并且想知道如何将 XYZ 坐标导出到 3D 可打印文件！我在从 XYZ 坐标到表面时遇到了麻烦（我应该这样做吗？）。此外，即使将 Lorenz Attractor 视为“表面”对我来说似乎也不正确，因为我的理解是它更像是一条“路径”。

到目前为止我的代码：

我之前试过这个：

但没有任何工作。错误说 z 必须是一个矩阵，我有点迷茫。

我正在（数字上）使用不同的方法解决洛伦兹系统。我正在使用 matplotlib 绘制它，但我想要一种更好地区分点的方法。

例如：

假设要绘制的点存储在a具有以下形式的数组中

现在这些代码行

生产：

不是很描述，是吗？所以我认为绘制离散点会更好。所以这些：

生产：

但它不像我想要的那样具有描述性。我想知道绘制 Lorenz 系统的最具描述性的方法是什么。

目前我使用以下代码生成洛伦兹级数

但nolds.lyap_r产生无效值（我假设有效值为 0.91）

import nolds l = generate([1, 0, 0]) nolds.lyap_r(l, tau=0.1, emb_dim=5) 1.0030932070169234

知道我在哪里犯了错误吗？

我正在尝试使用四阶龙格库塔方法解决洛伦兹系统，其中

由于该系统不明确依赖时间，因此可能会在迭代中忽略该部分，所以我只有 dX=F(x,y,z)

结果应该是 Lorenz atractors，但是我的代码在第五次迭代前后出现了分歧，这是因为我创建的内容是kconstants做的，但是我检查了，我很确定 runge kutta 的实现没有错......（至少我想）

编辑：

找到了一个类似的帖子，但无法弄清楚我做错了什么

我正在研究洛伦兹系统。我使用 'scipy.integrate.odeint' 内置函数按照 Wikipedia 1的建议进行集成。Lorenz 系统具有三个变量：x、y、z。当我比较具有相同初始条件、相同时间差（dt）但不同时间设定点的两个洛伦兹系统的 x 演化时，我得到了不同的集合。这两个系统在一段时间内发展相似，但后来出现了分歧2。是什么原因？唯一的区别是它们具有不同的时间集（但具有相同的时差）。时间点设置如何影响集成？

我还附上了我为上述代码获得的图表：x 在这两个系统中的演变。

1. 您好，我必须编写一个 python 函数来使用 Runge-Kutta 2cond 等级求解 Lorenz 微分方程

sigma=10, r=28 和 b=8/3

初始条件 (x,y,z)=(0,1,0)

这是我写的代码，但它给我一个错误，说在 中遇到溢出double_scalars，我看不出程序有什么问题

所以，我已经在 Mathematica 中看到了我的问题的编码解决方案，但对数学知之甚少，我还没有能够重现它。

这就是我想用 Python 做的事情：https ://mathematica.stackexchange.com/questions/159211/how-to-make-a-bifurcation-diagram-of-the-lorenz-system-under-a-变参数

我认为我的错误在于理解如何计算我正在寻找的内容以及如何调整我的可视化以使其看起来像链接中的那样，但欢迎任何想法。

到目前为止，我的代码如下所示：