分形插值曲线函数

%分形插值曲线函数 %x为插值点，五个插值点(1,2),(2,2.5),(3,4),(4,6),(5,3) % n为迭代次数 % example1; 来自于：周承新,陈慧琴;分形插值曲线的MATLAB实现[J];科技广场;2009,3:119-120. % 但有所改动，没有初始输入点（本人觉得不合适） % xy=[1 2;2 2.5;3 4;4 6;5 3]'; % n=4; % example2: 《分析与拟合》；作者：沙震，阮火军；出版社：浙江大学出版社；出版年月：2005年3月；129页，例5.1 % xy=[0 0;1/3 1;2/3 1;1 0]'; % d=ones(1:3

function [XY]=fractal_curves（xyn）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%分形插值曲线函数
%x为插值点，五个插值点（12）（22.5）（34）（46）（53）
% n为迭代次数

% example1;  来自于：周承新陈慧琴;分形插值曲线的MATLAB实现[J];科技广场;20093:119-120.
% 但有所改动，没有初始输入点（本人觉得不合适）
% xy=[1 2;2 2.5;3 4;4 6;5 3]‘;
% n=4;

% example2:  《分析与拟合》；作者：沙震，阮火军；出版社：浙江大学出版社；出版年月：2005年3月；129页，例5.1
% xy=[0 0;1/3 1;2/3 1;1 0]‘;
% d=ones（1:3）/3;
% n=3;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    
oldXY=xy;
 
col=length（xy）;
    
d=rand（1col-1）/2;


L=xy（1col）-xy（11）;
a=zeros（col-11）;  
e=zeros（col-11）; 
c=zeros（col-11）;
f=zeros（col-11）;
    

for j=1:col-1  
    a（j）=（xy（1j+1）-xy（1j））/L;
    e（j）=（xy（1col）