基于MatLab整体线性分形插值的研究.pdfVIP

2007年第3期 九江学院学报 (总第140期) No，3，2007 Journal of Jiujiang University (SumNO 14o) 基于MatLab整体线性分形插值的研究冰 袁利国 余荣忠2 凌和良。 聂笃宪 (1华南农业大学数学系 广东广州 510642；2九江学院理学院 江西九江 332005； 3南昌工程学院理学系 江西南昌 330099) 摘要：本文分析了整体线性分形插值原理，利用曲线上的有限个插值点，选取适当的垂 直比例因子，可以近似地构造出原曲线，在此逼近计算的过程中，给出了一种易实现的分形 插值误差计算方法，并基于MatLab说明了整体线性分形插值的实现过程。同时提出了粒子 群优化算法寻最优垂直比例因子的可能性。最后讨论了分形插值中的一些问题及发展状况。 关键词：分形插值；垂直比例因子；迭代函数系统 中图分类号：TP301．6；TP 391．41 文献标识码：A 文章编号：1673-4580(2007)03-0058一(04) 0 引言 数，其中N≥2，记K=I X R。记Ii=[x．．1，x．]，i= 传统的插值方法(Lagrange插值法、样条插值 1，2，…，N，令L．是I—I；的一压缩线性映射， 法等)对平稳数据和光滑曲线的插值拟合是有效 记：L．(x)=all+e；，且满足条件 的，但对剧烈震荡曲线与数据效果就不理想。1986 L_(】【0)=x．_1，Li(xN)=x； (1) 年，Barnsley在文【6 中引入分形插值函数，从而为 令F；是K—R上关于第二个变量是压缩的线 这类曲线与数据的插值拟合提出了一种新的思路。 性映射，记：Fi(x，Y)=CiX+dly+‘，即存qi∈[0， 分形插值是依据迭代函数系统(IFS)的吸引子而 1)，使得I F；(x，Y )一F；(x，Y2)I≤qi I Yl—Y2 I， 建立起来的，由于分形插值函数的盒维数可以介于 Vx∈I，Yl，Y2∈R，且满足条件 1到2之间的任意数值，因此，分形插值函数不仅可 )【0，y。) Y¨ (2) 以用来拟合光滑数据或曲线，更可以在震荡剧烈的 【 (xN，YN)=Yi 离散数据或曲线中显示其独特的优越性。这其中如 定义压缩仿射函数：wi：K—K： 何有效地估计其逼近的误差，如何分析垂直比例因 Wi(x )：fLi【x 1，i：1，2，…，N． (3) 子对插值拟合精度的影响，是本文所关心的问题。 、Fi(x，Y)， 1 分形插值原理 则{K：w；，i=1，2，…，N}构成一个迭代函数 1．1迭代函数系统(IFS)的随机算法¨ J 系统(简记为IFS)，其中a；，c；，d；三个数的绝对值 设{X；wl，W ，…，WN}是一个IFS，其中每一个 都小于1。下面定理1说明存在一个分形插值函数 n 与以上迭代函数系统对应，定理2进一步说明了分 W；具有一个概率Pio，且∑Pi=1。选取任一点 形插值函数的唯一性…。 xo∈X为初始点，然后递归地随机选取下述集合中 定理1 存在I上连续函数f，使得f的图像G 一 个点作为x ，于x ∈{W (x 一 )，…，W (xn一 )}， = Graph(f)={(x，f(x))I x∈Il是迭代函