fft1.m

基于MATLAB的fft电力谐波分析 matlab代码

%fs=1000;     %采样率为1000
%t=0:1/fs:1;  %时长1s
%f1=50;         %信号频率50Hz
fs=100;%设定采样频率
N=512;
n=0:N-1;
t=n/fs;
f1=20;
x=sin(2*pi*f1*t); %生成随机噪声的正弦波函数
x=x+randn(1,length(x));
subplot(211)
plot(x);
title('原正弦信号(含随机噪声)')
xlabel('序列(n)')
grid on

y=dfft(x);
n=0:length(y)-1;
f=fs*n/length(y);
subplot(212)
plot(f,abs(y));
title('原正弦信号(含随机噪声)的FFT')
xlabel('频率Hz')
grid on