trigintpoly:函数 trigintpoly 使用 fft 来求三角插值多项式-matlab开发

% 假设 f(t) 是区间 [0,2pi] 上的实数 2pi 周期函数% 并且 1*n 向量 x 是函数 f(t) 在 n 处的值% 等距点（n 必须是偶数） % t_j=(j-1)*2*pi/n, j=1,2,...,n。 ％ 功能% [y , yp , ypp] = trigintpoly (x,s) % 使用 fft 找到三角插值多项式% 在 n 个点 t_1,t_2,...,t_n 处对函数 f(t) 进行插值。 那么% 函数 trigintpoly 计算函数 f(t)、f'(t)、 % 和 f''(t) 在点 s（s 是一个 m*1 的点向量），即% y = f(s), yp=f'(s), ypp=f''(s) % % ％示例1： % n = 100; % t = 0:2*pi/n:2*pi-2*pi/n; % x = cos(2.*t).^3; % s = [-pi/4,0,p