%function?[allTmpliulei]=test1（IM）
%IM是輸入的原圖像
tmp=imread（‘vvvvv.jpg‘）;
tmp=tmp（::1）;
IM=double（tmp）;
IM=IM（1:1181:174）;
figure（1）;
imshow（uint8（IM））;
[sizeXsizeY]=size（IM）;

%顯示方向可調(diào)小波
jiao=-3.1415926/3;
pi=3.1415926;
for?f=1:6
????[x1y1]=meshgrid（-5:0.25:5-5:0.25:5）;
????[mn]=size（x1）;
????for?i=1:m
????????for?j=1:n
????????????x=x1（ij）;
????????????y=y1（1j）;
????????????
????????????z（ij）=-（x/（2*pi）*exp（-（x*x+y*y）/2））*cos（pi/6*（f-1））-（y/（2*pi）*exp（-（x*x+y*y）/2））*sin（pi/6*（f-1））;
????????end
????end
????figure（2）;
????subplot（23f）;
????mesh（x1y1z）;
????colormap（[0?0?0]）;
end

l=0;
%x方向變換
for?i=1:sizeX
????for?j=1:sizeY
????????tmp1=0;
????????for?m=-2:2
????????????for?n=-2:2
????????????????if（i-m）>0&&（i-m）<=sizeX&&（j-n）>0&&（j-n）<=sizeY
????????????????????tmp1=tmp1+IM（i-mj-n）*m*（-1/（2*pi）*exp（-（m*m+n*n）/2））;
????????????????end
????????????end
????????end
????????IMtmp1（ij）=tmp1;
????end
end

%y方向變換
for?i=1:sizeX
????for?j=1:sizeY
????????tmp1=0;
????????for?m=-2:2
????????????for?n=-2:2
????????????????if（i-m）>0&&（i-m）<=sizeX&&（j-n）>0&&（j-n）<=sizeY
????????????????????tmp1=tmp1+IM（i-mj-n）*n*（-1/（2*pi）*exp（-（m*m+n*n）/2））;
????????????????end
????????????end
????????end
????????IMtmp2（ij）=tmp1;
????end
end


allTmp=0;
%求6個(gè)方向的平均
for?pp=2:7
????jiao=pi/6*（pp-2）;
????
????
????IMtmp11=-IMtmp1*cos（jiao）;
????IMtmp22=-IMtmp2*sin（jiao）;
????IMtmp=IMtmp11+IMtmp22;
????
????
????IMtmpAbs=abs（IMtmp）;
????maxv=max（max（IMtmpAbs））;
????minv=min（min（IMtmpAbs））;
????minv1=minv*ones（sizeXsizeY）;
????
????
????IMtmpNom=（IMtmpAbs-minv1）*255/（maxv-minv）;
????liulei（::pp-1）=IMtmpNom;
????allTmp=allTmp+IMtmpNom;
end
%結(jié)果輸出
allTmp=allTmp/6;
figure（3）;
imshow（uint8（allTmp））;







????????????

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件??????70438??2010-04-14?13:00??實(shí)驗(yàn)3\p6.bmp

?????文件??????70438??2010-04-14?13:00??實(shí)驗(yàn)3\p8.bmp

?????文件???????1841??2010-04-14?12:30??實(shí)驗(yàn)3\test1.asv

?????文件???????1852??2010-04-14?12:59??實(shí)驗(yàn)3\test1.m

?????文件???????1872??2010-04-14?18:30??實(shí)驗(yàn)3\test8.asv

?????文件???????1871??2010-04-14?19:00??實(shí)驗(yàn)3\test8.m

?????文件?????157025??2009-06-03?08:22??實(shí)驗(yàn)3\vvvvv.JPG

?????文件?????255210??2010-06-15?12:02??實(shí)驗(yàn)3\多媒體實(shí)驗(yàn)報(bào)告3.docx

?????目錄??????????0??2010-06-15?12:02??實(shí)驗(yàn)3

-----------?---------??----------?-----??----

???????????????560547????????????????????9