per.m

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

        t4down1down_up(6*k+1,j)=t4down1down(k+1,j);
        t4down1down_up(6*k+2,j)=0;
        t4down1down_up(6*k+3,j)=0;
        t4down1down_up(6*k+4,j)=0;
        t4down1down_up(6*k+5,j)=0;
        t4down1down_up(6*k+6,j)=0;
        
        t4down2down_up(6*k+1,j)=t4down2down(k+1,j);
        t4down2down_up(6*k+2,j)=0;
        t4down2down_up(6*k+3,j)=0;
        t4down2down_up(6*k+4,j)=0;
        t4down2down_up(6*k+5,j)=0;
        t4down2down_up(6*k+6,j)=0;
        
        t4down3down_up(6*k+1,j)=t4down3down(k+1,j);
        t4down3down_up(6*k+2,j)=0;
        t4down3down_up(6*k+3,j)=0;
        t4down3down_up(6*k+4,j)=0;
        t4down3down_up(6*k+5,j)=0;
        t4down3down_up(6*k+6,j)=0;
        
        t4down4down_up(6*k+1,j)=t4down4down(k+1,j);
        t4down4down_up(6*k+2,j)=0;
        t4down4down_up(6*k+3,j)=0;
        t4down4down_up(6*k+4,j)=0;
        t4down4down_up(6*k+5,j)=0;
        t4down4down_up(6*k+6,j)=0;
        
        t4down5down_up(6*k+1,j)=t4down5down(k+1,j);
        t4down5down_up(6*k+2,j)=0;
        t4down5down_up(6*k+3,j)=0;
        t4down5down_up(6*k+4,j)=0;
        t4down5down_up(6*k+5,j)=0;
        t4down5down_up(6*k+6,j)=0;
        
        t4down6down_up(6*k+1,j)=t4down6down(k+1,j);
        t4down6down_up(6*k+2,j)=0;
        t4down6down_up(6*k+3,j)=0;
        t4down6down_up(6*k+4,j)=0;
        t4down6down_up(6*k+5,j)=0;
        t4down6down_up(6*k+6,j)=0;
        
        t4down7down_up(6*k+1,j)=t4down7down(k+1,j);
        t4down7down_up(6*k+2,j)=0;
        t4down7down_up(6*k+3,j)=0;
        t4down7down_up(6*k+4,j)=0;
        t4down7down_up(6*k+5,j)=0;
        t4down7down_up(6*k+6,j)=0;
        
        t4down8down_up(6*k+1,j)=t4down8down(k+1,j);
        t4down8down_up(6*k+2,j)=0;
        t4down8down_up(6*k+3,j)=0;
        t4down8down_up(6*k+4,j)=0;
        t4down8down_up(6*k+5,j)=0;
        t4down8down_up(6*k+6,j)=0;
        
        t5down1down_up(6*k+1,j)=t5down1down(k+1,j);
        t5down1down_up(6*k+2,j)=0;
        t5down1down_up(6*k+3,j)=0;
        t5down1down_up(6*k+4,j)=0;
        t5down1down_up(6*k+5,j)=0;
        t5down1down_up(6*k+6,j)=0;
        
        t5down2down_up(6*k+1,j)=t5down2down(k+1,j);
        t5down2down_up(6*k+2,j)=0;
        t5down2down_up(6*k+3,j)=0;
        t5down2down_up(6*k+4,j)=0;
        t5down2down_up(6*k+5,j)=0;
        t5down2down_up(6*k+6,j)=0;
        
        t5down3down_up(6*k+1,j)=t5down3down(k+1,j);
        t5down3down_up(6*k+2,j)=0;
        t5down3down_up(6*k+3,j)=0;
        t5down3down_up(6*k+4,j)=0;
        t5down3down_up(6*k+5,j)=0;
        t5down3down_up(6*k+6,j)=0;
        
        t5down4down_up(6*k+1,j)=t5down4down(k+1,j);
        t5down4down_up(6*k+2,j)=0;
        t5down4down_up(6*k+3,j)=0;
        t5down4down_up(6*k+4,j)=0;
        t5down4down_up(6*k+5,j)=0;
        t5down4down_up(6*k+6,j)=0;
        
        t5down5down_up(6*k+1,j)=t5down5down(k+1,j);
        t5down5down_up(6*k+2,j)=0;
        t5down5down_up(6*k+3,j)=0;
        t5down5down_up(6*k+4,j)=0;
        t5down5down_up(6*k+5,j)=0;
        t5down5down_up(6*k+6,j)=0;
        
        t5down6down_up(6*k+1,j)=t5down6down(k+1,j);
        t5down6down_up(6*k+2,j)=0;
        t5down6down_up(6*k+3,j)=0;
        t5down6down_up(6*k+4,j)=0;
        t5down6down_up(6*k+5,j)=0;
        t5down6down_up(6*k+6,j)=0;
        
        t5down7down_up(6*k+1,j)=t5down7down(k+1,j);
        t5down7down_up(6*k+2,j)=0;
        t5down7down_up(6*k+3,j)=0;
        t5down7down_up(6*k+4,j)=0;
        t5down7down_up(6*k+5,j)=0;
        t5down7down_up(6*k+6,j)=0;
        
        t5down8down_up(6*k+1,j)=t5down8down(k+1,j);
        t5down8down_up(6*k+2,j)=0;
        t5down8down_up(6*k+3,j)=0;
        t5down8down_up(6*k+4,j)=0;
        t5down8down_up(6*k+5,j)=0;
        t5down8down_up(6*k+6,j)=0;
        
        t6down1down_up(6*k+1,j)=t6down1down(k+1,j);
        t6down1down_up(6*k+2,j)=0;
        t6down1down_up(6*k+3,j)=0;
        t6down1down_up(6*k+4,j)=0;
        t6down1down_up(6*k+5,j)=0;
        t6down1down_up(6*k+6,j)=0;
        
        t6down2down_up(6*k+1,j)=t6down2down(k+1,j);
        t6down2down_up(6*k+2,j)=0;
        t6down2down_up(6*k+3,j)=0;
        t6down2down_up(6*k+4,j)=0;
        t6down2down_up(6*k+5,j)=0;
        t6down2down_up(6*k+6,j)=0;
        
        t6down3down_up(6*k+1,j)=t6down3down(k+1,j);
        t6down3down_up(6*k+2,j)=0;
        t6down3down_up(6*k+3,j)=0;
        t6down3down_up(6*k+4,j)=0;
        t6down3down_up(6*k+5,j)=0;
        t6down3down_up(6*k+6,j)=0;
        
        t6down4down_up(6*k+1,j)=t6down4down(k+1,j);
        t6down4down_up(6*k+2,j)=0;
        t6down4down_up(6*k+3,j)=0;
        t6down4down_up(6*k+4,j)=0;
        t6down4down_up(6*k+5,j)=0;
        t6down4down_up(6*k+6,j)=0;
        
        t6down5down_up(6*k+1,j)=t6down5down(k+1,j);
        t6down5down_up(6*k+2,j)=0;
        t6down5down_up(6*k+3,j)=0;
        t6down5down_up(6*k+4,j)=0;
        t6down5down_up(6*k+5,j)=0;
        t6down5down_up(6*k+6,j)=0;
        
        t6down6down_up(6*k+1,j)=t6down6down(k+1,j);
        t6down6down_up(6*k+2,j)=0;
        t6down6down_up(6*k+3,j)=0;
        t6down6down_up(6*k+4,j)=0;
        t6down6down_up(6*k+5,j)=0;
        t6down6down_up(6*k+6,j)=0;
        
        t6down7down_up(6*k+1,j)=t6down7down(k+1,j);
        t6down7down_up(6*k+2,j)=0;
        t6down7down_up(6*k+3,j)=0;
        t6down7down_up(6*k+4,j)=0;
        t6down7down_up(6*k+5,j)=0;
        t6down7down_up(6*k+6,j)=0;
        
        t6down8down_up(6*k+1,j)=t6down8down(k+1,j);
        t6down8down_up(6*k+2,j)=0;
        t6down8down_up(6*k+3,j)=0;
        t6down8down_up(6*k+4,j)=0;
        t6down8down_up(6*k+5,j)=0;
        t6down8down_up(6*k+6,j)=0;
        
        t7down1down_up(6*k+1,j)=t7down1down(k+1,j);
        t7down1down_up(6*k+2,j)=0;
        t7down1down_up(6*k+3,j)=0;
        t7down1down_up(6*k+4,j)=0;
        t7down1down_up(6*k+5,j)=0;
        t7down1down_up(6*k+6,j)=0;
        
        t7down2down_up(6*k+1,j)=t7down2down(k+1,j);
        t7down2down_up(6*k+2,j)=0;
        t7down2down_up(6*k+3,j)=0;
        t7down2down_up(6*k+4,j)=0;
        t7down2down_up(6*k+5,j)=0;
        t7down2down_up(6*k+6,j)=0;
        
        t7down3down_up(6*k+1,j)=t7down3down(k+1,j);
        t7down3down_up(6*k+2,j)=0;
        t7down3down_up(6*k+3,j)=0;
        t7down3down_up(6*k+4,j)=0;
        t7down3down_up(6*k+5,j)=0;
        t7down3down_up(6*k+6,j)=0;
        
        t7down4down_up(6*k+1,j)=t7down4down(k+1,j);
        t7down4down_up(6*k+2,j)=0;
        t7down4down_up(6*k+3,j)=0;
        t7down4down_up(6*k+4,j)=0;
        t7down4down_up(6*k+5,j)=0;
        t7down4down_up(6*k+6,j)=0;
        
        t7down5down_up(6*k+1,j)=t7down5down(k+1,j);
        t7down5down_up(6*k+2,j)=0;
        t7down5down_up(6*k+3,j)=0;
        t7down5down_up(6*k+4,j)=0;
        t7down5down_up(6*k+5,j)=0;
        t7down5down_up(6*k+6,j)=0;
        
        t7down6down_up(6*k+1,j)=t7down6down(k+1,j);
        t7down6down_up(6*k+2,j)=0;
        t7down6down_up(6*k+3,j)=0;
        t7down6down_up(6*k+4,j)=0;
        t7down6down_up(6*k+5,j)=0;
        t7down6down_up(6*k+6,j)=0;
        
        t7down7down_up(6*k+1,j)=t7down7down(k+1,j);
        t7down7down_up(6*k+2,j)=0;
        t7down7down_up(6*k+3,j)=0;
        t7down7down_up(6*k+4,j)=0;
        t7down7down_up(6*k+5,j)=0;
        t7down7down_up(6*k+6,j)=0;
        
        t7down8down_up(6*k+1,j)=t7down8down(k+1,j);
        t7down8down_up(6*k+2,j)=0;
        t7down8down_up(6*k+3,j)=0;
        t7down8down_up(6*k+4,j)=0;
        t7down8down_up(6*k+5,j)=0;
        t7down8down_up(6*k+6,j)=0;
        
        t8down1down_up(6*k+1,j)=t8down1down(k+1,j);
        t8down1down_up(6*k+2,j)=0;
        t8down1down_up(6*k+3,j)=0;
        t8down1down_up(6*k+4,j)=0;
        t8down1down_up(6*k+5,j)=0;
        t8down1down_up(6*k+6,j)=0;
        
        t8down2down_up(6*k+1,j)=t8down2down(k+1,j);
        t8down2down_up(6*k+2,j)=0;
        t8down2down_up(6*k+3,j)=0;
        t8down2down_up(6*k+4,j)=0;
        t8down2down_up(6*k+5,j)=0;
        t8down2down_up(6*k+6,j)=0;
        
        t8down3down_up(6*k+1,j)=t8down3down(k+1,j);
        t8down3down_up(6*k+2,j)=0;
        t8down3down_up(6*k+3,j)=0;
        t8down3down_up(6*k+4,j)=0;
        t8down3down_up(6*k+5,j)=0;
        t8down3down_up(6*k+6,j)=0;
        
        t8down4down_up(6*k+1,j)=t8down4down(k+1,j);
        t8down4down_up(6*k+2,j)=0;
        t8down4down_up(6*k+3,j)=0;
        t8down4down_up(6*k+4,j)=0;
        t8down4down_up(6*k+5,j)=0;
        t8down4down_up(6*k+6,j)=0;
        
        t8down5down_up(6*k+1,j)=t8down5down(k+1,j);
        t8down5down_up(6*k+2,j)=0;
        t8down5down_up(6*k+3,j)=0;
        t8down5down_up(6*k+4,j)=0;
        t8down5down_up(6*k+5,j)=0;
        t8down5down_up(6*k+6,j)=0;
        
        t8down6down_up(6*k+1,j)=t8down6down(k+1,j);
        t8down6down_up(6*k+2,j)=0;
        t8down6down_up(6*k+3,j)=0;
        t8down6down_up(6*k+4,j)=0;
        t8down6down_up(6*k+5,j)=0;
        t8down6down_up(6*k+6,j)=0;
        
        t8down7down_up(6*k+1,j)=t8down7down(k+1,j);
        t8down7down_up(6*k+2,j)=0;
        t8down7down_up(6*k+3,j)=0;
        t8down7down_up(6*k+4,j)=0;
        t8down7down_up(6*k+5,j)=0;
        t8down7down_up(6*k+6,j)=0;
        
        t8down8down_up(6*k+1,j)=t8down8down(k+1,j);
        t8down8down_up(6*k+2,j)=0;
        t8down8down_up(6*k+3,j)=0;
        t8down8down_up(6*k+4,j)=0;
        t8down8down_up(6*k+5,j)=0;
        t8down8down_up(6*k+6,j)=0;
    end
