分形插值

          [x]    [ai  0][x]    [bi]

     wi [y] = [ci si][y] + [di]

delx = x(n) - x(1);
s = [0 s];

a = [0];
b = [0];
c = [0];
d = [0];
for i = 2:n
   a(i) = (x(i) - x(i-1))/delx;
   b(i) = (x(n)*x(i-1) - x(1)*x(i))/delx;
   c(i) = (y(i) - y(i-1) - s(i)*(y(n) - y(1)))/delx;
   d(i) = (x(n)*y(i-1) - x(1)*y(i) - s(i)*(x(n)*y(1) - x(1)*y(n)))/delx;
end

插值维数 为D=1+log(|s1|+|s2|+|sN|)/logN

若s1=s2=sn,D与N s有关

不同的s即纵向变化因子，则不同的波形

知道s越大，波动越大。

s比1小