主要内容

interpft

一维插值(FFT法)

描述

例子

Y=内部英尺(X,N)插值中函数值的傅里叶变换X产生N等距的点。interpft对大小不等于1的第一个维度进行操作。

例子

Y=内部英尺(X,N,昏暗的)沿维度操作昏暗的。例如,如果X是一个矩阵interpft(X,n,2)对行进行操作X

例子

全部崩溃

用FFT方法插值一维数据并将结果可视化。

在间隔中生成一些采样点 [ 0 , 3. π ] 对于函数 F ( x ) = 2. ( x ) 余弦 ( x ) .使用间隔dx为确保数据间隔均匀,请绘制采样点。

dx=3*pi/30;x=0:dx:3*pi;f=sin(x)。^2.*cos(x);绘图(x,f,“o”)

图形包含一个轴。轴包含一个line类型的对象。

用FFT插值法求出200个查询点处的函数值。

N=200;y=内部英尺(f,N);

根据采样点的间距计算插值数据的间距dy = dx *长度(x) / N哪里N是插值点的数目。截断中的数据Y以匹配x2

dy=dx*length(x)/N;x2=0:dy:3*pi;y=y(1:length(x2));

绘制结果。

持有图(x2,y,'.')头衔(“周期函数的FFT插值”)

图中包含一个轴。标题为周期函数FFT插值的轴包含2个线型对象。

生成三个独立的正态分布随机数数据集。假设数据以正整数进行采样,1:N。将数据集存储为矩阵中的行。

一个= randn (20);x = 1:20;

在每个500个查询点处插入矩阵行。指定尺寸=2以便interpft在一排排的A.

N=500;y=interpft(A,N,2);

计算插值数据的间距dy.截断数据Y以匹配x2

dy=长度(x)/N;x2=1:dy:20;y=y(:,1:length(x2));

绘制结果。

子地块(3,1,1)图(x,A(1,:)',“o”);持有图(x2,y(1,:)','--')头衔(“第1排”)子地块(3,1,2)图(x,A(2,:)',“o”);持有图(x2,y(2,:)','--')头衔(“第2排”次要情节(3,1,3)情节(x, (: 3) ',“o”);持有图(x2,y(3,:)','--')头衔(“第3排”)

图中包含3个轴。第一行包含2个类型为line的对象。第2行包含2个类型为line的对象。第3行包含2个类型为line的对象。

输入参数

全部崩溃

输入数组,指定为向量、矩阵或多维数组。中的数据X假设是在自变量的均匀间隔时间内抽样的。interpft适用于周期性数据。

数据类型:单一的|双重的
复数支持:万博1manbetx

点数,指定为正整数标量。

数据类型:单一的|双重的

要沿其操作的维度,指定为正整数标量。如果未指定值,则默认值为大小不等于1的第一个数组维度。

  • interpft(X,n,1)插入的列X

  • interpft(X,n,2)插入的行X

例子:interpft(X,n,3)

输出参数

全部崩溃

插值点,作为向量、矩阵或多维数组返回。如果长度(X,dim)=m,X抽样间隔是dx,则为新的采样间隔Ydy=dx*m/n哪里n>m

如果昏暗的那么,是指定的吗interpft填充或截断X长度N在尺寸上昏暗的,以便尺寸(y,dim)=n

算法

这个interpft函数使用FFT方法。原始向量x使用fft,然后使用更多点将其转换回。

另见

|

之前介绍过的R2006a