Code: Selecteer alles
close all; clear all; clc;
fs = 200; % sampling frequentie
f0 = 50; % signaal frequentie
t = 0:1/fs:1; % tijdsas
x = sin(2*pi*f0*t); % signaal
figure; plot(t,x);
X = fft(x)/(length(x)/2); % neem FFT en herschaal
X = X(1:floor(length(X)/2)); % gooi de 2de helft van het spectrum weg
f = linspace(0,1,length(X))*fs/2; % stel de frequentie as in Hz op
figure; plot(f,abs(X));