我被弄糊涂了？

mjj_ljh

怎么回事，大家看看。

xiaongxp

应该是程序中的函数用多项式逼近而引起的误差累积吧。

mjj_ljh

2# xiaongxp


要是这样，这个算法要改，不知老外知不知道？

榕坚

按此说来，所有的数值运算在超过36次迭代后的结果就没有保证了。

xiaongxp

3# mjj_ljh
算法应该不是问题，因为计算机进行的函数数值运算都是用多项式逼近或三角逼近的（泰勒展开或傅里叶展开），只是保留的展开式的项数要增大，以满足迭代迭代的需要。

xiaongxp

难怪许多分式变换下的复变分形作不出来，看来不一定只是迭代次数达不到那么高，而是计算的误差所致。

柳烟

老外的此软件，精度最高到百万分位为止，就是造复分形，这精度都不够满足需要，老外能否在下个板本考虑考虑，帮我辈一个大忙？不过，发明老外的不懂汉语，也听不到我在这儿牢骚。

inRm

美国佬要是知道这儿有几位能用GSP把分形玩到这种程度，恐怕得跌跟斗了。

榕坚

就觉得奇怪，这样看来M集的扫描放大倍数10^16以内出现马赛克现象均为这种原因导致的结果（有些点提早收敛了，没有跑到最后的位置上去）。

mjj_ljh

这可能就是画板做分形的瓶井所在。歪打正着，这正好说明了混沌的初值敏感性。