per.m

用于多通道滤波器组实现完全重构,主要应用于图象处理

s4=t4down1down_up1+t4down2down_up2+t4down3down_up3+t4down4down_up4+t4down5down_up5+t4down6down_up6+t4down7down_up7+t4down8down_up8;
s5=t5down1down_up1+t5down2down_up2+t5down3down_up3+t5down4down_up4+t5down5down_up5+t5down6down_up6+t5down7down_up7+t5down8down_up8;
s6=t6down1down_up1+t6down2down_up2+t6down3down_up3+t6down4down_up4+t6down5down_up5+t6down6down_up6+t6down7down_up7+t6down8down_up8;
s7=t7down1down_up1+t7down2down_up2+t7down3down_up3+t7down4down_up4+t7down5down_up5+t7down6down_up6+t7down7down_up7+t7down8down_up8;
s8=t8down1down_up1+t8down2down_up2+t8down3down_up3+t8down4down_up4+t8down5down_up5+t8down6down_up6+t8down7down_up7+t8down8down_up8;




for j=1:539
    for k=0:87
        d1up(j,6*k+1)=s1(j,k+1);
        d1up(j,6*k+2)=0;
        d1up(j,6*k+3)=0;
        d1up(j,6*k+4)=0;
        d1up(j,6*k+5)=0;
        d1up(j,6*k+6)=0;
        
        d2up(j,6*k+1)=s2(j,k+1);
        d2up(j,6*k+2)=0;
        d2up(j,6*k+3)=0;
        d2up(j,6*k+4)=0;
        d2up(j,6*k+5)=0;
        d2up(j,6*k+6)=0;
        
        d3up(j,6*k+1)=s3(j,k+1);
        d3up(j,6*k+2)=0;
        d3up(j,6*k+3)=0;
        d3up(j,6*k+4)=0;
        d3up(j,6*k+5)=0;
        d3up(j,6*k+6)=0;
        
        d4up(j,6*k+1)=s4(j,k+1);
        d4up(j,6*k+2)=0;
        d4up(j,6*k+3)=0;
        d4up(j,6*k+4)=0;
        d4up(j,6*k+5)=0;
        d4up(j,6*k+6)=0;
        
        d5up(j,6*k+1)=s5(j,k+1);
        d5up(j,6*k+2)=0;
        d5up(j,6*k+3)=0;
        d5up(j,6*k+4)=0;
        d5up(j,6*k+5)=0;
        d5up(j,6*k+6)=0;
        
        d6up(j,6*k+1)=s6(j,k+1);
        d6up(j,6*k+2)=0;
        d6up(j,6*k+3)=0;
        d6up(j,6*k+4)=0;
        d6up(j,6*k+5)=0;
        d6up(j,6*k+6)=0;
        
        d7up(j,6*k+1)=s7(j,k+1);
        d7up(j,6*k+2)=0;
        d7up(j,6*k+3)=0;
        d7up(j,6*k+4)=0;
        d7up(j,6*k+5)=0;
        d7up(j,6*k+6)=0;
        
        d8up(j,6*k+1)=s8(j,k+1);
        d8up(j,6*k+2)=0;
        d8up(j,6*k+3)=0;
        d8up(j,6*k+4)=0;
        d8up(j,6*k+5)=0;
        d8up(j,6*k+6)=0;
    end
end                                             %539*528

for v=1:539
    d1up1(v,:)=conv(d1up(v,:),l1);
    d2up2(v,:)=conv(d2up(v,:),l2);
    d3up3(v,:)=conv(d3up(v,:),l3);
    d4up4(v,:)=conv(d4up(v,:),l4);
    d5up5(v,:)=conv(d5up(v,:),l5);
    d6up6(v,:)=conv(d6up(v,:),l6);
    d7up7(v,:)=conv(d7up(v,:),l7);
    d8up8(v,:)=conv(d8up(v,:),l8);
end                                             %539*539

d=d1up1+d2up2+d3up3+d4up4+d5up5+d6up6+d7up7+d8up8;
g=uint8(d);
for v=1:539
    m(1:512,v)=g(12:523,v);
end
for v=1:512
    u(v,1:512)=m(v,12:523);
end                                             %539*512
imshow(u)


%小图像

%l1=min(min(m1));
%e1=m1-l1;
%v1=uint8(e1);

l11=min(min(t1down1down));
e11=t1down1down-l11;
v11=uint8(e11);

l12=min(min(t1down2down));
e12=t1down2down-l12;
v12=uint8(e12);

l13=min(min(t1down3down));
e13=t1down3down-l13;
v13=uint8(e13);

l14=min(min(t1down4down));
e14=t1down4down-l14;
v14=uint8(e14);

l15=min(min(t1down5down));
e15=t1down5down-l15;
v15=uint8(e15);

l16=min(min(t1down6down));
e16=t1down6down-l16;
v16=uint8(e16);

l17=min(min(t1down7down));
e17=t1down7down-l17;
v17=uint8(e17);

l18=min(min(t1down8down));
e18=t1down8down-l18;
v18=uint8(e18);



l21=min(min(t2down1down));
e21=t2down1down-l21;
v21=uint8(e21);

