tianchong.m

本源代码主要实现谐波小波的分解及重构的matlab程序。

function out=tianchong(aa,j,t)
%aa为谐波小波分解系数
%j谐波小波分解层数
%t是时间
bb=zeros(length(j),length(t));
cc=length(aa);
%将aa填充到谐波小波时频网格图中
if cc==(2.^(j(length(j))+1)-1)+1
    for i=1:length(j)    
        if i==1
            bb(1,:)=aa(1);
        else
            k=2.^j(i);
            for s=1:length(t)
                bb(i,s)=aa(2.^j(i)+1+fix(k.*s./(length(t)+1)));
            end
        end
    end
else
    for i=1:length(j)-1    
        if i==1
            bb(1,:)=aa(1);
        else
            k=2.^j(i);
            for s=1:length(t)
                bb(i,s)=aa(2.^j(i)+1+fix((k.*s./length(t)+1)));
            end
        end
    end
    bb(length(j),:)=aa(length(a));
end
out=bb;