榕坚

40# changxde
会不会是起始点及初始导数值的问题，或者那个图形是否经过某种变换？

changxde

f(z)=(z(^3)+a)/(k(z^2)+b)的IFS法作Julia集

榕坚

这个图形还是一直没办法弄出来：

xiaongxp

下面的图是用所谓“point trapping（点陷阱或点诱捕）”法作的，其代码为：
float iterate( complex c )
{
    float dist = 1e20f;
    complex z = 0
    for( int i=0; i < maxIterations; i++ )
    {
        z = z*z + c;
        if( |z|>2.0f ) return 0.0f;
        dist = min( dist, |z-point| );
    }

    return dist;
}
不知point指什么？

下载 (48.39 KB)

2011-10-29 11:17

柳烟

下载 (8.75 KB)

2011-10-28 23:23

f(z)=(z3-0,99)/(0,01z3+0,999z2-0,999i), dargestellt auf [-8;8]x[-12;4].
这个图形，我用距离估计法，也弄不出来。

柳烟

上图，你按普通J 集造，造出的图相去太远

榕坚

45# 柳烟


形状大致一样，只是迭代次数要求较高。扫描速度非常慢。

榕坚

44# xiaongxp


是自己设定的陷阱点吗?这个代码看不懂。谁知道的解释一下。

榕坚

这个差不多了，不亚于扫一张大图啊。

dyk

45# 柳烟
我用距离估值或DEM法造J集时也很少成功，估计原因在导数上，由于迭代n次，一次一个函数，用那个函数的导数呢？我一直理解不了。