function?[]=graytwothree11（i）;
i=imread（‘11111111.bmp‘）;
figure（1）
imshow（i）
title（‘ori?image‘）
height=size（i1）;
width=size（i2）;
scale=1.5;
scale1=1.5;
H=zeros（uint16（scale*height）uint16（scale1*width））;
h1=1/scale;
w1=1/scale1;
i=double（i）;
for?h=3*h1:h1:height-2*h1
????for?w=3*w1:w1:width-2*w1
????????a=fix（w）;
????????b=fix（h）;
????????j=w-a;
????????k=h-b;
????????u1=-1/2*（1+j）^3?+?5/2*（1+j）^2?-?4*（1+j）+2;
????????u2=3/2*j^3?-?5/2*j^2+1;
????????u3=3/2*（1-j）^3?-?5/2*（1-j）^2+1;
????????u4=-1/2*（2-j）^3?+?5/2*（2-j）^2?-?4*（2-j）+2;