[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小