基于艾特肯加速抛物线算法分形图形的研究_第1页
基于艾特肯加速抛物线算法分形图形的研究_第2页
基于艾特肯加速抛物线算法分形图形的研究_第3页
基于艾特肯加速抛物线算法分形图形的研究_第4页
免费预览已结束,剩余1页可下载查看

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、    基于艾特肯加速抛物线算法分形图形的研究    姜卓睿段汕摘要:抛物线法和艾特肯法都是数值分析中在实数范围内方程求根的快速迭代算法。针对julla集的绘图原理。给出了 一种用抛物线法在复数范围内的迭代运算。然后用艾特肯法对其加速,最后将运算结果绘成分形图形的改进算法。关键词:分形抛物线法艾特肯法复数迭代逸逸时间算法tp311.11b:1002-2422(2007)03-0052-021抛物线法设已知方程f(x)=0的三个近似根xk,xk-1,xk-22,我们以这三点为节点构造二次插值多项式p(x),并适当选取p(x)的一个零点xk+1作为新的近似根

2、,这样确定的迭代过程称为抛物线法,亦称密勒(muller)法。在几何图形上,这种方法的基本思想是用抛物线y=p(x)与x轴的交点xk+1作为所求根的近似位置。抛物线法是超线形收敛的算法,其收敛阶p=1840,收敛速度比割线法更接近牛顿法。2艾特肯加速方法如果迭代序列收敛很慢,要达到要求的精度将使计算量很大,为此,需采用加速迭代收敛性的方法,艾特肯加速方法就是其中的一种,其迭代公式如下:3复平面上的艾特肯加速抛物线算法上述抛物线和艾特肯算法都是在实平面上实现的,计算一次,其结果只和初值有关,也不可能得到二维图形,本算法选择在复平面上实现,算法步骤如下:(1)获取图片大小信息:picx;pictu

3、re1.sealewidth:picy=picture1.scaleheight(picx和picy分别为图片横向和纵向的最大像素值),设定初始近似值的变化范围ss=1.5:ww.x=-ss:ww.y=ss:hh.x=-ss:hh.y=ss,同时设定初值row=0;(2)如果row< p>(3)col=0;(4)如果col< p>(5)预设抛物线迭代时的初值0为复数):上述算法中最大的问题是要先计算出f(z)=o的根值形算法最关键的技术之一就是迭代,但对任意一函数迭代的结果只有两种可能:收敛或逃逸(发散)如图2所示,假设有一个充分大的整数n,当未述逸区域m中的初始点a经

4、过小于n次迭代就达到未逃逸区域m的边界,甚至超出了边界,就认为a点逃逸出去了:而经过n次迭代后a的轨迹仍未达到m的边界,就认为a是收敛区域a上的点。用这种方法描绘出收敛区域a的边界图形,这便是逃逸时间算法的基本思想。运用这种思想,不必计算f(z)=o的根值,只要事先确定整数n和逃逸区域m的边界范围即可。则算法基本框架不变,只是将步骤(9)做如下变更:(9)如果m为事先设定的逃逸边界值,若满足条件,则根据此结果设计颜色在屏幕上打点,并终止迭代转步骤(11);否则执行步骤(10):采用此种方式绘制的分形图形无穷无尽。再也不必因为找不到f(z)=o的根值z*烦恼。4结束语上述分形图像有着非常明显的几何意义,为复平面上的函数提供了新的解释,同时也给我们带来了美的享受。分形图像有着非常广泛的用途,例如可以制作成各种尺寸的装饰画(用卡纸

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论