end                                                                     %528*88


for v=1:88
    t1down1down_up1(:,v)=conv(t1down1down_up(:,v),l1);
    t1down2down_up2(:,v)=conv(t1down2down_up(:,v),l2);
    t1down3down_up3(:,v)=conv(t1down3down_up(:,v),l3);
    t1down4down_up4(:,v)=conv(t1down4down_up(:,v),l4);
    t1down5down_up5(:,v)=conv(t1down5down_up(:,v),l5);
    t1down6down_up6(:,v)=conv(t1down6down_up(:,v),l6);
    t1down7down_up7(:,v)=conv(t1down7down_up(:,v),l7);
    t1down8down_up8(:,v)=conv(t1down8down_up(:,v),l8);
    
    t2down1down_up1(:,v)=conv(t2down1down_up(:,v),l1);
    t2down2down_up2(:,v)=conv(t2down2down_up(:,v),l2);
    t2down3down_up3(:,v)=conv(t2down3down_up(:,v),l3);
    t2down4down_up4(:,v)=conv(t2down4down_up(:,v),l4);
    t2down5down_up5(:,v)=conv(t2down5down_up(:,v),l5);
    t2down6down_up6(:,v)=conv(t2down6down_up(:,v),l6);
    t2down7down_up7(:,v)=conv(t2down7down_up(:,v),l7);
    t2down8down_up8(:,v)=conv(t2down8down_up(:,v),l8);
    
    t3down1down_up1(:,v)=conv(t3down1down_up(:,v),l1);
    t3down2down_up2(:,v)=conv(t3down2down_up(:,v),l2);
    t3down3down_up3(:,v)=conv(t3down3down_up(:,v),l3);
    t3down4down_up4(:,v)=conv(t3down4down_up(:,v),l4);
    t3down5down_up5(:,v)=conv(t3down5down_up(:,v),l5);
    t3down6down_up6(:,v)=conv(t3down6down_up(:,v),l6);
    t3down7down_up7(:,v)=conv(t3down7down_up(:,v),l7);
    t3down8down_up8(:,v)=conv(t3down8down_up(:,v),l8);
    
    t4down1down_up1(:,v)=conv(t4down1down_up(:,v),l1);
    t4down2down_up2(:,v)=conv(t4down2down_up(:,v),l2);
    t4down3down_up3(:,v)=conv(t4down3down_up(:,v),l3);
    t4down4down_up4(:,v)=conv(t4down4down_up(:,v),l4);
    t4down5down_up5(:,v)=conv(t4down5down_up(:,v),l5);
    t4down6down_up6(:,v)=conv(t4down6down_up(:,v),l6);
    t4down7down_up7(:,v)=conv(t4down7down_up(:,v),l7);
    t4down8down_up8(:,v)=conv(t4down8down_up(:,v),l8);
    
    t5down1down_up1(:,v)=conv(t5down1down_up(:,v),l1);
    t5down2down_up2(:,v)=conv(t5down2down_up(:,v),l2);
    t5down3down_up3(:,v)=conv(t5down3down_up(:,v),l3);
    t5down4down_up4(:,v)=conv(t5down4down_up(:,v),l4);
    t5down5down_up5(:,v)=conv(t5down5down_up(:,v),l5);
    t5down6down_up6(:,v)=conv(t5down6down_up(:,v),l6);
    t5down7down_up7(:,v)=conv(t5down7down_up(:,v),l7);
    t5down8down_up8(:,v)=conv(t5down8down_up(:,v),l8);
    
    t6down1down_up1(:,v)=conv(t6down1down_up(:,v),l1);
    t6down2down_up2(:,v)=conv(t6down2down_up(:,v),l2);
    t6down3down_up3(:,v)=conv(t6down3down_up(:,v),l3);
    t6down4down_up4(:,v)=conv(t6down4down_up(:,v),l4);
    t6down5down_up5(:,v)=conv(t6down5down_up(:,v),l5);
    t6down6down_up6(:,v)=conv(t6down6down_up(:,v),l6);
    t6down7down_up7(:,v)=conv(t6down7down_up(:,v),l7);
    t6down8down_up8(:,v)=conv(t6down8down_up(:,v),l8);
    
    t7down1down_up1(:,v)=conv(t7down1down_up(:,v),l1);
    t7down2down_up2(:,v)=conv(t7down2down_up(:,v),l2);
    t7down3down_up3(:,v)=conv(t7down3down_up(:,v),l3);
    t7down4down_up4(:,v)=conv(t7down4down_up(:,v),l4);
    t7down5down_up5(:,v)=conv(t7down5down_up(:,v),l5);
    t7down6down_up6(:,v)=conv(t7down6down_up(:,v),l6);
    t7down7down_up7(:,v)=conv(t7down7down_up(:,v),l7);
    t7down8down_up8(:,v)=conv(t7down8down_up(:,v),l8);
    
    t8down1down_up1(:,v)=conv(t8down1down_up(:,v),l1);
    t8down2down_up2(:,v)=conv(t8down2down_up(:,v),l2);
    t8down3down_up3(:,v)=conv(t8down3down_up(:,v),l3);
    t8down4down_up4(:,v)=conv(t8down4down_up(:,v),l4);
    t8down5down_up5(:,v)=conv(t8down5down_up(:,v),l5);
    t8down6down_up6(:,v)=conv(t8down6down_up(:,v),l6);
    t8down7down_up7(:,v)=conv(t8down7down_up(:,v),l7);
    t8down8down_up8(:,v)=conv(t8down8down_up(:,v),l8);
end                                                      %539*88

s1=t1down1down_up1+t1down2down_up2+t1down3down_up3+t1down4down_up4+t1down5down_up5+t1down6down_up6+t1down7down_up7+t1down8down_up8;
s2=t2down1down_up1+t2down2down_up2+t2down3down_up3+t2down4down_up4+t2down5down_up5+t2down6down_up6+t2down7down_up7+t2down8down_up8;
s3=t3down1down_up1+t3down2down_up2+t3down3down_up3+t3down4down_up4+t3down5down_up5+t3down6down_up6+t3down7down_up7+t3down8down_up8;