一个电力系统信号,其中频率为 200 和 300. 应用付立叶变换并求频谱图
clc;
clf;
clear all;
fs=1000;
t=0:1/fs:0.6;
f1=200;
f2=300;
x=sin(2*pi*f1*t)+sin(2*pi*f2*t);
subplot(411);
plot(x);
title('f1 ( 100Hz ) 2(300Hz) 的正弦信号,初相 0');
xlabel(' 序列( n)');
grid on;
number=1024;
y=fft(x,number);
n=0:length(y)-1;
f=fs*n/length(y);
subplot(412);
plot(f,abs(y)/max(abs(y)));
hold on;
plot(f,abs(fftshift(y))/max(abs(y)),'r');
title('f12 的正弦信号的 FFT ( 1024 点) ');
xlabel(' 频率 Hz');
grid on;
x=x+randn(1,length(x));
subplot(413);
plot(x);
title(' 原 f12 的正弦信号(含随机噪声) ');
xlabel(' 序列( n)');
grid on;
y=fft(x,number);
n=0:length(y)-1;
f=fs*n/length(y);
subplot(414);
plot(f,abs(y)/max(abs(y)));
title(' 原 f12 的正弦信号(含随机噪声)的 FFT ( 1024 点) ');
xlabel(' 频率 Hz');
grid on;