function f=biorh1(c,d);%實現二次樣條的一維重構，f是重構信號，c是分解后的低頻系數，d是分解后的高頻系數（喬世杰）
h1=sqrt(2)*[-0.0061 0.0869 -0.5798 -0.5798 0.0869 -0.0061];
g1=sqrt(2)*[0 -0.1250 0.3750 -0.3750 0.1250 0];
n1=size(c);
N=n1(2);%N是c的長度
n2=size(d);
L1=size(h1);
L2=size(g1);
if n1(2)==n2(2)%只針對低頻系數與高頻系數長度相等的情況
    for n=0:2*n2(2)-1
        t=0;
        s=0;
        for k=ceil((n-3)/2):floor(n/2+1)%mallat算法中，由h和g的支撐長度算出的參與計算的低頻系數和高頻系數的下標范圍
            if k<0
                tt=c(-k); %低頻左邊界關于0半樣本對稱延拓
                ss=-d(-k);
            elseif k>N-1 
                tt=c(2*N-k); %低頻右邊界關于N-1半樣本對稱延拓
                ss=-d(2*N-k);
            else tt=c(k+1);
                ss=d(k+1);
            end
            t=t+tt*h1(n-2*k+3) %低頻部分對應項相乘
            s=s+ss*g1(n-2*k+3) %高頻部分對應項相乘
        end
        f(n+1)=t+s
    end
end