fft 一维傅里叶变换使用（Matlab）

注：将一个带有噪音的时域信号转化为频域信号

clc,clear;
figure,
Hz=0.02;%周期；也可以说x轴的步长
Fs=1/Hz; %频率

%定义函数
x=0:Hz:10-1/50; 
y=5*sin(2*pi*15*x) + 10*sin(2*pi*20*x)+ 17*cos(2*pi*5*x);%此处可以添加噪音
y=y+10*randn(size(x));
%上面的参数定义好，下面的不需要修改

subplot(2,1,1);
plot(x,y);
title("6*sin(10*2*pi*t)");

Y = fft(y);
L=length(x); % x 轴数据的个数
P2 = abs(Y/L);  % Y 为傅里叶变换得到的结果  L为Y的长度
P1 = P2(1:L/2+1);   %从 P2 取出元素 的 1--(500/2)+1 长度的元素  即，取前一半+1 个元素
P1(1:end-1) = 2*P1(1:end-1);
f = Fs*(0:(L/2))/L;
subplot(2,1,2);
plot(f,P1) ;
title("fft");
xlabel('Hz');
ylabel('振幅');

注：上面一些函数说明：
randn（）函数：https://blog.csdn.net/minyeling/article/details/105741721.