快速傅立叶变换FFT的实现

最新推荐文章于 2021-03-21 01:10:39 发布

zhengqijun_

最新推荐文章于 2021-03-21 01:10:39 发布

阅读量1.5k

点赞数

分类专栏： Matlab 文章标签：快速傅立叶变换的实现 FFT的实现 FFT频谱图

本文链接：https://blog.csdn.net/zhengqijun_/article/details/53333000

版权

该博客通过Matlab程序展示了如何对连续信号x(t) = cos(200πt) + sin(100πt) + cos(50πt)进行快速傅立叶变换(FFT)，并得到了符合预期的频谱分析结果。

摘要由CSDN通过智能技术生成

使用Matlab软件来做频谱分析。

题目：已知连续信号x(t) = cos(200πt) + sin(100πt) + cos(50πt) ，用FFT对该信号做频谱分析。

Matlab参考程序如下：

clear; close all
fs = 400; T = 1 / fs;
Tp = 0.04; N = Tp * fs;
N1 = [N,4*N,8*N];
st = ['|X1(jf)|';'|X4(jf)|';'|X8(jf)|'];
for m = 1 : 3
    n = 1 : N1(m);
    xn = cos(200*pi*n*T) + sin(100*pi*n*T) + cos(50*pi*n*T) + cos(50*pi*n*T);
    Xk = fft(xn,4096);
    fk = [0:4095] / 4096 / T;
    subplot(3,2,2*m-1)
    plot(fk,abs(Xk)/max(abs(Xk)));ylabel(st(m,:))
    if m == 1
        title('矩形窗截取');
    end
end
for m = 1 : 3
    n = 1 : N1(m);
    wn = hamming(N1(m));
    xn = (cos(200*pi*n*T) + sin(100*pi*n*T) + cos(50*pi*n*T)).*wn';
    Xk = fft(xn,4096);
    fk