function out=selfdwt(in,J)
% 完成J尺度Haar小波變換
% in 是輸入信號(hào)
% out 是輸出的小波系數(shù)
%h0,h1是相應(yīng)的小波和尺度向量
h0=[1/sqrt(2),1/sqrt(2)]; 
h1=[-1/sqrt(2),1/sqrt(2)];
out=zeros(size(in));
m=size(in);%此處的m為一向量
x=length(in);
k=0;
l=0;
Wl=in;
for i=1:J
%此處對(duì)應(yīng)于高頻部分，結(jié)果中抽取偶數(shù)點(diǎn)
    Wh=conv(h1,Wl);
    for j=1:(length(h1)+length(Wl)-1)
        if mod(j,2)==0
            k=k+1;
            out(round(x/2)+k)=Wh(j);
        end
    end
 %下面對(duì)應(yīng)于低頻部分，結(jié)果中抽取偶數(shù)點(diǎn)
    Wl1=conv(h0,Wl);
    for j=1:(length(h0)+length(Wl)-1)
        if mod(j,2)==0
            l=l+1;
            out(l)=Wl1(j);
        end
    end
    
    x=round(x/2);
    m=round(m/2);
    Wl=zeros(m);
    for j=1:x
        Wl(j)=out(j);
    end
    k=0;
    l=0;
end