[MATLAB基础] MATLAB应用大全 书连载

[lili123]

2．一维快速傅立叶插值
在MATLAB中，一维快速傅立叶插值通过函数interpft()来实现，该函数利用傅立叶变换将输入数据变换到频域，然后用更多点的傅立叶逆变换，变换回时域，其结果是对数据进行增采样。调用格式如下。
 y=interpft(x, n)：对x进行傅立叶变换，然后采用n点傅立叶反变换，变回到时域。如果x是一个向量，数据x的长度为m，采样间隔为dx，则数据y的采样间隔是dx×m/n，其中n必须大于m。如果x是矩阵，该函数对矩阵x的列进行操作，其返回的结果y与x具有相同的列，行数为n。
 y=interpft(x, n, dim)：在dim指定的维度上进行操作。
【例5-13】 利用一维快速傅立叶插值实现数据增采样，其实现的MATLAB代码如下：

clear all;
x=0:1.1:8;
y=sin(x);
n=2*length(x); %增采样1倍
yi=interpft(y,n); %采用一维快速傅立叶插值
xi=0:0.55:8.3; %要插值的数据
figure; %画图显示
hold on;
plot(x,y,'ro',xi,yi,'b-');
legend('原始数据','插值后结果');

运行程序后，输出的结果如图5.7所示。