% 假设 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
2025-11-23 13:26:44
2KB
matlab
1