Fs = 1000; % 采样频率
T = 1/Fs; % 采样周期：0.001s
L = 1500; % 信号长度
t = (0:L-1)*T; % 时间向量. 时间向量从0开始递增，0s~1.499s
S = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t); % 模拟原信号
X = S + 2*randn(size(t)); % 加入干扰信号
subplot(221);
plot(t,S);
xlabel("时间（s）");
ylabel("幅值");
title("原信号波形");
subplot(222);
plot(t,X);
xlabel("时间（s）");
ylabel("幅值");
title("加入干扰信号波形");
Y=fft(X);
len = length(Y); % Len = L
P2=abs(Y/L); % 计算幅值并进行能量还原
% 751个，P1提取P2中左半部分的第1~751个
P1 = P2(1:L/2+1);
% 绘制单边谱，所以除了首尾元素，其余信号幅度翻倍
P1(2:end-1) = 2*P1(2:end-1);
%根据奈奎斯特定理，此处采样频率Fs为1000，因此频域分析的信号频率最大为500
f = Fs*(0:(L/2))/L; % 0 ~ 500Hz
subplot(223);
plot(f,P1);
xlabel("频率（Hz）");
ylabel("幅值");
title("含有干扰信号的傅里叶变换");运行效果：文章来源：https://www.toymoban.com/news/detail-697854.html
文章来源地址https://www.toymoban.com/news/detail-697854.html
到了这里，关于Matlab信号处理3：fft(快速傅里叶变换)标准使用方式的文章就介绍完了。如果您还想了解更多内容，请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章，希望大家以后多多支持TOY模板网！
