function y=IWT(x)
x=imread('lenna.gif');
subplot(2,3,1);imshow(x);title('Ô*ͼÏñ');
x=double(x);
y=x;
[hp,lp]=size(x);
hc=hp/2; lc=lp/2;
%----------------ÆæżÁÐÖØÅÅ
for n=1:lc
j=n*2-1;
y(:,n)=x(:,j);
y(:,lc+n)=x(:,j+1);
end
subplot(2,3,2);imshow(uint8(y)); title('ÆæżÁÐÖØÅÅ');
%-----------------ÐÐÖØÅÅ
for n=1:lc
y(:,lc+n)=y(:,lc+n)-y(:,n);
end
for n=1:lc
k=y(:,lc+n)/2-mod(y(:,lc+n)/2,1);
y(:,n)=y(:,n)+k;
end
subplot(2,3,3);imshow(uint8(y));title('ÐÐÖØÅÅ');
%---------------ÆæżÐÐÖØÅÅ
x=y;
for n=1:hc
j=n*2-1;
y(n,:)=x(j,:);
y(hc+n,:)=x(j+1,:);
end
subplot(2,3,4);imshow(uint8(y));title('ÆæżÐÐÖØÅÅ');
%---------------ÁÐÖØÅÅ
for n=1:hc
y(hc+n,:)=y(hc+n,:)-y(n,:);
end
for n=1:hc
k=y(hc+n,:)/2-mod(y(hc+n,:)/2,1);
y(n,:)=y(n,:)+k;
end
subplot(2,3,5);imshow(uint8(y));title('ÁÐÖØÅÅ');
%----------------Öظ´ÉÏÊö²½Ö裬½øÐжþ´ÎѹËõ
for n=1:lc
j=n*2-1;
y(:,n)=x(:,j);
y(:,lc+n)=x(:,j+1);
end
for n=1:lc
y(:,lc+n)=y(:,lc+n)-y(:,n);
end
for n=1:lc
k=y(:,lc+n)/2-mod(y(:,lc+n)/2,1);
y(:,n)=y(:,n)+k;
end
x=y;
for n=1:hc
j=n*2-1;
y(n,:)=x(j,:);
y(hc+n,:)=x(j+1,:);
end
for n=1:hc
y(hc+n,:)=y(hc+n,:)-y(n,:);
end
for n=1:hc
k=y(hc+n,:)/2-mod(y(hc+n,:)/2,1);
y(n,:)=y(n,:)+k;
end
%----------------ÔÙ½øÐÐÆæżÐÐÖØÅźÍÁÐÖØÅÅ
x=y;
for n=1:hc
j=n*2-1;
y(n,:)=x(j,:);
y(hc+n,:)=x(j+1,:);
end
for n=1:hc
y(hc+n,:)=y(hc+n,:)-y(n,:);
end
for n=1:hc
k=y(hc+n,:)/2-mod(y(hc+n,:)/2,1);
y(n,:)=y(n,:)+k;
end
subplot(2,3,6);imshow(uint8(y));title('¶þ´ÎѹËõ ');