如何向量化初始条件

时间:2019-10-07 17:26:39

标签: vectorization octave

我正在尝试对初始条件(x)(y)进行矢量化处理,以在不同的初始条件下绘制多个奇怪的吸引子,同时让每个初始条件都有单独的颜色图。但是我不确定如何对初始条件进行矢量化处理,以及如何使每个初始条件图成为单独的颜色。

我的想法是让 x y 作为向量x = [7.4,6.3,3.4,5.4,4.5] y = [。6,.3, [8,.3,.14]以及每个x和y初始条件都使用不同的颜色创建自己的图(但这不起作用,因为它希望一个值而不是多个值)。参见下面的x(k)和y(k)初始值图。

plot

以下工作代码未矢量化:

x(1)=[7.4];
y(1)=[.6];
a = -1.24458;
b = -1.25191;
c = -1.815908;
d = -1.90866;

maxiter = 10000;
z = linspace(1,maxiter+1,maxiter+1);
for k=1:maxiter

% Clifford Attractors

  x(k+1) = sin(a*y(k)) + c*cos(a*x(k));
  y(k+1) = sin(b*x(k)) + d*cos(b*y(k));
end

scatter(x,y,'.')

Ps:我正在使用类似于Matlab的64位Ubuntu Linux Octave 4.2.2

1 个答案:

答案 0 :(得分:0)

您可以只使用矩阵而不是向量,并将每一行视为一次迭代。为了稍微加快速度,我们可以在开始时分配矩阵。要绘制不同的运行,我们只需要遍历它们并每次使用不同的颜色即可:

b.txt
相关问题