function?Y?=?fuse_lap（M1M2ztflag）
%UNtitleD?Summary?of?this?function?goes?here
%???Detailed?explanation?goes?here

[z1s1]=size（M1）;
TM1?=?double（M1）;

[z2s2]=size（M2）;
TM2=double（M2）;

if（z1~=z2）||（s1~=s2）
error（‘Input?image?are?not?of?the?same?size!‘）;
return;
end

w=1/256*[1?4?6?4?1;
????????4?16?24?16?4;
????????6?24?36?24?6;
????????4?16?24?16?4;
????????1?4?6?4?1];
?w=double（w）;

E1=cell（1zt）;
E2=cell（1zt）;

for?i1=1:zt-1
????[zs]=size（TM1）;
????zl（i1）=z;
????sl（i1）=s;

????G1=conv2（TM1w‘same‘）;
????G2=conv2（TM2w‘same‘）;
????
????M1T?=?conv2（（undec2（dec2（G1）））（4*w）‘same‘）;
	M2T?=?conv2（（undec2（dec2（G2）））（4*w）‘same‘）;

????[m_t1n_t1]=size（TM1）;
????[m_t2n_t2]=size（M1T）;
????m_t1=?max（m_t1m_t2）;
????n_t1?=?max（n_t1n_t2）;???
????TM1（m_t1n_t1）=0;
????M1T（m_t1n_t1）=0;
????
????[m_t3n_t3]=size（M2T）;
????[m_t4n_t4]=size（TM2）;
????m_t3=?max（m_t3m_t4）;
????n_t3?=?max（n_t3n_t4）;???
????M2T（m_t3n_t3）=0;
????TM2（m_t3n_t3）=0;

???
????E1（i1）?=?{TM1?-?M1T};
????E2（i1）?=?{TM2?-?M2T};
????
????TM1?=?dec2（G1）;
????TM2?=?dec2（G2）;???

end
%循環重構
%????[pq]=s