matlab中如何编fft函数_Matlab中对时域函数作傅里叶变换 即FFT函数的使用

1.这是通常情况下，输入一个多频率成分的函数，然后得到傅里叶变换

clear;clc;
Fs = 10000;           % 采样频率，要大于时域信号中最大频率的两倍
T = 1/Fs;             % 采样周期，或者说是时域信号的时间步长，因为你处理的信号不会是连续的

L = 1000;             % 信号长度，或者说你选取的信号离散点的个数，就是总共1000个时间点 由此知，频率分辨率为 1hz
t = (0:L-1)*T;        % 时间相量，可以看成时域信号的横坐标，类似于1dt,2dt,3dt，直到（L-1）t
% 这意味着，你的Fs越大，或者T越小，得到的结果越准确。

%*****************************************************************
% 这一段是自己假设一个输入时域信号
f0=10;                % 频率值
y1 = sin(2*pi*f0*t);  % 时域函数
S = y1;               % 原函数
%*****************************************************************

Y = fft(S);           % 调用fft(x)，不用fft(x,n)，据说前者更准确

P2 = abs(Y/L);        % 每个量除以数列长度 L
                      % **********(为什么要除以L？)********
P1 = P2(1:L/2+1);     % 取交流部分，这个L/2+1是对的ÿ