通常我们在做处理信号的时候，会需要使用重采样来提高采样率是信号能够跟好的分析及处理。

具体在语音处理的领域，我们有专用的函数，但是在图像领域的插值算法同样适用于语音的重采样。

一般有三种常用的方法来进行插值：

1. 临近插值法

2. 线性插值法

3. 立方插值法

对于语音信号当然我们还可以使用平方插值法，与图像不同，分块的语音信号，如果在边沿过度不平滑时会导致听觉上的问题，比如有兹兹的噪声或者不连贯感。

先写上：三个插值函数

线性：

function signal_out=LinearResample(signal_in,length_desire,sample_last)

length_desire=length_desire+1; % For Padding

signal_append=[sample_last(end-1:end);signal_in];

signal_in=[sample_last(end);signal_in];

length_signal=length(signal_in);

size_step=(length_signal-1)./(length_desire-1);

signal_out=zeros(length_desire-1,1);

for i=2:length_desire

x(2,1)=ceil(round(((i-1).*size_step+1)*10)./10)+1;

x(1,1)=x(2,1)-1;

y=signal_append(x);

c=LinearFun(x,y);

a=(i-1).*size_step+2; % To suitable the Appended Function

a=[a,1];

b=a*c;

signal_out(i-1,1)=b;

end

end

function c=LinearFun(x,y)

x=[x,[1;1]];

c=inv(x)*y;

end平方：

function signal_out=QuadraticResample(signal_in,length_desire,sample_last)

length_desire=length_desire+1;

signal_append=[sample_last(end-2:end);signal_in];

signal_in=[sample_last(end);signal_in];

length_signal=length(signal_in);

size_step=(length_signal-1)./(length_desire-1);

signal_out=zeros(length_desire-1,1);

for i=2:length_desire

x(3,1)=ceil(round(((i-1).*size_step+1)*10)./10)+2;

x(2,1)=x(3,1)-1;

x(1,1)=x(3,1)-2;

y=signal_append(x);

c=QuadraticFun(x,y);

a=(i-1).*size_step+3; % To suitable the Appended Function

a=[a.^2,a,1];

b=a*c;

signal_out(i-1,1)=b;

end

end

function c=QuadraticFun(x,y)

x=[x.^2,x,[1;1;1]];

c=inv(x)*y;

end立方：

function signal_out=CubicResample(signal_in,length_desire,sample_last)

length_desire=length_desire+1;

signal_append=[sample_last(end-3:end);signal_in];

signal_in=[sample_last(end);signal_in];

length_signal=length(signal_in);

size_step=(length_signal-1)./(length_desire-1);

signal_out=zeros(length_desire-1,1);

for i=2:length_desire

x(4,1)=ceil(round(((i-1).*size_step+1)*10)./10)+3;

x(3,1)=x(4,1)-1;

x(2,1)=x(4,1)-2;

x(1,1)=x(4,1)-3;

y=signal_append(x);

c=CubicFun(x,y);

a=(i-1).*size_step+4; % To suitable the Appended Function

a=[a.^3,a.^2,a,1];

b=a*c;

signal_out(i-1,1)=b;

end

end

function c=CubicFun(x,y)

x=[x.^3,x.^2,x,[1;1;1;1]];

c=inv(x)*y;

end