clc;
clear?all;
X1=imread（‘1.tif‘）;
X1=double（X1）/256;??%這里轉化成double類型，否則使用小波變換輸出的會有大量大于1的存在，會導致圖像顯示有問題
figure（1）;
imshow（X1）title（‘左焦距‘）;
axis?square;
X2=imread（‘2.tif‘）;
X2=double（X2）/256;
imshow（X2）title（‘右焦距‘）;
axis?square;
[c1s1]=wavedec2（X13‘db3‘）;?%將x1進行2維，使用‘sym4’進行變換
sizec1=size（c1）;
for?I=1:sizec1（2）
????c1（I）=1.2*c1（I）;?%將分解后的值都擴大1.2倍
end
[c2s2]=wavedec2（X23‘db3‘）;
c=c1+c2;?%計算平均值
c=0.5*c;
s=s1+s2;
s=0.5*s;
xx=waverec2（cs‘db3‘）;?%進行重構
figure;
imshow（xx）title（‘融合后的‘）;?
axis?square;
%?[cH1cV1cD1]=detcoef2（‘all‘cs1）;%尺度1的所有方向的高頻系數
%?[cH2cV2cD2]=detcoef2（‘