对fft的一些理解.docxVIP

例子1：对50Hz的正弦波做fft分析，默认加矩形窗 close all; clc; clear; %% Fs = 1e3; % Sampling frequency T = 1/Fs; % Sample time L = 1000; % Length of signal t = (0:L-1)*T; % Time vector x = 1*sin(2*pi*50*t)+0*sin(2*pi*50.5*t); % Signal win = rectwin(L); % Windown x_data = x .* win; % plot Signal figure(Position,[65 278 560 420]); plot(t(1:L),x_data(1:L)) grid on; title(Time Domain Signal) xlabel(time (s)) % fft NFFT = 1*2^nextpow2(L); % Next power of 2 from length of x X = fft(x_data,NFFT)/L; f = (Fs/2)*linspace(0,1,NFF