我正在尝试canvas,并且试图修改这段代码,但不幸的是,我不了解其中的某些部分。
我的问题是-如何自定义上述代码,例如通过
f(z) = c^e(-z)
(公式取自一本带有分形示例的书)?
我知道我需要更改这部分代码:
function computeRow(task) {
var iter = 0;
var c_i = task.i;
var max_iter = task.max_iter;
var escape = task.escape * task.escape;
task.values = [];
for (var i = 0; i < task.width; i++) {
var c_r = task.r_min + (task.r_max - task.r_min) * i / task.width;
var z_r = 0, z_i = 0;
for (iter = 0; z_r*z_r + z_i*z_i < escape && iter < max_iter; iter++) {
// z -> z^2 + c
var tmp = z_r*z_r - z_i*z_i + c_r;
z_i = 2 * z_r * z_i + c_i;
z_r = tmp;
}
if (iter == max_iter) {
iter = -1;
}
task.values.push(iter);
}
return task;
}
但是z_i,z_r,c_i,c_r的真正含义以及我如何将它们绑定到上面的公式上并不能。
任何帮助将不胜感激。
复数有两个部分:实数,虚数。
因此z = a + b*i
,哪里a
是真实的部分,b*i
还是虚构的。
在提供的示例中z=z^2+c
,其中z=z_r+z_i*i
注: i*i = -1
所以z^2 = (z_r+z_i*i)*(z_r+z_i*i) = z_r*z_r+2*z_r*z_i*i + z_i*i*z_i*i = z_r*z_r+2*z_r*z_i*i - z_i*z_i
现在添加c
:z_r*z_r+2*z_r*z_i*i - z_i*z_i + c_r + c_i*i
组它
z_r*z_r+2*z_r*z_i*i - z_i*z_i + c_r + c_i*i = (z_r*z_r - z_i*z_i + c_r) + (2*z_r*z_i + c_i)*i
因此,我们tmp
从代码中获取var-是new的真正组成部分z
tmp = z_r*z_r - z_i*z_i + c_r
和虚部
2*z_r*z_i + c_i
由于z = z_r + z_i * i
,我们需要分配
z_r = z_r*z_r - z_i*z_i + c_r
z_i = 2*z_r*z_i + c_i
更新:用于f(z) = e^z - c
第一,一些复杂的形式:x = a+b*i = |x|(cos(p)+i*sin(p)) = |x|*e^(i*p)
其中,|x| = sqrt(a*a + b*b)
和p = b/a
在我们的例子:p=z_i/z_r
,|z| = sqrt(z_r*z_r+z_i*z_i)
e^z = e^(z_r+z_i*i) = e^z_r * (e^z_i*i) = e^z_r * (cos(p)+i*sin(p)) = (e^z_r * cos(p)) + i * (e^z_r * sin(p))
减去c
:
(e^z_r * cos(p)) + i * (e^z_r * sin(p)) - c_r - c_i*i = (e^z_r * cos(p) - c_r) + i * (e^z_r * sin(p) - c_i)
那么新的z
z_r = (e^z_r * cos(p) - c_r) = (e^z_r * cos(z_i/z_r) - c_r)
z_i = (e^z_r * sin(p) - c_i) = (e^z_r * sin(z_i/z_r) - c_i)
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句