l22=min(min(t2down2down));
e22=t1down2down-l22;
v22=uint8(e22);

l23=min(min(t2down3down));
e23=t2down3down-l23;
v23=uint8(e23);

l24=min(min(t2down4down));
e24=t2down4down-l24;
v24=uint8(e24);

l25=min(min(t2down5down));
e25=t2down5down-l25;
v25=uint8(e25);

l26=min(min(t2down6down));
e26=t2down6down-l26;
v26=uint8(e26);

l27=min(min(t2down7down));
e27=t2down7down-l27;
v27=uint8(e27);

l28=min(min(t2down8down));
e28=t2down8down-l28;
v28=uint8(e28);


l31=min(min(t3down1down));
e31=t3down1down-l31;
v31=uint8(e31);

l32=min(min(t3down2down));
e32=t3down2down-l32;
v32=uint8(e32);

l33=min(min(t3down3down));
e33=t3down3down-l33;
v33=uint8(e33);

l34=min(min(t3down4down));
e34=t3down4down-l34;
v34=uint8(e34);

l35=min(min(t3down5down));
e35=t3down5down-l35;
v35=uint8(e35);

l36=min(min(t3down6down));
e36=t3down6down-l36;
v36=uint8(e36);

l37=min(min(t3down7down));
e37=t3down7down-l37;
v37=uint8(e37);

l38=min(min(t3down8down));
e38=t3down8down-l38;
v38=uint8(e38);




l41=min(min(t4down1down));
e41=t4down1down-l41;
v41=uint8(e41);

l42=min(min(t4down2down));
e42=t4down2down-l42;
v42=uint8(e42);

l43=min(min(t4down3down));
e43=t4down3down-l43;
v43=uint8(e43);

l44=min(min(t4down4down));
e44=t4down4down-l44;
v44=uint8(e44);

l45=min(min(t4down5down));
e45=t4down5down-l45;
v45=uint8(e45);

l46=min(min(t4down6down));
e46=t4down6down-l46;
v46=uint8(e46);

l47=min(min(t4down7down));
e47=t4down7down-l47;
v47=uint8(e47);

l48=min(min(t4down8down));
e48=t4down8down-l48;
v48=uint8(e48);

l51=min(min(t5down1down));
e51=t5down1down-l51;
v51=uint8(e51);

l52=min(min(t5down2down));
e52=t5down2down-l52;
v52=uint8(e52);

l53=min(min(t5down3down));
e53=t5down3down-l53;
v53=uint8(e53);

l54=min(min(t5down4down));
e54=t5down4down-l54;
v54=uint8(e54);

l55=min(min(t5down5down));
e55=t5down5down-l55;
v55=uint8(e55);

l56=min(min(t5down6down));
e56=t5down6down-l56;
v56=uint8(e56);

l57=min(min(t5down7down));
e57=t5down7down-l57;
v57=uint8(e57);

l58=min(min(t5down8down));
e58=t5down8down-l58;
v58=uint8(e58);

l61=min(min(t6down1down));
e61=t6down1down-l61;
v61=uint8(e61);

l62=min(min(t6down2down));
e62=t6down2down-l62;
v62=uint8(e62);

l63=min(min(t6down3down));
e63=t6down3down-l63;
v63=uint8(e63);

l64=min(min(t6down4down));
e64=t6down4down-l64;
v64=uint8(e64);

l65=min(min(t6down5down));
e65=t6down5down-l65;
v65=uint8(e65);

l66=min(min(t6down6down));
e66=t6down6down-l66;
v66=uint8(e66);

l67=min(min(t6down7down));
e67=t6down7down-l67;
v67=uint8(e67);

l68=min(min(t6down8down));
e68=t6down8down-l68;
v68=uint8(e68);

l71=min(min(t7down1down));
e71=t7down1down-l71;
v71=uint8(e71);

l72=min(min(t7down2down));
e72=t7down2down-l72;
v72=uint8(e72);

l73=min(min(t7down3down));
e73=t7down3down-l73;
v73=uint8(e73);

l74=min(min(t7down4down));
e74=t7down4down-l74;
v74=uint8(e74);

l75=min(min(t7down5down));
e75=t7down5down-l75;
v75=uint8(e75);

l76=min(min(t7down6down));
e76=t7down6down-l76;
v76=uint8(e76);

l77=min(min(t7down7down));
e77=t7down7down-l77;
v77=uint8(e77);

l78=min(min(t7down8down));
e78=t7down8down-l78;
v78=uint8(e78);


l81=min(min(t8down1down));
e81=t8down1down-l81;
v81=uint8(e81);

l82=min(min(t8down2down));
e82=t8down2down-l82;
v82=uint8(e82);

l83=min(min(t8down3down));
e83=t8down3down-l83;
v83=uint8(e83);

l84=min(min(t8down4down));
e84=t8down4down-l84;
v84=uint8(e84);

l85=min(min(t8down5down));
e85=t8down5down-l85;
v85=uint8(e85);

l86=min(min(t8down6down));
e86=t8down6down-l86;
v86=uint8(e86);

l87=min(min(t8down7down));
e87=t8down7down-l87;
v87=uint8(e87);

l88=min(min(t8down8down));
e88=t8down8down-l88;
v88=uint8(e88);