per.m

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

clear t1 t2 t3 t4 t5 t6 t7 t8;
clear m1 m2 m3 m4 m5 m6 m7 m8 m m11 m22 m33 m44 m55 m66 m77 m88 k;
clear new1 new2 new3 new4 new5 new6 new7 new8;
clear new11 new22 new33 new44 new55 new66 new77 new88;

x=imread('lena1.jpg');
I=imread('water_t.jpg');
L=double(I);
t=double(x);
L1=L(:);
for v=1:512
    t1(v,:)=conv(t(v,:),P1);
    t2(v,:)=conv(t(v,:),P2);
    t3(v,:)=conv(t(v,:),P3);
    t4(v,:)=conv(t(v,:),P4);
    t5(v,:)=conv(t(v,:),P5);
    t6(v,:)=conv(t(v,:),P6);
    t7(v,:)=conv(t(v,:),P7);
    t8(v,:)=conv(t(v,:),P8);
end                                               %512*523
 
for j=1:512
    for k=0:87
        t1down(j,k+1)=t1(j,6*k+1);
        t2down(j,k+1)=t2(j,6*k+1);
        t3down(j,k+1)=t3(j,6*k+1);
        t4down(j,k+1)=t4(j,6*k+1);
        t5down(j,k+1)=t5(j,6*k+1);
        t6down(j,k+1)=t6(j,6*k+1);
        t7down(j,k+1)=t7(j,6*k+1);
        t8down(j,k+1)=t8(j,6*k+1);
    end
end                                              %512*88


for v=1:88
    t1down1(:,v)=conv(t1down(:,v),P1);
    t1down2(:,v)=conv(t1down(:,v),P2);
    t1down3(:,v)=conv(t1down(:,v),P3);
    t1down4(:,v)=conv(t1down(:,v),P4);
    t1down5(:,v)=conv(t1down(:,v),P5);
    t1down6(:,v)=conv(t1down(:,v),P6);
    t1down7(:,v)=conv(t1down(:,v),P7);
    t1down8(:,v)=conv(t1down(:,v),P8);
    
    t2down1(:,v)=conv(t2down(:,v),P1);
    t2down2(:,v)=conv(t2down(:,v),P2);
    t2down3(:,v)=conv(t2down(:,v),P3);
    t2down4(:,v)=conv(t2down(:,v),P4);
    t2down5(:,v)=conv(t2down(:,v),P5);
    t2down6(:,v)=conv(t2down(:,v),P6);
    t2down7(:,v)=conv(t2down(:,v),P7);
    t2down8(:,v)=conv(t2down(:,v),P8);
    
    t3down1(:,v)=conv(t3down(:,v),P1);
    t3down2(:,v)=conv(t3down(:,v),P2);
    t3down3(:,v)=conv(t3down(:,v),P3);
    t3down4(:,v)=conv(t3down(:,v),P4);
    t3down5(:,v)=conv(t3down(:,v),P5);
    t3down6(:,v)=conv(t3down(:,v),P6);
    t3down7(:,v)=conv(t3down(:,v),P7);
    t3down8(:,v)=conv(t3down(:,v),P8);
    
    t4down1(:,v)=conv(t4down(:,v),P1);
    t4down2(:,v)=conv(t4down(:,v),P2);
    t4down3(:,v)=conv(t4down(:,v),P3);
    t4down4(:,v)=conv(t4down(:,v),P4);
    t4down5(:,v)=conv(t4down(:,v),P5);
    t4down6(:,v)=conv(t4down(:,v),P6);
    t4down7(:,v)=conv(t4down(:,v),P7);
    t4down8(:,v)=conv(t4down(:,v),P8);
    
    t5down1(:,v)=conv(t5down(:,v),P1);
    t5down2(:,v)=conv(t5down(:,v),P2);
    t5down3(:,v)=conv(t5down(:,v),P3);
    t5down4(:,v)=conv(t5down(:,v),P4);
    t5down5(:,v)=conv(t5down(:,v),P5);
    t5down6(:,v)=conv(t5down(:,v),P6);
    t5down7(:,v)=conv(t5down(:,v),P7);
    t5down8(:,v)=conv(t5down(:,v),P8);
    
    t6down1(:,v)=conv(t6down(:,v),P1);
    t6down2(:,v)=conv(t6down(:,v),P2);
    t6down3(:,v)=conv(t6down(:,v),P3);
    t6down4(:,v)=conv(t6down(:,v),P4);
    t6down5(:,v)=conv(t6down(:,v),P5);
    t6down6(:,v)=conv(t6down(:,v),P6);
    t6down7(:,v)=conv(t6down(:,v),P7);
    t6down8(:,v)=conv(t6down(:,v),P8);
    
    t7down1(:,v)=conv(t7down(:,v),P1);
    t7down2(:,v)=conv(t7down(:,v),P2);
    t7down3(:,v)=conv(t7down(:,v),P3);
    t7down4(:,v)=conv(t7down(:,v),P4);
    t7down5(:,v)=conv(t7down(:,v),P5);
    t7down6(:,v)=conv(t7down(:,v),P6);
    t7down7(:,v)=conv(t7down(:,v),P7);
    t7down8(:,v)=conv(t7down(:,v),P8);
    
    t8down1(:,v)=conv(t8down(:,v),P1);
    t8down2(:,v)=conv(t8down(:,v),P2);
    t8down3(:,v)=conv(t8down(:,v),P3);
    t8down4(:,v)=conv(t8down(:,v),P4);
    t8down5(:,v)=conv(t8down(:,v),P5);
    t8down6(:,v)=conv(t8down(:,v),P6);
    t8down7(:,v)=conv(t8down(:,v),P7);
    t8down8(:,v)=conv(t8down(:,v),P8);
