立方内插的演示

由于AD的采样率一般而言是固定的，可能并不是最佳采样点。在这种情况下，有必要通过内插，得到最佳的采样点。

工程常用的内插函数，为立方采样。具体的定义、表达形式，可以参考《Digital communications: A Discrete-Time Approach》之8.4.2节。

首先是立方插值的函数

function interpolated_value=interpolation_cubic(x,x1,x2,x3,u)
%% 程序说明
% u: 小数间隔 [0 1]之间，位于[x1 x2]之间
% 按照时间先后 [x x1 x2 x3] 
% x1为当前采样点；x为前一个采样点；
% x2为x1的后一个采样点；x3为x1的后两个采样点

c3 = 1/6*u^3 - 1/6*u;          % 对应x3
c2 = -1/2*u^3+1/2*u^2+u;       % 对应x2 
c1 = 1/2*u^3-u^2-1/2*u+1;      % 对应x1
c0 = -1/6*u^3+1/2*u^2-1/3*u;   % 对应x

interpolated_value = x*c0 + x1*c1 + x2*c2 + x3*c3;

然后是利用这个function，来做测试。

clc;
clear;
close all;

fs = 1e3;
n = 0:999; %注意时间是从0开始的，后面for循环要注意统一时间点
t = n*1/fs;
signals = sin(2*pi*10*t);

u = 0.1; %x1和x2之间
sig =[]; %插值
tm = []; %记录真实的插值时间点
for index = 2:900
    x = signals(index - 1);