matlab帮助文档-fft（快速傅里叶变换）

matlab帮助文档-fft（快速傅里叶变换）

fft语法为：

• Y = fft(X)
• Y = fft(X,n)
• Y = fft(X,n,dim)

解释：

Y = fft(X) ，返回 X的离散傅里叶变换（DFT）

• 若X为列向量，对该列向量做DFT。
• 若X为矩阵，对X的每一列分别做DFT。
• 若X为多维数组，对X的每一页的每一列分别做DFT。

---------------------------------------------------------------------------------------------------------------

Y = fft(X,n) ，若X点的个数少于n，则在后面加0，凑到n。若X的点的个数多于n，则删除多余的数。其余与上相同。

---------------------------------------------------------------------------------------------------------------

Y = fft(X,n,dim) ，若dim=1，则按照列向量进行DFT。若dim=2，则按照行向量进行DFT。若dim>2,也可以计算，但是还不知是什么~

---------------------------------------------------------------------------------------------------------------

例子： 使用fft 计算一段波的 基本信息

同时我们可以给波增加小部分噪声  

计算波形的频率、幅值、相位分量。

原始函数为：S = 0.7*sin(2*pi*50*t+pi/4) + sin(2*pi*120*t-pi/2)

利用fft计算其中的 频率、幅值、相位。

clear;clc
Fs = 1000;            % 采样频率
T = 1/Fs;             % 采样周期
L = 1000;             % 信号长度   由此知，频率分辨率为 1hz
t = (0:L-1)*T;        % 时间相量
S = 0.7*sin(2*pi*50*t+pi/4) + sin(2*pi*150*t-pi/2);   % 原始函数

S = S + 2*randn(size(t)); %增加噪声 -------------------------------------------

Y = fft(S);

P2 = abs(Y/L);  % 每个量除以数列长度 L
P1 = P2(1:L/2+1);  % 取交流部分
P1(2:end-1) = 2*P1(2:end-1); % 交流部分模值乘以2
f = Fs*(0:(L/2))/L;
plot(f,P1)
title('Single-Sided Amplitude Spectrum of S(t)')
xlabel('f (Hz)')
ylabel('|P1(f)|')

f1=Y(1:L/2+1);
f1=atan2(imag(f1),real(f1))+pi/2;
%需要将fft的结果加上pi/2才是真实的值，实验出来的，原理未明。
theta = [f1(51) f1(151)]
%结果为弧度制

-------------------------计算结果---------------------------------------
theta =

    0.7854    4.7124（减去2pi，则为-pi/2）