step1.m

分形插值的matlab程序（以一组随机数组为例）

function step1
af=input('af==');
oo=input('数组维数==');
z=rand(oo)
[m,n]=size(z);

subplot(2,1,1);
mazz=max(max(z))*3;
meshz(z);
view(3);
axis([0 n-1 0 m-1 0 mazz]);

x=0:100:(n-1)*100;
xl=x(n)-x(1);
y=0:100:(m-1)*100;
yl=y(m)-y(1);
a=(x(2:n)-x(1:n-1))/xl;
b=(x(n)*x(1:n-1)-x(1)*x(2:n))/xl;
c=(y(2:m)-y(1:m-1))/yl;
d=(y(m)*y(1:m-1)-y(1)*y(2:m))/yl;

dn=ones(m-1,n-1);
dn=dn*af;
cc=(z(1:m-1,1:n-1)-z(1:m-1,2:n)-z(2:m,1:n-1)+z(2:m,2:n)-dn*(z(1,1)+z(m,n)-z(1,n)-z(m,1)))/(x(1)*y(1)+x(n)*y(m)-x(n)*y(1)-x(1)*y(m));
bb=(z(1:m-1,1:n-1)-z(1:m-1,2:n)-dn*(z(1,1)-z(1,n))-cc*(x(1)*y(1)-x(n)*y(1)))/(x(1)-x(n));
dd=(z(1:m-1,1:n-1)-z(2:m,1:n-1)-dn*(z(1,1)-z(m,1))-cc*(x(1)*y(1)-x(1)*y(m)))/(y(1)-y(m));
kk=z(2:m,2:n)-bb*x(n)-dd*y(m)-dn*z(m,n)-cc*x(n)*y(m);
for j=1:m-1;
    for j0=1:m;
        yv=c(j)*y(j0)+d(j);
        jj=(j-1)*(m-1)+j0;
        for i=1:n-1;
            for i0=1:n;
                ii=(i-1)*(n-1)+i0;
                xv=a(i)*x(i0)+b(i);
                zt=bb(j,i)*x(i0)+dd(j,i)*y(j0);
                zz(jj,ii)=zt+cc(j,i)*x(i0)*y(j0)+dn(j,i)*z(j0,i0)+kk(j,i);
            end;
        end;
    end;
end;
mm=(m-1)*(m-1)+1;
nn=(n-1)*(n-1)+1;
subplot(2,1,2);
meshz(zz);
view(3);
axis([0 nn-1 0 mm-1 0 mazz]);
axis off
end