end                                                       %523*88

for j=1:88
    for k=0:87
        t1down1down(k+1,j)=t1down1(6*k+1,j);
        t1down2down(k+1,j)=t1down2(6*k+1,j);
        t1down3down(k+1,j)=t1down3(6*k+1,j);
        t1down4down(k+1,j)=t1down4(6*k+1,j);
        t1down5down(k+1,j)=t1down5(6*k+1,j);
        t1down6down(k+1,j)=t1down6(6*k+1,j);
        t1down7down(k+1,j)=t1down7(6*k+1,j);
        t1down8down(k+1,j)=t1down8(6*k+1,j);
        
        t2down1down(k+1,j)=t2down1(6*k+1,j);
        t2down2down(k+1,j)=t2down2(6*k+1,j);
        t2down3down(k+1,j)=t2down3(6*k+1,j);
        t2down4down(k+1,j)=t2down4(6*k+1,j);
        t2down5down(k+1,j)=t2down5(6*k+1,j);
        t2down6down(k+1,j)=t2down6(6*k+1,j);
        t2down7down(k+1,j)=t2down7(6*k+1,j);
        t2down8down(k+1,j)=t2down8(6*k+1,j);
        
        t3down1down(k+1,j)=t3down1(6*k+1,j);
        t3down2down(k+1,j)=t3down2(6*k+1,j);
        t3down3down(k+1,j)=t3down3(6*k+1,j);
        t3down4down(k+1,j)=t3down4(6*k+1,j);
        t3down5down(k+1,j)=t3down5(6*k+1,j);
        t3down6down(k+1,j)=t3down6(6*k+1,j);
        t3down7down(k+1,j)=t3down7(6*k+1,j);
        t3down8down(k+1,j)=t3down8(6*k+1,j);
        
        t4down1down(k+1,j)=t4down1(6*k+1,j);
        t4down2down(k+1,j)=t4down2(6*k+1,j);
        t4down3down(k+1,j)=t4down3(6*k+1,j);
        t4down4down(k+1,j)=t4down4(6*k+1,j);
        t4down5down(k+1,j)=t4down5(6*k+1,j);
        t4down6down(k+1,j)=t4down6(6*k+1,j);
        t4down7down(k+1,j)=t4down7(6*k+1,j);
        t4down8down(k+1,j)=t4down8(6*k+1,j);
        
        t5down1down(k+1,j)=t5down1(6*k+1,j);
        t5down2down(k+1,j)=t5down2(6*k+1,j);
        t5down3down(k+1,j)=t5down3(6*k+1,j);
        t5down4down(k+1,j)=t5down4(6*k+1,j);
        t5down5down(k+1,j)=t5down5(6*k+1,j);
        t5down6down(k+1,j)=t5down6(6*k+1,j);
        t5down7down(k+1,j)=t5down7(6*k+1,j);
        t5down8down(k+1,j)=t5down8(6*k+1,j);
        
        t6down1down(k+1,j)=t6down1(6*k+1,j);
        t6down2down(k+1,j)=t6down2(6*k+1,j);
        t6down3down(k+1,j)=t6down3(6*k+1,j);
        t6down4down(k+1,j)=t6down4(6*k+1,j);
        t6down5down(k+1,j)=t6down5(6*k+1,j);
        t6down6down(k+1,j)=t6down6(6*k+1,j);
        t6down7down(k+1,j)=t6down7(6*k+1,j);
        t6down8down(k+1,j)=t6down8(6*k+1,j);
        
        t7down1down(k+1,j)=t7down1(6*k+1,j);
        t7down2down(k+1,j)=t7down2(6*k+1,j);
        t7down3down(k+1,j)=t7down3(6*k+1,j);
        t7down4down(k+1,j)=t7down4(6*k+1,j);
        t7down5down(k+1,j)=t7down5(6*k+1,j);
        t7down6down(k+1,j)=t7down6(6*k+1,j);
        t7down7down(k+1,j)=t7down7(6*k+1,j);
        t7down8down(k+1,j)=t7down8(6*k+1,j);
        
        t8down1down(k+1,j)=t8down1(6*k+1,j);
        t8down2down(k+1,j)=t8down2(6*k+1,j);
        t8down3down(k+1,j)=t8down3(6*k+1,j);
        t8down4down(k+1,j)=t8down4(6*k+1,j);
        t8down5down(k+1,j)=t8down5(6*k+1,j);
        t8down6down(k+1,j)=t8down6(6*k+1,j);
        t8down7down(k+1,j)=t8down7(6*k+1,j);
        t8down8down(k+1,j)=t8down8(6*k+1,j);
    end
