makesynthesistable.sci

小波分解源代码

function spkt = MakeSynthesisTable(atomic,pkt)
// MakeSynthesisTable -- Convert atom list to packet table
//  Usage
//    spkt = MakeSynthesisTable(atomic,pkt)
//  Inputs
//    atomic     atoms from wavelet packet or cosine packet table
//    pkt        packet table (only the shape is used)
//  Outputs
//    spkt       packet table filled with coefficients from atom list  
//
//  Copyright Aldo I Maalouf

    a=size(pkt);
    spkt = zeros(a(1),a(2));
    [n,L] = size(pkt);
    D = L-1;
//
//
    [natom,ncharm] = size(atomic);
    for j=1:natom,
        a = atomic(j,1); d = atomic(j,2); 
        b = atomic(j,3); k = atomic(j,4);
        if a ~= 0,
           spkt(pkt2ix(d,b,k,D,n),d+1) = a;
        end    
    end
    endfunction