add_save.m

来自「重叠保留法函数」· M 代码 · 共 23 行

function y=add_save(h,x);
L1=length(h);
L2=length(x);
if L1>=8 error('length of h should be shorter')
end;
if L2<=16 error('length of x should be longer')
end;
y=[];
m=(L2-mod(L2-8,8-L1+1)-8)/(8-L1+1);
h(8)=0;
H=fft(h);
for M=0:m
    temp=[];
    for n=1:8
        temp=[temp x(n+M*(8-L1+1))];
    end
    TEMP=fft(temp);
    HT=H.*TEMP;
    ht=ifft(HT);
    for k=L1:8
        y=[y ht(k)];
    end
end