end                                                            %88*88





for j=1:88
    for k=0:87
        t1down1down_up(6*k+1,j)=t1down1down(k+1,j);
        t1down1down_up(6*k+2,j)=0;
        t1down1down_up(6*k+3,j)=0;
        t1down1down_up(6*k+4,j)=0;
        t1down1down_up(6*k+5,j)=0;
        t1down1down_up(6*k+6,j)=0;
        
        t1down2down_up(6*k+1,j)=t1down2down(k+1,j);
        t1down2down_up(6*k+2,j)=0;
        t1down2down_up(6*k+3,j)=0;
        t1down2down_up(6*k+4,j)=0;
        t1down2down_up(6*k+5,j)=0;
        t1down2down_up(6*k+6,j)=0;
        
        t1down3down_up(6*k+1,j)=t1down3down(k+1,j);
        t1down3down_up(6*k+2,j)=0;
        t1down3down_up(6*k+3,j)=0;
        t1down3down_up(6*k+4,j)=0;
        t1down3down_up(6*k+5,j)=0;
        t1down3down_up(6*k+6,j)=0;
        
        t1down4down_up(6*k+1,j)=t1down4down(k+1,j);
        t1down4down_up(6*k+2,j)=0;
        t1down4down_up(6*k+3,j)=0;
        t1down4down_up(6*k+4,j)=0;
        t1down4down_up(6*k+5,j)=0;
        t1down4down_up(6*k+6,j)=0;
        
        t1down5down_up(6*k+1,j)=t1down5down(k+1,j);
        t1down5down_up(6*k+2,j)=0;
        t1down5down_up(6*k+3,j)=0;
        t1down5down_up(6*k+4,j)=0;
        t1down5down_up(6*k+5,j)=0;
        t1down5down_up(6*k+6,j)=0;
        
        t1down6down_up(6*k+1,j)=t1down6down(k+1,j);
        t1down6down_up(6*k+2,j)=0;
        t1down6down_up(6*k+3,j)=0;
        t1down6down_up(6*k+4,j)=0;
        t1down6down_up(6*k+5,j)=0;
        t1down6down_up(6*k+6,j)=0;
        
        t1down7down_up(6*k+1,j)=t1down7down(k+1,j);
        t1down7down_up(6*k+2,j)=0;
        t1down7down_up(6*k+3,j)=0;
        t1down7down_up(6*k+4,j)=0;
        t1down7down_up(6*k+5,j)=0;
        t1down7down_up(6*k+6,j)=0;
        
        t1down8down_up(6*k+1,j)=t1down8down(k+1,j);
        t1down8down_up(6*k+2,j)=0;
        t1down8down_up(6*k+3,j)=0;
        t1down8down_up(6*k+4,j)=0;
        t1down8down_up(6*k+5,j)=0;
        t1down8down_up(6*k+6,j)=0;
        
        
        t2down1down_up(6*k+1,j)=t2down1down(k+1,j);
        t2down1down_up(6*k+2,j)=0;
        t2down1down_up(6*k+3,j)=0;
        t2down1down_up(6*k+4,j)=0;
        t2down1down_up(6*k+5,j)=0;
        t2down1down_up(6*k+6,j)=0;
        
        t2down2down_up(6*k+1,j)=t2down2down(k+1,j);
        t2down2down_up(6*k+2,j)=0;
        t2down2down_up(6*k+3,j)=0;
        t2down2down_up(6*k+4,j)=0;
        t2down2down_up(6*k+5,j)=0;
        t2down2down_up(6*k+6,j)=0;
        
        t2down3down_up(6*k+1,j)=t2down3down(k+1,j);
        t2down3down_up(6*k+2,j)=0;
        t2down3down_up(6*k+3,j)=0;
        t2down3down_up(6*k+4,j)=0;
        t2down3down_up(6*k+5,j)=0;
        t2down3down_up(6*k+6,j)=0;
        
        t2down4down_up(6*k+1,j)=t2down4down(k+1,j);
        t2down4down_up(6*k+2,j)=0;
        t2down4down_up(6*k+3,j)=0;
        t2down4down_up(6*k+4,j)=0;
        t2down4down_up(6*k+5,j)=0;
        t2down4down_up(6*k+6,j)=0;
        
        t2down5down_up(6*k+1,j)=t2down5down(k+1,j);
        t2down5down_up(6*k+2,j)=0;
        t2down5down_up(6*k+3,j)=0;
        t2down5down_up(6*k+4,j)=0;
        t2down5down_up(6*k+5,j)=0;
        t2down5down_up(6*k+6,j)=0;
        
        t2down6down_up(6*k+1,j)=t2down6down(k+1,j);
        t2down6down_up(6*k+2,j)=0;
        t2down6down_up(6*k+3,j)=0;
        t2down6down_up(6*k+4,j)=0;
        t2down6down_up(6*k+5,j)=0;
        t2down6down_up(6*k+6,j)=0;
        
        t2down7down_up(6*k+1,j)=t2down7down(k+1,j);
        t2down7down_up(6*k+2,j)=0;
        t2down7down_up(6*k+3,j)=0;
        t2down7down_up(6*k+4,j)=0;
        t2down7down_up(6*k+5,j)=0;
        t2down7down_up(6*k+6,j)=0;
        
        t2down8down_up(6*k+1,j)=t2down8down(k+1,j);
        t2down8down_up(6*k+2,j)=0;
        t2down8down_up(6*k+3,j)=0;
        t2down8down_up(6*k+4,j)=0;
        t2down8down_up(6*k+5,j)=0;
        t2down8down_up(6*k+6,j)=0;
        
        t3down1down_up(6*k+1,j)=t3down1down(k+1,j);
        t3down1down_up(6*k+2,j)=0;
        t3down1down_up(6*k+3,j)=0;
        t3down1down_up(6*k+4,j)=0;
        t3down1down_up(6*k+5,j)=0;
        t3down1down_up(6*k+6,j)=0;
        
        t3down2down_up(6*k+1,j)=t3down2down(k+1,j);
        t3down2down_up(6*k+2,j)=0;
        t3down2down_up(6*k+3,j)=0;
        t3down2down_up(6*k+4,j)=0;
        t3down2down_up(6*k+5,j)=0;
        t3down2down_up(6*k+6,j)=0;
        
        t3down3down_up(6*k+1,j)=t3down3down(k+1,j);
        t3down3down_up(6*k+2,j)=0;
        t3down3down_up(6*k+3,j)=0;
        t3down3down_up(6*k+4,j)=0;
        t3down3down_up(6*k+5,j)=0;
        t3down3down_up(6*k+6,j)=0;
        
        t3down4down_up(6*k+1,j)=t3down4down(k+1,j);
        t3down4down_up(6*k+2,j)=0;
        t3down4down_up(6*k+3,j)=0;
        t3down4down_up(6*k+4,j)=0;
        t3down4down_up(6*k+5,j)=0;
        t3down4down_up(6*k+6,j)=0;
        
        t3down5down_up(6*k+1,j)=t3down5down(k+1,j);
        t3down5down_up(6*k+2,j)=0;
        t3down5down_up(6*k+3,j)=0;
        t3down5down_up(6*k+4,j)=0;
        t3down5down_up(6*k+5,j)=0;
        t3down5down_up(6*k+6,j)=0;
        
        t3down6down_up(6*k+1,j)=t3down6down(k+1,j);
        t3down6down_up(6*k+2,j)=0;
        t3down6down_up(6*k+3,j)=0;
        t3down6down_up(6*k+4,j)=0;
        t3down6down_up(6*k+5,j)=0;
        t3down6down_up(6*k+6,j)=0;
        
        t3down7down_up(6*k+1,j)=t3down7down(k+1,j);
        t3down7down_up(6*k+2,j)=0;
        t3down7down_up(6*k+3,j)=0;
        t3down7down_up(6*k+4,j)=0;
        t3down7down_up(6*k+5,j)=0;
        t3down7down_up(6*k+6,j)=0;
        
        t3down8down_up(6*k+1,j)=t3down8down(k+1,j);
        t3down8down_up(6*k+2,j)=0;
        t3down8down_up(6*k+3,j)=0;
        t3down8down_up(6*k+4,j)=0;
        t3down8down_up(6*k+5,j)=0;
        t3down8down_up(6*k+6,j)=0;