function?[LLHLLHHH]=improve97（x）
m=length（size（x））;
x=double（x）;
LL=[];HL=[];LH=[];HH=[];
m1=-1.5861343420693648;
m2=-0.0529801185718856;
m3=0.8829110755411875;
m4=0.4435068520511142;
m5=1.1496043988602418;
if?m>2
????error（‘this?is?not?a?2d?image!‘）;
else?if?m==2
????%先處理列，使列縮小為原來大小的一半，得到LH
???s1=x（:1:2:end-1）;
???s0=x（:2:2:end）;
???s01=conv2（[s1s1（:end）][m1m1]）;
???s01=s01（:2:1:end-1）;
???s0=s0+s01;
???s11=conv2（[s0（:1）s0][m2m2]）;
???s11=s11（:2:1:end-1）;
???s1=s1+s11;
???s02=conv2（[s1s1（:end）][m3m3]）;
???s02=s02（:2:1:end-1）;
???s0=s0+s02;
???s12=conv2（[s0（:1）s0][m4m4]）;
???s12=s12（:2:1:end-1）;
???s1=s1+s12;
???s1=s1.*m5;
???s0=s0./m5;
???%在處理行（L部分），使行縮小為原來大小的一半，得到LLLH
???s1=s1‘;
???s21=s1（:1:2:end-1）;
???s20=s1（:2:2:end）;
???s01=conv2（[s21s21（:end）][m1m1]）;
???s01=s01（:2:1:end-1）;
???s20=s20+s01;
???s11=conv2（[s20（:1）s20][m2m2]）;
???s11=s11（:2:1:end-1）;
???s21=s21+s11;
???s02=conv2（[s21s21（:end）][m3m3]）;
???s02=s02（:2:1:end-1）;
???s20=s20+s02;
???s12=conv2（[s20（:1）s20][m4m4]）;
???s12=s12（:2:1:end-1）;
???s21=s21+s12;
???s21=s21.*m5;
???s20=s20./m5;
???LL=s21‘;LH=s20‘;
???%處理行（H部分），得到HHHL
???s0=s0‘;
???s31=s0（:1:2:end-1）;
???s30=s0（:2:2:end）;
???s01=conv2（[s31s31（:end）][m1m1]）;
???s01=s01（:2:1:end-1）;
???s30=s30+s01;
???s11=conv2（[s30（:1）s30][m2m2]）;
???s11=s11（:2:1:end-1）;
???s31=s31+s11;
???s02=conv2（[s31s31（:end）][m3m3]）;
???s02=s02（:2:1:end-1）;
???s30=s30+s02;
???s12=conv2（[s30（:1）s30][m4m4]）;
???s12=s12（:2:1:end-1）;
???s31=s31+s12;
???s31=s31.*m5;
???s30=s30./m5;
???HH=s30‘;HL=s31‘;
else?
???s1=x（:1:2:end-1）;
???s0=x（:2:2:end）;
???s01=conv（[s1s1（:end）][m1m1]）;
???s01=s01（:2:1:end-1）;
???s0=s0+s01;
???s11=conv（[s0（:1）s0][m2m2]）;
???s11=s11（:2:1:end-1）;
???s1=s1+s11;
???s02=conv（[s1s1（:end）][m3m3]）;
???s02=s02（:2:1:end-1）;
???s0=s0+s02;
???s12=conv（[s0（:1）s0][m4m4]）;
???s12=s12（:2:1:end-1）;
???s1=s1+s12;
???s1=s1.*m5;
???s0=s0./m5;
???LL=s1;HH=s0;
????end
end
save?X?LL;
figure（1）;
subplot（231）;imshow（x[]）;
title?（‘原始圖像‘）;
subplot（232）;imshow（LL[]）;
title?（‘分解后的細節部分‘）;
subplot（233）;imshow（HL[]）;
title?（‘分解后的垂直部分‘）;
subplot（234）;imshow（LH[]）;
title?（‘分解后的水平部分‘）;
subplot（235）;imshow（HH[]）;
title?（‘分解后的對角線部分‘）;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2452??2015-05-29?10:35??improve97.m