我給你一個程式：

load wbarb

X1=X;

map1=map;

subplot(2,2,1);

image(X1);colormap(map1);

title("original image wbarb");

load woman;

X2=X;map2=map;

subplot(2,2,2);

image(X2);colormap(map2);

title("original image woman");

[c1,l1]=wavedec2(X1,2,"db3");

[c2,l2]=wavedec2(X2,2,"db3");

c=c1+c2; %fusing the two decomposition coefficients

XX=waverec2(c,l1,"db3");

subplot(2,2,3);

image(XX);

title("fusing image I");

Csize1=size(c1);

for i=1:Csize1(2)

c1(i)=1.2*c1(i);

end;

Csize2=size(c2);

for j=1:Csize1(2)

c2(i)=0.8*c2(j);

end;

c=0.5*(c1+c2);

XXX=waverec2(c,l2,"db3");

subplot(2,2,4);

image(XXX);

title("fusing image II");

這個可以用，直接執行就好