【MATLAB】利用FFT分析其频谱

clc; clear; %% 参数设置 Fs = 100; % Sampling frequency T = 1/Fs; % Sampling period f0 = 5; % 设置的频率 N = 100; % 整周期数 M = N/f0; x = 0:T:M; L = length(x); %% y = sin(2*pi*f0*x); figure(1); subplot(2,1,1); plot(x,y) title('100个整周期的频率为5的正弦波') xlabel('f (Hz)') f = Fs*(0:(L/2))/L; Y = fft(y); P2 = abs(Y/L); P1 = P2(1:L/2+1); P1(2:end-1) = 2*P1(2:end-1); subplot(2,1,2); plot(f,P1) xlabel('f (Hz)') title('100个整周期的频率为5的正弦波的频谱')