伸缩.m

用MATLAB实现了该算法

%function F2=renew(F1,N2,N2_1,L)
function F2=renew(F1,L,L_1);
%F1待处理音频
%F2恢复音频
%N2水印音频F变换前两个同步信号间的已知音频长度
%N2_1经DA/AD处理后的两个同步信号的音频长度
%L恢复音频F2长度
%L_1待处理音频长度
F2=zeros(1,L);

%a=N2_1/N2;
a=(L_1-1)/(L-1);

F2(1)=F1(1);
F2(L)=F1(L);

for i=2:L-1
    p=a*i;
    q=floor(a*i);
    b=p-q;
    F2(i)=(1-b)*F1(q)+b*F1(q+1);
end;