interp_dfs.m

matlab从入门到精通第6章书后源码对初学者很有帮助

function [xi,yi]=interp_dfs(T,x,Ws,ti)
% T样本间隔，x离散时间样本，Ws归一化频率，ti插值的时间序列
% 判断输入参数的个数和性质
if nargin<4
    ti=5;
end
if nargin<3 | Ws>1
    Ws=1;
end
% 确定数据采样点数据个数
N=length(x);
% 确定采样时间间隔
if length(ti)==1;
    ti=0:T/ti:(N-1)*T;
end
ks=ceil(Ws*N/2);
% 进行fft变换
yi=fft(x);
% 筛选时间数据序列
yi(ks+2:N-ks)=zeros(1,N-2*ks-1);
xi=zeros(1,length(ti));
for k=2:N/2
    xi=xi+yi(k)*exp(j*2*pi*(k-1)*ti/N/T);
end
% 给出根据公式计算的插值结果
xi=real(2*xi+yi(1)+yi(N/2+1)*cos(pi*ti/T))/N;