x=imread（‘lena.gif‘）;???%?載入圖像lena.gif，需要更改路徑到指定的文件
figure（1）;
imshow（x）;
title（‘原始圖像‘）;
x=double（x）;

%圖像分解
lpd=[1/2?1/2];hpd=[-1/4?2/3];???????????%?默認的低通、高通濾波器
[rowcol]=size（x）;??????????????????????%?讀取輸入矩陣的大小

%一級分解
for?j=1:row?????????????????????????????%?首先對輸入矩陣的每一行序列進行一維離散小波分解
????tmp1=x（j:）;
????[ca1cd1]=mydwt（tmp1lpdhpd1）;
????x（j:）=ca1;???????????????????%?將分解系數序列再存入矩陣x中，得到[L|H]
end
for?k=1:col?????????????????????????????%?再對輸入矩陣的每一列序列進行一維離散小波分解
????tmp2=x（:k）;
????[ca2cd2]=mydwt（tmp2lpdhpd1）;
????x（:k）=ca2;???????????????????%?將分解所得系數存入矩陣x中，得到[LLHl;LHHH]
end

%二級分解
for?j=1:row/2?????????????????????????????%?首先對輸入矩陣左上角四分之一矩陣的每一行序列進行一維離散小波分解
????tmp1=x（j1:col/2）;
????[ca1cd1]=mydwt（tmp1lpdhpd1）;
????x（j1:col/2）=ca1;???????????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????359424??2010-05-09?21:54??實驗二?圖像變換.doc

?????文件????????308??2010-05-08?14:15??downspl.m

?????文件???????1204??2010-05-08?14:17??myidwt.m

?????文件????????431??2010-05-08?14:18??upspl.m

?????文件????????846??2010-05-09?14:54??mydwt.m

?????文件???????4033??2010-05-09?21:54??2.m

?????文件???????4804??2010-05-09?21:19??wave2gray.m

?????文件???????1089??2004-12-16?13:36??wavecut.m

?????文件???????3528??2004-12-16?13:36??wavework.m

?????文件????????940??2004-12-16?13:36??wavecopy.m

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

???????????????376607????????????????????10